Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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 OSTI2023

Investigation on Tuning Power–Frequency Droop for Improved Grid–Forming Inverter and Synchronous Generator Transient Load Sharing

Poor transient load sharing has been observed during the parallel operation of synchronous generators (SGs) and droop-controlled grid-forming (GFM) inverter-based resources (IBRs) in islanded microgrids. This can result in overcurrent protections, causing the inverters to trip, which leads to the collapse of voltage and frequency. In this article, small-signal analysis supported by electromagnetic transients (EMTs) simulation is performed on a detailed model of a microgrid containing IBRs and an SG to identify underdamped modes and their sensitivities to the power–frequency (P-F) droop parameter. Further EMT simulation of a more complex microgrid model with multiple GFM IBRs and SGs is performed to determine whether trends identified from the modal analysis extend to cases where the number of SGs in operation may vary. The objective is to investigate the effectiveness of tuning this parameter to reduce the likelihood of overcurrents in response to load and generation disturbances. The work uncovers both fast and slow modes of concern. Primarily, the findings show that increasing the P-F droop value improves the damping of the slower mode while degrading the damping of the faster mode, leading to an upper limit on the damping improvement of the GFM IBR output current. In EMT simulations of the more complex microgrid, the greatest damping improvement still exhibits significant maximum overshoot. Furthermore, the optimal value of the P-F droop parameter is sensitive to the number of SGs in operation. The conclusion is that tuning the P-F droop parameter is not an effective means to avoid the overcurrent tripping of GFM IBRs during response to large load or generation disturbances when operating in parallel to SGs in islanded microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Hot Water Extraction: Short Rotation Willow, Mixed Hardwoods, and Process Considerations

Short rotation woody crops (SRWC) like shrub willow are highly productive biomass resources of interest for energy and fuel applications. Hot water extraction (HWE) as an upgrading tool to enable the use of willow biomass in pellet applications has been proposed, and is of increasing interest. This study treats willow and mixes of willow and conventional mixed hardwood feedstock with HWE in a tumbling laboratory reactor to elucidate the effects of time, temperature, feedstock mixes, and other process considerations (water:biomass ratio, presteaming, counter-current processing) on mass removals and other extraction outcomes (e.g., sugar, acetate, and furan yields). Results demonstrated alignment of extraction outcomes with P-factor from 155 °C to 175 °C, with a good compromise of removed mass and co-product potential in the range from 575–800 P-factor. The preferred condition was chosen as 575 P-factor. HWE of mixes of willow and hardwood feedstocks showed a linear response of extraction outcomes to willow:hardwood ratios. Testing of water:biomass ratios demonstrated that this is a significant consideration, with each outcome being affected somewhat differently, and indicating that HWE is more diffusion dependent than expected. Presteaming shows little to no effect on extraction outcomes, while multi-stage cooks simulating counter-current operation indicate a significant potential value in counter-current extraction.

09 BIOMASS FUELS↗
DOE OSTI2021

Dynamics of laser-generated magnetic fields using long laser pulses

In this work, we report on the experimental investigation of magnetic field generation with a half-loop gold sheet coil driven by long-duration (10 ns) and high-power (0.5 TW) laser pulses. The amplitude of the magnetic field was characterized experimentally using proton deflectometry. The field rises rapidly in the first 1 ns of laser irradiation, and then increases slowly and continuously up to 10 ns during further laser irradiation. The transient dynamics of current shape were investigated with a two-dimensional (2D) numerical simulation that included Ohmic heating of the coil and the resultant change of electrical resistivity determined by the coil material temperature. The numerical simulations show rapid heating at the coil edges by current initially localized at the edges. This current density then diffuses to the central part of the sheet coil in a way that depends both on normal current diffusion as well as temporal changes of the coil resistance induced by the Ohmic heating. The measured temporal evolution of the magnetic field is compared with a model that determines a solution to the coil current and voltage that is consistent with a plasma diode model of the drive region and a 2D simulation of current diffusion and dynamic resistance due to Ohmic heating in the laser coil.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dynamics of Laterally Propagating Flames in X-Ray Bursts. II. Realistic Burning and Rotation

In this work, we continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and thermal structure on the flame. We find that at lower rotation rates the flame becomes harder to ignite, whereas at higher rotation rates the nuclear burning is enhanced by increased confinement from the Coriolis force and the flame propagates steadily. At higher crustal temperatures, the flame moves more quickly and accelerates as it propagates through the atmosphere. If the temperature is too high, instead of a flame propagating across the surface the entire atmosphere burns uniformly. Our findings could have implications for the relationship between observed burst rise times and neutron star rotation and accretion rates. All of the software used for these simulations is freely available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Numerical study of coronal plasma jet formation

