Search NASA⌕ Search

SEARCH · Search NASA

Results for “MACROSCOPIC EQUATION”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Mechanisms of Zn removal from water by amorphous geopolymer: Molecular-level insights from X-ray absorption spectroscopy, isotope fractionation, and surface complexation modeling

Porous geopolymers have attracted widespread attention as promising heavy metal adsorbents that can be synthesized from aluminosilicate solid wastes. However, the precise microstructural evidence for adsorbed heavy metals on geopolymers remains unclear due to the insensitivity of conventional characterization techniques on minerals with amorphous structure and surface disorder. Batch adsorption and column experiments coupled with X-ray absorption spectroscopy (XAS), Zn stable isotope, and surface complexation model (SCM) were employed to reveal the Zn removal mechanisms with coal fly ash porous geopolymer (CFAPG) at a molecular scale. The macroscopic kinetic and isothermal adsorption of Zn on CFAPG were well described by the pseudo-second-order model and Bi_Langmuir equation, respectively, indicating the presence of abundant heterogeneous active sites on the CFAPG surface. Further, two types of active sites on the CFAPG surface were identified by XAS coupled with Zn isotopes in batch experiments at pH <= 6.0: one is pH-dependent and associated with tetrahedral zinc coordi-nation, and the other is pH-insensitive and associated with octahedral zinc coordination; these sites were confirmed by the bidentate SCM as the variable charge site (surface complexation, >S-OH) and the permanent negative charge site (cation exchange, >X - ), respectively. Furthermore, the important contribution of surface co -precipitation besides surface complexation and cation exchange to the Zn adsorption on CFAPG was identified by XAS coupled with SCM in a flow-through column experiment at pH >6.0. These investigations provide a systemic understanding of the Zn adsorption mechanisms on CFAPG and an SCM reference for the application and prediction of geopolymers in heavy metal-contaminated water remediations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Structural and Electronic Transitions in Liquid FeO Under High Pressure

Abstract FeO represents an important end‐member for planetary interiors mineralogy. However, its properties in the liquid state under high pressure are poorly constrained. Here, in situ high‐pressure and high‐temperature X‐ray diffraction experiments, ab initio simulations, and thermodynamic calculations are combined to study the local structure and density evolution of liquid FeO under extreme conditions. Our results highlight a strong shortening of the Fe‐Fe distance, particularly pronounced between ambient pressure and ∼40 GPa, possibly related with the insulator to metal transition occurring in solid FeO over a similar pressure range. Liquid density is smoothly evolving between 60 and 150 GPa from values calculated for magnetic liquid to those calculated for non‐magnetic liquid, compatibly with a continuous spin crossover in liquid FeO. The present findings support the potential decorrelation between insulator/metal transition and the high‐spin to low‐spin continuous transition, and relate the changes in the microscopic structure with macroscopic properties, such as the closure of the Fe‐FeO miscibility gap. Finally, these results are used to construct a parameterized thermal equation of state for liquid FeO providing densities up to pressure and temperature conditions expected at the Earth's core‐mantle boundary.

Morard, Guillaume↗
DOE OSTI2024

How macroscopic limits on neutron-star baryon loss yield microscopic limits on non-standard-model baryon decay

We investigate how our baryon-loss limits from anomalous binary-pulsar period lengthening can be interpreted microscopically to yield specific constraints on the particle physics of baryon number violation within a neutron star. We focus on the possibility of anomalous baryon disappearance via dark baryon processes and on scenarios in which the produced dark-sector particles do not survive to influence the response of the star to baryon-number-violating effects. We flesh out the conditions for which this may occur, as well as other key assumptions. We then turn to the analysis of particle processes in the dense nuclear medium found at the core of a neutron star, employing the techniques of relativistic mean-field theory. Using our study of in-medium effects and limits on macroscopic baryon number violation, we extract limits on in-vacuum baryon-number-violating processes, and we determine them for various equations of state. We conclude by noting the implications of our results for models of dark-sector-enabled baryogenesis. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Slowly rotating gravastars

