Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulators”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Finite time step and spatial grid effects in δf simulation of warm plasmas

This paper introduces a technique for analyzing time integration methods used with the particle weight equations in δf method particle-in-cell (PIC) schemes. The analysis applies to the simulation of warm, uniform, periodic or infinite plasmas in the linear regime and considers the collective behavior similar to the analysis performed by Langdon for full-f PIC schemes [1,2]. We perform both a time integration analysis and spatial grid analysis for a kinetic ion, adiabatic electron model of ion acoustic waves. An implicit time integration scheme is studied in detail for δf simulations using our weight equation analysis and for full-f simulations using the method of Langdon. It is found that the δf method exhibits a CFL-like stability condition for low temperature ions, which is independent of the parameter characterizing the implicitness of the scheme. The accuracy of the real frequency and damping rate due to the discrete time and spatial schemes is also derived using a perturbative method. The theoretical analysis of numerical error presented here may be useful for the verification of simulations and for providing intuition for the design of new implicit time integration schemes for the δf method, as well as understanding differences between δf and full-f approaches to plasma simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Predicting Reaction Products and Automating Reactive Trajectory Characterization in Molecular Simulations with Support Vector Machines

A machine learning-based methodology for the prediction of chemical reaction products, along with automated elucidation of mechanistic details via phase space analysis of reactive trajectories, is introduced using low-dimensional heuristic models and then applied to ab initio computer simulations of the photodissociation of acetaldehyde, an important chemical system in atmospheric chemistry. Our method is centered around training Support Vector Machines (SVMs) to identify optimal separatrices that delineate the regions of phase space that lead to different chemical reaction products. In contrast to the more common “black box” type machine learning methodologies for analyzing chemical simulation data, this SVM-based methodology allows for mechanistic insight to be gleaned from further analysis of the positioning of the phase space points used to train the SVM with respect to the separatrices. For example, a pair of phase space points that are in close proximity to each other but on opposite sides of a separatrix may be situated on opposite sides of a transition state, while phase space points occurring early in a simulation that are distant from a separatrix are likely to belong to trajectories that are highly biased toward the product state associated with the basin of attraction to which they belong. In addition to inferring mechanistic details about multiple-pathway chemical reactions, our method can also be used to increase reactive trajectory sampling efficiency in molecular simulations via rejection sampling, with trajectories leading to undesired product states being identified and terminated early in the simulation rather than being carried to completion.

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

Systematic and simulation-free coarse graining of homopolymer melts: A structure-based study

We propose a systematic and simulation-free strategy for coarse graining of homopolymer melts, where each chain of N{sub m} monomers is uniformly divided into N segments, with the spatial position of each segment corresponding to the center-of-mass of its monomers. We use integral-equation theories suitable for the study of equilibrium properties of polymers, instead of many-chain molecular simulations, to obtain the structural and thermodynamic properties of both original and coarse-grained (CG) systems, and quantitatively examine how the effective pair potentials between CG segments and the thermodynamic properties of CG systems vary with N. Our systematic and simulation-free strategy is much faster than those using many-chain simulations, thus effectively solving the transferability problem in coarse graining, and provides the quantitative basis for choosing the appropriate N-values. It also avoids the problems caused by finite-size effects and statistical uncertainties in many-chain simulations. Taking the simple hard-core Gaussian thread model [K. S. Schweizer and J. G. Curro, Chem. Phys. 149, 105 (1990)] as the original system, we demonstrate our strategy applied to structure-based coarse graining, which is quite general and versatile, and compare in detail the various integral-equation theories and closures for coarse graining. Our numerical results show that the effective CG potentials for various N and closures can be collapsed approximately onto the same curve, and that structure-based coarse graining cannot give thermodynamic consistency between original and CG systems at any N < N{sub m}.

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

Particle-in-cell simulations of electron beam control using an inductive current divider

Kinetic, time-dependent, electromagnetic, particle-in-cell simulations of the inductive current divider are presented. The inductive current divider is a passive method for controlling the trajectory of an intense, hollow electron beam using a vacuum structure that inductively splits the beam’s return current. The current divider concept was proposed and studied theoretically in a previous publication [Phys. Plasmas 22, 023107 (2015)] A central post carries a portion of the return current (I 1 ) while the outer conductor carries the remainder (I 2 ) with the injected beam current given by I b =I 1 +I 2 . The simulations are in agreement with the theory which predicts that the total force on the beam trajectory is proportional to (I 2 -I 1 ) and the force on the beam envelope is proportional to I b . For a fixed central post, the beam trajectory is controlled by varying the outer conductor radius which changes the inductance in the return-current path. The simulations show that the beam emittance is approximately constant as the beam propagates through the current divider to the target. As a result, independent control over both the current density and the beam angle at the target is possible by choosing the appropriate return-current geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

