Search NASA⌕ Search

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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Impact of 3D Structure on Magnetic Reconnection

Results from 2.5D and 3D studies of the onset and development of the tearing instability are presented, using high fidelity resistive MHD simulations. A limited parameter study of the strength of the reconnecting field (or shear angle) was performed. An initially simple 1D equilibrium was used, consisting of a modified force-free current sheet, with periodic boundary conditions in all directions. In all cases, the linear and nonlinear evolution led to a primary current sheet between two large flux ropes. The global reconnection rate during this later stage was analyzed in all simulations. It was found that in 2.5D the primary current sheet fragmented owing to plasmoids, and as expected, the global reconnection rate, calculated using multiple methods, increases with the strength of the reconnecting field owing to a stronger Alfvén speed. In 3D, the presence of interacting oblique modes of the tearing instability complicates the simple 2.5D picture, entangling the magnetic field of the inflow and introducing a negative effect on the reconnection rate. The two competing effects of stronger Alfvén speed and entangling, which both increase with the strength of the reconnecting field, resulted in a decrease in the reconnection rate with increasing reconnecting field. For all simulations, the 3D rates were less than in 2.5D but suggest that as one goes to weak reconnecting field (or strong guide field) the system becomes more 2.5D-like and the 2.5D and 3D rates converge. These results have relevance to situations like nanoflare heating and flare current sheets in the corona.

Lars K S Daldorff↗
DOE OSTI2021

Quasisymmetric Optimization of Nonaxisymmetry in Tokamaks

Predictive 3D optimization reveals a novel approach to modify a nonaxisymmetric magnetic perturbation to be entirely harmless for tokamaks, by essentially restoring quasisymmetry in perturbed particle orbits as much as possible. Such a quasisymmetric magnetic perturbation (QSMP) has been designed and successfully tested in the KSTAR and DIII-D tokamaks, demonstrating no performance degradation despite the large overall amplitudes of nonaxisymmetric fields and strong response otherwise expected in the tested plasmas. The results indicate that a quasisymmetric optimization is a robust path of error field correction across the resonant and nonresonant field spectrum in a tokamak, leveraging the prevailing concept of quasisymmetry for general 3D plasma confinement systems such as stellarators. Finally, the optimization becomes, in fact, a simple eigenvalue problem to the so-called torque response matrices if a perturbed equilibrium is calculated consistent with nonaxisymmetric neoclassical transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Versatile Approach to Stabilize Liquid–Liquid Interfaces using Surfactant Self-Assembly

Stabilizing liquid–liquid interfaces, whether between miscible or immiscible liquids, is crucial for a wide range of applications, including energy storage, microreactors, and biomimetic structures. Here, in this study, a versatile approach for stabilizing the water-oil interface is presented using the morphological transitions that occur during the self-assembly of anionic, cationic, and nonionic surfactants mixed with fatty acid oils. The morphological transitions underlying this approach are characterized and extensively studied through small-angle X-ray scattering (SAXS), rheometry, and microscopy techniques. Dissipative particle dynamics (DPD) as a simulation tool is adopted to investigate these morphological transitions both in the equilibrium ternary system as well as in the dynamic condition of the water-oil interface. Such a versatile strategy holds promise for enhancing applications such as liquid-in-liquid 3D printing. Moreover, it has the potential to revolutionize a wide range of fields where stabilizing liquid–liquid interfaces not only offers unprecedented opportunities for fine-tuning nanostructural morphologies but also imparts interesting practical features to the resulting liquid shapes. These features include perfusion capabilities, self-healing, and porosity, which could have significant implications for various industries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Integrated three-dimensional characterization of reactive phase formation and coarsening during isothermal annealing of metastable Zn–3Mg–4Al eutectic

