Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conservative variables”

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 451 records · Page 25

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗

Small particle transport across turbulent nonisothermal boundary layers

The interaction between turbulent diffusion, Brownian diffusion, and particle thermophoresis in the limit of vanishing particle inertial effects is quantitatively modeled for applications in gas turbines. The model is initiated with consideration of the particle phase mass conservation equation for a two-dimensional boundary layer, including the thermophoretic flux term directed toward the cold wall. A formalism of a turbulent flow near a flat plate in a heat transfer problem is adopted, and variable property effects are neglected. Attention is given to the limit of very large Schmidt numbers and the particle concentration depletion outside of the Brownian sublayer. It is concluded that, in the parameter range of interest, thermophoresis augments the high Schmidt number mass-transfer coefficient by a factor equal to the product of the outer sink and the thermophoretic suction.

Rosner, D. E.↗

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulence models. The difficulty lies in the fact that the reaction rate is in general an exponential function of the temperature, and the higher order correlations in the conventional moment closure models of the chemical source term can not be neglected, making the applications of such models impractical. The probability density function (pdf) method offers an attractive alternative: in a pdf model, the chemical source terms are closed and do not require additional models. A grid dependent Monte Carlo scheme was studied, since it is a logical alternative, wherein the number of computer operations increases only linearly with the increase of number of independent variables, as compared to the exponential increase in a conventional finite difference scheme. A new algorithm was devised that satisfies a restriction in the case of pure diffusion or uniform flow problems. Although for nonuniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗

Alternative strategies of nutrient acquisition and energy conservation map to the biogeography of marine ammonia-oxidizing archaea

Ammonia-oxidizing archaea (AOA) are among the most abundant and ubiquitous microorganisms in the ocean, exerting primary control on nitrification and nitrogen oxides emission. Although united by a common physiology of chemoautotrophic growth on ammonia, a corresponding high genomic and habitat variability suggests tremendous adaptive capacity. Here, we compared 44 diverse AOA genomes, 37 from species cultivated from samples collected across diverse geographic locations and seven assembled from metagenomic sequences from the mesopelagic to hadopelagic zones of the deep ocean. Comparative analysis identified seven major marine AOA genotypic groups having gene content correlated with their distinctive biogeographies. Phosphorus and ammonia availabilities as well as hydrostatic pressure were identified as selective forces driving marine AOA genotypic and gene content variability in different oceanic regions. Notably, AOA methylphosphonate biosynthetic genes span diverse oceanic provinces, reinforcing their importance for methane production in the ocean. Together, our combined comparative physiological, genomic, and metagenomic analyses provide a comprehensive view of the biogeography of globally abundant AOA and their adaptive radiation into a vast range of marine and terrestrial habitats.

59 BASIC BIOLOGICAL SCIENCES↗

Evolution of oxygen and stratification and their relationship in the North Pacific Ocean in CMIP6 Earth system models

Abstract. This study examines the linkages between the upper-ocean (0–200 m) oxygen (O2) content and stratification in the North Pacific Ocean using four Earth system models (ESMs), an ocean hindcast simulation, and an ocean reanalysis. The trends and variability in oceanic O2 content are driven by the imbalance between physical supply and biological demand. Physical supply is primarily controlled by ocean ventilation, which is responsible for the transport of O2-rich surface waters to the subsurface. Isopycnic potential vorticity (IPV), a quasi-conservative tracer proportional to density stratification that can be evaluated from temperature and salinity measurements, is used herein as a dynamical proxy for ocean ventilation. The predictability potential of the IPV field is evaluated through its information entropy. The results highlight a strong O2–IPV connection and somewhat higher (as compared to the rest of the basin) predictability potential for IPV across the tropical Pacific, where the El Niño–Southern Oscillation occurs. This pattern of higher predictability and strong anticorrelation between O2 and stratification is robust across multiple models and datasets. In contrast, IPV at mid-latitudes has low predictability potential and its center of action differs from that of O2. In addition, the locations of extreme events or hotspots may or may not differ between the two fields, with a strong model dependency, which persists in future projections. On the one hand, these results suggest that it may be possible to monitor ocean O2 in the tropical Pacific based on a few observational sites co-located with the more abundant IPV measurements; on the other, they lead us to question the robustness of the IPV–O2 relationship in the extratropics. The proposed framework helps to characterize and interpret O2 variability in relation to physical variability and may be especially useful in the analysis of new observation-based data products derived from the BGC-Argo float array in combination with the traditional but far more abundant Argo data.

