Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient flow dynamics”

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 325 records · Page 18

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗

Polar vortex dynamics

Recent work with high resolution, one-layer numerical models of fluid flows resembling those in the real stratosphere has suggested that: (1) the interiors of strong cyclonic vortices like the Antarctic polar vortex may be almost completely isolated laterally from their surroundings - perhaps even completely isolated, under some circumstances; (2) by contrast, material near the edge of such and isolated region can easily be eroded (or mixed one-sidedly) into the surrounding region; and (3) the erosion characteristically produces extremely steep gradients in isentropic distributions of potential vorticity (PV) and of other tracers, possibly down to horizontal length scales of a few kilometers only. Such length scales may occur both at the edge of the main polar vortex and in smaller features outside it, such as thin filamentary structures, produced by the erosion process.

Mcintyre, Michael↗

Wall-Modeled Large-Eddy Simulations for High-Lift Configurations using FUN3D

Wall-modeled large-eddy simulation (WMLES) capability has recently been implemented into FUN3D, an unstructured, node-centered, finite-volume solver developed at the NASA Langley Research Center. In this paper, WMLES is assessed for two configurations that are representative for high-lift applications. The first configuration is a nominal two-dimensional multielement airfoil that has been extensively studied in the literature. WMLES solutions are computed for four angles of attack and compared with previously reported solutions. Good agreement of integrated forces, surface pressures, and boundary-layer velocity profiles is shown with available experimental data especially at lower angles of attack. WMLES solutions are also computed for the NASA High-Lift Common Research Model over a large range of angles of attack. Forces, pitching moments, and pressure distributions are favorably compared with the experimental data up to the maximum lift, including the angle of attack where the maximum lift is obtained experimentally. Eddy visualization techniques of q-criterion and density-gradient magnitude illustrate the resolved content.

High-Lift flows↗

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗

Dynamic models of flux tubes in the interpretation of polarization measurements

Recent observations of Stokes parameter profiles indicate the presence of mass motions with large velocity gradients associated with small-scale magnetic elements. Dynamic models of flux tubes were used in order to interpret observations of unresolved elements. It is clear that the physical picture of the dynamic models will be quite different from the hydrostatic ones since there is a strong coupling between the magnetic and the velocity field. Polarization measurements have to be interpreted in terms of dynamic models. Two-D steady flow solutions in slender magnetic tubes have been worked out. It was found that the main properties of the intensity line profiles as well as the asymmetries of the V Stokes profiles can be explained best in terms of magnetic elements with moderate field strength.

Ribes-Nesme, E.↗

Experimental and numerical study of the British Experimental Rotor Programme blade

Wind-tunnel tests on the British Experimental Rotor Programme (BERP) tip are described, and the results are compared with computational fluid dynamics (CFD) results. The test model was molded using the Lynx-BERP blade tooling to provide a semispan, cantilever wing comprising the outboard 30 percent of the rotor blade. The tests included both surface-pressure measurements and flow visualization to obtain detailed information of the flow over the BERP tip for a range of angles of attack. It was observed that, outboard of the notch, favorable pressure gradients exist which ensure attached flow, and that the tip vortex also remains stable to large angles of attack. On the rotor, these features yield a very gradual break in control loads when the retreating-blade limit is eventually reached. Computational and experimental results were generally found to be in good agreement.

Brocklehurst, Alan↗

Sulfur gas exchange in Sphagnum-dominated wetlands

Sulfur gases are important components of the global cycle of S. They contribute to the acidity of precipitation and they influence global radiation balance and climate. The role of terrestrial sources of biogenic S and their effect on atmospheric chemistry remain as major unanswered questions in our understanding of the natural S cycle. The role of northern wetlands as sources and sinks of gaseous S was investigated by measuring rates of S gas exchange as a function of season, hydrologic conditions, and gradients in trophic status. The effects of inorganic S input on the production and emission of gaseous S were also investigated. Experiments were conducted in wetlands in New Hampshire, particularly a poor fen, fens within the Experimental Lakes Area (ELA) in Ontario, Canada and in freshwater and marine tundra. Emissions were determined using Teflon enclosures, gas cryotrapping methods, and gas chromatography (GC) with flame photometric detection. Dynamic (sweep flow) and static enclosures were employed. Dissolved gases were determined by gas stripping followed by GC.

