Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and Modeling Tools”

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 1,675 records · Page 93

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗

Data initiatives for ocean-driven melt of Antarctic ice shelves

Ocean-driven melt of Antarctic ice shelves is an important control on mass loss from the ice sheet, but is complex to study due to significant variability in melt rates both spatially and temporally. Here we assess the strengths and weakness of satellite and field-based observations as tools for testing models of ice-shelf melt. We discuss how the complementary use of field, satellite and model data can be a powerful but underutilised tool for studying melt processes. Finally, we identify some community initiatives working to collate and publish coordinated melt rate datasets, which can be used in future for validating satellite-derived maps of melt and evaluating processes in numerical simulations.

54 ENVIRONMENTAL SCIENCES↗

On the importance of excited state species in low pressure capacitively coupled plasma argon discharges

In the past three decades, first principles-based fully kinetic particle-in-cell Monte Carlo collision (PIC/MCC) simulations have been proven to be an important tool for the understanding of the physics of low pressure capacitive discharges. However, there is a long-standing issue that the plasma density determined by PIC/MCC simulations shows quantitative deviations from experimental measurements, even in argon discharges, indicating that certain physics may be missing in previous modeling of the low pressure radio frequency (rf) driven capacitive discharges. In this work, we report that the energetic electron-induced secondary electron emission (SEE) and excited state atoms play an important role in low pressure rf capacitive argon plasma discharges. The ion-induced secondary electrons are accelerated by the high sheath field to strike the opposite electrode and produce a considerable number of secondary electrons that lead to additional ionizing impacts and further increase of the plasma density. Importantly, the presence of excited state species even further enhances the plasma density via excited state neutral and resonant state photon-induced SEE on the electrode surface. The PIC/MCC simulation results show good agreement with the recent experimental measurements in the low pressure range (1–10 Pa) that is commonly used for etching in the semiconductor industry. At the highest pressure (20 Pa) and driving voltage amplitudes 250 and 350 V explored here, the plasma densities from PIC/MCC simulations considering excited state neutrals and resonant photon-induced SEE are quantitatively higher than observed in the experiments, requiring further investigation on high pressure discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas (Final Technical Report)

This document details the work accomplished during the period 04/15/2012-01/31/2020 for the collaborative project “Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas,” and serves as the final report on the University of California, San Diego part of this collaborative project for Department of Energy Office of Science - Fusion Energy Sciences Grant No. DE- SC0007783. This project focused on the development of analytic and computational tools, seeking to better understand the mutual interactions between small-scale temperature-gradient driven microturbulence and macroscopic magnetohydrodynamic (MHD) instabilities, particularly tearing modes, in magnetically confined plasmas relevant for the generation of fusion energy. This project advanced understanding through numerical simulations of different sets of model equations, corresponding to different microturbulence instabilities. Questions addressed included how the magnitudes, spatial structures, and scaling's of turbulent forcing's and drags impact magnetic island dynamics, as well as quantifying the “region of influence” of an island on the microturbulence dynamics. A key finding of the project is that the microturbulence fluctuations and associated cross-field particle and energy fluxes exhibit spatial asymmetries for sufficiently large islands. The results have improved our fundamental understanding of plasma dynamics, motivated new experimental investigations, and will help inform future fusion predictive modeling efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

Model Validation and Uncertainty Quantification on the KRUSTY Microreactor Design Using GRIFFIN Neutron Transport Code [Poster]

Argonne National Laboratory (ANL) and INL have developed a GRIFFIN steady state neutronics model for the multiphysics simulations of the Kilopower Reactor Using Sterling TechnologY (KRUSTY) microreactor in the Multiphysics Object Oriented Simulation Environment (MOOSE). The reliability of such deterministic neutronics models can be validated by comparing with computations from Monte Carlo codes (e.g. MCNP, SERPENT, OpenMC, Shift, etc). Furthermore, potential modeling/design improvements can be identified by incorporating uncertainty quantification (UQ), which can be performed by MOOSE’s Stochastic Tools Module (STM). KRUSTY is a prototype for a 5-kW thermal nuclear-powered space reactor. Its primary components consist of nuclear fuel, heat pipes, a control rod, a reflector, and the shielding. The fuel consists of 3 stacked U-7.65Mo cylinders with a hole in the center for the control rod. 8 liquid sodium heat pipes transfer fission energy from the solid fuel block to the Sterling power conversion system where the energy is extracted, and the cooled sodium flows back to the core via capillary action . The movable Boron Carbide control rod regulates the neutron population during startup or when a reactor temperature boost is needed . The beryllium oxide reflector is in 3 places in the reactor; it surrounds the core axially, it lies beneath the core on a platen, and it is present in the shim. The axial and lower reflectors rest on an adjustable stainless-steel platen that moves upward to cover the fuel and help the reactor reach criticality. Lastly, radial stainless steel surrounds the core offering protection from radiation exposure .

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Detecting Coastal Wetland Degradation by Combining Remote Sensing and Hydrologic Modeling

