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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Divide and conquer: Learning chaotic dynamical systems with multistep penalty neural ordinary differential equations

Forecasting high-dimensional dynamical systems is a fundamental challenge in various fields, such as geosciences and engineering. Neural Ordinary Differential Equations (NODEs), which combine the power of neural networks and numerical solvers, have emerged as a promising algorithm for forecasting complex nonlinear dynamical systems. However, classical techniques used for NODE training are ineffective for learning chaotic dynamical systems. In this work, we propose a novel NODE-training approach that allows for robust learning of chaotic dynamical systems. Here, our method addresses the challenges of non-convexity and exploding gradients associated with underlying chaotic dynamics. Training data trajectories from such systems are split into multiple, non-overlapping time windows. In addition to the deviation from the training data, the optimization loss term further penalizes the discontinuities of the predicted trajectory between the time windows. The window size is selected based on the fastest Lyapunov time scale of the system. Multi-step penalty(MP) method is first demonstrated on Lorenz equation, to illustrate how it improves the loss landscape and thereby accelerates the optimization convergence. MP method can optimize chaotic systems in a manner similar to least-squares shadowing with significantly lower computational costs. Our proposed algorithm, denoted the Multistep Penalty NODE, is applied to chaotic systems such as the Kuramoto-Sivashinsky equation, the two-dimensional Kolmogorov flow, and ERA5 reanalysis data for the atmosphere. It is observed that MP-NODE provide viable performance for such chaotic systems, not only for short-term trajectory predictions but also for invariant statistics that are hallmarks of the chaotic nature of these dynamics.

Chaotic dynamical systems↗

DOE OSTI2024

Extreme elevational migration spurred cryptic speciation in giant hummingbirds

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how “outlier” species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world’s largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

Science & Technology - Other Topics↗

DOE OSTI2022

FLOW Estimation and Rose Superposition (FLOWERS): an integral approach to engineering wake models

Abstract. Annual energy production (AEP) is often the objective function in wind plant layout optimization studies. The conventional method to compute AEP for a wind farm is to first evaluate power production for each discrete wind direction and speed using either computational fluid dynamics simulations or engineering wake models. The AEP is then calculated by weighted-averaging (based on the wind rose at the wind farm site) the power produced across all wind directions and speeds. We propose a novel formulation for time-averaged wake velocity that incorporates an analytical integral of a wake deficit model across every wind direction. This approach computes the average flow field more efficiently, and layout optimization is an obvious application to exploit this benefit. The clear advantage of this new approach is that the layout optimization produces solutions with comparable AEP performance yet is completed 2 orders of magnitude faster. The analytical integral and the use of a Fourier expansion to express the wind speed and wind direction frequency create a relatively smooth solution space for the gradient-based optimizer to excel in comparison to the existing weighted-averaging power calculation.

17 WIND ENERGY↗

DOE OSTI2025

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗

DOE OSTI2024

Inter-ELM pedestal turbulence dynamics dependence on q 95 and temperature gradient

A series of dedicated experiments from DIII-D tokamak provide spatially and temporally resolved measurements of electron density and temperature, and multiscale and multichannel fluctuations over a wide range of conditions. Measurements of long wavelength density fluctuations in the type-I ELMing H-mode pedestal routinely revealed a coexistence of multiple instabilities that exhibit dramatic different dynamic behaviors as q 95 and temperature gradient are varied, apparently responsible for limiting pedestal temperature profiles. Two distinct frequency bands of density fluctuations are modulated with ELM cycle with frequency above 200kHz propagating in the electron diamagnetic direction in the lab frame (electron mode) and below 200kHz propagating in the ion diamagnetic direction (ion mode). The electron mode amplitude peaks near the electron temperature gradient region and increases with q 95 which seems to be correlated with the increased χ e at higher q 95 , similar to the characteristics expected for Micro-tearing Mode (MTM). At higher q 95 , during the inter-ELM period, the ion mode decays at later phase of the ELM cycle. Consistently, the poloidal correlation length of the ion mode is also found to reduce which suggests the possible E×B flow shear suppression of ion mode at later phase of the ELM cycle as the Er well recovers. In contrast, the electron mode grows during the ELM cycle and reaches saturation at around 50-60% of ELM period. Linear gyrokinetic simulations find the MTMs as the most unstable mode in the pedestal electron temperature gradient region. The higher q 95 and lower magnetic shear destabilize MTMs. These observations provide key insights of the underlying physics of multifield properties and rich dataset of experimental 'fingerprints' that enable new tests of the theoretical pedestal models and lead to developing a predictive model for pedestal formation on ITER and future burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗

