Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Charge-exchange losses of beam ions in MAST Upgrade: experiments and modelling

Dedicated experiments were performed on MAST Upgrade to study beam-ion losses caused by charge exchange (CX) with edge neutrals. The fuelling was switched from the high-field side to the low-field side mid-discharge. Direct measurements suggest a strong increase in the neutral density around the plasma and a decrease in the beam-ion density, which is qualitatively explained by CX losses. Measurements by a resistive bolometer have suggested particle bombardment during neutral beam injection, providing a unique opportunity to separate CX from other loss mechanisms. To verify and quantify CX losses, the orbit-following code ASCOT, which accounts both for CX neutralization and reionization, was used to simulate beam-particle power loads on the bolometer. Simulations reproduce measured bolometer power loads during high-field-side fuelling, verifying CX losses of approximately 10% of the off-axis beam power. Toroidally symmetric simulations overestimate power loads on the bolometer during low-field-side fuelling, which is explained by toroidal asymmetry in the neutral density distribution, as is demonstrated by toroidally asymmetric simulations. Results suggest significantly higher CX losses during low-field-side fuelling, up to about 50% of off-axis beam power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Numerical Modeling of CO2 Sequestration within a Five-Spot Well Pattern in the Morrow B Sandstone of the Farnsworth Hydrocarbon Field: Comparison of the TOUGHREACT, STOMP-EOR, and GEM Simulators

The objectives of this study were (1) to assess the fate and impact of CO2 injected into the Morrow B Sandstone in the Farnsworth Unit (FWU) through numerical non-isothermal reactive transport modeling, and (2) to compare the performance of three major reactive solute transport simulators, TOUGHREACT, STOMP-EOR, and GEM, under the same input conditions. The models were based on a quarter of a five-spot well pattern where CO2 was injected on a water-alternating-gas schedule for the first 25 years of the 1000 year simulation. The reservoir pore fluid consisted of water with or without petroleum. The results of the models have numerous broad similarities, such as the pattern of reservoir cooling caused by the injected fluids, a large initial pH drop followed by gradual pH neutralization, the long-term persistence of an immiscible CO2 gas phase, the continuous dissolution of calcite, very small decreases in porosity, and the increasing importance over time of carbonate mineral CO2 sequestration. The models differed in their predicted fluid pressure evolutions; amounts of mineral precipitation and dissolution; and distribution of CO2 among immiscible gas, petroleum, formation water, and carbonate minerals. The results of the study show the usefulness of numerical simulations in identifying broad patterns of behavior associated with CO2 injection, but also point to significant uncertainties in the numerical values of many model output parameters.

STOMP↗
DOE OSTI2024

Downfolding from ab initio to interacting model Hamiltonians: comprehensive analysis and benchmarking of the DFT+cRPA approach

Abstract Model Hamiltonians are regularly derived from first principles to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited-state energies and wave functions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double-counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double-counting corrections diminish the quality. Background screening of the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians.

Chemistry↗
DOE OSTI2023

SMPS data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow, 1.0 LPM, is directed through TSI's 0.071 cm impactor attached to the classifier of the SMPS setup. The SMPS’s DMA and CPC are connected through a 20.3 cm length of 0.48 cm diameter tubing. Measured SMPS size distributions were used to calculate size-dependent particle losses for each SMPS scan. Particle losses from diffusion (based on Kesten, 1991 and Gormley, 1949) and inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977) were included in the loss calculation. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interaction energy and isosteric heat of adsorption between hydrogen and magnesium diboride

Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).

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

A phase-field study of stainless-steel oxidation from high-temperature carbon dioxide exposure

An electrochemical phase-field model has been developed to investigate the oxidation mechanisms of the 21-2N valve stainless steel alloy exposed to carbon dioxide at 973 K. Three oxide phases observed in oxidation experiments are included in the model: Mn 3 O 4 , Cr 2 O 3 , and MnCr 2 O 4 . Local charge neutrality and the conserved current condition are assumed to include the impact of the electric potential on oxidation, though it is found to be negligible when the system is electrically isolated. The sensitivity of the oxidation processes to the values of the diffusion mobilities is examined, and the oxidation rate is calibrated against experimental data by modifying the sensitive mobilities. We find that both inward oxygen and outward metal diffusion are important for oxidation. Furthermore, by investigating the impact of the order of the initial oxide layers, we find that Cr 2 O 3 serves as a barrier to outward Mn diffusion. Moreover, controlling the initial oxide layer order can be used to limit oxidation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Climatic Impact‐Driver Framework for Assessment of Risk‐Relevant Climate Information