Novi, Lyuba↗

Integrable higher-spin deformations of sigma models from auxiliary fields

We construct a new infinite family of integrable deformations of the principal chiral model (PCM) parametrized by an interaction function of several variables, which extends the formalism of [C. Ferko and L. Smith, An infinite family of integrable sigma models using auxiliary fields, .] and includes deformations of the PCM by functions of both the stress tensor and higher-spin conserved currents. We show in detail that every model in this class admits a Lax representation for its equations of motion, and that the Poisson bracket of the Lax connection takes the Maillet form, establishing the existence of an infinite set of Poisson-commuting conserved charges. We argue that the non-Abelian T-dual of any model in this family is classically integrable, and that T-duality “commutes” with a general deformation in this class, in a sense which we make precise. Finally, we demonstrate that these higher-spin auxiliary field deformations can be extended to accommodate the addition of a Wess-Zumino term, and we exhibit the Lax connection in this case. Published by the American Physical Society 2025

Bielli, Daniele (ORCID:0009000640034768)↗

Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone

We examined phylogeographic structure in gray fox (Urocyon cinereoargenteus) across the United States to identify the location of secondary contact zone(s) between eastern and western lineages and investigate the possibility of additional cryptic intraspecific divergences. We generated and analyzed complete mitochondrial genome sequence data from 75 samples and partial control region mitochondrial DNA sequences from 378 samples to investigate levels of genetic diversity and structure through population- and individual-based analyses including estimates of divergence (FST and SAMOVA), median joining networks, and phylogenies. We used complete mitochondrial genomes to infer phylogenetic relationships and date divergence times of major lineages of Urocyon in the United States. Despite broad-scale sampling, we did not recover additional major lineages of Urocyon within the United States, but identified a deep east-west split (~0.8 million years) with secondary contact at the Great Plains Suture Zone and confirmed the Channel Island fox (Urocyon littoralis) is nested within U. cinereoargenteus. Genetic diversity declined at northern latitudes in the eastern United States, a pattern concordant with post-glacial recolonization and range expansion. Beyond the east-west divergence, morphologically-based subspecies did not form monophyletic groups, though unique haplotypes were often geographically limited. Gray foxes in the United States displayed a deep, cryptic divergence suggesting taxonomic revision is needed. Secondary contact at a common phylogeographic break, the Great Plains Suture Zone, where environmental variables show a sharp cline, suggests ongoing evolutionary processes may reinforce this divergence. Follow-up study with nuclear markers should investigate whether hybridization is occurring along the suture zone and characterize contemporary population structure to help identify conservation units. Comparative work on other wide-ranging carnivores in the region should test whether similar evolutionary patterns and processes are occurring.

54 ENVIRONMENTAL SCIENCES↗

Attribution of the Australian bushfire risk to anthropogenic climate change

Abstract. Disastrous bushfires during the last months of 2019 and January 2020 affected Australia, raising the question to what extent the risk of these fires was exacerbated by anthropogenic climate change. To answer the question for southeastern Australia, where fires were particularly severe, affecting people and ecosystems, we use a physically based index of fire weather, the Fire Weather Index; long-term observations of heat and drought; and 11 large ensembles of state-of-the-art climate models. We find large trends in the Fire Weather Index in the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) Atmospheric Reanalysis (ERA5) since 1979 and a smaller but significant increase by at least 30 % in the models. Therefore, we find that climate change has induced a higher weather-induced risk of such an extreme fire season. This trend is mainly driven by the increase of temperature extremes. In agreement with previous analyses we find that heat extremes have become more likely by at least a factor of 2 due to the long-term warming trend. However, current climate models overestimate variability and tend to underestimate the long-term trend in these extremes, so the true change in the likelihood of extreme heat could be larger, suggesting that the attribution of the increased fire weather risk is a conservative estimate. We do not find an attributable trend in either extreme annual drought or the driest month of the fire season, September–February. The observations, however, show a weak drying trend in the annual mean. For the 2019/20 season more than half of the July–December drought was driven by record excursions of the Indian Ocean Dipole and Southern Annular Mode, factors which are included in the analysis here. The study reveals the complexity of the 2019/20 bushfire event, with some but not all drivers showing an imprint of anthropogenic climate change. Finally, the study concludes with a qualitative review of various vulnerability and exposure factors that each play a role, along with the hazard in increasing or decreasing the overall impact of the bushfires.