Sea-level rise and climate change stresses pose increasing threats to coastal wetlands that are vital to wildlife habitats, carbon sequestration, water supply, and other ecosystem services with global significance. However, existing studies are limited in individual sites, and large-scale mapping of coastal wetland degradation patterns over a long period is rare. Our study developed a new framework to detect spatial and temporal patterns of coastal wetland degradation by analyzing fine-scale, long-term remotely sensed Normalized Difference Vegetation Index (NDVI) data. Then, this framework was tested to track the degradation of coastal wetlands at the Alligator River National Wildlife Refuge (ARNWR) in North Carolina, United States, during the period from 1995 to 2019. We identified six types of coastal wetland degradation in the study area. Most of the detected degradation was located within 2 km from the shoreline and occurred in the past five years. Further, we used a state-of-the-art coastal hydrologic model, PIHM-Wetland, to investigate key hydrologic processes/variables that control the coastal wetland degradation. The temporal and spatial distributions of simulated coastal flooding and saltwater intrusion confirmed the location and timing of wetland degradation detected by remote sensing. The combined method also quantified the possible critical thresholds of water tables for wetland degradation. The remote sensing–hydrologic model integrated scheme proposed in this study provides a new tool for detecting and understanding coastal wetland degradation mechanisms. Our study approach can also be extended to other coastal wetland regions to understand how climate change and sea-level rise impact wetland transformations.

54 ENVIRONMENTAL SCIENCES↗

Modelling and Optimizing CSP-PV Hybrid Systems Using the Hybrid Optimization and Performance Platform (HOPP)

The hybridization of concentrating solar power (CSP) with thermal energy storage (TES), photovoltaics (PV), and electrochemical battery energy storage systems (BESS) has the potential to provide continuous, high-capacity-factor energy production at a lower cost than a PV-BESS or CSP with TES alone. Because of the system complexity of CSP technology, it is challenging to evaluate the technological and financial performance of a CSP-PV hybrid system without detailed modeling of annual operations. To address this challenge, we have developed a modeling framework for evaluating the performance and financial viability of CSP systems hybridized with PV and BESS technologies. This modeling effort incorporates CSP tower and trough configurations into an existing modeling tool recently developed by NREL, the Hybrid Optimization and Performance Platform (HOPP). This paper provides a brief overview of our methodology, as well as an example case study. CSP with TES hybridized with PV provides the best benefit-to-cost ratio compared to the other simulated technology combinations. However, for the conditions considered, this configuration only increases the benefit-to-cost ratio by about 1% compared to the CSP with TES configuration. The PV-BESS system provides the lowest benefit-to-cost ratio compared to the other configurations explored because of the relatively low capacity credit received by the system.

Hamilton, William (ORCID:0000000244654239)↗

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Splat Quenching as a High Throughput Rapid Solidification Testing Method for Developing Additive Manufacturing Alloys

Splat quenching as a high throughput rapid solidification screening tool is evaluated. Samples were made using known alloys (SS316, IN625, Ti-5553) to evaluate and establish the microstructures for a variety of alloy systems as well as the predicted cooling rates. The samples are inductively heated and levitated prior to being struck between two platens that produce a high contact pressure and a thin sample resulting in cooling rates on the order of 10 6 to 10 7 K/s. Process parameters were evaluated with analytical models in addition to numerical simulations to provide an effect of process variables such as substrate material, melt superheat, and platen velocity on the resulting solidification. The splat thickness was found to be controlled by platen velocity, melt superheat, as well as the feedstock volume. The sample thickness is the key controlling factor for varying the average cooling rate experienced by the splat quenched sample. The splat quenching techniques can reach regions of rapid solidification space that meet and exceed laser and electron beam techniques across the sample. In conclusion, the results of the study provide a useful foundation in understanding the splat quenching technique and its potential as a low effort tool to explore rapid solidification effects on alloys.

Alloy Design↗

Extended quantitative rescattering model for simulating high-order harmonic streaking spectra by synchronization of an intense IR laser and a time-delayed attosecond XUV pulse