NASA NTRS1992

Species and temperature measurement in H2/O2 rocket flow fields by means of Raman scattering diagnostics

Validation of CFD codes developed for prediction and evaluation of rocket performance is hampered by a lack of experimental data. Nonintrusive laser based diagnostics are needed to provide spatially and temporally resolved gas dynamic and fluid dynamic measurements. This paper reports the first nonintrusive temperature and species measurements in the plume of a 110 N gaseous hydrogen/oxygen thruster at and below ambient pressures, obtained with spontaneous Raman spectroscopy. Measurements at 10 mm downstream of the exit plane are compared with predictions from a numerical solution of the axisymmetric Navier-Stokes and species transport equations with chemical kinetics, which fully model the combustor-nozzle-plume flowfield. The experimentally determined oxygen number density at the centerline at 10 mm downstream of the exit plane is four times that predicted by the model. The experimental number density data fall between those numerically predicted for the exit and 10 mm downstream planes in both magnitude and radial gradient. The predicted temperature levels are within 10 to 15 percent of measured values.

De Groot, Wim A.↗

NASA NTRS1991

Interaction of two glancing, crossing shock waves with a turbulent boundary-layer at various Mach numbers

A preliminary experimental investigation was conducted to study two crossing, glancing shock waves of equal strengths, interacting with the boundary-layer developed on a supersonic wind tunnel wall. This study was performed at several Mach numbers between 2.5 and 4.0. The shock waves were created by fins (shock generators), spanning the tunnel test section, that were set at angles varying from 4 to 12 degrees. The data acquired are wall static pressure measurements, and qualitative information in the form of oil flow and schlieren visualizations. The principle aim is two-fold. First, a fundamental understanding of the physics underlying this flow phenomena is desired. Also, a comprehensive data set is needed for computational fluid dynamic code validation. Results indicate that for small shock generator angles, the boundary-layer remains attached throughout the flow field. However, with increasing shock strengths (increasing generator angles), boundary layer separation does occur and becomes progressively more severe as the generator angles are increased further. The location of the separation, which starts well downstream of the shock crossing point, moves upstream as shock strengths are increased. At the highest generator angles, the separation appears to begin coincident with the generator leading edges and engulfs most of the area between the generators. This phenomena occurs very near the 'unstart' limit for the generators. The wall pressures at the lower generator angles are nominally consistent with the flow geometries (i.e. shock patterns) although significantly affected by the boundary-layer upstream influence. As separation occurs, the wall pressures exhibit a gradient that is mainly axial in direction in the vicinity of the separation. At the limiting conditions the wall pressure gradients are primarily in the axial direction throughout.

Hingst, Warren R.↗

DOE OSTI2023

A Gradient Based Subgrid-Scale Parameterization for Ocean Mesoscale Eddies

Mesoscale eddies play an important role in transport of heat and biogeochemical tracers in the global ocean circulation. Resolving these energetic eddies, however, is challenging in ocean general circulation models (OGCM) because it requires a horizontal grid spacing of ≲1/8° that is computationally expensive. As a result, we are required to parameterize mesoscale eddy effects on large-scale ocean flows. In this work, we introduce a new subgrid-scale (SGS) model that is developed based on a Taylor series expansion of resolved variables to parameterize subgrid mesoscale eddy transports and momentum fluxes in OGCM. We have performed an a priori study to evaluate the performance of our new gradient model using high-resolution ocean simulations. Our results show that the gradient model well predicts the actual SGS thickness fluxes in the zonal and meridional directions in coarse-resolution simulations with the grid spacing ≳1/4°. The unresolved kinetic energy at the ocean surface is also skillfully estimated. More importantly, unlike current mesoscale eddy parameterizations, which are mainly developed based on an assumption of flat bottom ocean, our new SGS model can capture the structure of unresolved standing meanders at the ocean surface. We have also developed a dynamic procedure for setting in non-dimensional parameters in our new parameterization through a non-ad hoc and tuning-free method. Overall, this work suggests that implementing the gradient model in OGCM can improve the model accuracy with an affordable computational cost in eddy-permitting and non-eddying simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Parameter space for magnetization effects in high-energy-density plasmas

