Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dimensional changes”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

The effect of gamma ray irradiation on few layered MoSe2: A material for nuclear and space applications

In recent years, emerging two-dimensional (2D) materials, such as molybdenum diselenide (MoSe2), have been at the center of attention for many researchers. This is due to their unique and fascinating physicochemical properties that make them attractive in space and defense applications that include shielding harsh irradiation environments. In this study, we examined the effects of gamma (γ) rays at various doses on the structural, chemical, and optical properties of MoSe2 layers. After the samples were exposed to intense gamma radiation (from a 60Co source) with various exposure times to vary the total accumulated dosage (up to 100 kGy), Raman and photoluminescence spectroscopies were used to study and probe radiation-induced changes to the samples. When compared to pristine materials, very few changes in optical properties were typically observed, indicating good robustness with little sensitivity, even at relatively high doses of gamma radiation. The imaging using scanning electron microscopy revealed a number of nano-hillocks that were connected to substrate alterations. X-ray photoelectron spectroscopies revealed that Mo’s binding energies remained the same, but Se’s binding energies blueshifted. We associated this shift with the decrease in Se vacancies that occurred after irradiation as a result of Mo atoms creating adatoms next to Se atoms. When compared to pristine materials, very few changes in optical, chemical, and structural properties were typically observed. These findings highlight the inherent resilience of MoSe2 in hostile radioactive conditions, which spurs additional research into their optical, electrical, and structural characteristics as well as exploration for potential space, energy, and defense applications.

Materials Science↗
DOE OSTI2021

The Co-Optimization of Sustainable Aviation Fuel: Cost, Emissions, and Performance

The combustion of petroleum-based fuels contributes to increases in atmospheric CO2, contributing to climate change. As other sectors electrify, current battery technology makes this impractical for the aviation industry. The fastest pathway to reducing the carbon contributions of aviation is to create low or no net carbon-emitting drop-in petroleum alternatives. This research explores 3 potential sorghum derived jet fuel molecules from 4 different production routes, hydro-processed esters and fatty acids (HEFA), and Jet A, to identify how blends including these fuels could offer improved performance (MJ/kg, MJ/L) and improved emissions (gCO2/MJ) while minimizing costs. Different applications will value each target metric differently. This research utilizes the Jet Fuel Blend Optimizer (JudO) to create a 4-dimensional Pareto front across potential solutions. The molecules considered will not offer reduced cost compared to conventional jet fuel (Jet A), but they can be utilized to create blends with higher specific energy, greater energy density, and reduced greenhouse gas emissions. Considerations including carbon credits, the sale of byproducts, and the valuation of improved performance will make the proposed molecules more commercially viable. At the lowest region of GHG solutions, coupled with equivalent LCFS approximations, JudO determined blends that could achieve as high as 69% overall carbon reductions at a premium of $0.34/L.

Feldhausen, John J.↗
DOE OSTI2023

Determining the biogeochemical transformations of organic matter composition in rivers using molecular signatures

Inland waters are hotspots for biogeochemical activity, but the environmental and biological factors that govern the transformation of organic matter (OM) flowing through them are still poorly constrained. Here we evaluate data from a crowdsourced sampling campaign led by the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to investigate broad continental-scale trends in OM composition compared to localized events that influence biogeochemical transformations. Samples from two different OM compartments, sediments and surface water, were collected from 97 streams throughout the Northern Hemisphere and analyzed to identify differences in biogeochemical processes involved in OM transformations. By using dimensional reduction techniques, we identified that putative biogeochemical transformations and microbial respiration rates vary across sediment and surface water along river continua independent of latitude (18°N–68°N). In contrast, we reveal small- and large-scale patterns in OM composition related to local (sediment vs. water column) and reach (stream order, latitude) characteristics. These patterns lay the foundation to modeling the linkage between ecological processes and biogeochemical signals. We further showed how spatial, physical, and biogeochemical factors influence the reactivity of the two OM pools in local reaches yet find emergent broad-scale patterns between OM concentrations and stream order. OM processing will likely change as hydrologic flow regimes shift and vertical mixing occurs on different spatial and temporal scales. As our planet continues to warm and the timing and magnitude of surface and subsurface flows shift, understanding changes in OM cycling across hydrologic systems is critical, given the unknown broad-scale responses and consequences for riverine OM.