van Oldenborgh, Geert Jan (ORCID:0000000268989535)↗

Drag characteristics of unsteady, perturbed boundary layer flows

A series of time-dependent numerical computations have been performed for the flow of an unsteady boundary layer over an embedded cavity, which is aligned normal to the flow direction. The unsteady flow is the result of the imposition of a Stuart vortex, at the inflow boundary of the computational domain, onto an otherwise unperturbed two-dimensional laminar boundary layer flow. This produces a strongly vortical, but spatially monochromatic, motion at inflow which is allowed to evolve downstream through the time-dependent numerical solution of appropriate conservation equations. As the imposed vortex evolves over the solid boundary, a weaker induced vortical motion is produced and together these perturbing motions interact with the fluid motion in the embedded cavity. An analysis of the pertinent dynamic variables as well as the time-dependent drag characteristics for this type of non-planar geometry is performed; with particular emphasis on the relative contribution of both pressure drag and frictional drag at various times in the evolution process.

Gatski, T. B.↗

The alpha(3) Scheme - A Fourth-Order Neutrally Stable CESE Solver

The conservation element and solution element (CESE) development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To provide a solid foundation for a systematic CESE development of high order schemes, in this paper we describe a new 4th-order neutrally stable CESE solver of the advection equation Theta u/Theta + alpha Theta u/Theta x = 0. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables u(sup n) (sub j), (u(sub x))(sup n) (sub j) , and (uxz)(sup n) (sub j) (the numerical analogues of u, Theta u/Theta x, and Theta(exp 2)u/Theta x(exp 2), respectively) and four equations per mesh point, the new scheme is referred to as the alpha(3) scheme. As in the case of other similar CESE neutrally stable solvers, the alpha(3) scheme enforces conservation laws in space-time locally and globally, and it has the basic, forward marching, and backward marching forms. These forms are equivalent and satisfy a space-time inversion (STI) invariant property which is shared by the advection equation. Based on the concept of STI invariance, a set of algebraic relations is developed and used to prove that the alpha(3) scheme must be neutrally stable when it is stable. Moreover it is proved rigorously that all three amplification factors of the alpha(3) scheme are of unit magnitude for all phase angles if |v| <= 1/2 (v = alpha delta t/delta x). This theoretical result is consistent with the numerical stability condition |v| <= 1/2. Through numerical experiments, it is established that the alpha(3) scheme generally is (i) 4th-order accurate for the mesh variables u(sup n) (sub j) and (ux)(sup n) (sub j); and 2nd-order accurate for (uxx)(sup n) (sub j). However, in some exceptional cases, the scheme can achieve perfect accuracy aside from round-off errors.

Chang, Sin-Chung↗

Monitoring Avian Productivity and Survivorship at the Sandia Wetlands at Los Alamos National Laboratory (2014-2019)

Biologists at Los Alamos National Laboratory completed six years (2014–2019) of demographic monitoring of passerines (songbirds) during the breeding season in Los Alamos, New Mexico. We captured songbirds at a mist netting station located in the Sandia Wetlands in Sandia Canyon (Technical Areas 60 and 61). We identified, measured, and then banded captured birds with federal migratory bird bands. Banding operations took place between May and August of each year, with 10 mist netting sessions annually. We conducted this project as part of the implementation of the Biological Resources Management Plan. The project complies with the 2013 Memorandum of Understanding between the United States Fish and Wildlife Service (USFWS) and the United States Department of Energy regarding the implementation of Executive Order 13186. The overall number of birds and species captured was variable, and we saw no decreases over the course of six years in bird species diversity, abundance, or the percentage of breeding birds. Our results indicate that the Sandia Wetlands support numerous species of breeding birds, including species of conservation concern. We did not see a declining trend for the top 10 most-captured species or for any of the species of conservation concern over time; however, we need more years of data to make robust conclusions about population trends through time. Our results indicated that adult recapture rates at the Sandia Wetlands are at least twice as low when compared with estimated adult survival in migratory passerines from peer-reviewed literature. We also did not recapture any birds that were banded as juveniles in subsequent years throughout the entire six-year period. Low recapture rates for adults and the non-existent recapture rates for juveniles should be explored further. To accurately assess individual species survival rates for adults and juveniles using mark-recapture models, we need a longer dataset, which would allow for a more robust analysis.

59 BASIC BIOLOGICAL SCIENCES↗

The terminal area simulation system. Volume 1: Theoretical formulation

