Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Level Structure of {sup 31}S Via {sup 32}S({ital p,d}) {sup 31}S

Background: Properties of proton-unbound S-31 states determine the P-30(p, gamma) S-31 reaction rate, which has a significant impact on explosive hydrogen burning in classical novae and type-I x-ray bursts. Despite several previous studies, uncertainties still remain with respect to the nuclear structure of S-31 near the proton threshold. Purpose: The level structure of S-31 has been presently investigated via a charged-particle spectroscopy experiment using the S-32(p, d) S-31 reaction. Method: Deuterons corresponding to S-31 excited states with 3.285 <= E-x <= 10.8 MeV were momentum analyzed via an Enge split-pole spectrograph at six laboratory angles between 10 degrees and 62 degrees. Differential cross sections of the S-32(p, d)S-31 reaction were measured at E-p = 34.5 MeV. Distorted-wave Born approximation calculations were performed to constrain the spin-parity assignments of several of the observed levels. Results: We have detected 72 excited states of S-31, out of which 17 are within the astrophysical region of interest corresponding to the temperature range of 0.1-1.5 GK. We have resolved the discrepancy in the spin and parity of an excited state with E-x = 6542 keV, showing that is it not J(pi) = 3/2(-), and therefore the contribution of this state to the P-30(p, gamma) reaction rate is likely much less significant than previously thought owing to the larger angular-momentum transfer required to populate this excited state. Moreover, our measurement results help consolidate the spin-parity assignments for the 6377 and 6636 keV states in S-31. Conclusions: This work presents the most comprehensive spin-parity assignments to date from a single-neutron transfer reaction on S-32 to S-31 excited states in the region between 6 to 7 MeV excitation energy. This region is significant for the determination of the P-30(p, gamma)S-30 reaction rate over the temperatures characteristic of explosive hydrogen burning in novae.

Setoodehnia, K.↗
DOE OSTI2023

Pre-equilibrium evolution of conserved charges with initial conditions in the ICCING Monte Carlo event generator

Heavy-ion collisions can be well described through relativistic viscous hydrodynamics, but questions still remain when hydrodynamics is applicable because the initial state may begin very far from equilibrium. Thus, a pre-equilibrium evolution phase is used to bridge the gap between the initial state and hydrodynamics. KøMPøST is one such pre-equilibrium model that propagates the energy-momentum tensor by decomposing it into the background and fluctuations around that background, whose evolution is captured by Green's functions. We extend this formalism to include conserved charges and calculate the corresponding nonequilibrium Green's functions in the relaxation-time approximation. The ICCING algorithm initializes conserved charges in the initial state by sampling $g$ → $q$$\overline{q}$ splitting probabilities and is, thus, perfectly positioned to implement Green's functions for charge propagation. As a result, we show that this method alters the initial-state charge geometries and is applicable in central to mid-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Protecting Customer Privacy Through Distributed Energy Resource Anonymization

Due to their stochastic nature, the increase of Renewable Energy Resources (RERs) as a primary source of energy for power grids creates challenges regarding the reliability and resilience of the system. In order to combat these obstacles, expansion of Distributed Energy Resources (DERs) and their participation in Demand Response (DR) programs is necessary. Widespread participation requires prioritizing customer privacy and addressing concerns that may arise regarding communication between DERs and the Grid Service Provider (GSP). This paper discusses the use of flow reservation resources to split the operating cycles of DER load profiles into unique phases. The splitting of phases increases anonymization of the DERs by making it more difficult to determine the individual characteristics of the device. We discuss an example of this using simulated DER load profile data and examine the resulting effectiveness by using a machine learning algorithm for classification, called Support Vector Machine (SVM).

Distributed Energy Resource, Anonymization, Renewa↗
DOE OSTI2025

Design of a Structure for Assembly and Cooling the Magnet of the Next-Generation 45 GHz ECR Ion Source MARS-D

