SEARCH · Search NASA

Results for “3D equilibrium”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Neoclassical toroidal viscosity torque prediction via deep learning

GPECnet is a densely connected neural network that has been trained on GPEC data, to predict the plasma stability, neoclassical toroidal viscosity (NTV) torque, and optimized 3D coil current distributions for desired NTV torque profiles. Using NTV torque, driven by non-axisymmetric field perturbations in a tokamak, can be vital in optimizing pedestal performance by controlling the rotation profile in both the core, to ensure tearing stability, and the edge, to avoid edge localized modes (ELMs). The generalized perturbed equilibrium code (GPEC) software package can be used to calculate the plasma stability to 3D perturbations and the NTV torque profile generated by applied 3D magnetic fields. These calculations, however, involve complex integrations over space and energy distributions, which takes time to compute. Initially, GPECnet has been trained solely on data representative of the quiescent H-mode (QH) scenario, in which neutral beams are often balanced and toroidal rotation is low across the plasma profile. Lastly, this work provides the foundation for active control of the rotation shear using a combination of beams and 3D fields for robust and high performance QH mode operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

A comparison of the predictive capabilities of several turbulence models using upwind and central-difference computer codes

Four turbulence models are described and evaluated for transonic flows using the upwind code CFL3D and the central-difference code TLNS3D. In particular, the effects of recent modifications to the half-equation model of Johnson-King are explored in detail, and different versions of the model are compared. This model can obtain good results for both two-dimensional (2D) and three-dimensional (3D) separated flows. The one-equation models of Baldwin-Barth and Spalart-Allmaras perform well for separated airfoil flows, but can predict the shock too far forward at the outboard stations of a separated wing. The equilibrium model of Baldwin-Lomax predicts the shock location too far aft for both 2D and 3D separated flows, as expected. In general, all models perform well for attached or mildly separated flows.

Rumsey, Christopher L.↗

NASA NTRS1993

A Comparison of the Predictive Capabilities of Several Turbulence Models using Upwind and Central-Difference Computer Codes

Rumsey, Christopher L.↗

DOE OSTI2021

Direct measurements of the 3D plasma velocity in single-helical-axis RFP plasmas

The first local velocity measurements of helical equilibrium plasmas in the Reversed Field Pinch (RFP) Single Helical Axis (SHAx) state using a Charge Exchange Recombination Spectroscopy (CHERS) diagnostic are presented. Measurements show strong axisymmetric and non-axisymmetric flow, with n ¼ 5 components of flow related to the ðm; nÞ¼ð1; 5Þ dominant magnetic mode on the order of the axisymmetric flow in certain regions of the plasma, as well as significant n > 5 flow. Flow measurements are compared with NIMROD simulations of visco-resistive, single-fluid MHD in toroidal and cylindrical geometries with limited axial periodicity. Both measurements and the simulation with toroidal geometry show stronger inboard flows relative to the outboard flows, which is attributed to the toroidal geometry of the device. In the experiment, the n ¼ 5 component of flow is phase shifted from the reconnection-like flow pattern observed in the single-fluid simulations, possibly due to decoupling of the ion and electron fluids over much of the plasma. Finally, the strength of the helical angular flow shear relative to the critical shear necessary to disrupt nonlinear coupling between tearing modes is calculated around the helical magnetic axis. The shear in the measured flow is on the order of the theoretical critical threshold needed to nonlinearly decouple modes, but the measurement uncertainty in the gradient of the flow is large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Multiscale dynamics of charging and plating in graphite electrodes coupling operando microscopy and phase-field modelling

The phase separation dynamics in graphitic anodes significantly affects lithium plating propensity, which is the major degradation mechanism that impairs the safety and fast charge capabilities of automotive lithium-ion batteries. In this study, we present comprehensive investigation employing operando high-resolution optical microscopy combined with non-equilibrium thermodynamics implemented in a multi-dimensional (1D+1D to 3D) phase-field modeling framework to reveal the rate-dependent spatial dynamics of phase separation and plating in graphite electrodes. Here we visualize and provide mechanistic understanding of the multistage phase separation, plating, inter/intra-particle lithium exchange and plated lithium back-intercalation phenomena. A strong dependence of intra-particle lithiation heterogeneity on the particle size, shape, orientation, surface condition and C-rate at the particle level is observed, which leads to early onset of plating spatially resolved by a 3D image-based phase-field model. Moreover, we highlight the distinct relaxation processes at different state-of-charges (SOCs), wherein thermodynamically unstable graphite particles undergo a drastic intra-particle lithium redistribution and inter-particle lithium exchange at intermediate SOCs, whereas the electrode equilibrates much slower at low and high SOCs. These physics-based insights into the distinct SOC-dependent relaxation efficiency provide new perspective towards developing advanced fast charge protocols to suppress plating and shorten the constant voltage regime.

