Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Core-Pedestal Plasma Configurations in Advanced Tokamaks

Here, several configurations for the core and pedestal plasma are examined for a predefined tokamak design by implementing multiple heating/current drive (H/CD) sources to achieve an optimum configuration of high fusion power in a noninductive operation while maintaining an ideally magnetohydrodynamic (MHD) stable core plasma using the IPS-FASTRAN framework. IPS-FASTRAN is a component-based lightweight coupled simulation framework that is used to simulate magnetically confined plasma by integrating a set of high-fidelity codes to construct the plasma equilibrium (EFIT, TOQ, and CHEASE), calculate the turbulent heat and particle transport fluxes (TGLF), model various H/CD systems (TORIC, TORAY, GENRAY, and NUBEAM), model the pedestal pressure and width (EPED), and estimate the ideal MHD stability (DCON). The TGLF core transport model and EPED pedestal model are used to self-consistently predict plasma profiles consistent with ideal MHD stability and H/CD (and bootstrap) current sources. In order to evaluate the achievable and sustainable plasma beta, varying configurations are produced ranging from the no-wall stability to with-wall stability regimes, simultaneously subject to the self-consistent TGLF, EPED, and H/CD source profile predictions that optimize configuration performance. The pedestal density, plasma current, and total injected power are scanned to explore their impact on the target plasma configuration, fusion power, and confinement quality. A set of fully noninductive scenarios are achieved by employing ion-cyclotron, neutral beam injection, helicon, and lower-hybrid H/CDs to provide a broad profile for the total current drive in the core region for a predefined tokamak design. These noninductive scenarios are characterized by high fusion gain (Q ~ 4) and power (P fus ~ 600 MW), optimum confinement quality (H 98 ~ 1.1), and high bootstrap current fraction (f BS ~ 0.7) for Greenwald fraction below unity. The broad current profile configurations identified are stable to low-n kink modes either because the normalized pressure β N is below the no-wall limit or a wall is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Angular distribution of neutron production by proton and carbon-ion therapeutic beams

Carbon-ion beams are increasingly used in the clinical practice for external radiotherapy treatments of deep-seated tumors. At therapeutic energies, carbon ions yield significant secondary products, including neutrons, which may be of concern for the radiation protection of the patient and personnel. We simulated the neutron yield produced by proton and carbon-ion pencil beams impinging on a clinical phantom at three different angles: 15°, 45° and 90°, with respect to the beam axis. We validated the simulated results using the measured response of organic scintillation detectors. We compared the results obtained with FLUKA 2011.2 and MCNPX 2.7.0 based on three different physics models: Bertini, Isabel, and CEM. Over the different ions, energies, and angles, the FLUKA simulation results agree better with the measured data, compared to the MCNPX results. Simulations of carbon ions at low angles exhibit both the highest deviation from measured data and inter-model discrepancy, which is probably due to the different treatment of the pre-equilibrium stage. The reported neutron yield results could help in the comparison of carbon-ion and proton treatments in terms of secondary neutron production for radiation protection applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Dark radiation constraints on portal interactions with hidden sectors

Here, we update dark radiation constraints on millicharged particle (MCP) and gauged baryon-number-minus-lepton-number (B-L) extensions of the Standard Model (SM). In these models, a massive SM gauge singlet mediator couples the SM plasma to additional SM-singlet light degrees of freedom. In the early Universe, these new light particles are populated via the interaction of the SM with the MCP, or the new B-L gauge boson, and act as dark radiation. The presence of dark radiation in the early Universe is tightly constrained by current and upcoming cosmic microwave background (CMB) measurements. We update bounds on MCPs from current measurements of N eff and show that future CMB experiments will be able to rule out or discover the extended MCP model invoked to explain the EDGES anomaly. Our analysis of the gauged B-L model goes beyond previous studies by including quantum-statistical and out-of-equilibrium effects. Further, we account for the finite lifetime of the B-L gauge boson, which boosts the subsequent right-handed neutrino energy density. We also develop a number of approximations and techniques for simplifying and solving the relevant Boltzmann equations. We use our approximations to develop a lower bound on the radiation density in a generic hidden sector with a light relic that is insensitive to the details of the hidden sector, provided the mediator interacts more strongly with the hidden sector than with the SM.

