Structured and Unstructured Entropy-Stable High-Order Methods for Simulating High-Speed Compressible Turbulent Flows.
Abstract not provided.
Fisher, Travis C.↗
SEARCH · Engineering Papers
Results for “SIMULATION METHODS”
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.
Reduced-Order Atomistic Method for Simulating High Dose Irradiation in Metals.
Abstract not provided.
Chen, Elton Y↗
Assessment of Meshless Methods for Simulation of Empty Room Closure at the Waste Isolation Pilot Plant.
Abstract not provided.
Reedlunn, Benjamin↗
Comparison of Electromagnetic Transient and Phasor-based Simulation for the Stability of Grid-Forming-Inverter-based Microgrids
The increasing interest in studying feeder-level, inverter-based microgrids requires appropriate simulation methods that can effectively simulate the dynamics of distribution feeders that typically have thousands of nodes. To examine the accuracy of phasor-based simulation for the stability of feeder-level, inverter-based microgrids, a comparative study between electromagnetic transient (EMT) and phasor-based simulation was conducted. Case studies of a two-inverter small microgrid and a modified version of the IEEE 34-Node Test System that has seven grid-forming inverters show that phasor-based simulation agrees well with EMT simulation when the P-f droop gains of grid-forming inverter controllers are within normal ranges. However, under extremely high P-f droop gains, phasor-based simulation can result in over-optimistic stability results. This study provides insight into the applicability of phasor-based simulation for the stability of feeder-level, inverter-based microgrids.
Du, Wei↗
Open-top microfluidic devices and methods for simulating a function of a tissue
A device for simulating a function of a tissue includes a first structure, a second structure, and a membrane. The first structure defines a first chamber. The first chamber includes a matrix disposed therein and an opened region. The second structure defines a second chamber. The membrane is located at an interface region between the first chamber and the second chamber. The membrane includes a first side facing toward the first chamber and a second side facing toward the second chamber. The membrane separates the first chamber from the second chamber.
Varone, Antonio↗
Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods
Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Qualification of Numerical and Semi-Analytical Modeling Codes for Thermal Design of HLW Repository - 18455
Disposal in a geological repository is a viable option for safely isolating high-level radioactive waste. Thermal processes induced by decay heat from the radioactive waste may impact the integrity of the engineered barrier system, and may lead to fluid flow and transport; and the associated migration of radionuclides. This study looks at an analysis of thermal conduction for the emplacement of high-level waste in a repository located in bedded salt. The simulations used analytical and numerical methods to study conductive heat transport in a repository. This study concentrated on heat flow in a central waste package with contributions from adjacent waste packages and drifts. Of particular importance to this work is comparison of predictions of peak temperature at the drift wall and at the waste package surface using different solution methods. Simulation results of analytical and numerical modeling methods are compared. The results show that predictions of the codes were comparable under identical initial and boundary conditions. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Computational simulation of concurrent engineering for aerospace propulsion systems
Results are summarized of an investigation to assess the infrastructure available and the technology readiness in order to develop computational simulation methods/software for concurrent engineering. These results demonstrate that development of computational simulations methods for concurrent engineering is timely. Extensive infrastructure, in terms of multi-discipline simulation, component-specific simulation, system simulators, fabrication process simulation, and simulation of uncertainties - fundamental in developing such methods, is available. An approach is recommended which can be used to develop computational simulation methods for concurrent engineering for propulsion systems and systems in general. Benefits and facets needing early attention in the development are outlined.
Chamis, C. C.↗
Computational simulation for concurrent engineering of aerospace propulsion systems
Results are summarized for an investigation to assess the infrastructure available and the technology readiness in order to develop computational simulation methods/software for concurrent engineering. These results demonstrate that development of computational simulation methods for concurrent engineering is timely. Extensive infrastructure, in terms of multi-discipline simulation, component-specific simulation, system simulators, fabrication process simulation, and simulation of uncertainties--fundamental to develop such methods, is available. An approach is recommended which can be used to develop computational simulation methods for concurrent engineering of propulsion systems and systems in general. Benefits and issues needing early attention in the development are outlined.
Chamis, C. C.↗
Speeding up simulations by slowing down particles: Speed-limited particle-in-cell simulation
Based on the particle-in-cell (PIC) plasma simulation method, the speed-limited PIC (SLPIC) method delivers faster kinetic plasma simulation in cases where the particle distributions evolve slowly compared with the maximum stable PIC timestep. SLPIC thus offers more feasible, fully kinetic simulation in regimes that historically have required fluid approaches, such as magnetohydrodynamic (MHD), two-fluid, or Boltzmann electron treatments. In particular, SLPIC allows an explicit time advance with steps much larger than the inverse plasma frequency, avoiding the instability explicit PIC faces with large timesteps. SLPIC avoids this instability by slowing down fast particles (e.g., electrons) in a way that is rigorously underpinned by an approximate Vlasov equation; unlike MHD, two-fluid, and Boltzmann electron approaches, SLPIC does not fundamentally neglect any first-principles plasma physics, although the choices of the grid cell size, timestep, and number of macroparticles per cell naturally limit the physical phenomena that can be accurately represented. SLPIC can be implemented with minor modifications of a standard PIC code and does not require an implicit time advance. It enables large timesteps in first-principles kinetic plasma simulation of appropriately slow phenomena, and it can handle many of the same complications as PIC, such as boundary conditions and collisions. Finally, in an argon plasma sheath test problem, a SLPIC simulation achieved a speed-up of a factor of 160 over the corresponding PIC simulation, without loss of accuracy.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Shock melting method to determine melting curve by molecular dynamics: Cu, Pd, and Al
A melting simulation method, the shock melting (SM) method, is proposed and proved to be able to determine the melting curves of materials accurately and efficiently. The SM method, which is based on the multi-scale shock technique, determines melting curves by preheating and/or prepressurizing materials before shock. This strategy was extensively verified using both classical and ab initio molecular dynamics (MD). First, the SM method yielded the same satisfactory melting curve of Cu with only 360 atoms using classical MD, compared to the results from the Z-method and the two-phase coexistence method. Then, it also produced a satisfactory melting curve of Pd with only 756 atoms. Finally, the SM method combined with ab initio MD cheaply achieved a good melting curve of Al with only 180 atoms, which agrees well with the experimental data and the calculated results from other methods. It turned out that the SM method is an alternative efficient method for calculating the melting curves of materials.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
pnnl/mol_dgnn
Molecular dynamics simulations often rely on parameterized atomistic potentials, which can be inaccurate compared to ab initio simulation methods. However, dynamic ab initio simulations come with high computational cost. Deep learning techniques offer a method to obtain both high accuracy and low cost predictions of dynamics. Our method gives a prediction of the structural dynamics of a molecule accounting for the prior dynamics, giving high-accuracy predictions at low computational cost
Ashby, Michael↗
Hard radiation environment simulation
Test methods for simulating hard radiation environment of optical telescopes in earth orbit
Paul, F. W.↗
SU-G-IeP4-12: Performance of In-111 Coincident Gamma-Ray Counting: A Monte Carlo Simulation
Purpose: The decay of In-111 results in a non-isotropic gamma-ray cascade, which is normally imaged using a gamma camera. Creating images with a gamma camera using coincident gamma-rays from In-111 has not been previously studied. Our objective was to explore the feasibility of imaging this cascade as coincidence events and to determine the optimal timing resolution and source activity using Monte Carlo simulations. Methods: GEANT4 was used to simulate the decay of the In-111 nucleus and to model the gamma camera. Each photon emission was assigned a timestamp, and the time delay and angular separation for the second gamma-ray in the cascade was consistent with the known intermediate state half-life of 85ns. The gamma-rays are transported through a model of a Siemens dual head Symbia “S” gamma camera with a 5/8-inch thick crystal and medium energy collimators. A true coincident event was defined as a single 171keV gamma-ray followed by a single 245keV gamma-ray within a specified time window (or vice versa). Several source activities (ranging from 10uCi to 5mCi) with and without incorporation of background counts were then simulated. Each simulation was analyzed using varying time windows to assess random events. The noise equivalent count rate (NECR) was computed based on the number of true and random counts for each combination of activity and time window. No scatter events were assumed since sources were simulated in air. Results: As expected, increasing the timing window increased the total number of observed coincidences albeit at the expense of true coincidences. A timing window range of 200–500ns maximizes the NECR at clinically-used source activities. The background rate did not significantly alter the maximum NECR. Conclusion: This work suggests coincident measurements of In-111 gamma-ray decay can be performed with commercial gamma cameras at clinically-relevant activities. Work is ongoing to assess useful clinical applications.
60 APPLIED LIFE SCIENCES↗
Quantum Divide and Compute: Exploring the Effect of Different Noise Sources
Abstract Our recent work (Ayral et al. in Proceedings of IEEE computer society annual symposium on VLSI, ISVLSI, pp 138–140, 2020. 10.1109/ISVLSI49217.2020.00034 ) showed the first implementation of the Quantum Divide and Compute (QDC) method, which allows to break quantum circuits into smaller fragments with fewer qubits and shallower depth. This accommodates the limited number of qubits and short coherence times of quantum processors. This article investigates the impact of different noise sources—readout error, gate error and decoherence—on the success probability of the QDC procedure. We perform detailed noise modeling on the Atos Quantum Learning Machine, allowing us to understand tradeoffs and formulate recommendations about which hardware noise sources should be preferentially optimized. We also describe in detail the noise models we used to reproduce experimental runs on IBM’s Johannesburg processor. This article also includes a detailed derivation of the equations used in the QDC procedure to compute the output distribution of the original quantum circuit from the output distribution of its fragments. Finally, we analyze the computational complexity of the QDC method for the circuit under study via tensor-network considerations, and elaborate on the relation the QDC method with tensor-network simulation methods.
97 MATHEMATICS AND COMPUTING↗
Computational methods for digital simulation of continuous systems
Computational methods for digital simulation of continuous systems
Di Perna, R. A.↗
The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice
Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.
Geology↗
An Integrated Study of a Novel Thermal Coating for Nb-Based High Temperature Alloy
This report summarizes our recent works of ab initio density functional theory (DFT) method and molecular dynamics (MD) simulation on the interfaces between niobium substrate and coatings at atomic level. Potential oxidation barrier bond coat, Nb₂AlC and high entropy alloys, and top coat candidates were synthesized, characterized, and evaluated in our labs. The simulation methods, experimental validation techniques, achievements already reached, students and postdoc training, and future improvement are briefly introduced.
74 ATOMIC AND MOLECULAR PHYSICS↗