Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Critical Fluids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,459 records · Page 81

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

Effect of Fluid Flow on Zinc Electrodeposits from Acid Chloride Electrolytes

Zinc was deposited potentiostatically from acid chloride baths. Once bath chemistry and electrochemistry were controlled, the study was focused on convective mass transfer at horizontal electrodes and its effect on cell performance. A laser schlieren imaging technique allowed in situ observations of flow patterns and their correlation with current transients. Convection was turbulent and mass transfer as a function of Rayleigh number was well correlated by: Sh = 0.14 R to the 1/3 power. Similarly, convection initiation time was correlated by DT/d squared = 38 Ra to the -2/3 power. Time scale of fluctuations was about half the initiation time. Taking the boundary layer thickness as a characteristic length, a critical Rayleigh number for the onset of convection was deduced: Ra sub CR = 5000. Placing the anode on the top of the cathode completely changed the flow pattern but kept the I-t curves identical whereas the use of a cathode grid doubled the limiting current. A well defined plateau in the current voltage curves suggested that hydrogen evolution has been successfully inhibited. Finally, long time deposition showed that convection at horizontal electrodes increased the induction time for dentrite growth by at least a factor of 2 with respect to a vertical wire.

Abdelmassir, A. A.↗

The Local Electronic Structure of Supercritical CO 2 from X-ray Raman Spectroscopy and Atomistic-Scale Modeling

Supercritical CO 2 is encountered in several technical and natural systems related to biology, geophysics, and engineering. While the structure of gaseous CO 2 has been studied extensively, the properties of supercritical CO 2 , particularly close to the critical point, are not well-known. Here, in this work, we combine X-ray Raman spectroscopy, molecular dynamics simulations, and first-principles density functional theory (DFT) calculations to characterize the local electronic structure of supercritical CO 2 at conditions around the critical point. The X-ray Raman oxygen K-edge spectra manifest systematic trends associated with the phase change of CO 2 and the intermolecular distance. Extensive first-principles DFT calculations rationalize these observations on the basis of the 4sσ Rydberg state hybridization. X-ray Raman spectroscopy is found to be a sensitive tool for characterizing electronic properties of CO 2 under challenging experimental conditions and is demonstrated to be a unique probe for studying the electronic structure of supercritical fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗

Numerical simulation of fluid flow around a scramaccelerator projectile

Numerical simulations of the fluid motion and temperature distribution around a 'scramaccelerator' projectile are obtained for Mach numbers in the 5-10 range. A finite element method is used to solve the equations of motion for inviscid and viscous two-dimensional or axisymmetric compressible flow. The time-dependent equations are solved explicitly, using bilinear isoparametric quadrilateral elements, mass lumping, and a shock-capturing Petrov-Galerkin formulation. Computed results indicate that maintaining on-design performance for controlling and stabilizing oblique detonation waves is critically dependent on projectile shape and Mach number.

Pepper, Darrell W.↗

An experimental study of the fluid mechanics associated with porous walls

The fluid mechanics associated with the blowing phenomenon from porous walls is measured and characterized. The measurements indicate that the flow exiting a porous wall exhibits a lumpy velocity profile caused by the coalescence effects of smaller jets emerging from the surface. The velocity variations are spatially stable and prevail even at low flow rates. The intensity of this pseudoturbulence is found to be directly proportional to the filter rating of the porous wall and to increase linearly with the mean velocity. Beyond a critical mean velocity, the pseudoturbulence intensity shows a leveling trend with increase in the mean velocity. This critical velocity varies inversely as the filter rating and represents the onset of fully developed jetting action in the flow field. Based on the data, a more appropriate length scale for the flow field is proposed and a correlation is developed that can be used to predict the onset of fully developed jets in the flow emerging from a porous wall.

Ramachandran, N.↗

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗

Wetting layers and dispersion forces for a fluid in contact with a vertical wall

