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 253 records · Page 14

Validation of Simulation Models for Wall-Connected Scrape off Layer Currents during MHD Activity in the HBT-EP Tokamak (Final Report)

This is the final report for the DOE award DE-SC0021325, titled “Validation of Simulation Models for Wall-Connected Scrape off Layer Currents during MHD Activity in the HBT-EP Tokamak” for the period September 1, 2020 – May 31, 2023. This award is part of a broader collaborative project, which supported training of a PhD student at Columbia University under the joint supervision of Dr. Hansen and Columbia project members. The award supported research into computational models for plasma and 3D conducting structures in the vicinity relevant to tokamak disruptions and other mode activity (eg. RWMs) using the NIMROD, and PSI-Tet, codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A 3D simulation comparison of carrier transport in green and blue c-plane multi-quantum well nitride light emitting diodes

Until recently, the electrical efficiency of green nitride light-emitting diodes (LEDs) was considerably lower than that of blue LEDs. This is particularly surprising as one would expect a reduced forward voltage with increasing emission wavelength. In this paper, we theoretically investigated the impact of the number of quantum wells on the forward voltage of III-nitride LEDs with x = 0.15 (blue) and x = 0.24 (green) In x Ga1–xN QWs. The simulated dependence of current density (J) on applied diode bias (V) shows a significant increase of 1.9 V in the forward voltage between one and five quantum well (QW) c-plane green LED structures. Artificially turning off the polarization fields in the simulation does not entirely suppress this effect. Due to the large band offsets in the green LED multiple QW stack, simulations indicate a sequential band filling of the QW sequence. Furthermore, this mechanism should not be limited to c-plane LEDs and could also be present in nonpolar or semipolar devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Counter-Current Flow Limitation Studies in Complex Geometries Utilizing Interface Capturing Simulations Coupled with PID Flow Rate Controller

In nuclear thermal-hydraulic studies, counter-current flow limitation (CCFL) typically refers to steam rising at a fast rate such that it prevents coolant from draining down within a confined channel. CCFL is a crucial issue in nuclear reactor safety analysis. This study investigates CCFL in debris bed channels using high-resolution interface-capturing simulations. A novel proportional-integral-derivative flow rate controller is developed to efficiently achieve the CCFL conditions. Verification studies confirm that CCFL occurs under the same conditions with or without the controller, demonstrating that PID control ensures accurate prediction. Three debris bed channel geometries were examined: a cylindrical channel, a channel with small obstacles, and a channel with large obstacles. Results show that obstacles significantly impact flow behavior, interfacial shear, wall shear, and pressure gradients required for CCFL. Furthermore, the comparison with experimental data confirmed that simulations incorporating geometric complexities align more closely with experimental CCFL conditions. A pressure gradient correlation was also developed for CCFL prediction.

Counter-current flow limitation↗

Analyzing the Effects of the Interference between One- and Two-Body Currents on Neutrino-Nucleus Scattering in ACHILLES

Modeling improvements of neutrino–nucleus scattering are an essential aspect of reducing systematic uncertainties for current and future neutrino oscillation experiments. In particular, retaining quantum mechanical effects such as the interference between one-body (1p1h) and two-body (2p2h) currents in our cross-section models is an essential consideration. This work investigates the impact of this interference term on neutrino–argon scattering using the ACHILLES event generator, which directly implements the interference model developed by Lovato, Rocco, and Steinberg. Using fluxes relevant to SBND, the DUNE near detector, and MiniBooNE, charged- and neutral-current event samples are generated and predicted event rates are plotted against a variety of kinematic variables. For events simulated using the SBND flux, it is found that interference contributes approximately 10 \% to the total charged-current event rate, and serves to enhance the quasi-elastic peak. When considering events generated using the DUNE flux, resonance interactions contribute more to the predicted event rate, as the neutrino flux extends out to higher energies. The relative interference contribution subsequently decreases to roughly 7 \%. Neutral-current simulations show a higher rate of neutron knockout. Finally, a comparison of argon-to-carbon cross-section ratios is performed utilizing fluxes from SBND and MiniBooNE.