In this work, we theoretically investigate the modulated high-harmonic generation (HHG) driven by an intense few-cycle infrared (IR) laser field and a weak extreme-ultraviolet (XUV) pulse at a delayed time. We establish an extended quantitative rescattering (EQRS) model to simulate the HHG streaking spectra, with the ideas of correcting the IR ionization and the transition from the ground to continuum states in the strong-field approximation. The EQRS model has an accuracy comparable to that from “exactly” solving the time-dependent Schrödinger equation (TDSE). We reveal that the fringes in the streaking spectra are caused by the interference between the attosecond XUV pulse and harmonics resulting from different recombination pathways under the intense IR laser. We then demonstrate that the XUV pulse can be accurately retrieved by treating the single-atom TDSE results or macroscopic propagation results as the “input” data. This work provides with a tool for efficiently simulating and thoroughly analyzing the XUV-assisted HHG, which could also enhance its capability for tracing the electron dynamics involved in the strong-field phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Razorback – A reactor transient analysis code for large rapid reactivity additions in a natural circulation research reactor

Research reactors play an important role in higher education, scientific research, and medical radioisotope production around the world. It is thus important to ensure the safety of facility workers and the public. This work presents a new reactor transient analysis code referred to as Razorback. The code has been developed for the evaluation of large rapid reactivity addition in research reactors. Its initial focus is the Annular Core Research Reactor (ACRR) at Sandia National Laboratories. Results have been validated against ACRR pulse operations. Razorback models the reactor kinetics, fuel element heat transfer, fuel element thermal expansion, and natural circulation coolant channel thermal-hydraulic response. Simulation results for ACRR pulse operations are shown to agree very well with operational data obtained from the ACRR. Finally, Razorback is expected to be a valuable tool for ACRR pulse performance prediction and ACRR reactor safety analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Benchmarking collective effects of electron interactions in a wiggler with OPAL-FEL

OPAL-FEL is a recently developed tool for the modeling of particle accelerators containing wigglers or undulators. It extends the well established 3D electrostatic particle-tracking code OPAL, by merging it with the finite-difference time-domain electromagnetic solver MITHRA. We present results of two benchmark cases where OPAL-FEL simulations are compared to experimental results. Both experiments concern electron beamlines where the longitudinal phase space is modulated with a short magnetic wiggler. Good agreement was found in both the space charge and radiation dominated regimes.

43 PARTICLE ACCELERATORS↗

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗

Plasma physics in strong-field regimes: Theories and simulations

In strong electromagnetic fields, unique plasma phenomena and applications emerge whose description requires recently developed theories and simulations [Y. Shi, Ph.D. thesis, Princeton University (2018)]. In the classical regime, to quantify effects of strong magnetic fields on three-wave interactions, a convenient formula is derived by solving the fluid model to the second order in general geometry. As an application, magnetic resonances are exploited to mediate laser pulse compression, whereby higher intensity pulses can be produced in wider frequency ranges, as confirmed by particle-in-cell simulations. In even stronger fields, relativistic-quantum effects become important, and a plasma model based on scalar quantum electrodynamics (QED) is developed which unveils observable corrections to Faraday rotation and cyclotron absorption in strongly magnetized plasmas. Beyond the perturbative regime, lattice QED is extended as a numerical tool for plasma physics, using which the transition from wakefield acceleration to electron-positron pair production is captured when laser intensity exceeds the Schwinger threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗

Machine Learning-Enabled Image Classification for Automated Electron Microscopy

Abstract Traditionally, materials discovery has been driven more by evidence and intuition than by systematic design. However, the advent of “big data” and an exponential increase in computational power have reshaped the landscape. Today, we use simulations, artificial intelligence (AI), and machine learning (ML) to predict materials characteristics, which dramatically accelerates the discovery of novel materials. For instance, combinatorial megalibraries, where millions of distinct nanoparticles are created on a single chip, have spurred the need for automated characterization tools. This paper presents an ML model specifically developed to perform real-time binary classification of grayscale high-angle annular dark-field images of nanoparticles sourced from these megalibraries. Given the high costs associated with downstream processing errors, a primary requirement for our model was to minimize false positives while maintaining efficacy on unseen images. We elaborate on the computational challenges and our solutions, including managing memory constraints, optimizing training time, and utilizing Neural Architecture Search tools. The final model outperformed our expectations, achieving over 95% precision and a weighted F-score of more than 90% on our test data set. This paper discusses the development, challenges, and successful outcomes of this significant advancement in the application of AI and ML to materials discovery.

Materials Science↗