Abstract The climate science and applications communities need a broad and demand‐driven concept to assess physical climate conditions that are relevant for impacts on human and natural systems. Here, we augment the description of the “climatic impact‐driver” (CID) approach adopted in the Working Group I (WGI) contribution to the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report. CIDs are broadly defined as “physical climate system conditions (e.g., means, events, and extremes) that affect an element of society or ecosystems. Depending on system tolerance, CIDs and their changes can be detrimental, beneficial, neutral, or a mixture of each across interacting system elements and regions.” We give background information on the IPCC Report process that led to the development of the 7 CID types (heat and cold, wet and dry, wind, snow and ice, coastal, open ocean, and other) and 33 distinct CID categories, each of which may be evaluated using a variety of CID indices. This inventory of CIDs was co‐developed with WGII to provide a useful collaboration point between physical climate scientists and impacts/risk experts to assess the specific climatic phenomena driving sectoral responses and identify relevant CID indices within each sector. The CID Framework ensures that a comprehensive set of climatic conditions informs adaptation planning and risk management and may also help prioritize improvements in modeling sectoral dynamics that depend on climatic conditions. CIDs contribute to climate services by increasing coherence and neutrality when identifying and communicating relevant findings from physical climate research to risk assessment and planning activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Degradation and recovery of solid oxide fuel cell performance by control of cathode surface acidity: Case study – Impact of Cr followed by Ca infiltration

Solid oxide fuel cells (SOFC) have attracted attention as clean and efficient energy conversion devices with low emissions. However, several degradation mechanisms limit the electrochemical performance of current SOFCs, with cathode degradation due to Cr-poisoning from metal interconnects particularly problematic. The acidity/basicity of binary additives has been found to be a sensitive descriptor of the oxygen exchange kinetics, indicating that acidic Cr-species/basic Ca-species can be expected to deactivate/activate the cathode surface, respectively. Inspired by recent advances, the feasibility of relative acidity as a tool for reviving degraded SOFCs is demonstrated by neutralizing Cr-poisoned SOFCs by subsequent serial infiltration of Ca-species. Furthermore, a model mixed ionic and electronic conducting oxide, Pr 0.1 Ce 0.9 O 2-δ (PCO), is selected as the cathode material. Area-specific resistances (ASR) of symmetric cells obtained by electrochemical impedance spectroscopy show that Cr-infiltration results in a seven-fold increase in ASR, while subsequent infiltration of Ca-species leads to complete recovery. Performance degradation and recovery are attributed to depressed/enhanced redox properties at the PCO surface, as supported by XPS analysis. Experiments using anode-supported fuel cells show a reduction in peak power density by 26% upon Cr-infiltration, reversed following Ca-infiltration, after which no degradation is observed during subsequent operation for 100 h.

(Pr,Ce)O2-δ↗
DOE OSTI2022

Cryo-electron Microscopy of Adeno-associated Virus

Adeno-associated virus (AAV) has a single-stranded DNA genome encapsidated in a small icosahedrally symmetric protein shell with 60 subunits. AAV is the leading delivery vector in emerging gene therapy treatments for inherited disorders, so its structure and molecular interactions with human hosts are of intense interest. A wide array of electron microscopic approaches have been used to visualize the virus and its complexes, depending on the scientific question, technology available, and amenability of the sample. Approaches range from subvolume tomographic analyses of complexes with large and flexible host proteins to detailed analysis of atomic interactions within the virus and with small ligands at resolutions as high as 1.6 Å. Analyses have led to the reclassification of glycan receptors as attachment factors, to structures with a new-found receptor protein, to identification of the epitopes of antibodies, and a new understanding of possible neutralization mechanisms. AAV is now well-enough characterized that it has also become a model system for EM methods development. Heralding a new era, cryo-EM is now also being deployed as an analytic tool in the process development and production quality control of high value pharmaceutical biologics, namely AAV vectors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Stoichiometric Lanthanide Compounds with Diglycolamides: A Synthetic Approach toward Understanding Rare-Earth Speciation in Solution