Pando, Nathan [Fermilab]↗

Analyzing the Effects of the Interference between One- and Two-Body Currents on Neutrino-Nucleus Scattering in ACHILLES

Modeling improvements of neutrino–nucleus scattering are an essential aspect of reducing systematic uncertainties for current and future neutrino oscillation experiments. In particular, retaining quantum mechanical effects such as the interference between one-body (1p1h) and two-body (2p2h) currents in our cross-section models is an essential consideration. This work investigates the impact of this interference term on neutrino–argon scattering using the ACHILLES event generator, which directly implements the interference model developed by Lovato, Rocco, and Steinberg. Using fluxes relevant to SBND, the DUNE near detector, and MiniBooNE, charged- and neutral-current event samples are generated and predicted event rates are plotted against a variety of kinematic variables. For events simulated using the SBND flux, it is found that interference contributes approximately 10 \% to the total charged-current event rate, and serves to enhance the quasi-elastic peak. When considering events generated using the DUNE flux, resonance interactions contribute more to the predicted event rate, as the neutrino flux extends out to higher energies. The relative interference contribution subsequently decreases to roughly 7 \%. Neutral-current simulations show a higher rate of neutron knockout. Finally, a comparison of argon-to-carbon cross-section ratios is performed utilizing fluxes from SBND and MiniBooNE.

Pando, Nathan [Fermilab]↗

Constant-depth circuits for dynamic simulations of materials on quantum computers

Abstract Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a specific subset of one-dimensional (1D) materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. For an N -spin system, the constant-depth circuit contains only $\mathcal {O}(N^{2})$ O ( N 2 ) CNOT gates. Such compact circuits enable us to successfully execute long-time dynamic simulation of ubiquitous models, such as the transverse field Ising and XY models, on current quantum hardware for systems of up to 5 qubits without the need for complex error mitigation techniques. Aside from enabling long-time dynamic simulations with minimal Trotter error for a specific subset of 1D Hamiltonians, our constant-depth circuits can advance materials simulations on quantum computers more broadly in a number of indirect ways.

Materials simulation↗

Massively parallel phase-field simulations targeting exascale

The interface thickness in the phase-field (PF) method limits its simulation scales. Consequently, large-scale PF simulations become prohibitively expensive for resolving the extremely fine microstructures that typically form during rapid solidification processing. This challenge is significant in predicting microstructure evolution in metal additive manufacturing and has been identified by the United States Department of Energy’s Exascale Computing Project. Here, to address this, we develop a multi-GPU and MPI-based massively parallel simulation code, utilizing state-of-the-art algorithms, software, and libraries, for large-scale three-dimensional (3D) PF simulations. We report the first GPU-parallel PF simulations on Frontier (currently the second TOP500 exascale cluster) and Summit machines, taking dendritic growth as an example problem. We evaluate the parallel performance of our implementation using scaling studies with more than 24 000 GPUs (among the largest known computations to date) and the acceleration performance using large-scale simulations of dendritic growth in 3D. Finally, massively parallel GPUs in these supercomputers enabled the first coupled multiscale simulations of laser melting and subsequent dendritic solidification on the scale of a full melt-pool, demonstrating the feasibility of performing PF simulations with a point total over 2 billion grid points within an acceptable time.

Exascale↗

[Presentation] Long time scale Multiphysics simulation of spent nuclear fuel canister in MOOSE

