Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 523 records · Page 29

Predicting Wind Loading and Instability in Solar Tracking PV Arrays

Wind loading and the fluctuating pressure loads it creates on PV panel surfaces are associated with multiple degradation mechanisms and failures. Modest wind speeds create reversing loads that can initiate cell cracks and weather cracked cells. Stronger wind speeds and extreme weather events can lead to larger scale forces and the aerodynamic instability known as torsional galloping. All these effects are dependent on the complex coupling between wind speed, panel orientation, and a myriad of other hardware and site-specific factors. In this work, we present the latest developments from our work to build an open-source, high-performance computing (HPC) fluid dynamics solver to predict and mitigate these effects. This simulation package allows users to easily specify different array layouts, solar-tracking angles, panel geometries, and weather conditions before automatically generating a refined computational mesh and solving for the unsteady loading on each panel surface. Small domains (e.g., a single panel row in isolation) can be solved on a modern laptop, while larger domains or very high-fidelity studies can be solved on distributed or HPC resources with minimal modifications to the underlying problem specification. We present preliminary case studies obtained using this simulation package and highlight how increased wind speeds combined with sub-optimal tracking angles can exacerbate degradation drivers.

aerodynamics↗

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transmission-Distribution Dynamic Co-Simulation of Electric Vehicles Providing Grid Frequency Response

This paper investigates the impacts of electric vehicles (EVs) on power system frequency regulation based on an open-source transmission-and-distribution (T&D) dynamic co-simulation framework. The development of a EV dynamic model based on appropriate WECC dynamic model is introduced first, then the T&D dynamic co-simulation platform is described. The advantage of the overall platform is that DERs such as distributed photovoltaic (DPV) and EV are modeled explicitly in both transmission and distribution simulators for frequency dynamics and voltage, respectively. In the case studies, EVs frequency responses (i.e., primary and/or secondary) are simulated after the system is exposed to N-1 contingency such as a generation trip. Various EV frequency regulation participation strategies are investigated as well to study their impacts on system frequency response. The studies shows that EVs have potential capability to provide effective frequency regulation services.

electric vehicle↗

Transmission-Distribution Dynamic Co-Simulation of Electric Vehicles Providing Grid Frequency Response: Preprint

This paper investigates the impacts of electric vehicles (EVs) on power system frequency regulation based on an open-source transmission-and-distribution (T&D) dynamic co-simulation framework. The development of a EV dynamic model based on appropriate WECC dynamic model is introduced first, then the T&D dynamic co-simulation platform is described. The advantage of the overall platform is that DERs such as distributed photovoltaic (DPV) and EV are modeled explicitly in both transmission and distribution simulators for frequency dynamics and voltage, respectively. In the case studies, EVs frequency responses (i.e., primary and/or secondary) are simulated after the system is exposed to N-1 contingency such as a generation trip. Various EV frequency regulation participation strategies are investigated as well to study their impacts on system frequency response. The studies shows that EVs have potential capability to provide effective frequency regulation services.

electric vehicle↗

Impact of Floating PV (FPV) System on Water Evaporation: Case Study of an Open Loop Reservoir in North Carolina

Fresh water scarcity is becoming one of the important concerns worldwide. This situation of water scarcity is further exacerbated by the increase in water evaporation due to climate change. This manuscript will address this concern by considering the impact of floating photovoltaics (FPVs) systems on water evaporation. The installation of PV modules on the water bodies can help reduce the water evaporation rate from reservoirs by reducing the evaporating area of the water surface. This work evaluates the effect of FPV on water evaporation rates in the TuckerTown reservoir located in North Carolina state of United States. Meteorological data, such as mean temperature and dew-point, are used to calculate the rate of the water evaporation by Linacre method. The simulation results for the TuckerTown case study showed that the aggregated avoided evaporated volume for one year is 9 million gallons.

14 SOLAR ENERGY↗

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spent Crystalline Silicotitanate Storage Study—Post AP-105DF Processing

