Search NASA⌕ Search

SEARCH · Search NASA

Results for “Differential transport”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Directing Assembly of Mesoscale Multi‐Shell Morphologies of DNA Origami Crystals

Nature builds hierarchically ordered materials, such as seashells, wood, and bones, through spatially and temporally regulated growth. Mimicking such a level of control in synthetic systems remains challenging, particularly in achieving multiscale organizations with prescribed nanoscale arrangements and desired material morphologies. In this study, we introduce a DNA-based self-assembly strategy for constructing diverse multi-shell mesoscale morphologies from nanoscale lattices, enabling prescribed structural, and compositional 3D material patterns. Using DNA origami frames as modular monomers, we direct anisotropic epitaxial growth through addressable DNA frame binding motifs and encapsulate nanoparticles (NPs) in desired 3D patterns. Sequential monomer addition under thermodynamically favorable conditions enables shell growth through heterogeneous nucleation while minimizing unwanted homogeneous nucleation. Here, we demonstrate that DNA-encoded addressability enables epitaxial shell growth along specific lattice directions, yielding crystals with multilayered mesoscale organization, including tube-like (sushi roll) and plate-like (macaron) morphologies. Shell-specific NP configurations and compositions are achieved through addressable and differentiated placement of NPs within each shell, as validated by small-angle x-ray scattering and cross-sectional scanning transmission electron microscopy. We further demonstrate addressable NP release and reveal that shells modulate release kinetics. Together, these findings establish a platform for fabricating DNA origami crystals with programmable mesoscale morphologies, nanoscale structure, composition, and transport properties.

3D patterning↗
DOE OSTI2024

Cell‐type‐specific transcriptomics uncovers spatial regulatory networks in bioenergy sorghum stems

SUMMARY Bioenergy sorghum is a low‐input, drought‐resilient, deep‐rooting annual crop that has high biomass yield potential enabling the sustainable production of biofuels, biopower, and bioproducts. Bioenergy sorghum's 4–5 m stems account for ~80% of the harvested biomass. Stems accumulate high levels of sucrose that could be used to synthesize bioethanol and useful biopolymers if information about cell‐type gene expression and regulation in stems was available to enable engineering. To obtain this information, laser capture microdissection was used to isolate and collect transcriptome profiles from five major cell types that are present in stems of the sweet sorghum Wray. Transcriptome analysis identified genes with cell‐type‐specific and cell‐preferred expression patterns that reflect the distinct metabolic, transport, and regulatory functions of each cell type. Analysis of cell‐type‐specific gene regulatory networks (GRNs) revealed that unique transcription factor families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell‐specific transcriptome data was combined with known secondary cell wall (SCW) networks to identify the GRNs that differentially activate SCW formation in vascular sclerenchyma and epidermal cells. The spatial transcriptomic dataset provides a valuable source of information about the function of different sorghum cell types and GRNs that will enable the engineering of bioenergy sorghum stems, and an interactive web application developed during this project will allow easy access and exploration of the data ( https://mc‐lab.shinyapps.io/lcm‐dataset/ ).

09 BIOMASS FUELS↗
DOE OSTI2021

Towards a Unified theory of Fractional and Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green’s identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green’s identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. Here, the proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Unified Theory of Fractional Nonlocal and Weighted Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green's identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green's identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. The proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Genomics, Exometabolomics, and Metabolic Probing Reveal Conserved Proteolytic Metabolism of Thermoflexus hugenholtzii and Three Candidate Species From China and Japan

Thermoflexus hugenholtzii JAD2 T , the only cultured representative of the Chloroflexota order Thermoflexales, is abundant in Great Boiling Spring (GBS), NV, United States, and close relatives inhabit geothermal systems globally. However, no defined medium exists for T. hugenholtzii JAD2 T and no single carbon source is known to support its growth, leaving key knowledge gaps in its metabolism and nutritional needs. Here, we report comparative genomic analysis of the draft genome of T. hugenholtzii JAD2 T and eight closely related metagenome-assembled genomes (MAGs) from geothermal sites in China, Japan, and the United States, representing “Candidatus Thermoflexus japonica,” “Candidatus Thermoflexus tengchongensis,” and “Candidatus Thermoflexus sinensis.” Genomics was integrated with targeted exometabolomics and 13 C metabolic probing of T. hugenholtzii. The Thermoflexus genomes each code for complete central carbon metabolic pathways and an unusually high abundance and diversity of peptidases, particularly Metallo- and Serine peptidase families, along with ABC transporters for peptides and some amino acids. The T. hugenholtzii JAD2 T exometabolome provided evidence of extracellular proteolytic activity based on the accumulation of free amino acids. However, several neutral and polar amino acids appear not to be utilized, based on their accumulation in the medium and the lack of annotated transporters. Adenine and adenosine were scavenged, and thymine and nicotinic acid were released, suggesting interdependency with other organisms in situ. Metabolic probing of T. hugenholtzii JAD2 T using 13 C-labeled compounds provided evidence of oxidation of glucose, pyruvate, cysteine, and citrate, and functioning glycolytic, tricarboxylic acid (TCA), and oxidative pentose-phosphate pathways (PPPs). However, differential use of position-specific 13 C-labeled compounds showed that glycolysis and the TCA cycle were uncoupled. Thus, despite the high abundance of Thermoflexus in sediments of some geothermal systems, they appear to be highly focused on chemoorganotrophy, particularly protein degradation, and may interact extensively with other microorganisms in situ.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Grid-Integrated Production of Fischer-Tropsch Synfuels from Nuclear Power