Pebble-bed reactors are an important class of advanced reactors under consideration for various applications where their fuel would give a significant advantage in siting and high-quality heat production. However, the disposal of their fuel is not as thoroughly studied as other fuel forms. This article provides an example of evaluating advanced reactor spent nuclear fuel in MOOSE. In this study, pebble fuel is analyzed in a well-known spent fuel canister design to characterize the behavior of this fuel form over a long time scale. The results indicate that after approximately 100 years, decay heat is significantly reduced and the maximum temperature in the canister equalizes with the external temperature. The simulation goes on to an end time of one million years, demonstrating the capability of efficiently dealing with long time scales efficiently. We conclude that the canister temperatures seem manageable even with very aggressive loading times and while there are several improvements to be implemented in the future, MOOSE is currently capable of simulating the required scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Typical occupancy profiles and behaviors in residential buildings in the United States

The energy performance of a residential building is highly dependent on occupant's presence or non-presence in a building and their interactions with energy-consuming appliances. Typical occupancy schedules for residential buildings must be defined for applications such as building energy modeling as well as for assessing energy savings associated with the use of occupancy sensing technologies and occupancy-dependent controls. Currently, commonly used simulation programs assume a typical occupancy schedule, however, there is significant opportunity for improvement to these schedules as this is generally based on engineering judgement. This research uses 12 years of the American Time Use Survey (ATUS) data to develop typical occupancy schedules for a range of household types and occupant age ranges. This is compared to currently utilized residential occupancy schedules. In many cases the developed schedules exhibit similar patterns, however, differences are also found to be as high as 41% for certain periods of time. For occupancy sensing applications, the spatial-temporal distribution of occupants in residential buildings is also evaluated. These locations vary based on temporal factors as well as demographic factors such as age and number of occupants. Finally, the results of this research work towards improved occupancy schedule development can benefit both industry professionals and researchers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

UNR Grant Proposal Contribution

The work of Carrier et. al. (citation needed for Matt’s effect of surface roughness on phase transitions/ETI paper) demonstrates that 1D and 2D resistive magnetohydrodynamic (MHD) simulations can reliably model exploding aluminum rods driven by megaampere currents. Close agreement between simulation results and photonic Doppler velocimetry data from the Mykonos electrothermal instability (METI-II) campaign builds confidence in predictive modelling capabilities for pulsed-power HED experiments. Furthermore, 2D MHD simulations show how machined features and micro-scale surface roughness can seed the electrothermal instability (ETI). Surface roughness was observed to reduce the time of melt by 19%, with ETI growth driving the enhanced heating of the rod surface.

36 MATERIALS SCIENCE↗

Computing Free Energies with Fluctuation Relations on Quantum Computers

One of the most promising applications for quantum computers is the dynamic simulation of quantum materials. Current hardware, however, sets stringent limitations on how long such simulations can run before decoherence begins to corrupt results. The Jarzynski equality, a fluctuation theorem that allows for the computation of equilibrium free energy differences from an ensemble of short, non-equilibrium dynamics simulations, can make use of such short-time simulations on quantum computers. Here, we present a quantum algorithm based on the Jarzynski equality for computing free energies of quantum materials. We demonstrate our algorithm using the transverse field Ising model on both a quantum simulator and real quantum hardware. As the free energy is a central thermodynamic property that allows one to compute virtually any equilibrium property of a physical system, the ability to perform this algorithm for larger quantum systems in the future has implications for a wide range of applications including the construction of phase diagrams, prediction of transport properties and reaction constants, and computer-aided drug design.

Bassman, Lindsay↗

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SASHIMI-SIDM: semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales

We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that SASHIMI-SIDM agrees with SIDM subhalo populations from high-resolution cosmological zoom-in simulations in resolved regimes. Crucially, we predict that the fraction of core-collapsed subhalos peaks at a mass scale determined by the input SIDM cross section and decreases toward higher halo masses, consistent with the predictions of gravothermal models and cosmological simulations. For the first time, we also show that the core-collapsed fraction decreases toward lower halo masses. While the dependence of the collapse time on mass and concentration implies such behaviour, our semi-analytical approach allows us to quantify and illustrate this trend clearly across the full mass spectrum of subhalos, including for subhalo masses below the resolution limit of any current cosmological SIDM simulation. As a proof of principle, we apply SASHIMI-SIDM to predict the boost to the local dark matter density and annihilation rate from core-collapsed SIDM subhalos, which can be enhanced relative to CDM by an order of magnitude for viable SIDM models. Thus, SASHIMI-SIDM provides an efficient and reliable tool for scanning SIDM parameter space and testing it with astrophysical observations. The code is publicly available at https://github.com/shinichiroando/sashimi-si.

