The linear form of the parabolized linear stability equations is used in a variational approach to extend the previous body of results for the optimal, nonmodal disturbance growth in boundary-layer flows. This paper investigates the optimal growth characteristics in the hypersonic Mach number regime without any high-enthalpy effects. The influence of wall cooling is studied, with particular emphasis on the role of the initial disturbance location and the value of the spanwise wave number that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary-layer equations, mean flow solutions based on the full Navier-Stokes equations are used in select cases to help account for the viscousinviscid interaction near the leading edge of the plate and for the weak shock wave emanating from that region. Using the full Navier-Stokes mean flow is shown to result in further reduction with Mach number in the magnitude of optimal growth relative to the predictions based on the self-similar approximation to the base flow. Nomenclature A, B, C, D = linear matrix operators c = normalization coefficient E = total energy norm G = energy gain h 1 = streamwise metric factor h 3 = spanwise metric factor J = objective function K = kinetic energy norm K = bilinear concomitant L = flat-plate characteristic length L = Lagrangian function M = Mach number M E = energy weight matrix N = logarithmic amplitude ratio relative to neutral stability location of modal instability N η = number of discretization points along the wall-normal direction q = vector of base flow variables q = vector of perturbation variableŝ q = vector of amplitude variables R = local Reynolds number Re = Reynolds number T = temperature T ad = adiabatic wall temperature T w = wall temperature u; v; w = streamwise, wall-normal, and spanwise velocity components ν = kinematic viscosity x; y; z = Cartesian coordinates α = streamwise wave number β = spanwise wave number γ = heat capacity ratio δ = local similarity length scale of boundary layer δ 0.995 = 0.995 boundary-layer thickness ξ; η; ζ = streamwise, wall-normal, and spanwise coordinates in body-fitted coordinate system ρ = density Ω = domain of integration ω = angular frequency Subscripts ad = adiabatic wall condition r = reference value w = wall condition 0 = initial disturbance location 1 = final optimization location Superscripts H = conjugate transpose T = transpose = dimensional value † = adjoint