ON THE ROAD TO MORE REALISTIC GALAXY CLUSTER SIMULATIONS: THE EFFECTS OF RADIATIVE COOLING AND THERMAL FEEDBACK PRESCRIPTIONS ON THE OBSERVATIONAL PROPERTIES OF SIMULATED GALAXY CLUSTERS

Flux-limited X-ray surveys of galaxy clusters show that clusters come in two roughly equally proportioned varieties: 'cool core' clusters (CCs) and non-'cool core' clusters (NCCs). In previous work, we have demonstrated using cosmological N-body + Eulerian hydrodynamic simulations that NCCs are often consistent with early major merger events that destroy embryonic CCs. In this paper we extend those results and conduct a series of simulations using different methods of gas cooling and of energy and metal feedback from supernovae, where we attempt to produce a population of clusters with realistic central cooling times, entropies, and temperatures. We find that the use of metallicity-dependent gas cooling is essential to prevent early overcooling, and that adjusting the amount of energy and metal feedback can have a significant impact on observable X-ray quantities of the gas. We are able to produce clusters with more realistic central observable quantities than have previously been attained. However, there are still significant discrepancies between the simulated clusters and observations, which indicates that a different approach to simulating galaxies in clusters is needed. We conclude by looking toward a promising subgrid method of modeling galaxy feedback in clusters that may help to ameliorate the discrepancies between simulations and observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-D-BRF-06: A Brachytherapy Simulator with Realistic Haptic Force Feedback and Real-Time Ultrasounds Image Simulation for Training and Teaching

Purpose: Surgical procedures require dexterity, expertise and repetition to reach optimal patient outcomes. However, efficient training opportunities are usually limited. This work presents a simulator system with realistic haptic force-feedback and full, real-time ultrasounds image simulation. Methods: The simulator is composed of a custom-made Linear-DELTA force-feedback robotic platform. The needle tip is mounted on a force gauge at the end effector of the robot, which responds to needle insertion by providing reaction forces. 3D geometry of the tissue is using a tetrahedral finite element mesh (FEM) mimicking tissue properties. As the needle is inserted/retracted, tissue deformation is computed using a mass-tensor nonlinear visco-elastic FEM. The real-time deformation is fed to the L-DELTA to take into account the force imparted to the needle, providing feedback to the end-user when crossing tissue boundaries or needle bending. Real-time 2D US image is also generated synchronously showing anatomy, needle insertion and tissue deformation. The simulator is running on an Intel I7 6- core CPU at 3.26 MHz. 3D tissue rendering and ultrasound display are performed on a Windows 7 computer; the FEM computation and L-DELTA control are executed on a similar PC using the Neutrino real-time OS. Both machines communicate through an Ethernet link. Results: The system runs at 500 Hz for a 8333-tetrahedron tissue mesh and a 100-node angular spring needle model. This frame rate ensures a relatively smooth displacement of the needle when pushed or retracted (±20 N in all directions at speeds of up to 2 m/s). Unlike commercially-available haptic platforms, the oblong workspace of the L-DELTA robot complies with that required for brachytherapy needle displacements of 0.1m by 0.1m by 0.25m. Conclusion: We have demonstrated a real-life, realistic brachytherapy simulator developed for prostate implants (LDR/HDR). The platform could be adapted to other sites or training for other types of needle-based procedures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Simulated Polarimetric Fields of Ice Vapor Growth Using the Adaptive Habit Model. Part I: Large-Eddy Simulations

Abstract The bulk adaptive habit model (AHM) explicitly predicts ice particle aspect ratio, improving the representation of microphysical processes and properties, including ice–liquid-phase partitioning. With the unique ability to predict ice particle shape and density, the AHM is combined with an offline forward operator to produce fields of simulated polarimetric variables. An evaluation of AHM-forward-simulated dual-polarization radar signatures in an idealized Arctic mixed-phase cloud is presented. Interpretations of those signatures are provided through microphysical model output using the large-eddy simulation mode of the Weather Research and Forecasting Model. Vapor-grown ice properties are associated with distinct observable signatures in polarimetric radar variables, with clear sensitivities to the simulated ice particle properties, including ice number, size, and distribution shape. In contrast, the liquid droplet number has little influence on both polarimetric and microphysical variables in the case presented herein. Polarimetric quantities are sensitive to the dominating crystal habit type in a volume, with enhancements for aspect ratios much lower or higher than unity. This synthesis of a microphysical model and a polarimetric forward simulator is a first step in the evaluation of detailed AHM microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Simulations of the pairwise kinematic Sunyaev-Zel'dovich signal