rivers↗
DOE OSTI2023

Atomic motifs govern the decoration of grain boundaries by interstitial solutes

Abstract Grain boundaries, the two-dimensional defects between differently oriented crystals, tend to preferentially attract solutes for segregation. Solute segregation has a significant effect on the mechanical and transport properties of materials. At the atomic level, however, the interplay of structure and composition of grain boundaries remains elusive, especially with respect to light interstitial solutes like B and C. Here, we use Fe alloyed with B and C to exploit the strong interdependence of interface structure and chemistry via charge-density imaging and atom probe tomography methods. Direct imaging and quantifying of light interstitial solutes at grain boundaries provide insight into decoration tendencies governed by atomic motifs. We find that even a change in the inclination of the grain boundary plane with identical misorientation impacts grain boundary composition and atomic arrangement. Thus, it is the smallest structural hierarchical level, the atomic motifs, that controls the most important chemical properties of the grain boundaries. This insight not only closes a missing link between the structure and chemical composition of such defects but also enables the targeted design and passivation of the chemical state of grain boundaries to free them from their role as entry gates for corrosion, hydrogen embrittlement, or mechanical failure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Light-activated microtubule-based two-dimensional active nematic

We assess the ability of two light responsive kinesin motor clusters to drive dynamics of microtubule-based active nematics: opto-K401, a processive motor, and opto-K365, a non-processive motor. Measurements reveal an order of magnitude improvement in the contrast of nematic flow speeds between maximally- and minimally-illuminated states for opto-K365 motors when compared to opto-K401 construct. For opto-K365 nematics, we characterize both the steady-state flow and defect density as a function of applied light. We also examine the transient behavior as the system switches between steady-states upon changes in light intensities. Although nematic flows reach a steady state within tens of seconds, the defect density exhibits transient behavior for up to 10 minutes, showing a separation between small-scale active flows and system-scale structural states. Furthermore, our work establishes an experimental platform that can exploit spatiotemporally-heterogeneous patterns of activity to generate targeted dynamical states.

42 ENGINEERING↗
DOE OSTI2023

Designing the stripe-ordered cuprate phase diagram through uniaxial-stress

The ability to efficiently control charge and spin in the cuprate high-temperature superconductors is crucial for fundamental research and underpins technological development. Here, we explore the tunability of magnetism, superconductivity and crystal structure in the stripe phase of the cuprate La 2-x Ba x CuO 4 , with x = 0.115 and 0.135, by employing temperature-dependent (down to 400 mK) muon-spin rotation and AC susceptibility, as well as X-ray scattering experiments under compressive uniaxial stress in the CuO 2 plane. A sixfold increase of the 3-dimensional (3D) superconducting critical temperature T c and a full recovery of the 3D phase coherence is observed in both samples with the application of extremely low uniaxial stress of 0.1 GPa. This finding demonstrates the removal of the well-known 1/8-anomaly of cuprates by uniaxial stress. On the other hand, the spin-stripe order temperature as well as the magnetic fraction at 400 mK show only a modest decrease under stress. Moreover, the onset temperatures of 3D superconductivity and spin-stripe order are very similar in the large stress regime. However, a substantial decrease of the magnetic volume fraction and a full suppression of the low-temperature tetragonal structure is found at elevated temperatures, which is a necessary condition for the development of the 3D superconducting phase with optimal T c . Our results evidence a remarkable cooperation between the long-range static spin-stripe order and the underlying crystalline order with the three-dimensional fully coherent superconductivity. Overall, these results suggest that the stripe- and the SC order may have a common physical mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The SOUX AGN sample: optical/UV/X-ray SEDs and the nature of the disc