Here, the current Electron Cyclotron Resonance Ion Sources (ECRISs), constructed with Nb-Ti wires and the conventional racetrack-and-solenoid structure, have achieved operating frequencies up to 28 GHz and utilized about 90% of the critical current of the Nb-Ti wire. A Mixed Axial and Radial field System Demonstrator (MARS-D) is being developed at Lawrence Berkeley National Laboratory (LBNL). This system, which consists of an innovative hexagonal Closed-Loop Coil (CLC) and a set of solenoids, can generate higher magnetic fields (up to 150% ) while requiring only about 50% of the superconducting wire, enabling Nb-Ti wires to be used in the next-generation 45 GHz ECRIS. However, the assembly and cooling of such an efficient and compact magnet are particularly challenging due to the small radial gap between the CLC and solenoids, as well as the tight operating temperature margin. To address these challenges, a structure was developed that combines a three-section radially split solenoid mandrel with a series of shrink-fit reinforcement rings and cooling channels. This paper presents the detailed structure, manufacturing method, assembly procedure, impregnation method, mechanical Finite Element Analysis (FEA) comparison, and thermal FEA comparison.

ECR ion sources↗
DOE OSTI2024

Response of Graphite to Dynamic Loading and Hypervelocity Jet Impacts

The compressive strengths of three varieties of high purity graphite, PCEA, NBG-18, and NBG-25, as well as the depth of penetration of small-scale charges into these materials was experimentally determined. These grades are similar in density, ranging from 1.80 – 1.85 g/cc, and nominal apparent porosity, ranging from 18% to 20%, but provide a wide range in maximum grain or particle size from 10s to 1000s of µm. Two very different manufacturing methods are also represented; PCEA is extruded while NBG-18 and NBG-25 are iso-molded. The quasistatic and dynamic strengths of each grade were determined on a load frame and split-Hopkinson pressure bar, respectively. The depth of penetration (DOP) of two small-scale shaped charges, the Teledyne RP-1 and RP-4, was determined against graphite. The global response of the RP-4 impacts was markedly different as the PCEA samples remained intact while all the NBG-25 samples split into 2 or 3 pieces after the jet penetration had completed. However, for all tests, the trusted DOPs fell within 2 cm. Preliminary hydrocode modeling of the penetration events used existing models that were not designed for graphite. The results can be tuned to reasonably reproduce the DOP, but the wound channel geometry is not reproduced well. A model designed for graphite would need to represent graphite’s non-linear and energy dissipation characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Local element operations for curved simplex meshes

Summary Mesh optimization procedures are generally a combination of node smoothing and discrete operations which affect a small number of elements to improve the quality of the overall mesh. These procedures are useful as a post‐processing step in mesh generation procedures and in applications such as fluid simulations with severely deforming domains. In order to perform high‐order mesh optimization, these ingredients must also be extended to high‐order (curved) meshes. In this work, we present a method to perform local element operations on curved meshes. The mesh operations discussed in this work are edge/face swaps, edge collapses, and edge splitting (more generally refinement) for triangular and tetrahedral meshes. These local operations are performed by first identifying the patch of elements which contain the edge/face being acted on, performing the operation as a “straight‐sided one" by placing the high‐order nodes via an isoparametric mapping from the master element, and smoothing the high‐order nodes on the elements in the patch by minimizing a Jacobian‐based high‐order mesh distortion measure. Since the initial “straight‐sided guess” from the placement of the nodes via the isoparametric mapping frequently results in invalid elements, the distortion measure must be regularized which allows for mesh untangling for the optimization to succeed. We present several examples in 2D and 3D to demonstrate these local operations and how they can be combined with a high‐order node smoothing procedure to maintain mesh quality when faced with severe deformations.

Shi, A.↗
DOE OSTI2022

A new generation of effective core potentials from correlated and spin–orbit calculations: Selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with 4d, 5d, 6s, and 6p valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin–orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and SO splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin–orbit effects. The subsequent quantum Monte Carlo calculations point out the importance of accurate trial wave functions that, in some cases (mid-series transition elements), require treatment well beyond a single-reference.

Wang, Guangming (ORCID:0000000231327925)↗
DOE OSTI2025

Hydrogen and deuterium tunneling in niobium

We use density functional methods to identify the atomic configurations of H and D atoms trapped by O impurities embedded in bulk Nb. The O atoms are located at the octahedral position in the Nb body-centered cubic (BCC) lattice, and H (D) atoms tunnel between two degenerate tetrahedral sites separated by a mirror plane. Using nudged elastic band (NEB) methods, we calculate the double-well potential for O-H and O-D impurities and the wave functions and tunnel splittings for H and D atoms. Our results agree with those obtained from analysis of heat capacity and neutron scattering measurements on Nb with low concentrations of O-H and O-D.

density functional theory↗
DOE OSTI2025