We solve Einstein’s equations for slowly rotating gravitational condensate stars (gravastars) up to second order in the rotation by expanding about the spherically symmetric gravastar with de Sitter interior and Schwarzschild exterior matched at their common horizon. Requiring that the perturbations are finite on the null surface reduces the exterior geometry to that of a Kerr black hole, implying that a slowly rotating gravastar cannot be distinguished from a Kerr black hole by any measurement or observation restricted to the macroscopic spacetime exterior to the horizon. Here we determine the interior solution, the surface stress tensor, and the Komar mass and angular momentum localized on the slowly rotating horizon surface. With the interior equation of state fixed at p = –ρ, finite junction conditions on the null horizon surface necessarily lead to an interior solution with a singular core, where the perturbative expansion breaks down. Comparison to other models and implications for more rapidly rotating gravastars are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
DOE OSTI2020

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Tidal Love numbers of novel and admixed celestial objects

A sub-fraction of dark matter or new particles trapped inside celestial objects can significantly alter their macroscopic properties. Here, we investigate the new physics imprint on celestial objects by using a generic framework to solve the Tolman-Oppenheimer-Volkoff (TOV) equations for up to two fluids. We test the impact of populations of new particles on celestial objects, including the sensitivity to self-interaction sizes, new particle mass, and net population mass. Applying our setup to neutron stars and boson stars, we find rich phenomenology for a range of these parameters, including the creation of extended atmospheres. These atmospheres are detectable by their impact on the tidal love number, which can be measured at upcoming gravitational wave experiments such as Advanced LIGO, the Einstein Telescope, and LISA. We release our calculation framework as a publicly available code at this URL, allowing the TOV equations to be generically solved for arbitrary new physics models in novel and admixed celestial objects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Variance reduction techniques for Monte Carlo neutron noise simulations

The small fluctuations of the neutron flux caused by small perturbations of the macroscopic cross-sections take the name of neutron noise. Advanced Monte Carlo methods have been recently proposed in order to solve the neutron noise equations in the frequency domain, which allows establishing reference solutions to validate faster but approximate deterministic solvers. Due to the presence of particles carrying two statistical weights (for the real and imaginary components of the noise field), both of which may be positive or negative, Monte-Carlo simulations of neutron noise pose distinct challenges in terms of variance reduction. In this work we investigate two variance-reduction techniques, namely branchless collisions and weight cancellation, and probe their effectiveness for a benchmark con- figuration concerning the noise field induced by a pin with oscillating cross sections in a fuel assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Radiation-induced bowing of SiC/SiC composites under neutron flux gradients—integral experimental data for model validation

Here, the radiation-induced swelling of SiC and its composites, including strong dependencies on temperature and dose, can drive significant lateral bowing in the presence of temperature and/or dose gradients. In recent years, simulations have been performed to assess the extent of bowing in SiC composite light-water reactor (LWR) fuel cladding and boiling water reactor (BWR) channel boxes. However, to date, no integral experimental data exist to validate these models. This work provides the first experimental bowing evaluation of three ∼380 mm long SiC composite specimens irradiated under varying neutron dose gradients (∼50°C–60°C, 0.03–0.06 dpa): two tubes (∼9.8 mm diameter) and a miniature BWR channel box (∼30 mm square). The measured radiation-induced length swelling (∼0.3%–0.7% linear) was consistently 10%–21% higher than values obtained from 3D finite element structural analyses with inputs from 3D radiation transport calculations. This discrepancy could be at least partially explained by differences in dose rate (∼10 -8 dpa/s) compared to the literature data (∼10-6 dpa/s) used to establish the dose-to-swelling correlations in the model. Nevertheless, the modeled bowing magnitudes (<2 mm) obtained from finite element analyses and simple analytical equations were within the bounds of the experimental measurements for all specimens. With improved confidence in the ability to predict the structural response and measure the macroscopic deformations, future experiments will target transient bowing under neutron flux gradients at representative LWR temperatures and assess whether grid spacers can mitigate the tens of millimeters of bowing that would otherwise be expected in ∼4 m long LWR components.