Hines, Mark E.↗

Automated Fluid Feature Extraction from Transient Simulations

In the past, feature extraction and identification were interesting concepts, but not required to understand the underlying physics of a steady flow field. This is because the results of the more traditional tools like iso-surfaces, cuts and streamlines were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of much interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snap-shot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments like pV3). And methods must be developed to abstract the feature and display it in a manner that physically makes sense. The following is a list of the important physical phenomena found in transient (and steady-state) fluid flow: (1) Shocks, (2) Vortex cores, (3) Regions of recirculation, (4) Boundary layers, (5) Wakes. Three papers and an initial specification for the (The Fluid eXtraction tool kit) FX Programmer's guide were included. The papers, submitted to the AIAA Computational Fluid Dynamics Conference, are entitled : (1) Using Residence Time for the Extraction of Recirculation Regions, (2) Shock Detection from Computational Fluid Dynamics results and (3) On the Velocity Gradient Tensor and Fluid Feature Extraction.

Haimes, Robert↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Absolute wind measurements in the lower thermosphere of Venus using infrared heterodyne spectroscopy

The first absolute wind velocities above the Venusian cloud-tops were obtained using NASA/Goddard infrared heterodyne spectrometers at the NASA Infrared Telescope Facility (IRTF) and the McMath Solar Telescope. Beam-integrated Doppler displacements in the non-thermal emission core of (12)C(16)O2 10.33 micron R(8) sampled the line of sight projection of the lower thermospheric wind field (100 to 120 km). A field-usable Lamb-dip laser stabilization system, developed for spectrometer absolute frequency calibration to less than + or - 0.1 MHz, allowed S/N-limited line of sight velocity resolution at the 1 m/s level. The spectrometer's diffraction-limited beam (1.7 arc-second HPBW at McMath, 0.9 arc-second HPBW at IRTF), and 1 to 2 arc-second seeing, provided the spatial resolution necessary for circulation model discrimination. Qualitative analysis of beam-integrated winds provided definitive evidence of a dominant subsolar-antisolar circulation in the lower thermosphere. Beam-integrated winds were modelled with a 100x100 grid over the beam, incorporating beam spatial rolloff and across-the-beam gradients in non-thermal emission intensity, line of sight projection geometry, and horizontal wind velocity. Horizontal wind velocity was derived from a 2-parameter model wind field comprised of subsolar-antisolar and zonal components. Best-fit models indicated a dominant subsolar-antisolar flow with 120 m/s cross-terminator winds and a retrograde zonal component with a 25 m/s equatorial velocity. A review of all dynamical indicators above the cloud-tops allowed development of an integrated and self-consistent picture of circulation in the 70 to 200 km range.

Goldstein, Jeffrey J.↗

Review of the PDWA Concept for Combustion Enhancement in a Supersonic Air-Breathing Combustor Environment

This paper reviews the design of the Pulsed Detonation Wave Augmentor (PDWA) concept and the preliminary computational fluid dynamics studies that supported it. The PDWA relies on the rapid generation of detonation waves in a small tube, which are then injected into the supersonic stream of the main combustor. The blast waves thus generated are used to stimulate the mixing and combustion inside the main combustor. The mixing enhancement relies on various forms of the baroclinic interaction, where misaligned pressure and density gradients combine to produce vortical flow. By using unsteady shock waves, the concept also uses the Richtmyer-Meshkov effect to further increase the rate of mixing. By carefully designing the respective configurations of the combustor and the detonation tubes, one can also increase the penetration of the fuel into the supersonic air stream. The unsteady shocks produce lower stagnation pressure losses than steady shocks. Combustion enhancement can also be obtained through the transient shock-heating of the fuel-air interface, and the lowering of the ignition delay in these regions. The numerical simulations identify these processes, and show which configurations give the best results. Engineering considerations are also presented, and discuss the feasibility of the concept. Of primary importance are the enhancements in performance, the design simplicity, the minimization of the power, cost, and weight, and the methods to achieve very rapid cycling.

Canbier, Jean-Luc↗

