Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Doubling method”

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

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗
DOE OSTI2021

Vertical-Axis Wind Turbine Steady and Unsteady Aerodynamics for Curved Deforming Blades

Vertical-axis wind turbines’ simpler design and low center of gravity make them ideal for floating wind applications. However, efficient design optimization of floating systems requires fast and accurate models. Low-fidelity vertical-axis turbine aerodynamic models, including double multiple streamtube and actuator cylinder theory, were created during the 1980s. Commercial development of vertical-axis turbines all but ceased in the 1990s until around 2010 when interest resurged for floating applications. Despite the age of these models, the original assumptions (2-D, rigid, steady, straight bladed) have not been revisited in full. When the current low-fidelity formulations are applied to modern turbines in the unsteady domain, aerodynamic load errors nearing 50% are found, consistent with prior literature. However, a set of steady and unsteady modifications that remove the majority of error is identified, limiting it near 5%. This paper shows how to reformulate the steady models to allow for unsteady inputs including turbulence, deforming blades, and variable rotational speed. A new unsteady approximation that increases numerical speed by 5–10× is also presented. Combined, these modifications enable full-turbine unsteady simulations with accuracy comparable to higher-fidelity vortex methods, but over 5000× faster.

17 WIND ENERGY↗
DOE OSTI2024

Broadening of the Divertor Heat Flux Profile in High Confinement Tokamak Fusion Plasmas with Edge Pedestals Limited by Turbulence in DIII-D

Multimachine empirical scaling predicts an extremely narrow heat exhaust layer in future high magnetic field tokamaks, producing high power densities that require mitigation. Here, in the experiments presented, the width of this exhaust layer is nearly doubled using actuators to increase turbulent transport in the plasma edge. This is achieved in low collisionality, high confinement edge pedestals with their gradients limited by turbulent transport instead of large-scale, coherent instabilities. The exhaust heat flux profile width and divertor leg diffusive spreading both double as a high frequency band of turbulent fluctuations propagating in the electron diamagnetic direction doubles in amplitude. The results are quantitatively reproduced in electromagnetic XGC particle-in-cell simulations which show the heat flux carried by electrons emerges to broaden the heat flux profile, directly supported by Langmuir probe measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

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

An Assessment of the Laminar Hypersonic Double-Cone Experiments in the LENS-XX Tunnel

This is an investigation on two experimental datasets of laminar hypersonic flows, over a double-cone geometry, acquired in Calspan—University at Buffalo Research Center’s Large Energy National Shock (LENS)-XX expansion tunnel. These datasets have yet to be modeled accurately. A previous paper suggested that this could partly be due to mis-specified inlet conditions. The authors of this paper solved a Bayesian inverse problem to infer the inlet conditions of the LENS-XX test section and found that in one case they lay outside the uncertainty bounds specified in the experimental dataset. However, the inference was performed using approximate surrogate models. Here in this paper, the experimental datasets are revisited and inversions for the tunnel test-section inlet conditions are performed with a Navier–Stokes simulator. The inversion is deterministic and can provide uncertainty bounds on the inlet conditions under a Gaussian assumption. It was found that deterministic inversion yields inlet conditions that do not agree with what was stated in the experiments. An a posteriori method is also presented to check the validity of the Gaussian assumption for the posterior distribution. This paper contributes to ongoing work on the assessment of datasets from challenging experiments conducted in extreme environments, where the experimental apparatus is pushed to the margins of its design and performance envelopes.

42 ENGINEERING↗
DOE OSTI2020

Measurement of the Medium Energy NuMI Flux Using the Low-$\nu$ and High-$\nu$ Methods at MINERvA

As we continue to improve our experimental capabilities, the measurements we make of accelerator-based neutrino phenomena are becoming increasingly limited by systematic, rather than statistical uncertainties. Among these are uncertainties corresponding to our understanding of the neutrino flux and to the models that we use to simulate neutrino interactions. These already are significant for currently-running neutrino experiments, and will become more important for next-generation neutrino experiments such as DUNE. Dedicated efforts at experiments such as MINERvA are pursuing improved comprehensive modeling of neutrino cross sections and improving our understanding of how to precisely measure neutrino fluxes. This thesis discusses these themes in the context of MINERvA, a neutrino scattering experiment which was located on-axis in the NuMI beamline at Fermilab. A pair of \textit{in situ} measurements of the flux are presented which utilize Medium Energy $\nu_{\mu}$ data colle cted at MINERvA between 2013 and 2017. The low-$\nu$ technique is used to make a direct measurement of the flux from low-inelasticity events, and that flux is used to measure the total charged-current inclusive cross section. The novel ``high-$\nu$'' technique is used to make a second direct measurement of the flux from a complementary selection of high-inelasticity events, and that flux is used to measure a double-differential cross section of charged-current quasielastic-like events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Measurement of the central exclusive production of charged particle pairs in proton-proton collisions at $\sqrt{s}$ = 200 GeV with the STAR detector at RHIC