bowing↗
DOE OSTI2025

Supersaturation, Nucleation, and Phase Separation of Mesoscopic Systems

Supersaturation, nucleation, and phase separation are ubiquitous phenomena of great interest in both science and industry. However, a unified, quantitative understanding of these phenomena has yet to be achieved for mesoscopic systems. Here, we present a set of general equations that determine the monomer saturation degree, the size distribution, and the free energy of mesoscopic systems, as well as their phase-transition conditions. These equations reveal that, under supersaturation, the largest cluster size (LCS) is an important state variable; the supersaturation degree decreases with the LCS, approaching unity in the macroscopic limit. We identify the critical supersaturation, at which the nuclei undergo the phase transition to form large crystals. Below this critical supersaturation, the nucleus size distribution is either a unimodal function or a monotonically decreasing function of size, depending on the system and temperature. We also predict the most probable nucleus size and the direction of spontaneous changes of the LCS. Our theory provides a unified, quantitative explanation of the nucleus-size-distribution across six different systems, including nanoparticles and biological condensates. This work serves as a general theoretical framework useful for understanding and designing nucleation and phase transitions of mesoscopic systems.

Kang, Jingyu↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Grey Galaxies’ as an endpoint of the Kerr-AdS superradiant instability

Kerr-AdS d+1 black holes for d ≥ 3 suffer from classical superradiant instabilities over a range of masses above extremality. We conjecture that these instabilities settle down into Grey Galaxies (GGs) — a new class of coarse-grained solutions to Einstein’s equations which we construct in d = 3. Grey Galaxies are made up of a black hole with critical angular velocity ω = 1 in the ‘centre’ of AdS, surrounded by a large flat disk of thermal bulk gas that revolves around the centre of AdS at the speed of light. The gas carries a finite fraction of the total energy, as its parametrically low energy density and large radius are inversely related. GGs exist at masses that extend all the way down to the unitarity bound. Their thermodynamics is that of a weakly interacting mix of Kerr-AdS black holes and the bulk gas. Their boundary stress tensor is the sum of a smooth ‘black hole’ contribution and a peaked gas contribution that is delta function localized around the equator of the boundary sphere in the large N limit. We also construct another class of solutions with the same charges; ‘Revolving Black Holes (RBHs)’. RBHs are macroscopically charged SO(d, 2) descendants of AdS-Kerr solutions, and consist of ω = 1 black holes revolving around the centre of AdS at a fixed radial location but in a quantum wave function in the angular directions. RBH solutions are marginally entropically subdominant to GG solutions and do not constitute the endpoint of the superradiant instability. Nonetheless, we argue that supersymmetric versions of these solutions have interesting implications for the spectrum of supersymmetric states in, e.g. N = 4 Yang-Mills theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Deducing neutron star equation of state from telescope spectra with machine-learning-derived likelihoods

The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to perform maximum a posteriori estimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas (Final Technical Report)

This document details the work accomplished during the period 04/15/2012-01/31/2020 for the collaborative project “Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas,” and serves as the final report on the University of California, San Diego part of this collaborative project for Department of Energy Office of Science - Fusion Energy Sciences Grant No. DE- SC0007783. This project focused on the development of analytic and computational tools, seeking to better understand the mutual interactions between small-scale temperature-gradient driven microturbulence and macroscopic magnetohydrodynamic (MHD) instabilities, particularly tearing modes, in magnetically confined plasmas relevant for the generation of fusion energy. This project advanced understanding through numerical simulations of different sets of model equations, corresponding to different microturbulence instabilities. Questions addressed included how the magnitudes, spatial structures, and scaling's of turbulent forcing's and drags impact magnetic island dynamics, as well as quantifying the “region of influence” of an island on the microturbulence dynamics. A key finding of the project is that the microturbulence fluctuations and associated cross-field particle and energy fluxes exhibit spatial asymmetries for sufficiently large islands. The results have improved our fundamental understanding of plasma dynamics, motivated new experimental investigations, and will help inform future fusion predictive modeling efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