Magnetic fields are well known to affect the evolution of fluids via the J × B force, where J is the current density and B is the magnetic field. This force leads to the influence of magnetic fields on hydrodynamics (magnetohydrodynamics). Magnetic fields are often neglected in modeling of high-energy-density plasmas, since J × B is very small compared with the plasma pressure gradients. However, many experiments lie in a separate part of parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate magnetic fields. By comparing terms in the induction equation, we go on to estimate the regions of parameter space where these self-generated fields are strong enough to affect the hydrodynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1988

Thermal analysis of heat storage canisters for a solar dynamic, space power system

A thermal analysis was performed of a thermal energy storage canister of a type suggested for use in a solar receiver for an orbiting Brayton cycle power system. Energy storage for the eclipse portion of the cycle is provided by the latent heat of a eutectic mixture of LiF and CaF2 contained in the canister. The chief motivation for the study is the prediction of vapor void effects on temperature profiles and the identification of possible differences between ground test data and projected behavior in microgravity. The first phase of this study is based on a two-dimensional, cylindrical coordinates model using an interim procedure for describing void behavor in 1-g and microgravity. The thermal analysis includes the effects of solidification front behavior, conduction in liquid/solid salt and canister materials, void growth and shrinkage, radiant heat transfer across the void, and convection in the melt due to Marangoni-induced flow and, in 1-g, flow due to density gradients. A number of significant differences between 1-g and o-g behavior were found. This resulted from differences in void location relative to the maximum heat flux and a significantly smaller effective conductance in 0-g due to the absence of gravity-induced convection.

Wichner, R. P.↗

NASA NTRS1993

Dynamics of aging magnetic clouds

The dynamics of radially expanding magnetic clouds is rigorously analyzed within the framework of ideal MHD. The cloud is modelled as a cylindrically symmetric magnetic flux rope. In the force balance we include the gas pressure gradient and the Lorentz force. Interaction with the ambient solar wind due to expansion of the magnetic cloud is represented by a drag force proportional to the bulk velocity. We consider the self-similar expansion of a polytrope, and reduce the problem to an ordinary nonlinear differential equation for the evolution function. Analyzing the asymptotic behavior of the evolution function, we formulate theoretical expectations for the long-term behavior of cloud parameters. We focus on the temporal evolution of (1) the magnetic field strength; (2) the twist of the field lines; (3) the asymmetry of the total field profile; and (4) the bulk flow speed. We present data from two magnetic clouds observed at 1 AU and 2 AU, respectively, and find good agreement with theoretical expectations. For a peak magnetic field strength at 1 AU of 25 nT and a polytropic index of 0.5, we find that a magnetic cloud can be distinguished from the background interplanetary field up to a distance of about 5 AU. Taking larger magnetic fields and bigger polytropic indices this distance can double.

Osherovich, V. A.↗

NASA NTRS1999

Control of skeletal muscle perfusion at the onset of dynamic exercise

At the onset of exercise there is a rapid increase in skeletal muscle vascular conductance and blood flow. Several mechanisms involved in the regulation of muscle perfusion have been proposed to initiate this hyperemic response, including neural, metabolic, endothelial, myogenic, and muscle pump mechanisms. Investigators utilizing pharmacological blockade of cholinergic muscarinic receptors and sympathectomy have concluded that neither sympathetic cholinergic nor adrenergic neural mechanisms are involved in the initial hyperemia. Studies have also shown that the time course for vasoactive metabolite release, diffusion, accumulation, and action is too long to account for the rapid increase in vascular conductance at the initiation of exercise. Furthermore, there is little or no evidence to support an endothelium or myogenic mechanism as the initiating factor in the muscle hyperemia. Thus, the rise in muscle blood flow does not appear to be explained by known neural, metabolic, endothelial, or myogenic influences. However, the initial hyperemia is consistent with the mechanical effects of the muscle pump to increase the arteriovenous pressure gradient across muscle. Because skeletal muscle blood flow is regulated by multiple and redundant mechanisms, it is likely that neural, metabolic, and possibly endothelial factors become important modulators of mechanically induced exercise hyperemia following the first 5-10 s of exercise.

Review, Tutorial↗

NASA NTRSDate not supplied

Simulating Underexpanded Plumes in Martian and Lunar Environments