A fundamental understanding of coordination chemistry across the lanthanide series is essential for explaining the chemical behavior of rare-earth metals in complex liquid–liquid extraction processes. Probing the exact bonding between the extractant and the metal is sometimes done through the synthesis of solid-state compounds that can serve as models for metal speciation in solution. In the case of diglycolamide (DGA), a commonly used neutral diamide extractant, extensive studies identify the stepwise formation of 1:1 [Ln(DGA)(H 2 O) 6 ] 3+ , 1:2 [Ln(DGA) 2 (H 2 O) 3 ] 3+ , and 1:3 [Ln(DGA) 3 ] 3+ complexes in solution. The crystallographic reports, however, exclusively show a 1:3 [Ln(DGA) 3 ] 3+ moiety in the solid state, while the structures of 1:1 and 1:2 complexes are yet to be isolated and comprehensively studied. In this work, we report the synthesis and characterization of three new families of stoichiometric N,N,N′,N’ -tetramethyldiglycolamide (TMDGA) compounds: [Ln(TMDGA)(H 2 O) 5 Cl]Cl 2 ·2H 2 O, [Ln(TMDGA) 2 (H 2 O) 3 ]Cl 3 ·3H 2 O, and [Ln(TMDGA) 3 ]Cl 3 ·7H 2 O, where Ln = Nd, Eu, Gd, with Ln:TMDGA ratios of 1:1, 1:2, and 1:3, respectively. Further, the compounds have been analyzed using vibrational spectroscopy (both Raman and FT-IR), as well as variable temperature fluorescence spectroscopy. A spectrophotometric titration of Eu(III) was performed to confirm the presence of the [Eu(TMDGA) n ] 3+ ( n = 1, 2, and 3) species in solution and to compare the individual solid-state emission spectra to their respective analogues in solution.

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

Impact of toroidal rotation on the resistive ballooning modes in ASDEX Upgrade tokamak

Abstract In this work, we investigate the behavior of instabilities appearing between type-I edge localized modes (ELMs), with increasing neutral beam injection (NBI) power concomitant increase in toroidal rotation, and compare it to the modeling result of the linear magneto-hydrodynamic (MHD) code CASTOR3D. An injection of one NBI beam, increasing toroidal rotation, results in the mode slowing down from 12 kHz to 7 kHz, and its associated radial displacement decreases from 5 mm to 3.5 mm. In addition, modes shift radially outwards towards higher q , decreasing their poloidal mode numbers. The mode velocity is measured to be close to the E × B velocity with significant uncertainties. Through a set of CASTOR3D simulations with varying profiles, resistivity has been identified as the primary contributor to the growth rates. Only a small stabilizing effect due to toroidal rotation has been observed. While experimental results show a decrease of mode frequency with rotation, the opposite trend is observed in modeling. Reasons for discrepancies between modeling and experiment are discussed. Nevertheless, a main contributor to the mode frequency has been identified to be rotation velocity. CASTOR3D classifies modes as resistive ballooning modes as they do not appear unstable in ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Error Analysis on Numerical Integration Algorithms in a Hypoelasticity Framework

This report determines local truncation errors for common stress integration algorithms used in explicit finite element codes with hypoelastic material models. The hypoelastic integration algorithms in question utilize an operator splitting procedure in a rotation neutralized configuration, where the stress response is determined from de- coupling the total deformation into rotational and strain dependent components. This document analyzes the error in evolving the stress given a one-step time increment Δt and compares the errors associated with both the rotational and strain components of the operator splitting method. A slight modification to a traditional algorithm is proposed and studied, where the rate of deformation is appropriately rotated from the midstep configuration at t n+1/2 to the end step configuration at t n+1 before the constitutive evaluation. The proposed modification either completely eliminates the error associated with the rotation rate or is of the same order of magnitude as the original algorithm for the three test cases considered in this report. These cases consist of an unaxial stretch with a constant true strain rate with a rigid body rotation, an uniaxial stretch with a constant engineering strain rate with a rigid body rotation, and a simple shear deformation. All three cases are compared to a closed form solution, and in almost every test case the alternative algorithm yields the most accurate one-step local truncation error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Deterministic Transport: An Introduction [Slides]

