Numerical simulations focusing on convective heat transfer with complex wall temperature variations, as well as conjugate heat transfer problems involving solids are performed. Simulations are made using the GASP Computational Fluid Dynamics (CFD)code, which solves the Reynolds-Averaged Navier-Stokes equations (RANS). GASP has been modified by adding a solid heat conduction solver, enabling the coupling of RANS with the heat equation for an isotropic solid. Validation cases involving convective heat transfer are considered using both one-and two-equations turbulence models. Predicted Stanton numbers for low Mach number turbulent boundary layers closely agree with experimental results from Reynolds et al. 1, 2 for a constant wall temperature, a step and a double pulse in wall temperature. Simulations of supersonic boundary layers with a step in wall temperature are also performed. Good agreement is found for velocity and temperature profiles when compared with measurements from Debieve et al. 3 as well as for the skin friction and Stanton number. The flow field and wall temperature distribution inside a supersonic cooled nozzle is computed using a new conjugate heat transfer algorithm in GASP. Results agree well with measurements from Back et al. 4