We use the SOUX sample of ~700 active galactic nucleus (AGN) to form average optical-ultraviolet (UV)-X-rays spectral energy distributions (SEDs) on a two-dimensional (2D) grid of M BH and L 2500 . We compare these with the predictions of a new AGN SED model, qsosed, which includes prescriptions for both hot and warm Comptonization regions as well as an outer standard disc. This predicts the overall SED fairly well for 7.5 < log(M BH /M ⊙ ) < 9.0 over a wide range in L/L Edd , but at higher masses the outer disc spectra in the model are far too cool to match the data. We create optical-UV composites from the entire Sloan Digital Sky Survey sample and use these to show that the mismatch is due to there being no significant change in spectral shape of the optical-UV continuum across several decades of M BH at constant luminosity. We show for the first time that this cannot be matched by standard disc models with high black hole spin. These apparently fit, but are not self-consistent as they do not include the General Relativistic effects for the emission to reach the observer. At high spin, increased gravitational redshift compensates for almost all of the higher temperature emission from the smaller inner disc radii. The data do not match the predictions made by any current accretion flow model. Either the disc is completely covered by a warm Comptonization layer whose properties change systematically with L/L Edd , or the accretion flow structure is fundamentally different to that of the standard disc models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spectral properties of three-dimensional Anderson model

The three-dimensional Anderson model represents a paradigmatic model to understand the Anderson localization transition. In this work we first review some key results obtained for this model in the past 50 years, and then study its properties from the perspective of modern numerical approaches. Our main focus is on the quantitative comparison between the level sensitivity statistics and the level statistics. While the former studies the sensitivity of Hamiltonian eigenlevels upon inserting a magnetic flux, the latter studies the properties of unperturbed eigenlevels. We define two versions of dimensionless conductance, the first corresponding to the width of the level curvature distribution relative to the mean level spacing, and the second corresponding to the ratio of the Heisenberg time and the Thouless time obtained from the spectral form factor. We show that both conductances look remarkably similar around the localization transition, in particular, they predict a nearly identical critical point consistent with other well-established measures of the transition. We then study some further properties of those quantities: for level curvatures, we discuss particular similarities and differences between the width of the level curvature distribution and the characteristic energy studied by Edwards and Thouless (1972) in their pioneering work, in which the hopping at one lattice edge is changed from periodic to antiperiodic boundary conditions. In the context of the spectral form factor, we show that at the critical point it enters a broad time-independent regime, in which its value is consistent with the level compressibility obtained from the level variance. Finally, we test the scaling solution of the average level spacing ratio in the crossover regime using the cost function minimization approach introduced recently in Šuntajs et al. (2020). The latter approach seeks for the optimal scaling solution in the vicinity of the crossing point, while at the same time it allows for the crossing point to drift due to finite-size corrections. We find that the extracted transition point and the scaling coefficient agree with those from the literature to high numerical accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A numerical examination of the effect of sulfide dissolution on silicate weathering

Silicate weathering regulates climate on geological time scales as a net, climate-sensitive sink of atmospheric CO 2 . On the other hand, sulfuric acid produced through sulfide dissolution affects silicate weathering, diminishing its effectiveness as a climate regulator based on evidence from river chemistry. This study takes a theoretical approach to quantitatively examine the effect of sulfide dissolution, coupling a one-dimensional model of pyrite weathering to that of albite transformation to kaolinite. The coupled model reveals that when the reaction front of sulfide is deeper than that of silicates reacting with CO 2 , the silicate-weathering feedback on climate is not directly affected by sulfide dissolution, but only indirectly through oceanic processes such as reverse weathering and carbonate deposition. In turn, when sulfide dissolution occurs within zones of CO 2 -silicate reactions close to the surface, the feedback between climate and silicate weathering can be significantly weakened because sulfuric acid competes effectively with carbonic acid in mineral dissolution. As the depth of the sulfide reaction front can be changed by variations in water table depth, the latter feedback weakening could have occurred in concert with climate-change-driven water table fluctuations, amplifying atmospheric CO 2 variations during the Phanerozoic.

58 GEOSCIENCES↗
DOE OSTI2022

Impact of the electron density and temperature gradient on drift-wave turbulence in the Large Plasma Device