Advanced Brine Processing to Enable U.S. Lithium Independence (CRADA Report)

Current production of LiOH, which is needed to make Li-ion battery cathode active materials, utilizes a multistep process including solar evaporation, precipitation with Na 2 CO 3 and then conversion to LiOH using Ca(OH) 2 . This process requires a large amount of land area for solar evaporation, the right weather conditions, and chemicals for the conversion process that result in NaCl and CaCO 3 waste products. The production of Ca(OH) 2 is very energy intensive and evolves significant quantities of CO 2 . An alternative process flow utilizing direct lithium extraction techniques, followed by a chemical free conversion process can have benefits in reducing the needed land requirements and chemicals for traditional brine processing. There are many potential direct lithium extraction technologies that are currently being developed. The direct lithium extraction process from typical brine sources will produce a LiCl solution with some impurities including typically high concentrations of Na. This brine then needs to be converted to LiOH for use in battery cathode production. Ideally this conversion could occur without the use of additional chemicals. Electrochemistry can do this conversion either via electrolysis or bipolar membrane electrodialysis (BPED) to produce LiOH and HCl in solution. BPED utilized bipolar membranes to split water, which has a reduced potential as compared to splitting water at electrodes into hydrogen and oxygen gas. This reduced potential required results in a significant energy savings for BPED over electrolysis methods. This CRADA project aimed to develop such an integrated process using direct lithium extraction followed by BPED to produce a LiOH solution. That solution can then be crystallized into battery grade LiOH. In particular, Albemarle utilized a direct lithium extraction process to produce a concentrated LiCl solution that could be used for the BPED process. The BPED process was first tested using various LiCl solutions with impurity ions added at bench scale to understand the effects of impurities and determine processing parameters. Then testing was performed using the direct lithium extracted brine at the bench scale before scaling the process up. After the process was scaled up a long duration test was carried out to estimate the lifetime of the membranes, which is key to the economics of the BPED process.

25 ENERGY STORAGE↗
DOE OSTI2026

Double-real corrections to color singlet decay in a parton-shower inspired scheme

We introduce a local infrared subtraction method for next-to-next-to-leading order QCD calculations in color singlet decays, with counterterms based on scalar radiators and pure splitting functions. Overlapping singularities in the multipole radiation pattern are disentangled by partial fractioning, and the kinematics mapping corresponds to iterated next-to-leading order kinematics. We verify that the double-real remainder to $e^+e^-\to\;q\bar{q}$ is rendered finite in the single and double unresolved limits and investigate the numerical convergence of the Monte-Carlo integral. We compute the phase-space integrals of the scalar counterterms in the back-to-back configuration, both analytically and with the help of numerical techniques based on sector decomposition.

Campbell, John M. [Fermilab]↗
DOE OSTI2024

Dynamic population balance in molecular-level simulations of hypersonic flows

This report summarizes the work towards developing stochastic weighted particle methods (SWPM) for future application in hypersonic flows. Extensive changes to Sandia’s direct simulation Monte Carlo (DSMC) solver, SPARTA (Stochastic Particle Real Time Analyzer), were made to enable the necessary particle splitting and reduction capabilities for SWPM. The results from one-dimensional Couette and Fourier flows suggest that SWPM can reproduce the correct transport for a large range of Knudsen numbers with adequate accuracy. The associated velocity and temperature profiles are in good agreement with DSMC. An issue with particle placement during particle number reduction, is identified, to which, a simple but effective solution based on minimizing the center of mass error is proposed. High Mach wheel flows are simulated using the SWPM and DSMC methods. SWPM is capable of providing nearly an order of magnitude increase in efficiency over DSMC while retaining high accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hard x-ray – optical four-wave mixing using a split-and-delay line