A New View of the Dynamics of Reynolds Stress Generation in Turbulent Boundary Layers

The structure of a numerically simulated turbulent boundary layer over a flat plate at Re(theta) = 670 was studied using the invariants of the velocity gradient tensor (Q and R) and a related scalar quantity, the cubic discriminant (D = 27R(exp 2)/4 + Q(exp 3)). These invariants have previously been used to study the properties of the small-scale motions responsible for the dissipation of turbulent kinetic energy. In addition, these scalar quantities allow the local flow patterns to be unambiguously classified according to the terminology proposed by Chong et al. The use of the discriminant as a marker of coherent motions reveals complex, large-scale flow structures that are shown to be associated with the generation of Reynolds shear stress -u'v'(bar). These motions are characterized by high spatial gradients of the discriminant and are believed to be an important part of the mechanism that sustains turbulence in the near-wall region.

Cantwell, Brian J.↗

A numerical investigation of the finite element method in compressible primitive variable Navier-Stokes flow

The results of a comprehensive numerical investigation of the basic capabilities of the finite element method (FEM) for numerical solution of compressible flow problems governed by the two-dimensional and axis-symmetric Navier-Stokes equations in primitive variables are presented. The strong and weak points of the method as a tool for computational fluid dynamics are considered. The relation of the linear element finite element method to finite difference methods (FDM) is explored. The calculation of free shear layer and separated flows over aircraft boattail afterbodies with plume simulators indicate the strongest assets of the method are its capabilities for reliable and accurate calculation employing variable grids which readily approximate complex geometry and capably adapt to the presence of diverse regions of large solution gradients without the necessity of domain transformation.

Cook, C. H.↗

A shallow water model of the winter stratosphere

The shallow water equations are used to investigate the interaction of planetary wave-breaking with a 'diabatic' forcing. The numerical integrations demonstrate the formation of a sharp gradient in potential vorticity at the edge of the polar vortex, despite 'diabatic relaxation' toward a smooth gradient. Successive cycles of Rossby wave-breaking accentuate the distinction between the vortex and the surf-zone. Persistent easterly acceleration in the tropics eventually produces a flow regime of little relevance to the atmosphere. For the period when the flow has some qualitative resemblance to the middle stratosphere diffusion coefficients are calculated from the meridional fluxes of potential vorticity and of passive tracers. It is found that the unique dynamical nature of potential vorticity among tracers is not reflected in its diffusion coefficient.

Juckes, Martin↗

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗

SSME thrust chamber simulation using Navier-Stokes equations

The capability of the PHOENICS fluid dynamics code in predicting two-dimensional, compressible, and reacting flow in the combustion chamber and nozzle of the space shuttle main engine (SSME) was evaluated. A non-orthogonal body fitted coordinate system was used to represent the nozzle geometry. The Navier-Stokes equations were solved for the entire nozzle with a turbulence model. The wall boundary conditions were calculated based on the wall functions which account for pressure gradients. Results of the demonstration test case reveal all expected features of the transonic nozzle flows. Of particular interest are the locations of normal and barrel shocks, and regions of highest temperature gradients. Calculated performance (global) parameters such as thrust chamber flow rate, thrust, and specific impulse are also in good agreement with available data.

Przekwas, A. J.↗

Modeling the topographic influence on aboveground biomass using a coupled model of hillslope hydrology and ecosystem dynamics

Abstract. Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling.

54 ENVIRONMENTAL SCIENCES↗

Canonical solutions for unsteady flow fields

The initial value problem of one-dimensional gas-dynamics involving discontinuous, nonuniform initial data is discussed. Canonical solutions which are valid in a small x, t region aroung a discontinuity, and which include the first order effects of nonuniformities in the data, are derived explicitly. The theory is derived by considering a group of elementary piston problems. Solutions with a shock or with a centered expansion wave are worked out individually in order to relate initial flow properties and their gradients to the speed and acceleration of the discontinuity waves. They are then combined to represent the solution of a general initial value problem by regarding the piston path as a contact line. In addition, problems with chemical reaction are discussed in terms of elementary piston problems which involve strong detonation waves, Chapman-Jouguet detonation waves, and deflagration waves.

Liu, G. C.↗