In this paper we present an experimental study of edge turbulence in the Large Plasma Device at UCLA. We utilize a scan of discharge power and prefill pressure (neutral density) to show experimentally that turbulent density fluctuations decrease with decreasing density gradient, as predicted for resistive drift-wave turbulence (RDWT). As expected for RDWT, we observe that the cross-phase between the density and potential fluctuations is close to 0. Moreover, the addition of an electron temperature gradient leads to a reduction in the amplitude of the density fluctuations, as expected for RDWT. However, counter to theoretical expectations, we find that the potential fluctuations do not follow the same trends as the density fluctuations for changes either in density gradients or the addition of a temperature gradient. The disconnect between the density and potential fluctuations is connected to changes in the parallel flows as a result of differences in the prefill pressure, i.e. neutral density. Further analysis of the density and potential fluctuation spectra show that the electron temperature gradient reduces the low frequency fluctuations up to $10 \,{\rm kHz}$ and the introduction of a temperature gradient leads to an unexpected ${\sim }{\rm \pi}$ shift of the density–potential cross-phase at ${\sim }10\,{\rm kHz}$ , while maintaining the typical resistive drift-wave cross-phase at lower frequencies. These experiments partly confirm existing knowledge on resistive drift-wave turbulence, but also introduce new observations that indicate a need for dedicated nonlinear three-dimensional turbulence simulations that include neutrals.

Physics↗
DOE OSTI2021

Selective Biosorption of Valuable Rare Earth Elements Among Co-Occurring Lanthanides

To meet the increasing demand for rare earth elements (REEs) essential for low carbon-intensity energy technologies, new methods are needed for selective REE extraction from unconventional resources. Only a few REEs have significant economic value, but isolation of target rare earths from co-occurring lanthanides is challenging due to their similar chemical behavior. We present a novel approach to enhance separation of specific REEs from lanthanide mixtures. Escherichia coli cells engineered with lanthanide binding tags (LBTs) were immobilized in nonadsorbing, permeable polyethylene glycol diacrylate beads and packed into continuous flow, fixed-bed columns. Breakthrough of 15 rare earths in the +3 oxidation state resulted in notable differences in adsorption selectivity, with greatest separation between europium (Eu) and lanthanum (La) due to competitive displacement. REE adsorption onto fixed-bed columns was predicted by coupling a surface complexation model to a calibrated one-dimensional dual porosity transport model that accounts for interbead advective and intrabead diffusive transport. A tradeoff between high-abundance, low-affinity native carboxyl sites and low-abundance, high-affinity engineered LBT sites dictates process recovery efficiency and selectivity. Key chemical and operational parameters are identified to maximize selective extraction of high-value lanthanides, achieving a threefold enhancement of Eu recovery relative to La in a mixed REE solution.

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

Visualizing Energy Transfer at Buried Interfaces in Layered Materials Using Picosecond X-Rays

Understanding the fundamentals of nanoscale heat propagation is crucial for next-generation electronics. For instance, weak van der Waals bonds of layered materials are known to limit their thermal boundary conductance (TBC), presenting a heat dissipation bottleneck. In this paper we present a new non-destructive method to probe heat transport in nanoscale crystalline materials using time-resolved x-ray measurements of photo-induced thermal strain. This technique directly monitors time-dependent temperature changes in the crystal and the subsequent relaxation across buried interfaces by measuring changes in the c-axis lattice spacing after optical excitation. We investigate films of five different layered transition metal dichalcogenides MoX 2 [X= S, Se, Te] and WX 2 [X= S, Se] as well as graphite and a W-doped alloy of MoTe 2 , on c-plane sapphire substrates, finding TBC values in the range 10-30 MWm -2 K -1 at room temperature. In conjunction with molecular dynamics simulations, we show that the high thermal resistances are a consequence of weak interfacial van der Waals bonding and low phonon irradiance. This work paves the way for an improved understanding of thermal bottlenecks in emerging three-dimensional heterogeneously integrated technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microphysical Piggybacking in the Weather Research and Forecasting Model

