Engineering topics

Li, Fei

Publications and source records attributed to Li, Fei.

At least 19 records

DOE OSTI2024

A Scalable, High-Efficiency, Low-Energy-Spread Laser Wakefield Accelerator Using a Tri-Plateau Plasma Channel

The emergence of multi-petawatt laser facilities is expected to push forward the maximum energy gain that can be achieved in a single stage of a laser wakefield acceleration (LWFA) to tens of giga-electron volts, which begs the question—is it likely to impact particle physics by providing a truly compact particle collider? Colliders have very stringent requirements on beam energy, acceleration efficiency, and beam quality. In this article, we propose an LWFA scheme that can for the first time simultaneously achieve hitherto unrealized acceleration efficiency from the laser to the electron beam of >20% and a sub-1% energy spread using a stepwise plasma structure and a nonlinearly chirped laser pulse. Three-dimensional high-fidelity simulations show that the nonlinear chirp can effectively mitigate the laser waveform distortion and lengthen the acceleration distance. This, combined with an interstage rephasing process in the stepwise plasma, can triple the beam energy gain compared to that in a uniform plasma for a fixed laser energy, thereby dramatically increasing the efficiency. A dynamic beam loading effect can almost perfectly cancel the energy chirp that arises during the acceleration, leading to the sub-percent energy spread. This scheme is highly scalable and can be applied to petawatt LWFA scenarios. Scaling laws are obtained, which suggest that electron beams with parameters relevant for a Higgs factory could be reached with the proposed high-efficiency, low-energy-spread scheme.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Linking NH$^+_4$ motion to magnetism in molecular multiferroic (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]: A neutron vibrational spectroscopy study

Here, we present a neutron vibrational spectroscopy study to investigate the influence of NH$^+_4$ motion on the magnetism in [(NH 4 ) 1–x K x ] 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]. The parent compounds, (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)] (x = 0) and K 2 ⁡[FeCl 5 ⁢(H 2 ⁢O)]⁢(x = 1) are isostructural at room temperature, yet displaying drastically different magnetic and multiferroic behavior. K 2 ⁡[FeCl 5 ⁢(H 2 ⁢O)] is nonmultiferroic with type-A collinear antiferromagnetic structure below T N ≈ 14.06 K, whereas (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)] is a type-II multiferroic with incommensurate cycloidal spin structure below T FE ≈ 6.8 K. A recent study of the dielectric, structure, and magnetic properties in the mixed [(NH 4 ) 1–x ⁢K x ] 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)] shows that a small amount of potassium substitution to replace NH$^+_4$ transforms the spin structure from incommensurate cycloidal (x ≤ 0.06) into commensurate collinear antiferromagnetic (x ≥ 0.15), indicating NH$^+_4$ is essential to the emergent phenomena observed in this molecular multiferroic compound. Our vibrational spectroscopy study reveals that NH$^+_4$ libration and torsion motion exhibit substantial temperature dependence at low temperatures. The intensity of NH$^+_4$ libration and torsion modes increases slightly at 5 K in comparison with data at 25 K behaving like a magnon, indicating that they are coupled to the magnetism in (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]. Comparing data of x=0, 0.06, 0.09, and 0.15 samples further illustrates that the strength of the increased signal in NH$^+_4$ libration mode is very sensitive to potassium concentration. The signal diminishes quickly with increasing potassium concentration and vanishes in the x = 0.15 sample corresponding to the magnetic structure change for x ≥ 0.15. The results directly link the anomalous behavior in NH$^+_4$ libration motion to the magnetism in [(NH 4 ) 1–x ⁢K x ] 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)], providing new insights into the crucial role NH$^+_4$ plays in the coupled phenomena in (NH 4 ) 2 ⁢[FeCl 5 ⁢(H 2 ⁢O)]. The unique information opens a new door to go through in searching for new multifunctional materials by incorporation of NH 4 via a material-by-design approach.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Emittance preservation for the electron arm in a single PWFA-LC stage using quasi-adiabatic plasma density ramp matching sections