Microstructural instabilities associated with a metastable lamellar eutectic when exposed to an elevated temperature were analyzed in a Zn-3Mg-4Al alloy as a model system. X-ray diffraction and scanning electron microscopy showed that annealing at a temperature 43 K below the equilibrium eutectic temperature converted the metastable ternary MgZn 2 /η-Zn/β-ZnAl eutectic to a stable mixture of Mg 2 Zn 11 /η-Zn/α-Al. 3D EBSD suggested that supersaturation of η-Zn with Mg and Al eased difficulties in nucleation of Mg 2 Zn 11 and α-Al at the MgZn 2 /η-Zn phase interfaces. Quasi-in-situ X-ray nanotomography revealed that a reactive diffusion process prevailed in an early stage of annealing leading to rapid consumption of MgZn 2 for growth of Mg 2 Zn 11 and concomitant thinning and pinch-off of η-Zn lamellae. The reactive diffusion process also brought about significant and quantifiable changes in the topological characteristics of the η-Zn phase at remarkably shorter time scales compared to Rayleigh instability and Ostwald ripening mechanisms. Overall, the present results provide fresh insights on the short-term annealing effects on metastable eutectics, with relevance to structural evolution in additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the frequency bifurcations of the MHD startup modes in NSTX

The observed bifurcations of the low frequency (<50 kHz) and low toroidal periodicity (n < 5) magnetohydrodynamic (MHD) activity often present in the initial part of the National Spherical Tokamak Experiment (NSTX) discharges can be explained by the evolution of the radial profile of the safety factor (q=rB φ /RB θ ) crossing multiple rational surfaces in the core. Important performance limiting instability mechanisms in the NSTX spherical tokamak are often linked to low frequency and low-n MHD activity. They are quite common in long-pulse NSTX plasmas. They can be present at the beginning of the plasma current flat-top, at the end of the discharge or during the whole duration, and they have been observed to deleteriously impact performance over a wide range of q 95 . An interesting feature observed in some NSTX discharges is the presence of a bifurcation in the frequency of the low n modes, as low as n = 1, that have frequencies comparable to the plasma core rotation divided by n. Equilibrium reconstructions constrained by magnetic diagnostics data and motional stark effect pitch angle radial profiles suggest that the observed bifurcations are linked to a fast evolving minimum value of q. Finally, 3D non-linear resistive MHD simulations show that these modes are ideal and exist as non-resonant before the correspondent rational surface enters the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A versatile approach to stabilize liquid-liquid interfaces using surfactant self-assembly

Stabilizing liquid–liquid interfaces, whether between miscible or immiscible liquids, is crucial for a wide range of applications, including energy storage, microreactors, and biomimetic structures. In this study, a versatile approach for stabilizing the water-oil interface is presented using the morphological transitions that occur during the self-assembly of anionic, cationic, and nonionic surfactants mixed with fatty acid oils. The morphological transitions underlying this approach are characterized and extensively studied through small-angle X-ray scattering (SAXS), rheometry, and microscopy techniques. Dissipative particle dynamics (DPD) as a simulation tool is adopted to investigate these morphological transitions both in the equilibrium ternary system as well as in the dynamic condition of the water-oil interface. Such a versatile strategy holds promise for enhancing applications such as liquid-in-liquid 3D printing. Moreover, it has the potential to revolutionize a wide range of fields where stabilizing liquid–liquid interfaces not only offers unprecedented opportunities for fine-tuning nanostructural morphologies but also imparts interesting practical features to the resulting liquid shapes. These features include perfusion capabilities, self-healing, and porosity, which could have significant implications for various industries.

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

Systematic analysis of magnetic equilibrium reconstruction with eddy currents on LTX- β

Magnetic equilibrium reconstruction in the Lithium Tokamak Experiment-Beta (LTX-β) is complicated by strong eddy currents and toroidal asymmetries arising from its segmented, close-fitting conducting shell. In the earlier experiment LTX, these three-dimensional (3D) conductor effects significantly distorted diagnostic signals and challenged conventional axisymmetric reconstruction methods. In this work, we demonstrate that accurate plasma equilibria can be recovered on LTX-β by incorporating a small number of dominant eddy current modes derived from realistic conductor models. Using the open-source $\tt{TokaMaker}$ Grad–Shafranov solver, we reconstruct equilibria across a systematically selected set of LTX-β discharges and validate them against the legacy $\tt{PSI-Tri}$ hybrid 2D-3D code. Our results show that fully 2D $\tt{TokaMaker}$ reconstructions achieve significantly improved agreement with flux loop and Mirnov probe measurements, reducing total chi-squared fitting errors, especially during startup. Finally, we develop new hybrid 2D-3D $\tt{TokaMaker}$ reconstructions by integrating the $\tt{ThinCurr}$ 3D eddy current model, which we optimize for significant further reduction of chi-squared in most scenarios. These findings underscore the importance of realistic wall-current modeling in short-pulse tokamaks, and establish a physics-based reconstruction framework that is extensible to devices with complex passive structures and 3D wall interactions.