We report on the measurement of the Central Exclusive Production of charged particle pairs h+h- (h = π, K, p) with the STAR detector at RHIC in proton-proton collisions at $\sqrt{s}$ = 200 GeV. The charged particle pairs produced in the reaction pp → p' + h + h - + p' are reconstructed from the tracks in the central detector and identified using the specific energy loss and the time of flight method, while the forward-scattered protons are measured in the Roman Pot system. Exclusivity of the event is guaranteed by requiring the transverse momentum balance of all four final-state particles. Differential cross sections are measured as functions of observables related to the central hadronic final state and to the forward-scattered protons. They are measured in a fiducial region corresponding to the acceptance of the STAR detector and determined by the central particles’ transverse momenta and pseudorapidities as well as by the forward-scattered protons’ momenta. This fiducial region roughly corresponds to the square of the four-momentum transfers at the proton vertices in the range 0.04 GeV 2 < -t 1 , -t 2 < 0.2 GeV 2 , invariant masses of the charged particle pairs up to a few GeV and pseudorapidities of the centrally-produced hadrons in the range |η| < 0.7. The measured cross sections are compared to phenomenological predictions based on the Double Pomeron Exchange (DPE) model. Structures observed in the mass spectra of π + π - and K + K - pairs are consistent with the DPE model, while angular distributions of pions suggest a dominant spin-0 contribution to π + π - production. For π + π - production, the fiducial cross section is extrapolated to the Lorentz-invariant region, which allows decomposition of the invariant mass spectrum into continuum and resonant contributions. The extrapolated cross section is well described by the continuum production and at least three resonances, the f 0 (980), f 2 (1270) and f 0 (1500), with a possible small contribution from the f 0 (1370). Fits to the extrapolated differential cross section as a function of t 1 and t 2 enable extraction of the exponential slope parameters in several bins of the invariant mass of π + π - pairs. These parameters are sensitive to the size of the interaction region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Second-generation downscaled earth system model data using generative machine learning

The second-generation Sup3rCC dataset provides high-resolution meteorological data generated through the downscaling of multiple earth system models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). This downscaling is performed through application of a generative machine learning approach called Super-Resolution for Renewable Resource Data (sup3r). This dataset builds on the first-generation Sup3rCC data by applying improved bias correction methods and adding downscaled precipitation to the output variables. As with the first Sup3rCC version, the data still include temperature, wind speed and direction at multiple heights, pressure, three components of downwelling solar radiation, and relative humidity—all at 4-kilometer (km) hourly resolution over the contiguous United States. This is a 25x spatial enhancement and 24x temporal enhancement of the source 100-km daily-average ESM data. This extension of the Sup3rCC dataset includes data from six ESMs from two shared socioeconomic pathways (SSPs) totaling 400 years of data with multiple future projections of changing meteorological conditions. The scenario selection was based on a structured evaluation of historical ESM skill and comprehensive representation of possible trajectories of future climate change in temperature, humidity, precipitation, solar irradiance, and near-surface wind speeds. The inclusion of multiple future projections is intended to enable users to assess key drivers of un 36 certainty and variability. All data are double-bias corrected, resulting in a product that can be used out-of-the-box for energy system analysis with minimal historical bias. The potential applications of Sup3rCC data extend to various topics in renewable energy resource assessment, energy systems modeling, and grid resilience studies. High-resolution future meteorological projections are critical for evaluating the effects of changing meteorological conditions on renewable energy generation, energy demand, and for optimizing energy storage and grid infrastructure. The 4-km hourly resolution of the downscaled data enables understanding of spatial and temporal variability at the scales necessary for energy system operational planning. In addition, the dataset can support risk assessments by providing detailed information on possible future extreme weather events and long-term meteorological variability at scales relevant to energy infrastructure. By offering an enhanced representation of possible future meteorological conditions, the second-generation Sup3rCC dataset enables more precise modeling of energy resilience and adaptation strategies in response to changing meteorological conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

An Orbital Basis Set for Double Photoionization of Atoms and Molecules