Abstract This paper presents incorporation of the microphysical piggybacking into the Weather Research and Forecasting (WRF) model. Microphysical piggybacking is to run a single simulation applying two microphysical schemes, the first scheme driving the simulation and the second piggybacking this simulated flow. “Driving the simulation” means that the simulated microphysical processes, affect the cloud buoyancy and thus force the simulated flow. In contrast, the piggybacking variables are advected by the simulated flow and undergo microphysical transformation, but they do not affect the simulated flow (like in prescribed flow—kinematic—simulations). The two sets of variables (driver and piggybacker) include temperature, water vapor mixing ratio, and all microphysical variables. We provide details of implementing piggybacking into the WRF model, illustrate its applications, and demonstrate the benefits of this methodology in two idealized three‐dimensional cases: (a) a squall line case applying two microphysics schemes, the Thompson bulk microphysics scheme and the University of Pécs/NCAR bin (UPNB) scheme. The piggybacking simulations revealed that the microphysics‐dynamics interaction plays a more important role than the pure microphysical size sorting effect in the transition zone formation. (b) A case of daytime shallow‐to‐deep convective development over land. This case uses the UPNB scheme and contrasts convection developing in environments with either pristine or polluted cloud condensation nuclei (CCN). The piggybacking results indicated that the increase of cloud cover and decrease of supersaturation are mainly associated with the microphysical effect of increasing CCN while the change of precipitation on the ground is also influenced by microphysics‐dynamics interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hot-spot model for inertial confinement fusion implosions with an applied magnetic field

Imposing a magnetic field on inertial confinement fusion implosions magnetizes the electrons in the compressed fuel; this suppresses thermal losses, which increases temperature and fusion yield. Indirect-drive experiments at the National Ignition Facility with 12 and 26 T applied magnetic fields demonstrate up to 40% increase in temperature, 3× increase in fusion yield, and indicate that magnetization alters the radial temperature profile [Moody et al., Phys. Rev. Lett. 129, 195002 (2022); Lahmann et al., APS DPP (2022)]. In this work, we develop a semi-analytic hot-spot model, which accounts for the two-dimensional (2D) Braginskii anisotropic heat flow due to an applied axial magnetic field. First, we show that hot-spot magnetization alters the radial temperature profile, increasing the central peakedness, which is most pronounced for moderately magnetized implosions (with 8–14 T applied field), compared to both unmagnetized (with no applied field) and highly magnetized (with 26 T or higher applied field) implosions. This model explains the trend in the experimental data, which finds a similarly altered temperature profile in the 12 T experiment. Next, we derive the hot-spot model for gas-filled (Symcap) implosions, accounting for the effects of magnetization on the thermal conduction and in changing the radial temperature (and density) profiles. Using this model, we compute predicted central temperature amplification and yield enhancement scaling with the applied magnetic field. The central temperature fits the experimental data accurately, and the discrepancy in the yield suggests a systematic (independent of applied field) degradation, such as mix, and additional degradation in the reference unmagnetized shot, such as reduced laser drive, increased implosion asymmetry, or the magnetic field suppressing ablator mixing into the hot-spot.

Alpha particles↗
DOE OSTI2021

Steady-State Mechanical Analysis for Target Assembly in the Material Plasma Exposure eXperiment Facility

The Material Plasma Exposure eXperiment (MPEX) steady-state linear plasma facility is currently under design at Oak Ridge National Laboratory to expose target specimens to fusion divertor regimes. The neutron-irradiated target is actively cooled and remote handled in the MPEX facility for conducting plasma-material–interaction (PMI) experiments. Here, the steady-state stresses in the target and target assembly system are investigated using two-dimensional (2-D) and three-dimensional (3-D) models to provide expected stresses/strains under the heat loads to which various system components would be exposed during MPEX operation. The calculated temperatures from the 2-D axisymmetric mechanical model were found to be in excellent agreement with those from the full 3-D thermohydraulic model, providing a strong model validation. Numerical simulation results for the steady-state mechanical model indicate nonuniform distributions for the temperature, stress, and deformation within the critical components. For the initial design, the deformation results indicate possible gap openings between contacting surfaces below the plasma-facing materials. To reduce the possibility of interfacial gap opening, the target assembly was slightly changed and evaluated using the 2-D stress model. Numerical simulation results indicate that the interfacial gap openings can be minimized without drastically changing the entire target assembly. The stress-strain conditions for the target will be further used to assess the appropriate operation during MPEX experiments and gain insight into materials science phenomena during PMI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

The multi-region relaxed magnetohydrodynamics (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The stepped pressure equilibrium code (SPEC) was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad–Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Lastly, we provide an example equilibrium with a 3D helical field structure as the favored lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