Grad–Shafranov solver↗
DOE OSTI2022

Non-symmetric pinning of topological defects in living liquid crystals

Topological defects, such as vortices and disclinations, play a crucial role in spatiotemporal organization of equilibrium and non-equilibrium systems. The defect immobilization or pinning is a formidable challenge in the context of the out-of-equilibrium system, like a living liquid crystal, a suspension of swimming bacteria in lyotropic liquid crystal. Here we control the emerged topological defects in a living liquid crystal by arrays of 3D-printed microscopic obstacles (pillars). Our studies show that while –1/2 defects may be easily immobilized by the pillars, +1/2 defects remain motile. Due to attraction between oppositely charged defects, positive defects remain in the vicinity of pinned negative defects, and the diffusivity of positive defects is significantly reduced. Experimental findings are rationalized by computational modeling of living liquid crystals. Our results provide insight into the engineering of active systems via targeted immobilization of topological defects.

42 ENGINEERING↗
NASA NTRS2017

3D Material Response Analysis of PICA Pyrolysis Experiments

Primarily interested in improving ablation modeling for use in inverse reconstruction of flight environments on ablative heat shields. Ablation model is essentially a component of the heat flux sensor, so model uncertainties lead to measurement uncertainties. Non-equilibrium processes have been known to be significant in low density ablators for a long time, but increased accuracy requirements of the reconstruction process necessitates incorporating this physical effect. Attempting to develop a pyrolysis model for implementation in material response based on the PICA data produced by Bessire and Minton. Pyrolysis gas species molar yields as a function of temperature and heating rate. Several problems encountered while trying to fit Arrhenius models to the data led to further investigation of the experimental setup.

Oliver, Brandon A.↗
DOE OSTI2023

Measuring Multicomponent Adsorption of Tracer Gases on Natural Zeolites

A natural clinoptilolite sample near the Nevada National Security Site was obtained to study adsorption and retardation on gas transport. Of interest is understanding the competition for adsorption sites that may reduce tracer gas adsorption relative to single-component measurements, which may be affected by the multi-scale pore structure of clinoptilolite. Clinoptilolite has three distinct domains of pore size distributions ranging from nanometers to micrometers: micropores with 0.4–0.7 nm diameters, measured on powders by CO2 adsorption at 273 K, representing the zeolite cages; mesopores with 4–200 nm diameters, observed using liquid nitrogen adsorption at 77 K; and macropores with 300–1000 nm diameters, measured by mercury injection on rock chips (~ 100 mesh), likely representing the microfractures. These pore size distributions are consistent with X-ray computed tomography (CT) and focused ion beam scanning electron microscope (FIB-SEM) images, which are used to construct the three-dimensional (3D) pore network to be used in future gas transport modeling. To quantify tracer gas adsorption in this multi-scale pore structure and multicomponent gas species environment, natural zeolite samples initially in equilibrium in air were exposed to a mixture of tracer gases. As the tracer gases diffuse and adsorb in the sample, the remaining tracer gases outside the sample fractionate. Using a quadrupole mass spectrometer to quantify this fractionation, the degree of adsorption of tracer gases in the multicomponent gas environment and multi-scale pore structure is assessed. The major finding is that Kr reaches equilibrium much faster than Xe in the presence of ambient air, which leads to more Kr uptake than Xe over limited exposure periods. When the clinoptilolite chips were exposed to humid air, the adsorption capability decreases significantly for both Xe and Kr with relative humidity (RH) as low as 3%. Both Xe and Kr reaches equilibrium faster at higher RH. The different, unexpected, adsorption behavior for Xe and Kr is due to their kinetic diameters similar to the micropores in clinoptilolite which makes it harder for Xe to access compared to Kr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1973