Ellipsometric measurements of the thickness of the wetting layer of sulfur hexafluoride that forms on a fused-silica surface partially immersed in the liquid are presented. The components of the experiment included a He-Ne laser (0.6328 micron), polarizer, and quarter-wave plate with 45-deg azimuth. An optical-grade fused-silica prism was polished to wavelength/20, and strain-induced birefringence was accounted for. It was found that the layer thickness decreases from 40 to 20 nm with a height increase from 0.7 to 3 mm above the liquid-vapor interface. The thickness is independent of the temperature between 10 and 0.1 K below the critical temperature (319 K). The results are compared in detail to the theory of Dzyaloshinskii et al. (1961), and are noted to coincide closely.

Kayser, R. F.↗

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗

Critical mass flux through short Borda type inlets of various cross sections

Mass flux measurements associated with choked flows through four Borda-type inlet geometries: circular, square, triangular, and rectangular (two-dimensional) and two sharp-edged geometries are discussed for a wide range of inlet stagnation conditions. The results obtained indicate that the mass flux is independent of the inlet cross-section geometry, while it is dependent on the inlet stagnation conditions. The results also suggest that parallel surfaces as found in seals, dampers, bearings, and heat exchanger tubes of rectangular cross section with Borda-type inlet configurations are subject to large forces. It is noted that the reduced mass flux is independent of working fluid.

Hendricks, R. C.↗

Static Feed Water Electrolysis Subsystem Testing and Component Development

A program was carried out to develop and test advanced electrochemical cells/modules and critical electromechanical components for a static feed (alkaline electrolyte) water electrolysis oxygen generation subsystem. The accomplishments were refurbishment of a previously developed subsystem and successful demonstration for a total of 2980 hours of normal operation; achievement of sustained one-person level oxygen generation performance with state-of-the-art cell voltages averaging 1.61 V at 191 ASF for an operating temperature of 128F (equivalent to 1.51V when normalized to 180F); endurance testing and demonstration of reliable performance of the three-fluid pressure controller for 8650 hours; design and development of a fluid control assembly for this subsystem and demonstration of its performance; development and demonstration at the single cell and module levels of a unitized core composite cell that provides expanded differential pressure tolerance capability; fabrication and evaluation of a feed water electrolyte elimination five-cell module; and successful demonstration of an electrolysis module pressurization technique that can be used in place of nitrogen gas during the standby mode of operation to maintain system pressure and differential pressures.

Koszenski, E. P.↗

Development of Methods for In-line Monitoring of Yield Stress During the Transfer of Radioactive Waste - 20435

Radioactive waste slurries at the Savannah River Site (SRS) and Hanford often behave as non-Newtonian fluids with a yield stress. The yield stress increases the size of equipment needed to transport or mix these slurries and can increase the risk of process upsets. This property is often measured by collection of samples and transportation to a laboratory, which can increase the risk of personnel exposure and even change the rheology. In addition, slurry rheology has been observed to change with time. Real-time in-line monitoring (RTIM) of rheological properties of slurries can have substantial benefits to critical processes which are currently in operation or will be part of future plans in Department of Energy (DOE) complex processes. The need for an appropriate RTIM technology that can potentially minimize or eliminate numerous sampling operations was recently identified. This research at Florida International University focuses on addressing the need for a reliable RTIM technology at US Department of Energy sites. The objective of this work is to investigate and develop appropriate and reliable technologies that can perform accurate and easy measurements of yield stress with minimal changes to the flow conditions at US Department of Energy sites. Use of ultrasound in longitudinal and shear modes was considered as a method for monitoring of yield stress variations with no disturbance to the flow. Measurement of shear wave velocity (in shear mode) and monitoring of sound velocity variations with start of flow from rest (longitudinal mode) are the underlying principles. Measurement of pressure loss and liquid rise (due to static differential pressure) during flow of slurries were two other methods that were considered and investigated. In addition, analysis of stability depth of a penetrometer in slurry medium (penetrometer method) revealed a potential suitability. This paper presents efforts associated with investigation of different methods mentioned earlier for kaolin-water mixtures at different concentrations. Theoretical analysis showed that all methods under investigation were capable of monitoring yield stress in a desirable range of 0 to 70 pascal. A bench scale test setup was constructed for investigation of pressure loss and liquid rise methods. Simulants at different concentrations were created for initial testing campaigns. For ultrasonic testing, effects of yield stress on sound propagation speed in longitudinal and shear modes were investigated for simulants placed inside PVC pipes. Parameters such as simulant concentration and pipe diameter were varied to find optimal values. It was found that spacing (or medium thickness) of 1/8 inch or less was necessary to observe echoes with sufficient strength. In addition, a significant effect of yield stress on sound propagation velocity was observed. For the liquid rise method, kaolin-water simulant at 1.17 specific gravity was pumped through a 1/4 inch pipe and liquid height was measured in a vertical riser branching off a three-way connection. Static gauge pressure was varied within [0.3 to 1.1] psi range by increase of pumping pressure and flow. Results showed consistent yield stress data within 10% deviation from values reported by a rheometer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Optimization of a Welding Procedure for Making Critical Aluminum Welds on the LBNF Absorber Core Block

The LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. A weld quality optimization was performed to produce National Aeronautical Standard 1514 Class I quality welds on the aluminum core blocks. This was not successful in all cases. An existing Gas Tungsten Arc Welding Procedure Specification was fine tuned to minimize, in most cases, and eliminate detect-able tungsten inclusions in the welds. All the weld coupons, however passed welding inspection as per the piping code: ASME B31.3 Normal Fluid Service. Tungsten electrode diameter, type, and manufacturer were varied. Some of the samples were pre-heated and others were not. It was observed that larger diameter electrodes, 5/32 in., with pre-heated joints resulted in welds with the least number of tungsten inclusions. It is hypothesized that thinner electrodes breakdown easily and get lodged into the weld pool during the welding process. This breakdown is further enhanced by the large temperature differential between the un-preheated sample and the hot electrode.

43 PARTICLE ACCELERATORS↗

Evaporation Sub-model Development for Volume of Fluid (eVOF) Method Applicable to Spray-Wall Interaction Including Film Characteristics with Validation at High Pressure and Temperature Conditions (Final Report)

Internal combustion engines have seen a great evolution over the last several decades through application of high pressure direct injection, multiple injections, and other technologies to reduced fuel consumption, NOx, and PM. Although combustion systems with advanced injection strategies have been studied extensively, there exists a significant fundamental knowledge gap on the fuel-spray interactions with the piston surface and chamber walls. Advanced computational codes validated with experimental techniques have to be developed for accurate representation of the drop impingement, fuel film formation, and vaporization. Current engine CFD spray models utilize a Lagrangian framework for modeling which lacks critical considerations of the physics pertaining to these interactions and thus requiring extensive parameterization, tuning and validation. The team from Michigan Technological University, University of Massachusetts Dartmouth, and Argonne National Laboratory is composed of experts in sprays, combustion, engines and CFD with a wide spectrum of knowledge including specific expertise in the area under consideration. In the proposed work, a VOF modeling approach has been adopted for the spray-wall interaction, film formation and spreading, and vaporization. With the inclusion of a vaporization submodel, a more predictive and accurate simulation of the spray-film was performed without extensive need of parameterization and tuning. Extensive experimentation of the spray-wall interaction under the range of conditions matching the thermodynamic and surface temperatures that occur in diesel and gasoline engines were conducted to validate the SWI submodels and for development of the evaporation sub-model, which has been implemented in the Converge software.

42 ENGINEERING↗

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Differential transcriptomic alterations in nasal versus lung tissue of acrolein-exposed rats