While the Loci/Chem Computational fluid dynamics (CFD) application has been validated and used extensively for simulating launch environments in atmospheric conditions, use of this tool for simulating supersonic plumes at much lower ambient pressure requires further validation. Simulations of underexpanded plumes have been performed using Loci/Chem at Martian pressure to predict several metrics for laminar, turbulent, and impinging plumes. The CFD results are compared with an experimental data set for low Reynolds number plumes at these conditions in order to evaluate the current capability of the Loci/Chem tool for these types of environments. The CFD validation results to date show reasonable agreement with the experimental data across all the metrics of interest for the configurations considered. The Loci/Chem-Boltzmann tool is a hybrid continuum/rarefied flow solver which extends modeling capabilities to very low-pressure environments such as those on the Moon. The Loci/Chem-Boltzmann solver uses a gradient-based continuum breakdown criterion to restrict solution of the computationally expensive Boltzmann equation to only a subset of the domain, while using the Navier-Stokes equations elsewhere. This CFD application is under active development, and current efforts toward establishing production ready capability for evaluating Lunar plume surface interactions are well under way. Initial simulations of an Apollo Lunar excursion module indicate that several regions of the flow require solving the Boltzmann equations due to extreme rarefaction. Early simulations are promising, indicating reasonable overall computational time for a full 3D human scale lander simulation.

Plume Surface Interaction↗

DOE OSTI2023

Physics-Driven Convolutional Autoencoder Approach for CFD Data Compressions: Preprint

With the growing size and complexity of turbulent flow models, data compression approaches are of the utmost importance to analyze, visualize, or restart the simulations. Recently, in-situ autoencoder-based compression approaches have been proposed and shown to be effective at producing reduced representations of turbulent flow data. However, these approaches focus solely on training the model using point-wise sample reconstruction losses that do not take advantage of the physical properties of turbulent flows. In this paper, we show that training autoencoders with additional physics-informed regularizations, e.g., enforcing incompressibility and preserving enstrophy, improves the compression model in three ways: (i) the compressed data better conform to known physics for homogeneous isotropic turbulence without negatively impacting point-wise reconstruction quality, (ii) inspection of the gradients of the trained model uncovers changes to the learned compression mapping that can facilitate the use of explainability techniques, and (iii) as a performance byproduct, training losses are shown to converge up to 12x faster than the baseline model.

auto-encoders↗

DOE OSTI2022

Plankton energy flows using a global size-structured and trait-based model

Plankton community models are critical tools for understanding the processes that shape marine plankton communities, how plankton communities impact biogeochemical cycles, and the feedbacks between community structure and function. Here, using the flexible Marine Biogeochemistry Library (MARBL), we present the Size-based Plankton ECological TRAits (MARBL-SPECTRA) model, which is designed to represent a diverse plankton community while remaining computationally tractable. MARBL-SPECTRA is composed of nine phytoplankton and six zooplankton size classes represented using allometric scaling relationships for physiological traits and interactions within multiple functional types. MARBL-SPECTRA is embedded within the global ocean component of the Community Earth System Model (CESM) and simulates large-scale, emergent patterns in phytoplankton growth limitation, plankton phenology, plankton generation time, and trophic transfer efficiency. The model qualitatively reproduces observed global patterns of surface nutrients, chlorophyll biomass, net primary production, and the biogeographies of a range of plankton size classes. In addition, the model simulates how predator:prey dynamics and trophic efficiency vary across gradients in total ecosystem productivity. Shorter food chains that export proportionally more carbon from the surface to the ocean interior occur in productive, eutrophic regions, whereas in oligotrophic regions, the food chains are relatively long and export less organic matter from the surface. The union of functional type modeling with size-resolved, trait-based modeling approaches allows MARBL-SPECTRA to capture both large-scale elemental cycles and the structure of planktonic food webs affecting trophic transfer efficiency.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Water Level Data from Wells PLM1 and PLM6 for the East River Watershed, Colorado

This dataset (Williams et al., 2020) contains the original un-QA/QC-ed water level data for PLM1 and PLM6 and has been obsoleted. The data contained within this dataset is not to be used. Refer to Faybishenko et al., 2022 (DOI: 10.15485/1866836) for the latest QA/QC-ed data available via ESS-DIVE.This data set contains water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location. These wells used to monitor subsurface water and carbon inventories and fluxes at the East River Watershed, Colorado, USA. Complete metadata information on the PLM1 and PLM6 wells are available in the related data package reference Varadharajan C, et al (2020). https://doi.org/10.15485/1660962.Data are reported in .csv files per well. The latitude and longitude of each location are given in a file called locations.csv. These data are used for determining the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry can be used to estimate rates of solute export from the hillslope to the floodplain and river.These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales.

54 ENVIRONMENTAL SCIENCES↗