Theoretical study of H2/+/ spectroscopic properties. II, III

Description of the theoretical spectroscopic properties of the 2p pi/sub u/ and 3d sigma/sub g/ excited states of the H2/+/ hydrogen molecular ion. Numerical integration of the Schrodinger equation is used to determine vibration-rotation eigenvalues. Dunham power series expansions are used to determine the equilibrium separation, potential coefficients, and spectroscopic constants. The eigenvalues are used to determine delta-G, Bv, Dv, and Hv.

Beckel, C. L.↗
NASA NTRS1980

Atomic data for S IV and solar observations of the 3x/2/3p /2/P-3s3p/2/ /4/P multiplet

High resolution observations of intersystem lines of S IV near 1400 A are available from Skylab. These lines are potentially useful as density diagnostics for the solar atmosphere. Energy levels, transition probabilities and collision strengths have, therefore, been calculated for S IV, including the configurations 3x(2)3p, 3s3p(2), and 3s(2)3d. Line intensities and level populations have been calculated as a function of electron density. The calculated population of the 3s3p(2) (4)P(5/2) level is found to reach a pseudo-Boltzmann equilibrium at a density which is four times higher than is inferred from solar spectra and level population calculations of lighter ions such as O IV.

Bhatia, A. K.↗
NASA NTRS1983

Conical similarity of shock/boundary layer interactions generated by swept fins

A parametric experimental study has been made of the class of 3D shock wave/turbulent boundary layer interactions generated by swept-leading-edge fins. The fin sweepback angles ranged from 0 to 65 deg at angles of attack of 5, 9, and 15 deg. Two equilibrium 2D turbulent boundary layers with a free-stream Mach number of 2.95 and a Reynolds number of 6.3 x 10 to the 7th/m were used as incoming flow conditions. All the resulting interactions were found to possess conical symmetry of surface pressures and skin friction lines beyond an initial inception zone. Further, these interactions revealed a simple similarity based on inviscid shock strength irrespective of fin sweepback or angle of attack.

Lu, F. K.↗
DOE OSTI2021

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Vidyut3d: A Non-Equilibrium Plasma Modeling Tool [SWR-24-101]

Vidyut3d is a massively-parallel plasma-fluid solver for low-temperature plasmas (LTPs) that supports both local field (LFA) and local mean energy (LMEA) approximations, as well as complex gas and surface-phase chemistry. The solver supports 2D and 3D domains, and uses AMReX's adaptive mesh refinement capabilities to increase the grid resolution around complex structures (e.g. streamer heads and sheaths) while maintaining a tractable problem size. Vidyut specializes in simulating various types of gas-phase discharges, as well as plasma/surface interactions and surface chemistry (e.g. for plasma-mediated catalysis applications). The solver also supports hybrid CPU/GPU parallelization strategies, and has demonstrated excellent scaling on various HPC architectures for problem sizes consisting of O(100 M) control volumes.

Sitaraman, Hariswaran↗
NASA NTRS2016

The Structure of the Local Hot Bubble

Diffuse X-rays from the Local Galaxy (DXL) is a sounding rocket mission designed to quantify and characterize the contribution of Solar Wind Charge eXchange (SWCX) to the Diffuse X-ray Background and study the properties of the Local Hot Bubble (LHB). Based on the results from the DXL mission, we quantified and removed the contribution of SWCX to the diffuse X-ray background measured by the ROSAT All Sky Survey. The cleaned maps were used to investigate the physical properties of the LHB. Assuming thermal ionization equilibrium, we measured a highly uniform temperature distributed around kT = 0.097 keV +/- 0.013 keV (FWHM) +/- 0.006 keV(systematic). We also generated a thermal emission measure map and used it to characterize the three-dimensional (3D) structure of the LHB, which we found to be in good agreement with the structure of the local cavity measured from dust and gas.

X-rays↗