cosmological neutrinos↗
DOE OSTI2022

Prospects for H-mode inhibition in negative triangularity tokamak reactor plasmas

Instability to high toroidal mode number ($\textit{n}$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity (NT) (δ) shaping, preventing access to strong H-mode regimes when δ is sufficiently negative. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$\textit{n}$ ballooning stability as a function of internal profiles and equilibrium shape using a combination of the CHEASE and BALOO codes. While the critical δ required for avoiding 2nd stability to high-$\textit{n}$ modes is observed to depend in a complicated way on various shaping parameters, including the equilibrium inverse aspect ratio, elongation and squareness, equilibria with NT are robustly prohibited from accessing the 2nd stability region, offering the prediction that NT reactors should maintain L-mode-like operation. In order to access high-$\textit{n}$ 2nd stability, the local shear over the entire bad curvature region must be sufficiently negative to overcome curvature destabilization on the low field side. Finally, scalings of the ballooning-limited pedestal height are provided as a function of plasma parameters to aid future scenario design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Excluded-volume model for quarkyonic matter: Three-flavor baryon-quark mixture

The single-flavor excluded-volume model based on the effective size of baryons reproduces the hard-soft density evolution of the equation of state (EoS) required by the recent studies of GW170817. This phenomenological model basically realizes the concept of quarkyonic matter which is introduced from large-$N_c$ gauge theory for dense matter. Further, enhanced nucleon interactions and dynamically generated quark degrees of freedom can reproduce the hard-soft evolution of the EoS. In this paper, we extend the excluded-volume model to a three-flavor system by considering electromagnetic charge and possible weak equilibrium in order to obtain a proper description for the hard-soft behavior of the EoS inferred from the gravitational waves observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Transverse momentum fluctuations and their correlation with elliptic flow in nuclear collision

In this article, we propose observables v 0 and v 0 (p T ) which quantify the relative fluctuations in the total transverse momentum at fixed multiplicity. We first study the factorization of the fixed multiplicity momentum dependent two particle correlation function into a product of v 0 ($p^{a}_{T}$) and v 0 ($p^{b}_{T}$) within realistic hydrodynamic simulations. Then we present computations of v 0 (p T ) for different particle types. We determine the relation between the integrated v 0 and previously measured observables, and compare results from a hybrid hydrodynamics based model to experimental data. The effects of bulk viscosity and an initial pre-equilibrium stage on the results are quantified. We find that v 0 is strongly correlated with the initial state entropy per elliptic area, S/A. Using this result, we explain how the observed correlations between the elliptic flow and the transverse momentum (both in simulations and experiment) reflect the initial state correlations between 1/A and ellipticity ε 2 at fixed multiplicity. Furthermore, we argue that the systematic experimental study of v 0 , with the same sophistication as used for the other v n , can contribute significantly to our understanding of quark gluon plasma properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Nonequilibrium effects on stability of hybrid stars with first-order phase transitions

The stability of hybrid stars with first-order phase transitions as determined by calculating fundamental radial oscillation modes is known to differ from the predictions of the widely used Bardeen-Thorne-Meltzer criterion. We consider the effects of out-of-chemical-equilibrium physics on the radial modes and hence stability of these objects. For a barotropic equation of state, this is done by allowing the adiabatic sound speed to differ from the equilibrium sound speed. We show that doing so extends the stable branches of stellar models, allowing stars with rapid phase transitions to support stable higher-order stellar multiplets similarly to stars with multiple slow phase transitions. We also derive a new junction condition to impose on the oscillation modes at the phase transition. Termed the reactive condition, it is physically motivated, consistent with the generalized junction conditions between two phases, and has the common rapid and slow conditions as limiting cases. Unlike the two common cases, it can only be applied to nonbarotropic stars. Here, we apply this junction condition to hybrid stellar models generated using a two-phase equation of state consisting of nuclear matter with unpaired quark matter at high densities joined by a first-order phase transition and show that, like in the slow limiting case, stars that are classically unstable are stabilized by a finite chemical reaction speed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Wettability Effects on Falling Film Flow and Heat Transfer Over Horizontal Tubes in Jet Flow Mode

The performance of a falling-film heat exchanger is strongly linked to the surface characteristics and the heat transfer processes that take place over the tubes. The primary aim of this numerical study is to characterize the influence of surface wettability on the film flow behavior and its associated surface heat transfer in the jet-flow mode. Volume of fluid (VOF) based simulations are carried out for horizontal tubes with different surface wettabilities. The wettability of the tube surfaces is represented using the Kistler's dynamic contact angle model. Surface wettability effects ranging from superhydrophilic to superhydrophobic are studied by varying the equilibrium contact angle from 2 deg to 175 deg. Two different liquid mass flow rates of 0.06 and 0.18 kg/m-s corresponding to the inline and staggered jet flow modes are studied. Results are presented in terms of the liquid film thickness, the contact areas between the different phases (solid–liquid and liquid–air), and the heat transfer coefficient or Nusselt number. The resistance imposed by the increasing contact angles inhibits the extent of the liquid spreading over the tube surface, and this, in turn, influences the liquid film thickness, and the wetted area of the tube surface. A significant decrement in the heat transfer rate from the tube surfaces was observed as the equilibrium contact angle increased from 2 deg to 175 deg. The local distributions of the Nusselt number over the tube surface are strongly influenced by the flow recirculation in the liquid bulk.

Engineering↗
DOE OSTI2025

Dynamic control of molecular transport MXene transistor membranes

Controlled spatial confinement and surface properties of lamellar 2D nanomaterial membranes could enhance many precision separation processes. Traditionally, researchers view channel dimensions, surface properties, and permeation rates of these membranes as intrinsic properties that cannot be modulated in operando. We report that gate voltage applied to the conducting laminar MXene membrane can modulate the permeation rate of ions and neutral solutes, as well as its effective size rejection. In operando wide-angle x-ray scattering measurements reveal that these changes are not driven by electrically induced variations in the d spacing of the MXene layers. Instead, experimental data and continuum electrokinetic modeling reveal that ion transport through the MXene channels is primarily affected by Donnan equilibrium at the membrane-solution interface. We also report a strong increase in the permeation rates through the membrane under a low-frequency ac voltage gating regime that we attribute to diffusioosmotic flow oscillations induced in the membrane. Overall, MXene “transistor” membranes provide a previously unidentified approach to dynamic control of molecular separations.

Pendse, Aaditya [Lawrence Livermore National Labor↗
DOE OSTI2023

Topological kinetic crossover in a nanomagnet array

Ergodic kinetics, which are critical to equilibrium thermodynamics, can be constrained by a system’s topology. Here, we studied a model nanomagnetic array in which such constraints visibly affect the behavior of the magnetic moments. In this system, magnetic excitations connect into thermally active one-dimensional strings whose motion can be imaged in real time. At high temperatures, our data showed the merging, breaking, and reconnecting of strings, resulting in the system transitioning between topologically distinct configurations. Below a crossover temperature, the string motion is dominated by simple changes in length and shape. In this low-temperature regime, the system is energetically stable because of its inability to explore all possible topological configurations. This kinetic crossover suggests a generalizable conception of topologically broken ergodicity and limited equilibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

LLEPE

LLEPE (liquid-liquid equilibrium parameter estimator) is a tool for estimating thermodynamic parameters like standard enthalpies, standard molar volumes and Pitzer coefficients, by combining Gibbs minimization with regression to multiphase equilibrium isotherm data. It also includes visualization tools to access data quality and prediction model performance.

ILOEJE, CHUKWUNWIKE↗
DOE OSTI2020

Artificial Intelligence for Dynamic Phase Transitions in Materials [PowerPoint]

The presentation covers the following areas: emergent phenomena and non-equilibrium transitions in granular materials; materials science applications, including defects, interfaces and emergent phenomena; modeling with molecular dynamics simulations; principal component analysis (PCA); previous works; PCA applied to nonequilibrium transitions; PCA for varied quenched disorder; PCA incorporating particle velocity; and ending with a summary and impacts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Benchmark Specification of Advanced Burner Test Reactor

As an effort to assess the Argonne Reactor Computation (ARC) suite of fast reactor analysis codes, a numerical benchmark problem was developed using the reference 250 MWt Advanced Burner Test Reactor (ABTR) metallic core fueled with beginning of equilibrium cycle compositions (Chang et al., 2006). Material thermal expansion at operating condition was modeled by adjusting the hexagonal pitch, axial meshes, and the fuel and structure material densities appropriately. Irradiation swelling of metal fuel was considered, and the bond sodium was displaced into the lower part of fission gas plenum.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Safety Analysis for Accident-Tolerant Fuels with Increased Enrichment and Extended Burnup

The ERP research and development efforts in FY22 focused on the safety analysis of ATF with increased enrichment and extended burnup to provide scientific knowledge of the ATF fuel performance, failure mechanism and source term analysis during severe accident of the NPP. FeCrAl clad ATF was selected as basis of the reactor fuel. An optimized equilibrium cycle was developed for 24-months extended burnup operation. A 4-loop PWR model under loss-of-coolant accident (LOCA) was assessed and the behavior of the major source term from failed ATF has been analyzed. The result showed that grace to enhanced mechanical characteristics of ATF, the fuel failure starts up to 200~500 seconds later than conventional Zr clad fuel. This time gives sufficient buffer until low pressure safety injection (LPSI) of the emergency core cooling system (ECCS) restarted by the operator which reduces fuel fragmentation, relocation and dispersal (FFRD). As a consequence, the total amount of major radioactive materials released from FeCrAl clad fuel into the reactor containment was found up to three times less compare to Zr clad fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

1D Simulation of Avalanche to Streamer to Spark Transition of Plasma Discharge in Ammonia-Air Combustion

The backward problem of plasma assisted combustion emphasizes evaluating the effect of the evolving thermochemical state on the plasma discharge. This paper investigates the dependence of avalanche to streamer to spark formation dynamics and kinetics on the gas composition and temperature at different points in an ammonia-air premixed laminar flame using a self-consistent multigrid-based 1D plasma solver. Different values of ..alpha.., the coefficient for effective ionization events per unit length, have been reported for electron avalanches in air and stoichiometric NH3-air mixtures. The streamer inception has been shown to obey the Meek's criterion. An exponential reduction in streamer and spark formation time has been observed from plasma simulations at different points in the unburnt, pre-heat zone, reaction zone and the fully burnt regions of the premixed flame. While the enhancement of the reduced electric field with increasing temperature affects effective ionization, there exists a minimum breakdown field for streamer formation, which does not vary proportionally with the changing number density of the gas. The change in the mixture from reactants (NH3, O2, N2) to products of complete combustion of ammonia in air (N2, H2O) has also been shown to affect the streamer and spark formation. Finally, the major pathways during the streamer and spark phases which are responsible for producing important radicals used in combustion of NH3 are also discussed.

ADVANCED PROPULSION SYSTEMS,INORGANIC, ORGANIC, PH↗
DOE OSTI2024

Applications of Thermochemical Modeling in Molten Salt Reactors

The extensively evaluated and consistent thermodynamic database, the Molten Salt Thermal Properties Database—Thermochemical (MSTDB-TC), was used along with additional thermodynamic values from other sources as examples of ways to examine molten salt reactor (MSR) fuel behavior. Relative stability with respect to halide potential and temperature for likely fuel and fission product components were mapped in Ellingham diagrams for the chloride and fluoride systems. The Ellingham diagrams provide a rich, visual means for identifying halide-forming components in proposed fuel/solvent salt systems. Thermochemical models and values from MSTDB-TC and ancillary sources were used in global equilibrium calculations to provide compositions for a close analysis of the behavior of a possible Molten Chloride Salt Fast Reactor and a Molten Salt Reactor Experiment-type system at high burnup (100 GWd/t). The results illustrated the oxidative nature of burnup in MSRs and provided information about redox behavior and possible control.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