Idaho National Laboratory (INL) investigates the relative economic profitability of an integrated energy system (IES) coupling an NPP with a synfuel production process at selected case study locations across the United States. In the synfuel IES, a high-temperature steam electrolysis (HTSE) plant is thermally and electrically coupled with an NPP to produce zero-carbon hydrogen. The synthetic fuel is produced from this H 2 combined with a CO 2 supply using the reverse water gas shift process followed by the Fischer-Tropsch (FT) reaction. This analysis considers a system in which the CO 2 is sourced from regional CO 2 emitters via the construction and operation of pipeline-based CO 2 supply networks. Locating the FT plant at the same site as the NPP and HTSE plants enables the NPP to provide zero-carbon heat and power to the HTSE plant and zero-carbon power to the FT plant as well as avoid the requirement for long-distance H 2 product transport from the HTSE plant to the FT plant. Hydrogen storage is used to enable the NPP to dispatch power to the electrical grid (instead of the HTSE plant) when grid demand increases, thus enabling the FT plant to continue to operate in a steady-state production mode. The ability to cease hydrogen production for several hours within each day enables the NPP to provide power to the grid to balance the electricity market during peak periods and maximize revenues for the nuclear synfuel IES. The FT process design considered has a 99% carbon conversion efficiency. The use of nuclear energy and nuclear energy-derived hydrogen enables synfuel production to achieve this high level of carbon utilization. Additionally, the life-cycle carbon emissions of the nuclear-based synfuel production process are very low, with WTW emissions of approximately 25 gCO 2 e/MJ, including steam credits (generated from FT process excess heat), and approximately 7 gCO 2 e/MJ, if steam credits are excluded. This compares favorably with the WTW emissions of 90.5 gCO 2 e/MJ for a compression-ignition, direct injection (CIDI) vehicle with a fuel economy of 31.6 miles per gallon gasoline equivalent (MPGGE), using low-sulfur diesel produced using conventional petroleum production and refining processes. Several NPPs in various regions of the U.S. are considered as case study analyses. Supply locations and transportation via pipeline of the CO 2 feedstock to the NPP site are analyzed through the National Energy Technology Laboratory (NETL) CO 2 Transport Cost model. The team finds that the amount of CO 2 generated by different sectors is sufficient for the synfuel production process at all locations considered. The CO 2 transportation costs are functions of the distance of the source to the NPP location, the CO 2 capture cost at the source, and the quantity of CO 2 transported. Historical electricity prices for the NPP case study locations are collected and analyzed. Monthly average prices, price range, and duration of negative-price periods vary among these locations. For each location, an auto-regressive moving average (ARMA) model is trained on historical electricity price data. ARMA validation is done to ensure the synthetic price distributions represent one of historical prices with high fidelity. Synthetic time series from these ARMA models are used in a coupled dispatch and system optimization in the Holistic Energy Resource Optimization Network (HERON) to compute the differential net present value (NPV) of the IES. The team finds that this econometric is positive, ranging from $14M–1.3bn (2020) depending on the location. The optimal synfuel IES configuration to obtain this increase in NPV often maximizes the size of the synfuel production process with regards to the size of the NPP. However, the team shows that the NPP still plays a stabilizing role for the grid: In periods of high prices and high loads, more electricity from the NPP is sent to the grid. A high variability of electricity prices and extreme maximum prices tend to drive up electricity production. While it requires significant investment, the synfuel IES could increase the economic profitability for the existing fleet of LWRs across the country while still maintaining the grid stabilizer role of NPPs. During its lifetime, the main costs for the nuclear synfuel IES are the carbon feedstock transportation costs, followed by the capital expenses (CAPEX) and operation and maintenance (O&M) costs while the revenue comes first from the IRA H 2 production tax credit (PTC) and then from the sales of synfuel products. The profitability of the synfuel IES is most sensitive to the value of the hydrogen PTC and the synfuel products as well as the cost of the carbon feedstock, highlighting the importance of governmental incentives regarding hydrogen, carbon emissions, and synfuel in driving the deployment of future nuclear synfuel IESs.