Plasma-based acceleration (PBA) is being considered for a next generation linear collider (LC). In some PBA-LC designs for the electron arm, the extreme beam parameters are expected to trigger background ion motion within the witness beam, which can lead to longitudinally varying nonlinear focusing forces and result in an unacceptable emittance growth of the beam. To mitigate this, we propose to use quasi-adiabatic plasma density ramps as matching sections at the entrance and exit of each stage. We match the witness electron beam to the low density plasma entrance, where the beam initially has a large matched spot size so the ion motion effects are relatively small. As the beam propagates in the plasma density upramp, it is quasi-adiabatically focused, and its distribution maintains a non-Gaussian equilibrium distribution in each longitudinal slice throughout the process, even when severe ion collapse has occurred. This only causes small amounts of slice emittance growth. The phase mixing between slices with different betatron frequencies leads to additional projected emittance growth within the acceleration stage. A density downramp at the exit of an acceleration section can eliminate much of the slice and projected emittance growth as the beam and ion motion adiabatically defocuses and decreases, respectively. Simulation results from QuickPIC with Azimuthal Decomposition show that within a single acceleration stage with a 25 GeV energy gain, this concept can limit the projected emittance growth to only ∼2% for a 25 GeV, 100 nm emittance witness beam and ∼20% for a 100 GeV, 100 nm normalized emittance witness beam. The trade-off between the adiabaticity of the plasma density ramp and the initial ion motion at the entrance for a given length of the plasma density ramp is also discussed.

Physics↗

NASA NTRS2019

Effects of Wall Cooling on Supersonic Modes in High-Enthalpy Hypersonic Boundary Layers over a Cone

High-enthalpy hypersonic boundary-layer transition plays an important role in many entry/descent vehicles. Aerothermodynamic performance of these vehicles strongly depends on the transition location on the surface. However, detailed transition flow physics in these chemically reacting boundary layers are poorly understood and transition estimates during the design phase rely heavily on empirically-derived transition criteria such as 𝑹𝒆(𝜽)/𝑴(𝒆). One of the most intriguing characteristics in hypersonic boundary layers is the presence of unstable supersonic modes, first identified in the 1990s. Due to a recent surge in hypersonic applications, there has been renewed interest in studying the flow physics of supersonic modes using either theory or direct numerical simulations. This paper investigates the rise of supersonic modes in a high-enthalpy hypersonic flow over a 𝟓° half-angle cone at various wall temperatures using both quasiparallel linear stability theory and linear PSE. It was found that supersonic modes exist in all wall temperature conditions including the adiabatic wall case. Cooler wall temperature causes the second Mack mode to become an unstable supersonic mode naturally downstream of the upper-branch neutral location when the wall is sufficiently cooled. In terms of the integrated growth, the second mode is still the dominant mode. Nonetheless, supersonic modes can cause additional series of relatively weaker growth than that of the second mode beyond the peak amplitude location. According to the present linear PSE results, contrary to what was speculated in the literature, supersonic-mode pressure disturbance structure radiated into the freestream is nonacoustic in nature and the formation of unstable supersonic modes is mainly associated with the synchronization of phase speed between the instability and acoustic waves in a nonp ayer, not due to nonlinear modal interaction as suggested in the literature.

Chang, Chau-Lyan↗

NASA NTRS2019

Stabilization of a Swept-Wing Boundary Layer by Discrete Roughness Elements at High Reynolds Numbers

