Titanium alloy Ti-6242 (Ti-6Al-2Sn-4Zr-2Mo) is frequently used in the high-pressure compressor of aero engines due to its excellent resistance to fatigue and creep failure at high temperature. While exhibiting high strength at elevated temperatures, it is susceptible to dwell fatigue at temperatures below 473 K due in part to the presence of microtextured regions (MTRs), also known as macrozones. MTRs are clusters of similarly orientated alpha particles, which form during alpha/beta processing and remain stable even after large deformation. The major objective of this dissertation is to quantify the evolution of MTRs under different thermomechanical processing parameters, and predict the optimal processing parameters to eliminate the MTRs.Idealized MTRs with pure initial orientation are first employed as the benchmark case to investigate the loading direction effect on its breakdown efficiency. Three high-temperature compression processes are simulated with different loading directions using crystal plasticity finite element method, and the results are validated against high-temperature compression experiments and EBSD measurement. The evolution of equivalent plastic strain, accumulated shear strain, and misorientation distribution is analyzed in detail to reveal the relationship between loading direction and MTR breakdown efficiency. Lastly, the reorientation velocity divergence of arbitrary loading direction is expressed in the Rodrigues\u27 space in order to predict the optimal processing parameters for MTR elimination. The MTR breakdown efficiency also depends on the morphology and its position within the specimen. Two different length scales have to be analyzed in order to consider both factors, which present great challenge to the numerical simulation. In this dissertation, a high-efficient FE-FFT multiscale modeling framework is derived and developed to overcome this challenge. The Fourier-Galerkin method is utilized to solve the microscale unit cell problem, while total Lagrangian nite element is used to solve the macroscopic boundary value problems. Several numerical improvements are derived and implemented to further improve its numerical efficiency, including consistent linearization, consistent homogenized tangent stiffness, and inexact Newton method. A series of numerical studies is conducted to investigate the accuracy, efficiency, and robustness of this algorithm