The ab initio theoretical treatment of one-photon double photoionization processes has been limited to atoms and diatomic molecules by the challenges posed by large grid-based representations of the double ionized continuum wave function. To provide a path for extensions to polyatomics, an energy-adapted orbital basis approach is demonstrated that reduces the dimensions of such representations and simultaneously allows larger time steps in time-dependent computational descriptions of double ionization. Additionally, an algorithm that exploits the diagonal nature of the two-electron integrals in the grid basis and dramatically accelerates the transformation between grid and orbital representations is presented. Excellent agreement between the present results and benchmark theoretical calculations is found for H – and Be atoms, as well as the hydrogen molecule, including for the triply differential cross sections that relate the angular distribution and energy sharing of all of the particles in the molecular frame.

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

Pressure cells for in situ neutron total scattering: time and real-space resolution during deuterium absorption

In situ gas-loading sample holders for two-dimensionally arranged detectors in time-of-flight neutron total scattering experiments have been developed to investigate atomic arrangements during deuterium absorption using time and real-space resolution. A single-crystal sapphire container was developed that allows conditions of 473 K and 10 MPa hydrogen gas pressure. High-resolution transient measurements detected deuterium absorption by palladium that proceeded within a few seconds. A double-layered container with thick- and thin-walled vanadium allowed conditions of 423 K and 10 MPa hydrogen gas pressure. The deuterium occupation sites of a lanthanum–nickel–aluminium alloy are discussed in detail on the basis of real-space high-resolution data obtained from in situ neutron scattering measurements and reverse Monte Carlo structural modeling.

08 HYDROGEN↗
DOE OSTI2025

Revisiting the Relationship Between Induced Polarization and Surface Conductivity: Ratios From Laboratory to Field

Abstract Among the subsurface geophysical methods used in the critical zone investigations, induced polarization (IP) shows great vitality thanks to its unique ability to assess porosity via bulk conduction and estimate permeability through surface conduction and/or polarization. However, such an advantageous separation between bulk and surface is mostly implemented by multi‐salinity experiments in the laboratory, which is incredibly difficult to realize in the field. One promising approach to address such an obstinate issue is to gauge the surface conductivity ( σ s ) from the quadrature conductivity ( σ ″) or normalized chargeability ( M n ) with the ratios between ( l = σ ″/ σ s , l mn = M n / σ s ). While these ratios are known not to be universal, the underlying principles are not fully understood and relevant theoretical studies are rare, which makes quantitative IP applications difficult. Here we scrutinize the conduction and polarization mechanisms of geomaterials and pinpoint that the two ratios are inherently functions of salinities and frequencies rather than only determined by the properties of the electrical double layer (EDL), hence representative samples from the investigated field must be calibrated in the laboratory and a characteristic frequency should be chosen for their usage. Besides the macro‐scale ratios l and l mn , we define two micro‐scale ratios χ and χ mn directly from the EDL, such that the new ratios exclude the effect of salinity and frequency and offer the opportunity to characterize and monitor changes of the EDL. Our study demonstrates that the existing macro‐scale ratios converge toward the values of novel micro‐scale ratios at high water salinity.