Direct numerical simulations (DNS) are performed to study potential stabilizing ef- fect of spanwise periodic discrete roughness elements (DREs) on cross ow instabilities in a spatially developing three-dimensional boundary layer over an in nite-swept natural- laminar- ow wing at a freestream Mach number of 0:75 and a chord Reynolds number of approximately 25 million. In the DNS, both the spanwise periodic DREs and distributed roughness in the leading-edge region are implemented to simulate a typical experimen- tal scenario in which multiple steady cross ow modes including the most unstable mode (i.e., the \target" mode) emerge because of the presence of naturally distributed surface roughness in the leading edge region and spanwise periodic control cylinders of subcritical wavelength are used to force small-wavelength disturbances (i.e., the control mode) for damping the target mode. The DNS results show that the e ectiveness of DRE control is sensitive to roughness diameter, height, and chordwise placement. For the DRE parame- ters considered in this study, the stabilizing e ect on the target mode is small within the computational domain that ended at about 35% of the chord.

Nicholson, Gary↗

NASA NTRS2018

Direct Numerical Simulation of Receptivity to Roughness in a Swept-Wing Boundary Layer at High Reynolds Numbers

Direct numerical simulations (DNS) are performed to examine the receptivity to roughness in a spatially developing three-dimensional boundary layer over an in finite-swept natural-laminar-flow wing at a free stream Mach number of 0:75 and a chord Reynolds number of approximately 25 million based on the long, swept chord. Stationary cross ow disturbances are excited by applying either critically spaced discrete cylinders of micron size or naturally occurring distributed roughness in the leading-edge region. The DNS data show that the spanwise spectral content of the excited cross ow disturbances is highly dependent upon the shape of roughness elements, and the initial growth of the cross ow structures is a nonlinear function of the element height. The linear growth rate of the excited cross ow disturbances predicted by DNS shows good agreement with linear parabolized stability equations. The receptivity study lays the foundation for investigating the stabilization of the naturally most unstable steady cross ow mode by using spanwise periodic DREs.

Nicholson, Gary↗

NASA NTRS2018

Numerical Investigation of Roughness Effects on Transition on Spherical Capsules

The state of the boundary layer on space re-entry vehicles significantly affects the design of the thermal protection system. However, the physical mechanism that leads to the laminar-turbulent boundary-layer transition on blunt spherical capsules remains an open question in literature. This work numerically assesses the potential of roughness-induced non-modal disturbance growth on re-entry capsules with a spherical-section forebody by optimal transient-growth theory and direct numerical simulation. Two different sets of wind-tunnel experiments are considered. Optimal transient-growth studies have been performed for the blunt capsule experiments at Mach 5.9 in the Hypersonic Ludwieg tube Braunschweig (HLB) of the Technische Universitat Braunschweig. In some of these measurements, the capsule model was equipped with a specifically designed patch of distributed micron-sized surface roughness. The transient-growth results for the HLB capsule are compared to corresponding numerical data for a Mach 6 blunt capsule experiment in the Adjustable Contour Expansion (ACE) facility of the Texas A&M University (TAMU) at lower Reynolds number. Similar trends are observed for both configurations. In particular, a rather low maximum energy gain is noted for the surface temperature conditions of the experiments. It is shown that the surface temperature dependence of the optimal transient-growth results is very similar for both capsule configurations. Moreover, the generation of stationary disturbances by well-defined roughness patches on the capsule surface is studied for the conditions of the HLB experiment using direct numerical simulations (DNS). To help explain the observed laminar-turbulent transition downstream of the roughness patch in some of the HLB capsule experiments, additional simulations were carried out to study the evolution of unsteady perturbations within the steady disturbance flow field due to the roughness patch. However, the DNS did not provide any indication of modal or non-modal disturbance growth in the wake of the roughness patch, and hence, the physical mechanism underlying the observed onset of transition remains unknown.

Hein, Stefan↗

NASA NTRS2016

Direct Numerical Simulation of Transition Due to Traveling Crossflow Vortices

Previous simulations of laminar breakdown mechanisms associated with stationary crossflow instability over a realistic swept-wing configuration are extended to investigate the alternate scenario of transition due to secondary instability of traveling crossflow modes. Earlier analyses based on secondary instability theory and parabolized stability equations have shown that this alternate scenario is viable when the initial amplitude of the most amplified mode of the traveling crossflow instability is greater than approximately 0.03 times the initial amplitude of the most amplified stationary mode. The linear growth predictions based on the secondary instability theory and parabolized stability equations agree well with the direct numerical simulation. Nonlinear effects are initially stabilizing but subsequently lead to a rapid growth followed by the onset of transition when the amplitude of the secondary disturbance exceeds a threshold value. Similar to the breakdown of stationary vortices, the transition zone is rather short and the boundary layer becomes completely turbulent across a distance of less than 15 times the boundary layer thickness at the completion of transition.

Li, Fei↗

NASA NTRS2016

Stationary Crossflow Breakdown due to Mixed Mode Spectra of Secondary Instabilities

Numerical simulations are used to study laminar breakdown characteristics associated with stationary crossflow instability in the boundary-layer flow over a subsonic swept-wing configuration. Previous work involving the linear and nonlinear development of individual, fundamental modes of secondary instability waves is extended by considering the role of more complex, yet controlled, spectra of the secondary instability modes. Direct numerical simulations target a mixed mode transition scenario involving the simultaneous presence of Y and Z modes of secondary instability. For the initial amplitudes investigated in this paper, the Y modes are found to play an insignificant role during the onset of transition, in spite of achieving rather large, O(5%), amplitudes of RMS velocity fluctuation prior to transition. Analysis of the numerical simulations shows that this rather surprising finding can be attributed to the fact that the Y modes are concentrated near the top of the crossflow vortex and exert relatively small influence on the Z modes that reside closer to the surface and can lead to transition via nonlinear spreading that does not involve interactions with the Y mode. Finally, secondary instability calculations reveal that subharmonic modes of secondary instability have substantially lower growth rates than those of the fundamental modes, and hence, are less likely to play an important role during the breakdown process involving complex initial spectra.

Li, Fei↗

NASA NTRS2015

Computations of Disturbance Amplification Behind Isolated Roughness Elements and Comparison with Measurements

Computations are performed to study laminar-turbulent transition due to isolated roughness elements in boundary layers at Mach 3.5 and 5.95, with an emphasis on flow configurations for which experimental measurements from low disturbance wind tunnels are available. The Mach 3.5 case corresponds to a roughness element with right-triangle planform with hypotenuse that is inclined at 45 degrees with respect to the oncoming stream, presenting an obstacle with spanwise asymmetry. The Mach 5.95 case corresponds to a circular roughness element along the nozzle wall of the Purdue BAMQT wind tunnel facility. In both cases, the mean flow distortion due to the roughness element is characterized by long-lived streamwise streaks in the roughness wake, which can support instability modes that did not exist in the absence of the roughness element. The linear amplification characteristics of the wake flow are examined towards the eventual goal of developing linear growth correlations for the onset of transition.

Choudhari, Meelan↗

NASA NTRS2015

Discrete-Roughness-Element-Enhanced Swept-Wing Natural Laminar Flow at High Reynolds Numbers

Nonlinear parabolized stability equations and secondary-instability analyses are used to provide a computational assessment of the potential use of the discrete-roughness-element technology for extending swept-wing natural laminar flow at chord Reynolds numbers relevant to transport aircraft. Computations performed for the boundary layer on a natural-laminar-flow airfoil with a leading-edge sweep angle of 34.6 deg, freestream Mach number of 0.75, and chord Reynolds numbers of 17 × 10(exp 6), 24 × 10(exp 6), and 30 × 10(exp 6) suggest that discrete roughness elements could delay laminar-turbulent transition by about 20% when transition is caused by stationary crossflow disturbances. Computations show that the introduction of small-wavelength stationary crossflow disturbances (i.e., discrete roughness element) also suppresses the growth of most amplified traveling crossflow disturbances.

Malik, Mujeeb↗

NASA NTRS2014

An Experimental Study of Roughness-Induced Instabilities in a Supersonic Boundary Layer

Progress on an experimental study of laminar-to-turbulent transition induced by an isolated roughness element in a supersonic laminar boundary layer is reported in this paper. Here, the primary focus is on the effects of roughness planform shape on the instability and transition characteristics. Four different roughness planform shapes were considered (a diamond, a circle, a right triangle, and a 45 degree fence) and the height and width of each one was held fixed so that a consistent frontal area was presented to the oncoming boundary layer. The nominal roughness Reynolds number was 462 and the ratio of the roughness height to the boundary layer thickness was 0.48. Detailed flow- field surveys in the wake of each geometry were performed via hot-wire anemometry. High- and low-speed streaks were observed in the wake of each roughness geometry, and the modified mean flow associated with these streak structures was found to support a single dominant convective instability mode. For the symmetric planform shapes - the diamond and circular planforms - the instability characteristics (mode shapes, growth rates, and frequencies) were found to be similar. For the asymmetric planform shapes - the right-triangle and 45 degree fence planforms - the mode shapes were asymmetrically distributed about the roughness-wake centerline. The instability growth rates for the asymmetric planforms were lower than those for the symmetric planforms and therefore, transition onset was delayed relative to the symmetric planforms.

Kegerise, Michael A.↗

NASA NTRS2014

Nonlinear Development and Secondary Instability of Traveling Crossflow Vortices

Transition research under NASA's Aeronautical Sciences Project seeks to develop a validated set of variable fidelity prediction tools with known strengths and limitations, so as to enable "sufficiently" accurate transition prediction and practical transition control for future vehicle concepts. This paper builds upon prior effort targeting the laminar breakdown mechanisms associated with stationary crossflow instability over a swept-wing configuration relevant to subsonic aircraft with laminar flow technology. Specifically, transition via secondary instability of traveling crossflow modes is investigated as an alternate scenario for transition. Results show that, for the parameter range investigated herein, secondary instability of traveling crossflow modes becomes insignificant in relation to the secondary instability of the stationary modes when the relative initial amplitudes of the traveling crossflow instability are lower than those of the stationary modes by approximately two orders of magnitudes or more. Linear growth predictions based on the secondary instability theory are found to agree well with those based on PSE and DNS, with the most significant discrepancies being limited to spatial regions of relatively weak secondary growth, i.e., regions where the primary disturbance amplitudes are smaller in comparison to its peak amplitude. Nonlinear effects on secondary instability evolution is also investigated and found to be initially stabilizing, prior to breakdown.

Li, Fei↗

NASA NTRS2013

Direct Numerical Simulation of Transition in a Swept-Wing Boundary Layer

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing derived from a two-dimensional, partial-differential-equation based eigenvalue computation; and the mode selected for forcing corresponds to the most amplified secondary instability mode which, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. Both the growth of the secondary instability wave and the resulting onset of laminar-turbulent transition are captured within the DNS computations. The growth of the secondary instability wave in the DNS solution compares well with linear secondary instability theory when the amplitude is small; the linear growth is followed by a region of reduced growth resulting from nonlinear effects before an explosive onset of laminar breakdown to turbulence. The peak fluctuations are concentrated near the boundary layer edge during the initial stage of transition, but rapidly propagates towards the surface during the process of laminar breakdown. Both time-averaged statistics and flow visualization based on the DNS reveal a sawtooth transition pattern that is analogous to previously documented surface flow visualizations of transition due to stationary crossflow instability. The memory of the stationary crossflow vortex is found to persist through the transition zone and well beyond the location of the maximum skin friction.

Duan, Lian↗

NASA NTRS2013

Towards Bridging the Gaps in Holistic Transition Prediction via Numerical Simulations

The economic and environmental benefits of laminar flow technology via reduced fuel burn of subsonic and supersonic aircraft cannot be realized without minimizing the uncertainty in drag prediction in general and transition prediction in particular. Transition research under NASA's Aeronautical Sciences Project seeks to develop a validated set of variable fidelity prediction tools with known strengths and limitations, so as to enable "sufficiently" accurate transition prediction and practical transition control for future vehicle concepts. This paper provides a summary of selected research activities targeting the current gaps in high-fidelity transition prediction, specifically those related to the receptivity and laminar breakdown phases of crossflow induced transition in a subsonic swept-wing boundary layer. The results of direct numerical simulations are used to obtain an enhanced understanding of the laminar breakdown region as well as to validate reduced order prediction methods.

Choudhari, Meelan M.↗

NASA NTRS2013

Wake Instabilities Behind Discrete Roughness Elements in High Speed Boundary Layers

Computations are performed to study the flow past an isolated, spanwise symmetric roughness element in zero pressure gradient boundary layers at Mach 3.5 and 5.9, with an emphasis on roughness heights of less than 55 percent of the local boundary layer thickness. The Mach 5.9 cases include flow conditions that are relevant to both ground facility experiments and high altitude flight ("cold wall" case). Regardless of the Mach number, the mean flow distortion due to the roughness element is characterized by long-lived streamwise streaks in the roughness wake, which can support instability modes that did not exist in the absence of the roughness element. The higher Mach number cases reveal a variety of instability mode shapes with velocity fluctuations concentrated in different localized regions of high base flow shear. The high shear regions vary from the top of a mushroom shaped structure characterizing the centerline streak to regions that are concentrated on the sides of the mushroom. Unlike the Mach 3.5 case with nearly same values of scaled roughness height k/delta and roughness height Reynolds number Re(sub kk), the odd wake modes in both Mach 5.9 cases are significantly more unstable than the even modes of instability. Additional computations for a Mach 3.5 boundary layer indicate that the presence of a roughness element can also enhance the amplification of first mode instabilities incident from upstream. Interactions between multiple roughness elements aligned along the flow direction are also explored.

Choudhari, Meelan↗

NASA NTRS2012

Computational Investigation of Supersonic Boundary Layer Transition Over Canonical Fuselage Nose Configurations

Boundary layer transition over axisymmetric bodies at non-zero angle of attack in supersonic flow is numerically investigated as part of joint research between the National Aeronautics and Space Administration (NASA) and Japan Aerospace Exploration Agency (JAXA). Transition over four axisymmetric bodies (namely, Sears-Haack body, semi-Sears-Haack body, 5-degree straight cone and flared cone) with different axial pressure gradients has been studied at Mach 2 in order to understand the effect of axial pressure gradient on instability amplification along the leeward symmetry plane and in the region of nonzero crossflow away from it. Comparisons are made with measured transition data in Mach 2 facilities as well as with predicted and measured transition characteristics for a 5-degree straight cone in a Mach 3.5 low disturbance tunnel. Limitations of using linear stability correlations for predicting transition over axisymmetric bodies at angle of attack are pointed out.

Choudhari, Meelan M.↗

NASA NTRS2012

Transition Within Leeward Plane of Axisymmetric Bodies at Incidence in Supersonic Flow

Boundary layer transition along the leeward symmetry plane of axisymmetric bodies at nonzero angle of attack in supersonic flow was investigated experimentally and numerically as part of joint research between the Japan Aerospace Exploration Agency (JAXA) and National Aeronautics and Space Administration (NASA). Transition over four axisymmetric bodies (namely, Sears-Haack body, semi-Sears-Haack body, straight cone and flared cone) with different axial pressure gradients was measured in two different facilities with different unit Reynolds numbers. The semi-Sears-Haack body and flared cone were designed at JAXA to broaden the range of axial pressure distributions. For a body shape with an adverse pressure gradient (i.e., flared cone), the experimentally measured transition patterns show an earlier transition location along the leeward symmetry plane in comparison with the neighboring azimuthal locations. For nearly zero pressure gradient (i.e.,straight cone), this feature is only observed at the larger unit Reynolds number. Later transition along the leeward plane was observed for the remaining two body shapes with a favorable pressure gradient. The observed transition patterns are only partially consistent with the numerical predictions based on linear stability analysis. Additional measurements are used in conjunction with the stability computations to explore the phenomenon of leeward line transition and the underlying transition mechanism in further detail.

Tokugawa, Naoko↗