New, hard x-ray free electron lasers (FEL) produce intense femtosecond-to-attosecond pulses at angstrom wavelengths, giving access to the fundamental spatial and temporal scales of matter. These revolutionary light sources open the door to applying the suite of nonlinear, optical spectroscopy methods at hard x-ray photon energies. Nonlinear spectroscopy with hard x-rays can allow for measuring the coherence properties of short wavelength excitations with atomic specificity and for understanding how high energy excitations couple to other degrees of freedom in atomic, molecular or condensed-phase systems. As a step in this direction, here we present hard x-ray, optical four-wave mixing (4WM) measurements done at 9.8 keV at the split-and-delay line at the x-ray correlation spectroscopy (XCS) hutch of the Linac Coherent Light Source (LCLS). In this work, we create an x-ray transient grating (TG) from a pair of crossing x-ray beams and diffract optical laser pulses at 400 nm from the TG. The key technical advance here is being able to independently vary the delays of the x-ray pulses. Measurements were made in 3 different solid samples: bismuth germinate (BGO), zinc oxide (ZnO) and yttrium aluminum garnet (YAG). The resulting phase-matched, 4WM signal is measured in two different ways: by varying the x-ray, x-ray pulse delay which can reveal both material and light source coherence properties and also by varying the optical laser delay with respect to the x-ray TG to study how the x-ray excitation couples to the optical properties. Although no coherent 4WM signal was seen in these measurements, the absence of this signal gives important information on experimental requirements for detecting this in future work. Also, our laser-delay scans, although not a new measurement, were applied to different materials than in past work and reveal new examples x-ray induced lattice dynamics in solids. This work represents a key step towards extending nonlinear optics and time-resolved spectroscopy into the hard x-ray regime.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Connected Traffic Signal Coordination Optimization Framework through Network-Wide Adaptive Linear Quadratic Regulator–Based Control Strategy

Traffic congestion in metropolitan areas causes several significant challenges, such as longer travel times, decreased productivity, increased fuel consumption and vehicle emissions, and even severe injuries during crashes. Traffic signal control is a management approach to reduce traffic congestion and allocate the appropriate right of way for safety and mobility efficiency, both in temporal and spatial domains. Here, this study proposes a network-wide adaptive signal control coordination optimization framework based on the linear quadratic regulator algorithm. The traffic flow conditions driven by signal control inputs are formulated based on their network-wide state-space representation. After modeling traffic control regulation constraints, an adaptive linear quadratic regulator algorithm is designed to maximize the network-wide total throughput under the current conditions. Optimal signal control split time durations for multiple intersections in the network are derived by solving the algebraic Riccati equation. Furthermore, the recursive least square parameter estimation method is employed to quantify dynamic traffic condition changes. To verify the effectiveness of this proposed signal control framework, both simulation and real-world experimental tests are conducted for multiple intersections in downtown Chattanooga, Tennessee, United States. In preparation for real-world experimental tests, pipelines for real-time data processing implementation and historical traffic flow data analysis are conducted. The test results demonstrate that the proposed control framework achieves a decrease in travel time by up to 19.4%, total time spent (TTS) by up to 11.9%, and relative queue balance (RQB) by up to 15.6%. The research findings indicate that the proposed signal control framework can be generalized to handle large scale signal control optimization network-wide.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Analysis of the Fe 2p XPS for hematite α Fe2O3: Consequences of covalent bonding and orbital splittings on multiplet splittings

The origins of the complex Fe 2p X-Ray Photoelectron Spectra (XPS) of hematite (α-Fe2O3) are analyzed and related to the character of the bonding in this compound. This analysis provides a new and novel view of the reasons for XPS binding energies (BEs) and BE shifts, which deepens the current understanding and interpretation of the physical and chemical significance of the XPS. In particular, many-body effects are considered for the initial and the final, 2p-hole configuration wavefunctions. It is shown that a one-body or one-configuration analysis is not sufficient and that the many-body, many-determinantal, and many-configurational character of the wavefunctions must be taken into account to describe and understand why the XPS intensity is spread over an extremely large number of final 2p-hole multiplets. The focus is on the consequences of angular momentum coupling of the core and valence open shell electrons, the ligand field splittings of the valence shell orbitals, and the degree of covalent mixing of the Fe(3d) electrons with the O(2p) electrons. Novel theoretical methods are used to estimate the importance of these various terms. An important consequence of covalency is a reduction in the energy separation of the multiplets. Although shake satellites are not considered explicitly, the total losses of intensity from the angular momentum multiplets to shake satellites is determined and related to the covalent character of the Fe-O interaction. The losses are found to be the same for Fe 2p1/2 and 2p3/2 ionization.

Bagus, Paul S. (ORCID:0000000349315217)↗
DOE OSTI2020

Development and validation of a 30-day mortality index based on pre-existing medical administrative data from 13,323 COVID-19 patients: The Veterans Health Administration COVID-19 (VACO) Index