In this lecture, we cover deterministic particle transport. We focus on neutral particles (neutrons, gamma-rays, thermal photons). Our goal is to predictively and accurately model how free particles interact with their environments and properties of systems involving those particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

First Sterile Neutrino Search at NOvA Experiment Using both Neutrino and Anti-Neutrino Data

NOvA is a long-baseline experiment with two functionally identical detectors with a Near Detector placed 1 km from the neutrino source at Fermilab and a Far Detector 810 km away at Ash River, Minnesota. NOvA uses a very high-intensity ($\sim$900 kW) beam of neutrinos and antineutrinos from Fermilab's NuMI beamline, together with two functionally identical detectors placed 14 mrad off the beam axis. NOvA is built to investigate the complex properties of neutrinos with an emphasis on neutrino oscillations. NOvA not only probes active neutrino mixing but also explores exotic oscillations, including sterile neutrino searches. This analysis uses the joint $\nu_{\mu}$ and neutral current (NC) disappearance channels to probe active-sterile mixing in the 3+1 model. Furthermore, this analysis leverages our statistics of 27$\times$10$^{20}$ protons on target (POT) in neutrino mode and 12.5$\times$10$^{20}$ protons on target (POT) in antineutrino mode. In this poster, we will present the first dual-baseline search for sterile neutrinos using both the neutrino and antineutrino datasets.

Chaudhary, S. [Indian Inst. Tech., Guwahati]↗
DOE OSTI2026

Final Report: A Multi-Channel Fusion Product

The goal of this project was to measure charged fusion products from the d(d,p)t reaction in MAST-U plasmas as a function of time and position with good energy resolution using a system of up to six charged particle detectors. The data from this new diagnostic will make it possible to determine the neutral beam ion density profile as a function of R, z, and t with reduced model dependency and contribute new information to a global analysis of fast ion diagnostic data needed for the determination of the fast ion distribution function (velocity space tomography).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Regulation of the central safety factor and normalized beta under low NBI torque in DIII-D

An algorithm has been designed to simultaneously control the central safety factor (q 0 ) and normalized beta (β N ) while ensuring near-zero torque from the neutral beam injection in DIII-D. Feedback control of q 0 and β N in tokamaks can be beneficial due to the close relationship that these variables have with plasma performance and magneto-hydrodynamic stability. In addition, low neutral-beam-torque conditions are of special interest in present devices because future burning-plasma tokamaks such as ITER will most likely operate at very low plasma rotation. The control synthesis of the algorithm presented in this work is based on a linearized, one-dimensional (1D) model of the current-profile dynamics coupled with a zero-dimensional (0D) plasma-energy balance. The actuators considered are neutral beam injection and electron-cyclotron heating and current drive, and discrete logic determines the neutral-beam injection powers that deliver near-zero torque. Here, the algorithm has been tested in nonlinear, 1D simulations using COTSIM (Control-Oriented Transport SIMulator) and in DIII-D experiments, demonstrating satisfactory performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Exploring the potential of hydrogen in decarbonizing China's light-duty vehicle market

The Chinese government has pledged to achieve overall carbon neutrality by 2060. Currently, the transportation sector contributes to about 10% of total greenhouse gas (GHG) emissions in China. Hence, China has created a well-defined energy vehicle development strategy to reduce GHG emissions from the transportation sector, further expanding into hydrogen vehicle technologies. In this study, the Transportation Energy Analysis Model (TEAM) investigates the potential of hydrogen internal combustion engine vehicles (H2- ICEVs) and fuel cell vehicles (FCEVs) as a reliable pathway towards the government's aspiration of carbon neutrality in the transportation sector. According to TEAM, by adopting FCEVs and H2-ICEVs in the vehicle market, hydrogen demand could reach 25% of the total light-duty transportation energy demand in 2050. Consequently, this will lead to an annual reduction of more than 35 million tons GHG compared to only counting on the electrification pathway in the decarbonization task. Besides, FCEVs would take longer to penetrate the light-duty vehicle market compared to H2-ICEVs, as the current fuel cell technology still requires much improvement to attain a competitive vehicle cost of production.

Life cycle analysis↗