Search NASA⌕ Search

Engineering topics

Lian Duan

Publications and source records attributed to Lian Duan.

At least 19 records

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary layer transition↗

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary-layer transition↗

Simulation and Modeling of Cold-Wall Hypersonic Turbulent Boundary Layers on Flat Plate

Direct numerical simulations (DNS) of flat-plate, zero-pressure-gradient turbulent boundary layers are presented for nominal freestream Mach numbers of 11 and 14 and a highly cooled wall (wall-to-recovery temperature of approximately 0.2). The flow conditions of the DNS are representative of the experimental data for a Mach 11.1 turbulent boundary layer on a flat plate that was tested at Calspan–University of Buffalo Research Center (CUBRC) and the operational conditions of the AEDC Hypervelocity Tunnel No. 9 at Mach 14. The wall shear stress and turbulent heat flux predicted by DNS show good comparisons with those measured at CUBRC and those modeled with a Baldwin-Lomax model. The generated DNS database extends the range of Reynolds numbers from those typical of previous DNS studies for hypersonic boundary-layer flows and confirms the validity of Morkovin’s hypothesis within the relatively broader range of Reynolds numbers considered in this study. A priori assessment of turbulent heat flux models with DNS data shows that the algebraic energy flux (AEF) formulation gives better predictions of the streamwise turbulent heat flux than the commonly used constant turbulent Prandtl number model while both models give similarly good predictions of the wall-normal turbulent heat flux. A posteriori assessment shows that the AEF model improves the prediction of the temperature profiles while keeping a good prediction of the velocity field.

Junji Huang↗

Effect of 3D Roughness Patch on Instability Amplification in a Supersonic Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of high-speed vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, sinusoidal shape roughness elements on a Mach 3.5 flat plate boundary layer for flow conditions corresponding to the planned conditions of an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at the NASA Langley Research Center. Analysis of convective instabilities in the wake of the roughness patch was reported in a previous paper and the current work extends that analysis to instability amplification across the length of the roughness patch. Quasi-parallel stability analysis of the modified boundary layer flow over the patch indicates two dominant families of unstable disturbances, namely, a group of high frequency modes that amplify in localized regions along the roughness patch and another group of low frequency modes that have smaller peak amplification rates but amplify steadily both above the roughness patch and in the wake region behind it. The results suggest that the amplification factors associated with the high-frequency modes are sufficiently low, at least for the roughness patches considered in this paper, so that these modes are unlikely to have a major influence on the transition process. The amplification of the low-frequency modes within the region of the roughness patch is further quantified via direct numerical simulations. Results confirm the strongly destabilizing influence of the roughness patch on the first mode instabilities, yielding an N-factor increment of ΔN ≈ 3.6 for a roughness patch length of eight wavelengths.

Laminar turbulent transition↗

Nonlinear Evolution and Breakdown of Azimuthally Compact Crossflow Vortex Pattern over a Yawed Cone

Hypersonic boundary-layer flows over a circular cone at moderate angle of incidence can support strong crossflow instability in between the windward and leeward rays on the plane of symmetry. Due to the more efficient excitation of stationary crossflow vortices by surface roughness, a possible path to transition in such flows corresponds to rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. In the present paper, the previous analyses of crossflow instability over a 7- degree half-angle, yawed circular cone in a Mach 6 free stream have been extended to the nonlinear evolution of azimuthally localized crossflow vortex packets and the amplification characteristics and nonlinear breakdown of high-frequency secondary instabilities associated with those packets. A comparison between plane marching PSE and direct Navier-Stokes simulations (DNS) reveals favorable agreement in regard to mode shapes, most amplified disturbance frequencies, and N-factor evolution. In contrast, the quasi-parallel predictions are found to result in severe underprediction of the N-factors. The direct numerical simulations also indicate that the breakdown of secondary instabilities in a 3D hypersonic boundary layer shares certain common features with the previous computations of crossflow transition over subsonic swept wings.

Vortex structure↗

Pressure Fluctuations Induced by a Hypersonic Turbulent Boundary Layer