Qi, Youzheng [Earth and Environmental Sciences Are↗
DOE OSTI2021

Mixed precision s –step Lanczos and conjugate gradient algorithms

Compared to the classical Lanczos algorithm, the s-step Lanczos variant has the potential to improve performance by asymptotically decreasing the synchronization cost per iteration. However, this comes at a price; despite being mathematically equivalent, the s-step variant may behave quite differently in finite precision, potentially exhibiting greater loss of accuracy and slower convergence relative to the classical algorithm. It has previously been shown that the errors in the s-step version follow the same structure as the errors in the classical algorithm, but are amplified by a factor depending on the square of the condition number of the O(s)-dimensional Krylov bases computed in each outer loop. As the condition number of these s-step bases grows (in some cases very quickly) with s, this limits the s values that can be chosen and thus can limit the attainable performance. In this work, we show that if a select few computations in s-step Lanczos are performed in double the working precision, the error terms then depend only linearly on the conditioning of the s-step bases. This has the potential for drastically improving the numerical behavior of the algorithm with little impact on per-iteration performance. Our numerical experiments demonstrate the improved numerical behavior possible with the mixed precision approach, and also show that this improved behavior extends to mixed precision s-step CG. Here, we present preliminary performance results on NVIDIA V100 GPUs that show that the overhead of extra precision is minimal if one uses precisions implemented in hardware.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High dimensional binary classification under label shift: phase transition and regularization

Label Shift has been widely believed to be harmful to the generalization performance of machine learning models. Researchers have proposed many approaches to mitigate the impact of the label shift, e.g., balancing the training data. However, these methods often consider the underparametrized regime, where the sample size is much larger than the data dimension. The research under the overparametrized regime is very limited. Here, to bridge this gap, we propose a new asymptotic analysis of the Fisher Linear Discriminant classifier for binary classification with label shift. Specifically, we prove that there exists a phase transition phenomenon: Under certain overparametrized regime, the classifier trained using imbalanced data outperforms the counterpart with reduced balanced data. Moreover, we investigate the impact of regularization to the label shift: The aforementioned phase transition vanishes as the regularization becomes strong.

binary classification↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Full-Length ASFV B646L Gene Sequencing by Nanopore Offers a Simple and Rapid Approach for Identifying ASFV Genotypes

African swine fever (ASF) is an acute, highly hemorrhagic viral disease in domestic pigs and wild boars. The disease is caused by African swine fever virus, a double stranded DNA virus of the Asfarviridae family. ASF can be classified into 25 different genotypes, based on a 478 bp fragment corresponding to the C-terminal sequence of the B646L gene, which is highly conserved among strains and encodes the major capsid protein p72. The C-terminal end of p72 has been used as a PCR target for quick diagnosis of ASF, and its characterization remains the first approach for epidemiological tracking and identification of the origin of ASF in outbreak investigations. Recently, a new classification of ASF, based on the complete sequence of p72, reduced the 25 genotypes into only six genotypes; therefore, it is necessary to have the capability to sequence the full-length B646L gene (p72) in a rapid manner for quick genotype characterization. Here, we evaluate the use of an amplicon approach targeting the whole B646L gene, coupled with nanopore sequencing in a multiplex format using Flongle flow cells, as an easy, low cost, and rapid method for the characterization and genotyping of ASF in real-time.

Virology↗
DOE OSTI2024

Code-to-code comparison between FLASH and HYDRA in gas-puff Z-pinch modeling

The numerical modeling of gas-puff Z pinches involves the nonlinear coupling of a broad range of complex, multi-physics phenomena that makes such simulations challenging. The challenge is further compounded by nonlinear instabilities that can impact the dynamics of imploding gas-puff Z pinches, such as the magneto Rayleigh–Taylor instability (MRTI). If the growth rate and amplitude of the latter is comparable to the relevant timescales and properties of the imploding plasma, the MRTI can dramatically alter implosion dynamics, dictate pinch stability, and govern the plasma properties achievable in pulsed-power-driven laboratory experiments. National Laboratories and academic teams have developed numerical tools that can accurately model Z-pinch configurations and provide reliable design capabilities that can guide experimental choices and assist in interpreting experimental results. Most such tools, however, are not broadly available. Here, we present newly developed Z-pinch simulation capabilities of the publicly available FLASH code, applied in the study of MRTI growth and dynamical effects in gas-puff implosions. To verify the new implementations, we perform a comparison of FLASH gas-puff implosion simulations with previously published calculations with the HYDRA code from Lawrence Livermore National Laboratory, which have been validated with experimental data from the CESZAR pulsed-power driver at the University of California, San Diego. The experiments involved double- and triple-nozzle configurations, in an experimental attempt to stabilize the pinch to the MRTI. The code-to-code comparison shows similar results between the FLASH and HYDRA simulations, supporting the use of FLASH in the modeling of future gas-puff Z-pinch experiments at CESZAR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Proton-detected solid-state NMR spectroscopy of spin-1/2 nuclei with large chemical shift anisotropy

Constant-time (CT) dipolar heteronuclear multiple quantum coherence (D-HMQC) has previously been demonstrated as a method for proton detection of high-resolution wideline NMR spectra of spin-1/2 nuclei with large chemical shift anisotropy (CSA). However, 1 H transverse relaxation and t 1 -noise often reduce the sensitivity of D-HMQC experiments, preventing the theoretical gains in sensitivity provided by 1 H detection from being realized. In this paper, we demonstrate a series of improved pulse sequences for 1 H detection of spin-1/2 nuclei under fast MAS, with 195 Pt SSNMR experiments on cisplatin as an example. First, a t 1 -incrementation protocol for D-HMQC dubbed Arbitrary Indirect Dwell (AID) is demonstrated. AID allows the use of arbitrary, rotor asynchronous t 1 -increments, but removes the constant time period from CT D-HMQC, resulting in improved sensitivity by reducing transverse relaxation losses. Next, we show that short high-power adiabatic pulses (SHAPs), which efficiently invert broad MAS sideband manifolds, can be effectively incorporated into 1 H detected symmetry-based resonance echo double resonance (S-REDOR) and t 1 -noise eliminated (TONE) D-HMQC experiments. The S-REDOR experiments with SHAPs provide approximately double the dipolar dephasing, as compared to experiments with rectangular inversion pulses. We lastly show that sensitivity and resolution can be further enhanced with the use of swept excitation pulses as well as adiabatic magic angle turning (aMAT).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