cosmological simulations↗

Intelligently Partitioned Phasor-EMT Hybrid Simulations of Large-Scale, High-IBR Power Systems

As the penetration level of power electronics-interfaced renewables such as photovoltaics (PV) and wind has surged in modern electric grids, new operational risks caused by the dynamics of those inverter-based resources (IBRs) are emerging in parallel. Lessons learned from various grid events include that the impact of IBRs on system-level grid stability will become prominent along with the increase of renewables and that the short-timescale dynamic impacts of IBRs on grid stability are not fully captured by current commercial dynamic simulation tools [1] [2]. For example, IBRs can be controlled to mitigate those destabilizing interactions, but conventional phasor-domain tools (e.g. PSS/E, PSLF) often cannot capture that; likewise, the existing electromagnetic transient (EMT) simulation tools (e.g. PSCAD, EMTP) can simulate detailed IBR controls, but for large power systems with many IBRs, slow simulation speeds severely impede the ability to study dynamic events [3] [4]. Massively paralleling simulations using high-performance computing (HPC) can help address this, especially now that cloud-based HPC capability is widely available, but today s EMT tools are not HPC-compatible, and parallelization of dynamic simulation solvers is not trivial because each region can dynamically affect the others. Thus, dynamic simulation of grids with very large numbers of IBRs potentially poses a barrier to the ongoing energy transition.

24 POWER TRANSMISSION AND DISTRIBUTION↗

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↗

Learning effective physical laws for generating cosmological hydrodynamics with Lagrangian deep learning

The goal of generative models is to learn the intricate relations between the data to create new simulated data, but current approaches fail in very high dimensions. When the true data-generating process is based on physical processes, these impose symmetries and constraints, and the generative model can be created by learning an effective description of the underlying physics, which enables scaling of the generative model to very high dimensions. In this work, we propose Lagrangian deep learning (LDL) for this purpose, applying it to learn outputs of cosmological hydrodynamical simulations. The model uses layers of Lagrangian displacements of particles describing the observables to learn the effective physical laws. The displacements are modeled as the gradient of an effective potential, which explicitly satisfies the translational and rotational invariance. The total number of learned parameters is only of order 10, and they can be viewed as effective theory parameters. Further, we combine N-body solver fast particle mesh (FastPM) with LDL and apply it to a wide range of cosmological outputs, from the dark matter to the stellar maps, gas density, and temperature. The computational cost of LDL is nearly four orders of magnitude lower than that of the full hydrodynamical simulations, yet it outperforms them at the same resolution. We achieve this with only of order 10 layers from the initial conditions to the final output, in contrast to typical cosmological simulations with thousands of time steps. This opens up the possibility of analyzing cosmological observations entirely within this framework, without the need for large dark-matter simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Power Module Precursors and Prognostics

The ongoing work involves 1) the identification and evaluation of various electric signals as indicators to module failure and precursors to aging, and 2) the development of corresponding detection techniques. In the current phase, 1) the increase of on-state resistance is selected as the precursor to the power cycle aging of silicon carbide power MOSFET modules based on review and rest data. 2) A detection method leveraging the gate voltage fluctuations magnitude during commutation transients to inform the on-state resistance increase is investigated which has the benefit of being able to be integrated into the gate driver. A dedicated detection circuit for implementation is designed and simulated and is currently in fabrication.

ADVANCED PROPULSION SYSTEMS↗