A three-dimensional numerical cloud model was developed for the general purpose of studying convective phenomena. The model utilizes a time splitting integration procedure in the numerical solution of the compressible nonhydrostatic primitive equations. Turbulence closure is achieved by a conventional first-order diagnostic approximation. Open lateral boundaries are incorporated which minimize wave reflection and which do not induce domain-wide mass trends. Microphysical processes are governed by prognostic equations for potential temperature water vapor, cloud droplets, ice crystals, rain, snow, and hail. Microphysical interactions are computed by numerous Orville-type parameterizations. A diagnostic surface boundary layer is parameterized assuming Monin-Obukhov similarity theory. The governing equation set is approximated on a staggered three-dimensional grid with quadratic-conservative central space differencing. Time differencing is approximated by the second-order Adams-Bashforth method. The vertical grid spacing may be either linear or stretched. The model domain may translate along with a convective cell, even at variable speeds.

Proctor, F. H.↗

Modeling of Thermal Performance of Multiphase Nuclear Fuel Cell Under Variable Gravity Conditions

A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.

Ding, Z.↗

Modeling Variability of Plasma Conditions in the Io Torus

Telescopic observations an in situ measurements of the Io plasma torus show the density, temperature and composition to vary over time, sometimes up to a factor of 2. While previous models of the physical and chemical processes in the Io plasma torus have reasonably modeled the conditions of the Voyager 1 era, their authors have not addressed the observed variability nor explored the sensitivity of torus conditions to input parameters. In this paper we present a homogeneous torus model parameterized by five variables (transport timescale, neutral source strength, ratio of oxygen sulfur to atoms in the source, fraction of superthermal electrons, temperature of these hot electrons). The model incorporates the most recent data for ionization, recombination, charge exchange and radiative energy losses for the major torus species (S, S(sup +), S(sup ++), S(sup +++), O, O(sup +), O(sup ++). We solve equations of conservation of mass and energy to find equilibrium conditions for a set of input parameters. We compare model plasma conditions with those observed by Voyager 1 Voyager 2, and Cassini. Furthermore, we explore the sensitivity of torus conditions to each parameter. We find that (1) torus conditions are distinctly different for the Voyager 1, Voyager 2 and Cassini eras, (2) unique torus input parameters for any given era are poorly constrained given the wide range of solution space that is consistent with the range of observed torus conditions, (3) ion composition is highly sensitive to the specification of a non-thermal electron distribution, (4) neutral O/S source ratio is highly variable with model values ranging between 1.7 for Cassini to 4.0 for Voyager 1 conditions, (5) transport times range between 23 days for Voyager 2 to 50 days for Voyager 1 and Cassini, (6) neutral source strengths range between 7 to 30 x 10(sup -4) cm (sup -3) s(sup -1) which corresponds to a net production of 0.4 to 1.3 tons/s for a torus volume of 1.4 x 10(sup 31) cm(sup 3), or 38 R(sub j)(sup 3).

Delamere, P. A.↗

Agent-based simulations of shared automated vehicle operations: reflecting travel-party size, season and day-of-week demand variations

Here, this paper explores the effects of day of week and season of year demand variations for shared rides, along with realistic travel party sizes, on shared autonomous vehicle (SAV) services across the Austin, Texas region. Using the agent-based POLARIS program, synthetic person-trips that reflect travel-party size (from one to four persons) and demand variations over days and months, as evident in the National Household Travel Survey data were simulated in each scenario over a 24 h travel day. Results show that realistic party sizes can bring considerable changes to SAV fleet performance, including up to 8.5% higher service rates (number of requests accepted within 15 min), 5 min shorter journey times (wait time + travel time), 28% higher vehicle occupancies on weekends, and roughly 4% lower empty fleet VMT. Weekend travel is most impacted by season of year, with weekday travel patterns looking more uniform (thanks to work and school trips). Various performance metrics for the Austin network, like total and empty VMT, change by up to 30% when considering realistic variations in party size and time of year. This paper underscores the value of recognizing day-to-day and month-to-month variations in travel demand, and the importance of agent-based model equations to reflect travel-party size. Such realism can help quantify SAV seat occupancies more accurately, highlighting the importance of shared mobility. However, it also creates demand and supply issues for operators that now need more information on party size to manage dynamic ride-sharing, or those that may wish to shift their fleet vehicles to other regions for special events to protect profits while offering reasonable wait times to customers throughout the year.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