25 ENERGY STORAGE↗

DOE OSTI2025

Equilibrium expectations for non-Gaussian fluctuations near a QCD critical point

With the highly anticipated results from the Beam Energy Scan II program at RHIC being recently revealed, an understanding of particle-number fluctuations and their significance as a potential signature of a possible QCD critical point is crucial. Early works that embarked on this endeavor sought to estimate the fluctuations due to the presence of a critical point assuming they stay in equilibrium. From these results came the proposal to focus efforts on higher, non-Gaussian, moments of the event-by-event distributions, in particular of the number of protons. These non-Gaussian moments are especially sensitive to critical fluctuations, as their magnitudes are proportional to high powers of the critical correlation length. As the equation of state provides key input for hydrodynamical simulations of heavy-ion collisions, we estimate equilibrium fluctuations from the BEST equation of state (EoS) that includes critical features from the 3D Ising Model. In particular, the proton factorial cumulants and their dependence on non-universal mapping parameters is investigated within the BEST EoS. Furthermore, the correlation length, as a central quantity for the assessment of fluctuations in the vicinity of a critical point, is also calculated in a consistent manner with the scaling equation of state. An understanding of the equilibrium estimates of proton factorial cumulants will be useful for further comparison to estimates of out-of-equilibrium fluctuations in order to determine the magnitude of the observable fluctuations to be expected in heavyion collision experiments, in which the time spent near a critical point is short.