The pairwise kinematic Sunyaev–Zel’dovich (kSZ) signal from galaxy clusters is a probe of their line of sight momenta, and thus a potentially valuable source of cosmological information. In addition to the momenta, the amplitude of the measured signal depends on the properties of the intracluster gas and observational limitations such as errors in determining cluster centers and redshifts. In this work, we simulate the pairwise kSZ signal of clusters at z<1, using the output from a cosmological N-body simulation and including the properties of the intracluster gas via a model that can be varied in post-processing. We find that modifications to the gas profile due to star formation and feedback reduce the pairwise kSZ amplitude of clusters by ∼50%, relative to the naive “gas traces mass” assumption. We demonstrate that miscentering can reduce the overall amplitude of the pairwise kSZ signal by up to 10%, while redshift errors can lead to an almost complete suppression of the signal at small separations. We confirm that a high-significance detection is expected from the combination of data from current generation, high-resolution cosmic microwave background experiments, such as the South Pole Telescope, and cluster samples from optical photometric surveys, such as the Dark Energy Survey. Furthermore, we forecast that future experiments such as Advanced ACTPol in conjunction with data from the Dark Energy Spectroscopic Instrument will yield detection significances of at least 20σ, and up to 57σ in an optimistic scenario. Our simulated maps are publicly available at http://www.hep.anl.gov/cosmology/ksz.html.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A numerical simulation of a dry-out process for CO 2 sequestration in heterogeneous deep saline aquifers: Original Research Article: Dry-out problem: Multifactorial numerical simulation

Salt precipitation can be generated near the injection well when dry supercritical CO 2 (scCO 2 ) is injected into deep saline aquifers. Traditional laboratory experiments and numerical simulations focus on the generation of salt precipitation and its impact on core permeability. Current laboratory experimental studies have shown that the heterogeneity of the core and water content on the surface of the particles plays an important role in the formation of salt precipitation. The aim of this study is to explore the effects of brine back-flow, heterogeneity, and low water content on salt precipitation. Numerical simulations were performed using core-scale and site-scale models. Simulation results showed the following: 1. Back-flow plays an important role in the salt accumulation near the well – an homogeneous laboratory experiment used with a small core may not provide a full picture of the back-flow phenomenon. 2. The continuous dry-out process reduces the water content near the injection well to zero. As a result, the full range of saturation flow equations should be used to simulate the dry-out process accurately. 3. The heterogeneity of the matrix not only affects the distribution of salt precipitation but also increases the amount of salt precipitation near the well. Homogeneous simulation underestimates the amount of salt precipitation. 4. The reduction in porosity due to salt precipitation could extend the dry-out process in the low permeability zone, further expanding the salt precipitation area.

58 GEOSCIENCES↗
DOE OSTI2019

Assessing the mechanisms governing the daytime evolution of marine stratocumulus using large‐eddy simulation

Large‐eddy simulation (LES) output for a case of thin stratocumulus off the coast of California is examined in a mixed‐layer analysis framework to identify the specific mechanisms responsible for governing the evolution of the cloud system. An equation for cloud‐base height tendency isolates the individual cloud‐modulating mechanisms that control the evolution of boundary‐layer liquid‐water static energy ( S l ) and total water mixing ratio ( q T ). With a suitable spin‐up procedure, the control simulation performs admirably compared with observed estimates of liquid water content, vertical velocity variance, and radiative fluxes sampled during an aircraft field campaign. Investigation of the cloud response to various environmental forcing scenarios was addressed through a suite of sensitivity simulations, including variations in subsidence velocity, surface fluxes, wind shear near the inversion, and radiative forcing. In the control simulation, rising cloud‐base tendencies are associated with entrainment warming/drying and short‐wave absorption, whereas lowering cloud‐base tendencies are driven by long‐wave cooling. Even in the presence of substantial afternoon solar heating, entrainment fluxes remained active. The thin cloud demonstrated unexpected resiliency, with mixed‐layer analysis indicating that, as the short‐wave flux decreases later in the afternoon, the relative contribution of long‐wave cooling often becomes large enough to offset entrainment warming/drying and result in a reversal of cloud‐base tendency. The evolution of cloud‐base tendency is found to be insensitive to the net radiative flux divergence for most of the simulations (liquid water path ranging from ~10–50 g/m 2 ). Error analysis in comparison with LES S l and q T budgets suggests that our method of entrainment flux calculation could be improved by a more complete understanding of entrainment‐layer physics.

McMichael, Lucas A.↗
DOE OSTI2022

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulation of the BEAVRS benchmark using VERA

The Virtual Environment for Reactor Applications (VERA) is being developed by the Consortium for the Advanced Simulation of Light Water Reactors (CASL). Here, VERA is used to perform the Benchmark for Evaluation and Validation of Reactor Simulations (BEAVRS), which provides two cycles’ worth of operating power history, along with a full, detailed description of the geometry, and measured data. Cycle 1 and 2 are simulated with VERA and the results are compared to the measured zero power physics tests, critical boron concentration, and flux maps.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