Introduction: Acrolein is a significant component of anthropogenic and wildfire emissions, as well as cigarette smoke. Although acrolein primarily deposits in the upper respiratory tract upon inhalation, patterns of site-specific injury in nasal versus pulmonary tissues are not well characterized. This assessment is critical in the design of in vitro and in vivo studies performed for assessing health risk of irritant air pollutants. Methods: In this study, male and female Wistar-Kyoto rats were exposed nose-only to air or acrolein. Rats in the acrolein exposure group were exposed to incremental concentrations of acrolein (0, 0.1, 0.316, 1 ppm) for the first 30 min, followed by a 3.5 h exposure at 3.16 ppm. In the first cohort of male and female rats, nasal and bronchoalveolar lavage fluids were analyzed for markers of inflammation, and in a second cohort of males, nasal airway and left lung tissues were used for mRNA sequencing. Results: Protein leakage in nasal airways of acrolein-exposed rats was similar in both sexes; however, inflammatory cells and cytokine increases were more pronounced in males when compared to females. No consistent changes were noted in bronchoalveolar lavage fluid of males or females except for increases in total cells and IL-6. Acrolein-exposed male rats had 452 differentially expressed genes (DEGs) in nasal tissue versus only 95 in the lung. Pathway analysis of DEGs in the nose indicated acute phase response signaling, Nrf2-mediated oxidative stress, unfolded protein response, and other inflammatory pathways, whereas in the lung, xenobiotic metabolism pathways were changed. Genes associated with glucocorticoid and GPCR signaling were also changed in the nose but not in the lung. Discussion: These data provide insights into inhaled acrolein-mediated sex-specific injury/inflammation in the nasal and pulmonary airways. The transcriptional response in the nose reflects acrolein-induced acute oxidative and cytokine signaling changes, which might have implications for upper airway inflammatory disease susceptibility.

Alewel, Devin I.↗

Oscillatory/Chaotic Thermocapillary Flow Induced by Radiant Heating

There is a continuing need to understand the fluid physics occurring under low gravity conditions in processes such as crystal growth, materials processing, and the movement of bubbles or droplets. The fluid flow in such situations is often caused by a gradient in interfacial tension. If a temperature gradient is created due to a heat source, the resulting flow is called thermocapillary flow, a special case of Marangoni Convection. In this study, an experimental investigation was conducted using silicone oil in cylindrical containers with a laser heat source at the free surface. It was desired to determine the conditions under which steady, axisymmetrical thermocapillary flow becomes unstable and oscillatory three-dimensional flow states develop. The critical Marangoni number for each observed oscillatory state was measured as a function of the container aspect ratio and the dynamic Bond number, a measure of buoyant force versus ii thermocapillary force. Various oscillatory modes were observed during three- dimensional convection, and chaotic flow was reached in one test condition. The critical Marangoni numbers are compared with those measured in previous studies, and the power spectra and phase trajectories of the instantaneous surface temperature distributions are used to characterize the routes of transitions to the chaotic flow state. Results show that only superharmonic modes appear in the routes to chaos while infinite number of subharmonic modes occur in flow transitions for pure Rayleigh convection.

DeWitt, Kenneth J.↗

Mesoscale Variability in SUCCESS Data

Analysis of meteorological, chemical, and microphysical data from the airborne SUCCESS mission is reported. Careful analysis of the complex DC-8 flight pattern of May 2, 1996 reveals 19 linear or nearly linear flight segments within six main geographical areas, which we have analyzed. Significant mountain wave activity is revealed in the data from the MMS and MTP instruments on the DC-8, which resembles previous observations of mountain wave structures near Boulder, CO. Strong mountain-wave-induced upwelling downwind of the Rockies is noted. Turbulence is also noted in regions of the mountain wave consistent with overturning near the tropopause. Zonal winds recorded on the ER-2 are shown to consistent with mountain wave breaking at or near critical levels in the stratosphere, consistent with the strong turbulence reported by the pilot during the ER-2 flight. Those observations have been supported with spectral analyses and modeling studies. "Postcasts" of mountain wave activity on May 2, 1996, using the Naval Research Laboratory Mountain Wave Forecast Model (NRL/MWFM) predicts both strong mountain wave activity near the tropopause (as measured by the DC-8) and strong mountain-wave-induced turbulence in the stratosphere (as encountered by the ER-2). Two-dimensional simulations of fluid flow over topography reveal similar isentropic structures to observations.

Eckermann, Stephen D.↗