A Tank-Side Cesium Removal (TSCR) system is under development by Washington River Protection Solutions to support initial production of immobilized low-activity waste (LAW) by feeding Hanford tank supernate from tank farms to the Hanford Waste Treatment and Immobilization Plant (WTP) LAW Facility. Tank waste supernate will be filtered to remove suspended solids and then Cs will be removed by processing it through crystalline silicotitanate (CST) ion exchange media manufactured by Honeywell UOP, LLC. The Cs-loaded CST columns will be stored indefinitely, with a goal of eventual CST removal and treatment. Thus, the spent CST needs to be recoverable after undetermined storage time. Previous testing with AP-105 simulant showed that rinsing the CST bed with 3 bed volumes (1.4 apparatus volumes [AVs]) of 0.1 M NaOH resulted in a dried bed that maintained flow characteristics indicative of ease of recovery. This study explored the intermediate conditions (between feed dried in place and the 1.4 AVs of 0.1 M NaOH rinse) to evaluate CST bed properties after: 1) stoppage with feed in place; 2) stoppage after draining feed; 3) stoppage after 0.7 AV of 0.1 M NaOH rinse through column; 4) stoppage after 1.5 AVs of 0.1 M NaOH rinse through column. Post processing, each column was heated at 50 °C for 19 days under pseudo-storage conditions to simulate the expected dried and stored CST bed conditions. Testing was conducted at the small scale (12-mL bed volume); actual, Cs-depleted, AP-105 tank waste was used as the feed. Post-dried CST bed physical properties (angle of repose and penetration depth) were measured to evaluate how CST moved and flowed. All process stop-conditions resulted in a solidified CST bed except for the final condition, 1.5 AVs of 0.1 M NaOH rinse. At this small scale, the three CST beds presented an issue for retrievability after the short storage period (19 days at 50 °C). The latter case confirmed the results from simulant testing. The testing was intended to provide a preliminary assessment of issues that may arise from desiccation of CST during storage with the indicated salt solutions in place. Since these were small-scale tests, the processing system did not scale to full scale conditions exactly; however, the tests did provide insight into the impact on the dried and stored CST bed after stopping processing at an earlier step (upset condition) than normal. These results indicate that if an upset condition occurs at TSCR, a dilute hydroxide rinse should be considered before the CST dries from internal heating.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Enhancements to Overset Methods for Achieving Higher Accuracy and Solution Convergence: Preprint

Overset methods, when utilized within a multi-solver framework, provide a very flexible and powerful toolset for analyzing moving body problem. However, many overset framework suffer from several deficiencies such as (1) convergence degradation of the non-linear iterative schemes (2) difficulty in achieving higher-order accuracy and (3) high computational costs over long simulation times. In this work, we explore approaches that mitigate these deficiencies and provide a pathway for improved accuracy and fidelity of overset methods. A high-order, discontinuous Galerkin solver is integrated into an existing overset framework and applied to an unsteady ROBIN fuselage case. For convergence acceleration, an Overset-GMRES linear solver with a timestep controller is applied to a DLR-F6 aircraft steady simulation and an unsteady, notional rotor-fuselage. Finally the same rotor-fuselage case as well as the JVX rotor/wing case is simulated using a reduced-order aerodynamic model.

aerospace↗

Simulating and benchmarking neutron total scattering instrumentation from inception of events to reduced and fitted data

In the design and realization of modern neutron scattering instrumentation, particularly when designing beamline concepts from the ground up, it is desirable to fully benchmark against realistically simulated data. This is especially true for total scattering beamlines, where the future deliverable data is to be analysed in both reciprocal- and real-space representations, and needs must be carefully balanced to ensure sufficient range, resolution and flux of the instrument. An approach to optimize the design of neutron scattering instrumentation via a workflow including ray-tracing simulations, event-based data reduction, heuristic analysis and fitting against realistically simulated spectra is demonstrated here. The case of the DISCOVER beamline concept at the Spallation Neutron Source is used as an example. The results of the calculations are benchmarked through simulation of existing instrumentation and subsequent direct comparison with measured data. Finally, on the basis of the validated models, the ability to explore design characteristics for future beamline concepts or future instrument improvements is demonstrated through the examples of detector tube size and detector layout.

47 OTHER INSTRUMENTATION↗

On the Parameterization of Convective Downdrafts for Marine Stratocumulus Clouds

The role of nonlocal transport on the development and maintenance of marine stratocumulus (Sc) clouds in coarse-resolution models is investigated, with a special emphasis on the downdraft contribution. A new parameterization of cloud-top-triggered downdrafts is proposed and validated against large-eddy simulation (LES) for two Sc cases. The applied nonlocal mass-flux scheme is part of the stochastic multiplume eddy-diffusivity/mass-flux (EDMF) framework decomposing the turbulent transport into local and nonlocal contributions. The complementary local turbulent transport is represented with the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme. This EDMF version has been implemented in the Weather Research and Forecasting (WRF) single-column model (SCM) and tested for three model versions: without mass flux, with updrafts only, and with both updrafts and downdrafts. In the LES, the downdraft and updraft contributions to the total heat and moisture transport are comparable and significant. The WRF SCM results show a good agreement between the parameterized downdraft turbulent transport and LES. While including updrafts greatly improves the modeling of Sc clouds over the simulation without mass flux, the addition of downdrafts is less significant, although it helps improve the moisture profile in the planetary boundary layer.

Wu, Elynn↗

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗

Validation of the imaging neutral particle analyzer in nearly MHD quiescent plasmas using injected beam ions on DIII-D

The imaging neutral particle analyzer (INPA) is a scintillator-based diagnostic that provides energy and radially resolved measurements of confined fast ions. To verify operation, understanding, and modeling of the diagnostic performance, distributions of deuterium beam ions in low density, nearly MHD-quiescent plasmas are measured on the DIII-D tokamak using the INPA. Here, the distribution is altered by changing the pitch-angle scattering rate with electron cyclotron heating (ECH). Experimental images are mostly reproduced with simulation of core fast ions charge exchanging with injected beam and edge cold neutrals in classical cases. Image comparisons between the experiment and simulation show agreement with less than 25% discrepancies near the beam injection energy of 50 keV. Other strong passive signals, observed in a localized region on the image, are suggested to come from charge exchange between cold neutrals and fast ions around the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutrino heating in 1D, 2D, and 3D core-collapse supernovae: characterizing the explosion of high-compactness stars