Background Available COVID-19 mortality indices are limited to acute inpatient data. Using nationwide medical administrative data available prior to SARS-CoV-2 infection from the US Veterans Health Administration (VA), we developed the VA COVID-19 (VACO) 30-day mortality index and validated the index in two independent, prospective samples. Methods and findings We reviewed SARS-CoV-2 testing results within the VA between February 8 and August 18, 2020. The sample was split into a development cohort (test positive between March 2 and April 15, 2020), an early validation cohort (test positive between April 16 and May 18, 2020), and a late validation cohort (test positive between May 19 and July 19, 2020). Our logistic regression model in the development cohort considered demographics (age, sex, race/ethnicity), and pre-existing medical conditions and the Charlson Comorbidity Index (CCI) derived from ICD-10 diagnosis codes. Weights were fixed to create the VACO Index that was then validated by comparing area under receiver operating characteristic curves (AUC) in the early and late validation cohorts and among important validation cohort subgroups defined by sex, race/ethnicity, and geographic region. We also evaluated calibration curves and the range of predictions generated within age categories. 13,323 individuals tested positive for SARS-CoV-2 (median age: 63 years; 91% male; 42% non-Hispanic Black). We observed 480/3,681 (13%) deaths in development, 253/2,151 (12%) deaths in the early validation cohort, and 403/7,491 (5%) deaths in the late validation cohort. Age, multimorbidity described with CCI, and a history of myocardial infarction or peripheral vascular disease were independently associated with mortality–no other individual comorbid diagnosis provided additional information. The VACO Index discriminated mortality in development (AUC = 0.79, 95% CI: 0.77–0.81), and in early (AUC = 0.81 95% CI: 0.78–0.83) and late (AUC = 0.84, 95% CI: 0.78–0.86) validation. The VACO Index allows personalized estimates of 30-day mortality after COVID-19 infection. For example, among those aged 60–64 years, overall mortality was estimated at 9% (95% CI: 6–11%). The Index further discriminated risk in this age stratum from 4% (95% CI: 3–7%) to 21% (95% CI: 12–31%), depending on sex and comorbid disease. Conclusion Prior to infection, demographics and comorbid conditions can discriminate COVID-19 mortality risk overall and within age strata. The VACO Index reproducibly identified individuals at substantial risk of COVID-19 mortality who might consider continuing social distancing, despite relaxed state and local guidelines.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Electrolyzers in the System Advisor Model (SAM): A Techno-Economic Potential Study

A technoeconomic analysis of grid to low temperature electrolysis (Grid-LTE), photovoltaic to low temperature electrolysis (PV-LTE), concentrating solar power to high temperature electrolysis (CSP-HTSE) and a concentrating solar power with PV to high temperature electrolysis (CSP-PV-HTSE) centralized hydrogen production systems are analyzed to assess the economics of system and to provide a baseline for comparing these technologies against hydrogen production cost targets. A framework integrating the system advisor model (SAM) and US Department of Energy hydrogen production models (H2A) is developed to assess these systems. The hydrogen levelized cost given current and future assumptions for technology cost and performance is evaluated at optimal system configurations. The framework described in this report integrates SAM with H2A electrolyzer technologies and provides analysts a detailed technoeconomic method to analyze concentrating and photovoltaic solar technologies to produce energy that are directly coupled to LTE and HTSEs that use that energy to split water into hydrogen and oxygen. The baseline hydrogen levelized cost (HLCs) for the GRID-LTE, PV-LTE, CSP-HTSE, and CSP-PV-HTSE systems in Daggett, CA are 2.82, 3.86, 3.68, and 2.90 $\$$USD 2016/kg H2 and 2.50, 2.13, 2.84, 2.15 $\$$USD 2016/kg H2 in the 2020 and 2050 scenarios respectively. To achieve the $\$$2/kg H2 target in locations with excellent solar resources, cost parameters values aligned with aggressive R&D targets will need to be achieved for all the systems configurations. In Daggett, PV costs of $\$$0.68 /Wac or moderate ATB PV CAPEX projections result in HLCs of $\$$2 /kg H2. Similarly for PV-MSALT-HTSE systems, $\$$0.60 /Wac result in $\$$2/kg H2. For the MSALT-HTSE systems, better than aggressive 2050 ATB salt tower CAPEX projections would be needed to reach $\$$2/kg H2. Molten salt tower capital costs of $\$$2400/kW would enable $\$$2/kg H2 in 2050.

08 HYDROGEN↗