A new scenario for solar flare eruption in the coronal holes is analyzed by using MHD stability concepts for a spheromak configuration. The stability properties of a spheromak partially embedded into a conducting surface are studied using three dimensional MHD simulations. In agreement with the analytical theory, a large degree of line-tying stabilizes the spheromak's tilt instability, while the elongation has a destabilizing effect. High-resolution nonlinear simulations also demonstrate current sheet formation at the upper surface of the spheromak, where the tilted magnetic field of the spheromak reconnects with the background magnetic field. The calculated stability threshold and the observed magnetic reconnection support a model of coronal jet eruptions where a dome-like magnetic structure grows through flux emergence on the solar surface, tilts, reconnects, and erupts. Countering the effect from elongation, line-tying strongly stabilizes a spheromak growing from a flux-emergence process, suggesting that to accelerate the onset of eruptive coronal jets, there must be magnetic reconnection at the bottom of the spheromak to detach the structure from the solar surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

PDMATLAB2D: A Peridynamics MATLAB Two-dimensional Code

PDMATLAB2D is a meshfree peridynamics implementation in MATLAB suitable for simulation of two-dimensional fracture problems. The purpose of this code is twofold. First, it provides an entry-level peridynamics computational tool for educational and training purposes. Second, it serves as an accessible and easily modifiable computational tool for peridynamics researchers who would like to adapt the code for a multitude of peridynamics simulation scenarios. The current version of the code implements a bond-based brittle elastic peridynamic model and a critical stretch criterion for bond breaking. However, the code is designed to be extendable for other peridynamic models and computational features. In this paper, we provide an overview of the code structure and functions with illustrative examples. In conclusion, due to the integrated computation and postprocessing MATLAB capabilities, PDMATLAB2D can serve as an effective testbed for testing new constitutive models and advanced numerical features for peridynamics computations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Cross-scale modeling and experimental integration for advancing cathode electrolyte interphase studies in high energy density lithium-ion batteries

Electrochemical interfaces are critical to the performance and durability of lithium-ion batteries (LIBs). The solid electrode-electrolyte interphase (SEI and CEI) structures that form during cycling can passivate reactive surfaces, ensuring safe operation, but also may contribute to performance degradation. Understanding the microscopic factors influencing interphase formation, growth, and evolution is essential for balanced battery design. While significant research has focused on the anode-electrolyte interphase (SEI), the cathode-electrolyte interphase (CEI) remains less explored, despite its importance in high-voltage and advanced battery technologies. Challenges in conducting in-situ or operando experiments arise from the occluded nature of these interfaces and the long timescales involved, often leading to biased interpretations. A validated multi-scale, multi-physics modeling approach, integrated with advanced characterization techniques, can effectively elucidate the intrinsic stability of electrolyte and cathode surfaces, the impact of chemical heterogeneity, and the role of microstructural features on CEI performance. In conclusion, this article reviews current modeling and simulation strategies for studying CEI in advanced LIBs and highlights opportunities for future methodological advancements and experimental integration.

Cathode-electrolyte interphase↗
DOE OSTI2025

Cyclotron resonance accelerators for industrial applications

Here, this paper describes novel configurations for cyclotron resonance acceleration of electrons and ions that have several attractive features including: a compact robust room-temperature single-cell RF cavity as the accelerator structure; and continuous high current accelerated un-bunched beam output with self-scanning, obviating need for a separate beam scanner. An electron accelerator version, the electron Cyclotron Resonance Accelerator (eCRA), is under development to be an efficient source for high power electron and x-ray beams for medical, research, sterilization, and National Security applications, so as to replace radioactive materials. An ion accelerator version, the ion Cyclotron Auto-Resonance Accelerator (iCARA) is described here, suggesting its potential to produce, as an example, a high-current multi-MeV beam of deuterons which could be highly competitive with that produced either with linacs or cyclotrons. Such a deuteron beam could produce a high flux of fast neutrons via deuteron stripping, for applications including the transmutation of used nuclear fuel, material studies relevant for a fusion reactor inner wall, tritium breeding and medical isotope production. For the high-current, high efficiency simulated performance for eCRA and iCARA as described in this paper, the particle beams produced may not exhibit the low emittance values that are important for most discovery research. Rather, the beams could be useful for industrial applications where higher emittance and some energy spread can be tolerated, in favor of high beam power.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Nonperturbative heavy-flavor transport approach for hot QCD matter

The heavy charm and bottom quarks are unique probes of the transport properties of the quark-gluon plasma (QGP) and its hadronization in high-energy nuclear collisions. A key challenge in this context is to embed the interactions of the heavy quarks in the expanding medium compatible with the strong-coupling nature of the QGP, and thus to unravel the underlying microscopic mechanisms. In the present work we progress toward this goal by combining recent T -matrix interactions for elastic scattering with an effective transport implementation of gluon radiation, and apply these in a Langevin framework in a viscous hydrodynamic evolution. Hadronization of heavy quarks is evaluated using a modern recombination model with 4-momentum conservation, supplemented with fragmentation constrained by data in proton-proton collisions. Deploying this approach to charm-hadron observables in Pb-Pb collisions at the LHC yields fair agreement with experiment while also identifying areas of further systematic improvement of the simulations and its current input.

Heavy-flavor transport↗
DOE OSTI2022

Towards fast and accurate predictions of radio frequency power deposition and current profile via data-driven modelling: applications to lower hybrid current drive