Karthein, Jamie M. [Massachusetts Institute of Tec↗

NASA NTRS2014

Thermal Non-equilibrium Consistent with Widespread Cooling

Time correlation analysis has been used to show widespread cooling in the solar corona; this cooling has been interpreted as a result of impulsive (nanoflare) heating. In this work, we investigate wide-spread cooling using a 3D model for a solar active region which has been heated with highly stratified heating. This type of heating drives thermal non-equilibrium solutions, meaning that though the heating is effectively steady, the density and temperature in the solution are not. We simulate the expected observations in narrowband EUV images and apply the time correlation analysis. We find that the results of this analysis are qualitatively similar to the observed data. We discuss additional diagnostics that may be applied to differentiate between these two heating scenarios.

Winebarger, A.↗

NASA NTRS2014

Thermal Non-equilibrium Consistent with Widespread Cooling

Winebarger, A.↗

DOE OSTI2020

In-Operando Study of Shape Retention and Microstructure Development in a Hydrolyzing Sol–Gel Ink during 3D-Printing

3D printing of amorphous and crystalline ceramics is of paramount importance for the fabrication of a wide range of devices with applications across different technology fields. Printed ceramics are remarkably enabled by the sol-gel synthesis method in conjunction with continuous filament direct ink writing. During printing, multiple processes contribute to the evolution of inks including shape retention, chemical conversion, solidification, and microstructure formation. Traditionally, depending on the ink composition and printing environment, several mechanisms have been associated with the shape retention and solidification of 3D printed structures: gelation, rapid solvent evaporation, energy-driven phase transformation, and chemical-driven phase transformation. Understanding the fundamental differences between these mechanisms becomes key since they strongly influence the spatiotemporal evolution of the materials, as the out-of-equilibrium processes inherent to the extrusion, relaxation, and solidification of printed materials have significant effects on the materials properties. Herein, we investigate the shape retention mechanism and the hydrolysis-induced material conversion and microstructure formation during the 3D printing of a water reactive sol-gel ink that transforms into titanium dioxide-based ceramic. This study aims at identifying characteristic mechanisms associated with the material transformation, establishing connections between the microstructure development and the timescales associated with solidification under operando 3D-printing conditions. The investigation of this material's out-of-equilibrium pathways under processing conditions is enabled by time-resolved coherent X-ray scattering, providing simultaneous access to temporospatially resolved microstructural and dynamics information. Furthermore, we explore X-ray speckle tracking as a tool to resolve deformations of the microstructure in a printed filament associated with the deposition of consecutive filaments. Through this work, we aim at providing a fundamental understanding of the relationships behind these transformative processes in 3D printing and their timescales as the basis for achieving unprecedented control over printed materials microstructure.

36 MATERIALS SCIENCE↗

NASA NTRS1989

Size distribution of dust grains: A problem of self-similarity

Distribution functions describing the results of natural processes frequently show the shape of power laws, e.g., mass functions of stars and molecular clouds, velocity spectrum of turbulence, size distributions of asteroids, micrometeorites and also interstellar dust grains. It is an open question whether this behavior is a result simply coming about by the chosen mathematical representation of the observational data or reflects a deep-seated principle of nature. The authors suppose the latter being the case. Using a dust model consisting of silicate and graphite grains Mathis et al. (1977) showed that the interstellar extinction curve can be represented by taking a grain radii distribution of power law type n(a) varies as a(exp -p) with 3.3 less than or equal to p less than or equal to 3.6 (example 1) as a basis. A different approach to understanding power laws like that in example 1 becomes possible by the theory of self-similar processes (scale invariance). The beta model of turbulence (Frisch et al., 1978) leads in an elementary way to the concept of the self-similarity dimension D, a special case of Mandelbrot's (1977) fractal dimension. In the frame of this beta model, it is supposed that on each stage of a cascade the system decays to N clumps and that only the portion beta N remains active further on. An important feature of this model is that the active eddies become less and less space-filling. In the following, the authors assume that grain-grain collisions are such a scale-invarient process and that the remaining grains are the inactive (frozen) clumps of the cascade. In this way, a size distribution n(a) da varies as a(exp -(D+1))da (example 2) results. It seems to be highly probable that the power law character of the size distribution of interstellar dust grains is the result of a self-similarity process. We can, however, not exclude that the process leading to the interstellar grain size distribution is not fragmentation at all. It could be, e.g., diffusion-limited growth discussed by Sander (1986), who applied the theory of fractal geometry to the classification of non-equilibrium growth processes. He received D=2.4 for diffusion-limited aggregation in 3d-space.

Henning, TH.↗

DOE OSTI2020

Magnetic reconnection in 3D fusion devices: non-linear reduced equations and linear current-driven instabilities

Abstract Magnetic reconnection in 3D fusion devices is investigated. With the use of Boozer co-ordinates, we reduce the non-linear resistive magnetohydrodynamic equations in the limit of large aspect ratio and finite pressure fluctuations, to obtain a set of non-linear equations suitable for magnetic reconnection studies in stellarators. Magnetic flux unfreezing due to a finite electron mass is also considered. Equations that govern the linear regime and some of their general properties are given. We emphasise the role of magnetic geometry and identify how some aspects of stellarator optimisation could have an impact on reconnecting instabilities, in particular by exacerbating those enabled by electron inertia. The effect of 3D coupling on the linear reconnection rates and the mode structure is quantitatively addressed in the case in which the equilibrium rotational transform has one specific resonant location for which one mode can reconnect while coupled to an arbitrary number of non-resonant harmonics. The full problem is rigorously reduced to an equivalent cylindrical one, by introducing some geometrically modified plasma inertial and dissipative scales. The 3D scalings for the growth rates of reconnection instabilities and their destabilisation criteria are given.

Physics↗

DOE OSTI2025

Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂⁡(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1991

Comparison of hypersonic experiments and PNS predictions. I - Aerothermodynamics. II - Aerodynamics

A 3D parabolized Navier-Stokes (PNS) technique is examined in terms of applications to equilibrium-air and perfect-gas hypersonic flows at different Mach numbers and angles of attack. The method renders sublayer approximation unnecessary, and treats axis-normal coupling effects first followed by a shock solution and a corrector step. Predictions of hypersonic flowfields over several blunt conical configurations are found to agree with established wind-tunnel and flight data, although those for the Dhawan-Narasimha (1958) transition model vary from the experimental data. The PNS technique is also applied to predictions of wall-pressure distributions and force and moment data, and the results agree with experimental values. The PNS method allows the calculation of a wide variety of configurations with only a slight increase over traditional computing times.

Bhutta, Bilal A.↗

DOE OSTI2024

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗

DOE OSTI2023

Loss of energetic ions due to n = 1 internal kink instability in HL-2M

Effects of three-dimensional (3D) perturbations due to an unstable n = 1 (n is the toroidal mode number) internal kink (IK) on the energetic particles (EPs) are systematically investigated for the HL-2M tokamak, utilizing the MARS-F/K code and a recently developed test particle tracing module. A high-beta sawteething HL-2M scenario, simulated by the TRANSP code, is chosen for this study. In general, the 3D perturbation associated with an unstable IK is found to affect the EP drift orbit, confinement, and loss in HL-2M. The instability-induced EP loss fraction is found to be typically less than 10%, without counting for the prompt orbit loss associated with the 2D equilibrium field for counter-current particles. The latter reaches about 16% in HL-2M. For co-current EPs, a 100 G 3D magnetic field (inside the plasma) due to the IK does not induce any EP loss assuming a static perturbation. A sawtooth-like time-varying perturbation field, with the peak amplitude reaching 1000 G, can however produce about 30% loss for the co-current EPs in HL-2M. The majority of lost EPs tend to strike the lower divertor region, with a small fraction of particles striking the low-field side mid-plane region of the limiting surface.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Development of New Reactor Core Configuration for Power Uprate - Fuel Reload & Heat Processing Analyses, Core Design, System Safety Assessments, and Fuel Performance Analyses

With the passage of the Infrastructure Investment and Jobs Act in 2021 and the Inflation Reduction Act (IRA) in 2022, the United States stands at a critical juncture for the future of nuclear power. These landmark policies provide significant support for clean energy initiatives, positioning nuclear power as a key component of the nation’s strategy to reduce carbon emissions and achieve energy security. This growing emphasis on nuclear energy is driven by the need for reliable, low-carbon power sources as the country transitions away from fossil fuels. Federal policy, along with increasing state-level support, is encouraging investment in nuclear technology advancements to meet these demands. Building new nuclear power plants (NPPs), however, presents significant challenges due to high costs and long construction timelines. As a result, increasing the power output of existing NPPs through power uprates has emerged as a more feasible and cost-effective strategy. One key area of advancement is the development of accident-tolerant fuel (ATF), such as chromium-coated zirconium alloy cladding, which offers enhanced material performance, enabling power uprates in light water reactors (LWRs). Given the growing demand for nuclear energy fueled by federal policies and state initiatives, it is essential to evaluate the feasibility and benefits of significant power uprates in existing pressurized water reactors (PWRs) using advanced fuel technologies. The introduction of ATF concepts opens new opportunities for safely and economically achieving these power increases. Assessing whether these innovations can support substantial power uprates while maintaining operational safety is crucial to maximizing the potential of the nation’s existing nuclear infrastructure. This project aims to explore how power uprates can be achieved by boosting reactor thermal power output and optimizing reactor core design, while ensuring the safety and economic viability of NPPs. Specifically, it will focus on demonstrating the technical and economic feasibility of power uprates in a PWR using low 5-10% enrichment uranium (LEU+) high burnup (HBU) fuel combined with ATF concepts. In fiscal year 2024 (FY24), the research and development focus on building foundational models and conducting multi-physics performance and safety analyses to support the power uprate. The findings of the study would be shared through LWRS Seasonal Meetings, conferences and workshops with utility companies and researchers. These also serve as a basis for further study of fuel reloading optimization with ATF claddings.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Kinetic equilibrium reconstructions of plasmas in the MAST database and preparation for reconstruction of the first plasmas in MAST upgrade

Reconstructions of plasma equilibria using magnetic sensors were routine during operation of the Mega Ampere Spherical Tokamak (MAST) device, but reconstructions using kinetic profiles were not. These are necessary for stability and disruption analysis of the MAST database, as well as for operation in the upgrade to the device, MAST-U. The three-dimensional (3D) code VALEN is used to determine eddy currents in the 3D vessel structures for vacuum coil test shots, which are then mapped to effective resistances in the two-dimensional vessel groupings in the EFIT equilibrium reconstruction code to be used in conjunction with nearby loop voltage measurements for estimated currents in the structures during reconstruction. Kinetic equilibrium reconstructions with EFIT, using all available magnetic sensors as well as Thomson scattering measurements of electron temperature and density, charge exchange recombination spectroscopy measurements of ion temperature, and internal magnetic field pitch angle measurements from a motional Stark effect (MSE) diagnostic are performed for a large database of MAST discharges. Excellent convergence errors are obtained for the portions of the discharges where the stored energy was not too low, and it is found that reconstructions performed with temperature and density measurements but without MSE data usually already match the pitch angle measurements well. A database of 275 kinetic equilibria is used to test the ideal MHD stability calculation capability for MAST. In conclusion, the necessary changes to conducting structure in VALEN, and diagnostic setup in EFIT have been completed for the upgrade from MAST to MAST-U, enabling kinetic reconstructions to commence from the first plasma discharges of the upgraded device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

A three-dimensional upwind PNS code for chemically reacting scramjet flowfields

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional (3D) chemically reacting flow in scramjet (supersonic combustion ramjet) engines. The code is a modification of the 3D upwind PNS (UPS) airflow code which has been extended in the present study to permit internal flow calculations with hydrogen-air chemistry. With these additions, the new code has the capability of computing aerodynamic and propulsive flowfields simultaneously. The algorithm solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that has been modified to account for 'real gas' effects. The fluid medium is assumed to be a chemically reacting mixture of thermally perfect (but calorically imperfect) gases in thermal equilibrium. The new code has been applied to two test cases. These include the Burrows-Kurkov supersonic combustion experiment and a generic 3D scramjet flowfield. The computed results compare favorably with the available experimental data.

Wadawadigi, Ganesh↗