Direct numerical simulations (DNS) are used to examine the pressure fluctuations generated by a spatially-developed Mach 5.86 turbulent boundary layer. The unsteady pressure field is analyzed at multiple wall-normal locations, including those at the wall, within the boundary layer (including inner layer, the log layer, and the outer layer), and in the free stream. The statistical and structural variations of pressure fluctuations as a function of wall-normal distance are highlighted. Computational predictions for mean velocity pro les and surface pressure spectrum are in good agreement with experimental measurements, providing a first ever comparison of this type at hypersonic Mach numbers. The simulation shows that the dominant frequency of boundary-layer-induced pressure fluctuations shifts to lower frequencies as the location of interest moves away from the wall. The pressure wave propagates with a speed nearly equal to the local mean velocity within the boundary layer (except in the immediate vicinity of the wall) while the propagation speed deviates from the Taylor's hypothesis in the free stream. Compared with the surface pressure fluctuations, which are primarily vortical, the acoustic pressure fluctuations in the free stream exhibit a significantly lower dominant frequency, a greater spatial extent, and a smaller bulk propagation speed. The freestream pressure structures are found to have similar Lagrangian time and spatial scales as the acoustic sources near the wall. As the Mach number increases, the freestream acoustic fluctuations exhibit increased radiation intensity, enhanced energy content at high frequencies, shallower orientation of wave fronts with respect to the flow direction, and larger propagation velocity.

Lian Duan↗

Secondary Instability of Stationary Crossflow Vortices in Mach 6 Boundary Layer Over a Circular Cone

Hypersonic boundary layer flows over a circular cone at moderate incidence can support strong crossflow instability. Due to more efficient excitation of stationary crossflow vortices by surface roughness, such boundary layer flows may transition to turbulence via rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. The amplification characteristics of these secondary instabilities are investigated for crossflow vortices generated by an azimuthally periodic array of roughness elements over a 7-degree half-angle circular cone in a Mach 6 free stream. Depending on the local amplitude of the stationary crossflow mode, the most unstable secondary disturbances either originate from the second (i.e., Mack) mode instabilities of the unperturbed boundary layer or correspond to genuine secondary instabilities that reduce to stable disturbances at sufficiently small amplitudes of the stationary crossflow vortex. The predicted frequencies of dominant secondary disturbances are similar to those measured during wind tunnel experiments at Purdue University and the Technical University of Braunschweig, Germany.

Boundary layer flow↗

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models.The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar convexwall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University. The favorable pressure gradient is induced bythe streamwise surface curvature. The DNS results will be compared with experimental measurements and RANS computations for the same flow conditions and wall geometry. The DNS datasets will also be analyzed to provide an improved physical understanding of hypersonic boundary-layer turbulence subject to mechanical nonequilibrium, including the validity of Morkovin’s hypothesis and the various available compressibility transformations, as well as to obtain turbulence statistics that are relevant to RANS modeling, including the distribution ofeddy viscosity, the budgets of turbulent kinetic energy and Reynolds stresses, and the velocity-temperature correlations. The RANS predictions will be obtained by using both the commonly used one- and two-equation models that assume a constant turbulent Prandtl number and an algebraic, variable Prandtl number model for the turbulent energy flux.

turbulent boundary layers↗

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.

turbulent boundary layers↗

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Adverse Pressure Gradients due to Concave Streamline Curvature

Direct numerical simulations (DNS) of adverse-pressure-gradient turbulent boundary layers over a planar concave wall are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar concave wall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University (TAMU). The DNS was validated against the experimental results of TAMU for the same flow conditions and wall geometry. An analysis of the DNS datasets was also conducted to provide an assessment of the validity of Morkovin’s hypothesis and the strong Reynolds analog for turbulence subject to mechanical nonequilibrium. In addition to the DNS results, RANS predictions are obtained by using the Baldwin-Lomax (BL), Spalart-Allmaras (SA), and the k - w SST turbulence models. The comparisons between RANS and DNS showed little impact of an adverse pressure gradient on the accuracy of these models, at least up to an incompressible Clauser pressure gradient parameter of beta(sub inc) 1.22. While the Boussinesq assumption provided reasonable predictions for the Reynolds shear stress, it failed to adequately predict the normal components of the Reynolds stress.

turbulent boundary layers↗