Three machine learning techniques (multilayer perceptron, random forest and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ( $n_{e0}$ , $T_{e0}$ , $I_p$ , $B_t$ , $R_0$ , $n_{\|}$ , $Z_{{\rm eff}}$ , $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$ ) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗
DOE OSTI2021

Precision measurement of electron-electron scattering in GaAs/AlGaAs using transverse magnetic focusing

Electron-electron (e-e) interactions assume a cardinal role in solid-state physics. Quantifying the e-e scattering length is hence critical. In this paper we show that the mesoscopic phenomenon of transverse magnetic focusing (TMF) in two-dimensional electron systems forms a precise and sensitive technique to measure this length scale. Conversely we quantitatively demonstrate that e-e scattering is the predominant effect limiting TMF amplitudes in high-mobility materials. Using high-resolution kinetic simulations, we show that the TMF amplitude at a maximum decays exponentially as a function of the e-e scattering length, which leads to a ready approach to extract this length from the measured TMF amplitudes. The approach is applied to measure the temperature-dependent e-e scattering length in high-mobility GaAs/AlGaAs heterostructures. The simulations further reveal current vortices that accompany the cyclotron orbits - a collective phenomenon counterintuitive to the ballistic transport underlying a TMF setting.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Accelerating End-to-End Deep Learning for Particle Reconstruction using CMS open data

Machine learning algorithms are gaining ground in high energy physics for applications in particle and event identification, physics analysis, detector reconstruction, simulation and trigger. Currently, most data-analysis tasks at LHC experiments benefit from the use of machine learning. Incorporating these computational tools in the experimental framework presents new challenges. This paper reports on the implementation of the end-to-end deep learning with the CMS software framework and the scaling of the end-to-end deep learning with multiple GPUs. The end-to-end deep learning technique combines deep learning algorithms and low-level detector representation for particle and event identification. We demonstrate the end-to-end implementation on a top quark benchmark and perform studies with various hardware architectures including single and multiple GPUs and Google TPU.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

O’Hare Airport Short-Term Ground Transportation Modal Demand Forecast Using Gaussian Processes

Here, the principal objective of this study is to analyze the spatial and temporal variation of ground transportation airport demand and provide demand forecast to inform planning capability and explore alternatives for investments to accommodate airport growth. Because of its good adaptability and strong generalization ability for dealing with high-dimensional input, small-sample, and nonlinear spatial data, Gaussian process (GP) regression is used to provide forecast estimates using data from transportation network company (TNC) trips and urban rail passengers at Chicago's O'Hare International Airport. TNC airport trips differ significantly, with three times more distance, more than twice the travel time, and half of the share requests compared with nonairport trips. This highlights the need for separate demand models. Hourly analysis of the rail service indicates that this is likely heavily used by airport workers, whereas TNC services focus on travelers because of variations in the peak demand hours. Heteroscedastic GP regression is implemented because of differences in trip variance between night and day hours. Estimates are given for weekdays and weekend trips, and the 95% confidence intervals are calculated. The introduction of flight schedule information into the models shows marginal improvements in their performance. However, fitting a GP regression becomes computationally expensive with increased sample size and the introduction of spatial components. Transportation planners and policymakers can use the results and methods implemented in this study to optimize transportation assets and provide long-range simulations of the current and future conditions in the area.

42 ENGINEERING↗
DOE OSTI2021

SOLPS-ITER validation with TCV L-mode discharges

This work presents a quantitative test of SOLPS-ITER simulations against tokamak a configuration variable (TCV) L-mode experiments. These simulations account for drifts, currents, kinetic neutrals, and carbon impurities providing the most complete edge transport simulations for TCV to date. The comparison is performed on nominally identical discharges carried out to assess the effectiveness of TCV’s divertor baffles in the framework of the European Plasma Exhaust program and employs numerous edge diagnostics providing a detailed code-experiment benchmark for TCV. The simulations show a qualitative consistency, but the quantitative differences remain, which are assessed herein. It is found that, for a given separatrix density, the simulations most notably yield a colder, and denser, divertor state with a higher divertor neutral pressure than measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Global gyrokinetic simulations of ASDEX Upgrade up to the transport timescale with GENE–Tango

An accurate description of turbulence up to the transport timescale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport timescales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collisions, realistic geometries, toroidal rotation and electromagnetic effects for the first time. This approach overcomes gyrokinetic codes' limitations and enables high-fidelity profile calculations in experimentally relevant plasma conditions, significantly reducing the computational cost. We present numerical results of GENE–Tango for two ASDEX Upgrade discharges, one of which exhibits a pronounced peaking of the ion temperature profile not reproduced by TGLF–ASTRA. We show that GENE–Tango can correctly capture the ion temperature peaking observed in the experiment. By retaining different physical effects in the GENE simulations, e.g., collisions, toroidal rotation and electromagnetic effects, we show that the ion temperature profile's peaking can be linked to electromagnetic effects of submarginal (stable) KBM modes. Furthermore, based on these results, the expected GENE–Tango speedup for the ITER standard scenario is larger than two orders of magnitude compared to a single gyrokinetic simulation up to the transport timescale, possibly making first-principles ITER simulations feasible on current computing resources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