08 HYDROGEN↗
DOE OSTI2025

Emergent Nanostructure and Ion Transport in Polyzwitterion/Polyanion Blends

We investigated blends of poly(1-(3-sulfonatopropyl)-2-vinylpyridinium) (P2VPPS) and poly(lithium (trifluoromethane)sulfonimide methacrylate) (poly(MTFSI)Li) at varying molar ratios to gain a mechanistic understanding of ionic conductivity in a miscible polyzwitterion/polyanion system. This dataset contains the raw numerical data corresponding to the figures in the manuscript. The data files include the following information: (1) Experimental Data – includes X-ray and neutron scattering measurements, broadband dielectric spectroscopy (BDS) data, extracted DC conductivity values, differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) results, extracted glass transition temperatures, etc. (2) CGMD Data – includes molecular dynamics (MD) trajectory files and computed structural correlations. All data files are organized/named according to the figure numbers in the manuscript.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Setup and verification of a SCALE/KENO platform for generic FHR benchmark calculations

The work presented in this article is preliminary to downstream analysis of a generic fluoride salt-cooled high-temperature reactor (gFHR) core performed by the University of Tennessee in collaboration with Kairos Power (KP). A Monte Carlo transport model of the publicly available gFHR equilibrium core is developed in SCALE/KENO with multigroup energy treatment. Several output quantities of interest are used to verify the simulations against a benchmark model developed by KP using the continuous energy Monte Carlo code Serpent 2. Good agreement is seen in flux and fission rate profiles with a maximum relative difference of 1.4% and 2.8% respectively. Furthermore, an effective multiplication factor bias of 44 pcm was observed between the two simulations. The fuel temperature reactivity coefficient calculated with SCALE is within uncertainty to the reference model. This verification acts as a publicly reproducible benchmark for the gFHR in SCALE/KENO. A simplified depletion model is also presented where a single fuel pebble is depleted to discharge burnup through the equilibrium core while the equilibrium core is assumed to be invariant. This method produces results that intercept the equilibrium core concentrations in every case, however, an interesting artifact of this particular depletion model is uncovered. The phenomenon is shown to be a fundamental feature of the differential rate equations and inspires questions about how this system behaves when the time evolution of the equilibrium core is considered. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Sensitivity of Au + Au collisions to the symmetric nuclear matter equation of state at 2–5 nuclear saturation densities

We demonstrate that proton and pion flow measurements in heavy-ion collisions at incident energies ranging from 1 to 20 GeV per nucleon in the fixed target frame can be used for an accurate determination of the symmetric nuclear matter equation of state at baryon densities equal 2–4 times nuclear saturation density n 0 . We simulate Au + Au collisions at these energies using a hadronic transport model with an adjustable vector mean-field potential dependent on baryon density n B . Here, we show that the mean field can be parametrized to reproduce a given density dependence of the speed of sound at zero temperature $c$$^{2}_{s}$ (n B , T = 0), which we vary independently in multiple density intervals to probe the differential sensitivity of heavy-ion observables to the equation of state at these specific densities. Recent flow data from the STAR experiment at the center-of-mass energies √ s NN = {3.0, 4.5} GeV can be described by our model, and a Bayesian analysis of these data indicates a hard equation of state at n B ϵ (2, 3)n 0 and a possible phase transition at n B ϵ (3, 4)n 0 . More data at √ s NN = 2–5 GeV, as well as a more thorough analysis of the model systematic uncertainties will be necessary for a more precise conclusion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Deinococcus radiodurans UWO298 Dependence on Background Radiation for Optimal Growth

Ionizing radiation is a major environmental variable for cells on Earth, and so organisms have adapted to either prevent or to repair damages caused by it, primarily from the appearance and accumulation of reactive oxygen species (ROS). In this study, we measured the differential gene expression in Deinococcus radiodurans UWO298 cultures deprived of background ionizing radiation (IR) while growing 605 m underground at the Waste Isolation Pilot Plant (WIPP), reducing the dose rate from 72.1 to 0.9 nGy h –1 from control to treatment, respectively. This reduction in IR dose rate delayed the entry into the exponential phase of the IR-shielded cultures, resulting in a lower biomass accumulation for the duration of the experiment. The RNASeq-based transcriptome analysis showed the differential expression of 0.2 and 2.7% of the D. radiodurans genome after 24 and 34 h of growth in liquid culture, respectively. Gene expression regulation after 34 h was characterized by the downregulation of genes involved in folding newly synthesized and denatured/misfolded proteins, in the assimilation of nitrogen for amino acid synthesis and in the control of copper transport and homeostasis to prevent oxidative stress. We also observed the upregulation of genes coding for proteins with transport and cell wall assembly roles. These results show that D. radiodurans is sensitive to the absence of background levels of ionizing radiation and suggest that its transcriptional response is insufficient to maintain optimal growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics

We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Small Signaling Peptides in Sorghum bicolor : Integrating Phylogeny and Gene Expression to Characterize Roles in Stem Development

Small signaling peptides (SSPs) are important regulators of plant growth, development, and responses to biotic and abiotic stress, yet their role in the C4 grass Sorghum bicolor is largely uncharacterized. To help fill this knowledge gap, 219 sorghum genes that encode SSPs were identified based on SSP sequences previously identified in Arabidopsis thaliana, Zea mays, Oryza sativa, Triticum aestivum , and Brachypodium distachyon . The 219 sorghum SSP-encoding genes were assigned to 19 gene families, analyzed for the presence of motifs, and aligned with genes that encode SSPs in other plants using phylogenetic analysis. Sorghum genes in 12 of the 19 SSP gene families had not been previously characterized. Expression of the 219 SSP-encoding genes in sorghum organs, during stem development, and in stem tissues and cell types revealed distinct spatial, temporal, and developmental patterns of expression. Genes associated with the SbCEP and SbRGF families were preferentially expressed in roots, whereas SbEPF genes were expressed in stem epidermal and pith parenchyma cells and panicles. The expression of genes during bioenergy sorghum stem growth and development was investigated because stems account for ~80% of harvested biomass and serve as conduits for water and nutrient transport between leaves and roots. During stem development, 28 SSP genes in several families ( CLE, EPF, CEP, GASS, PSY, ES, PSK, CAPE, POE ) were expressed at higher levels in zones of cell proliferation. For example, the TDIF homologs SbCLE41 and SbCLE42 were expressed at high levels in nascent stem nodes where they may regulate vascular bundle cambial activity and cell differentiation. A different set of 15 genes in the CIF, POE, CAPE, PSY, CEP, RALF , and CLE families were expressed at higher levels in zones of stem tissue differentiation highlighted by elevated expression of five SbRALFR s in the stem nodal plexus. Cell type–specific expression of many sorghum genes that encode SSPs was observed in fully elongated internodes indicating gene expression is regulated with high spatial resolution. Overall, the results provide a foundation of information for analysis of SSP function in sorghum that can be integrated with knowledge of sorghum gene regulatory networks to modulate traits important for production of sorghum crops.

bioenergy sorghum↗
DOE OSTI2023

High-Performance Electron Sources: Numerical Methods and Beam Dynamics at the Precision Frontier

Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Nanodevices engineering and spin transport properties of MnBi2Te4 monolayer

Abstract Two-dimensional (2D) magnetic materials are essential for the development of the next-generation spintronic technologies. Recently, layered van der Waals (vdW) compound MnBi 2 Te 4 (MBT) has attracted great interest, and its 2D structure has been reported to host coexisting magnetism and topology. Here, we design several conceptual nanodevices based on MBT monolayer (MBT-ML) and reveal their spin-dependent transport properties by means of the first-principles calculations. The pn -junction diodes and sub-3-nm pin -junction field-effect transistors (FETs) show a strong rectifying effect and a spin filtering effect, with an ideality factor n close to 1 even at a reasonably high temperature. In addition, the pip - and nin -junction FETs give an interesting negative differential resistive (NDR) effect. The gate voltages can tune currents through these FETs in a large range. Furthermore, the MBT-ML has a strong response to light. Our results uncover the multifunctional nature of MBT-ML, pave the road for its applications in diverse next-generation semiconductor spin electric devices.

Chemistry↗
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Results from helioseismology indicate that the radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical models using the mean-field approach have been successful in explaining this property of the NSSL, while global direct or large-eddy magnetoconvection models have so far been unable to reproduce this. We investigate the reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under conditions mimicking those in the solar NSSL. Simulations with as few ingredients as possible to generate mean flows were studied. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. The parameters of the simulations were chosen such that they matched the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch was varied by changing the rotation rate such that the resulting Coriolis numbers covered the same range as in the NSSL. We measured the turbulent transport coefficient relevant for the nondiffusive (Λ-effect) and diffusive (turbulent viscosity) parts of the Reynolds stress and compared them with predictions of current mean-field theories. A negative radial gradient of the mean flow is generated only at the equator where meridional flows are absent. At other latitudes, the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that the meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional flows and stresses are neglected. However, the rotational dependence of the turbulent viscosity in our simulations agrees well with theoretical predictions. Moreover, our results agree qualitatively with global convection simulations in that an NSSL can only be obtained near the equator.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