Massive stars can end their lives with a successful supernova explosion (leaving behind a neutron star or, more rarely, a black hole), or a failed explosion that leaves behind a black hole. The density structure of the pre-collapse progenitor star already encodes much of the information regarding the outcome and properties of the explosion. However, the complexity of the collapse and subsequent shock expansion phases prevents drawing a straightforward connection between the pre-collapse and post-explosion properties. In order to derive such a connection several explodability studies have been performed in recent years. However, different studies can predict different explosion outcomes. In this article, we show how compactness, which is related to the average density of the star’s core, has an important role in determining the efficiency of neutrino heating, and therefore the outcome of the explosion. Commonly, high-compactness progenitors are assumed to yield failed explosions, due to their large mass accretion rates, preventing the shock from expanding. We show by analysing ~150 2D flash and F ornax simulations and 20 3D F ornax simulations that this is not the case. Instead, due to the rapid increase of neutrino heating with compactness, high-compactness progenitors lead to successful shock revival. We also show that 1D+ simulations that include v-driven convection using a mixing-length theory approach correctly reproduce this trend. Finally, we compare 1D+ models, which we show can reproduce some aspects of multi-D simulations with reasonable accuracy, with other widely used 1D models in the literature.

neutrinos↗

Quantum Computing for High-Energy Physics: State of the Art and Challenges

Quantum computers offer an intriguing path for a paradigmatic change of computing in the natural sciences and beyond, with the potential for achieving a so-called quantum advantage—namely, a significant (in some cases exponential) speedup of numerical simulations. The rapid development of hardware devices with various realizations of qubits enables the execution of small-scale but representative applications on quantum computers. In particular, the high-energy physics community plays a pivotal role in accessing the power of quantum computing, since the field is a driving source for challenging computational problems. This concerns, on the theoretical side, the exploration of models that are very hard or even impossible to address with classical techniques and, on the experimental side, the enormous data challenge of newly emerging experiments, such as the upgrade of the Large Hadron Collider. In this Roadmap paper, led by CERN, DESY, and IBM, we provide the status of high-energy physics quantum computations and give examples of theoretical and experimental target benchmark applications, which can be addressed in the near future. Having in mind hardware with about 100 qubits capable of executing several thousand two-qubit gates, where possible, we also provide resource estimates for the examples given using error-mitigated quantum computing. The ultimate declared goal of this task force is therefore to trigger further research in the high-energy physics community to develop interesting use cases for demonstrations on near-term quantum computers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large-Eddy Simulations of Convection Initiation over Heterogeneous, Low Terrain

Abstract Large-eddy simulations are conducted to investigate and physically interpret the impacts of heterogeneous, low terrain on deep-convection initiation (CI). The simulations are based on a case of shallow-to-deep convective transition over the Amazon River basin, and use idealized terrains with varying levels of ruggedness. The terrain is designed by specifying its power-spectral shape in wavenumber space, inverting to physical space assuming random phases for all wave modes, and scaling the terrain to have a peak height of 200 m. For the case in question, these modest terrain fields expedite CI by up to 2–3 h, largely due to the impacts of the terrain on the size of, and subcloud support for, incipient cumuli. Terrain-induced circulations enhance subcloud kinetic energy on the mesoscale, which is realized as wider and longer-lived subcloud circulations. When the updraft branches of these circulations breach the level of free convection, they initiate wider and more persistent cumuli that subsequently undergo less entrainment-induced cloud dilution and detrainment-induced mass loss. As a result, the clouds become more vigorous and penetrate deeper into the troposphere. Larger-scale terrains are more effective than smaller-scale terrains in promoting CI because they induce larger enhancements in both the width and the persistence of subcloud updrafts.

54 ENVIRONMENTAL SCIENCES↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An attempt to simulate laser-induced all-optical spin switching in a crystalline ferrimagnet

Interest in all-optical spin switching (AOS) is growing rapidly. The recent discovery of AOS in Mn 2 RuGa provides a much needed clean case of crystalline ferrimagnets for theoretical simulations. Here, we attempt to simulate it using the state-of-the-art first-principles method combined with the Heisenberg exchange model. We first compute the spin moments at two inequivalent manganese sites and then feed them into our model Hamiltonian. We employ an ultrafast laser pulse to switch the spins. We find that there is a similar optimal laser field amplitude to switch spins. However, we find that the exchange interaction has a significant effect on the system switchability. Weakening the exchange interaction could make the system unswitchable. This provides a crucial insight into the switching mechanism in ferrimagnets.

36 MATERIALS SCIENCE↗