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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2013

Defect formation by pristine indenter at the initial stage of nanoindentation

Nano-indentation is a sophisticated method to characterize mechanical properties of materials. This method samples a very small amount of material during each indentation. Therefore, this method is extremely useful to measure mechanical properties of nano-materials. The measurements using nanoindentation is very sensitive to the surface topology of the indenter and the indenting surfaces. The mechanisms involved in the entire process of nanoindentation require an atomic level understanding of the interplay between the indenter and the substrate. In this paper, we have used atomistic simulation methods with empirical potentials to investigate the effect of various types of pristine indenter on the defect nucleation and growth. Using molecular dynamics simulations, we have predicted the load-depth curve for conical, vickers, and sperical tip. The results are analyzed based on the coherency between the indenter tip and substrate surface for a fixed depth of 20 Å. The depth of defect nucleation and growth is observed to be dependent on the tip geometry. A tip with larger apex angle nucleates defects at a shallower depth. However, the type of defect generated is dependent on the crystalline orientation of the tip and substrate. For coherent systems, prismatic loops were generated, which released into the substrate along the close-packed directions with continued indentation. For incoherent systems, pyramidal shaped dislocation junctions formed in the FCC systems and disordered atomic clusters formed in the BCC systems. These defect nucleation and growth process provide the atomistic mechanisms responsible for the observed load-depth response during nanoindentation.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Impact of post-rainfall evaporation from porous roof tiles on building cooling load in subtropical China

Rainfall occurs frequently in subtropical regions of China, with the subsequent water evaporation from building roofs impacting the thermal performance and the energy consumption of buildings. In this paper, we proposed a novel simulation method using actual meteorological data to evaluate this impact. New features were developed in EnergyPlus to enable the simulation: (1) an evaporation latent heat flux source term was added to the heat balance equation of the external surface and (2) algorithms for the evaporative cooling module (ECM) were developed and implemented into EnergyPlus. The ECM experimental results showed good agreement with the simulated results. The ECM was used to assess the impact of evaporation from porous roof tiles on the cooling load of a one-floor building in subtropical China. The results show that the evaporation process decreased the maximal values of the external and internal roof surface temperatures by up to 6.4 °C and 3.2 °C, respectively, while the lower internal surface temperature decreased the room accumulated cooling load by up to 14.8% during the hot summer period. Lastly, the enhanced EnergyPlus capability can be used to evaluate the evaporative cooling performance of roofs with water-storage mediums, as well as to quantify their impact on building cooling loads.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

The magnetic effects on the energetic landscape of Fe-Cu alloy: A model Hamiltonian approach

Using ab initio calculations and the magnetic cluster expansion formula, a model with the magnetic effects on the energetic landscape of Fe–Cu was developed. The (partially) disordered local moment picture and Monte Carlo simulation method were combined to study the magnetically disordering effects on the thermodynamic properties in Fe–Cu alloys. The simulation results show nonlinear enhancements of Fe magnetism with increasing Cu concentration under varying magnetic ordering/disordering states, and the nonlinearity of the enhancement is more significant under the ferromagnetic state. It is also found that magnetism has subtle effects on the clustering of Cu because both the substitutional and binding energies could be affected but in different ways depending on the magnetic states. The results suggest that the alloying-induced enhancement of the Fe magnetism is a fundamental feature of Fe–Cu, which introduces difference on the concentration dependence of mixing, depending on the magnetic ordering/disordering states. In conjunction with the concentration dependence of short-range ordering, the enhancement of Fe magnetism also strongly affects the Curie temperatures in Fe–Cu.

Ab Initio Calculations, Fe-Cu, ultrasonic backscat↗

DOE OSTI2023

Regional infrasonic observations from surface explosions – influence of atmospheric variations and realistic terrain

SUMMARY A pair of 1-ton, conventional surface explosions were conducted at the Nevada National Security Site in the fall of 2020 producing seismoacoustic signatures observable hundreds of kilometres from the source location. Regional infrasonic observations include tropospheric ducting at large distances to the south, a wide stratospheric waveguide with signals observed more than 700 km to the east, and anomalous arrivals in the stratospheric shadow zone. Notable differences in propagation between the events are identified despite the explosions being conducted just two days apart due to a sharp temporal shift in the tropospheric winds as well as structural changes in the stratospheric winds. Propagation simulations of the two events have been completed using a combination of ray tracing and parabolic equation (PE) methods. Simulations have been conducted to quantify the impact of the temporal variations in the atmosphere as well as the influence of terrain on propagation. Temporal variations in reflection locations are found to produce notable changes in downrange propagation due to spatially varying terrain features. Finite frequency effects modelled by the PE are found to predict ensonification not included in corresponding 2D ray tracing simulations. Notable variations in predicted signal amplitude are found due to focusing by along-path and cross-path terrain gradients; though, the later of these is only modelled using fully 3D ray tracing analysis.

58 GEOSCIENCES↗

DOE OSTI2022

Towards string order melting of spin-1 particle chains in superconducting transmons using optimal control

Utilizing optimal control to simulate a model Hamiltonian is an emerging strategy that leverages the intrinsic physics of a device with digital quantum simulation methods. Here we evaluate optimal control for probing the nonequilibrium properties of symmetry-protected topological (SPT) states simulated with superconducting hardware. Assuming a tunable transmon architecture, we cast the evolution of these SPT states as a series of one- and two-site pulse optimization problems that are solved in the presence of leakage constraints. From the generated pulses, we classically simulate the time-dependent melting of the perturbed SPT string order across a six-site model with an average state infidelity of 10 -3 . The feasibility of these pulses as well as their efficient application indicate that high-fidelity simulations of string order melting are within reach of current quantum computing systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2011

Assessment of Computational Fluid Dynamics (CFD) Models for Shock Boundary-Layer Interaction

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and in general it was difficult to discern clear trends in the data. For the Reynolds Averaged Navier-Stokes methods the choice of turbulence model appeared to be the largest factor in solution accuracy. Large-eddy simulation methods produced error levels similar to RANS methods but provided superior predictions of normal stresses.

DeBonis, James R.↗

DOE OSTI2022

Phase-field simulations to inform nuclear fuel performance modeling

Software tools to simulate nuclear fuel performance at the engineering scale, such as Idaho National Laboratory (INL)’s BISON code, are increasingly relied upon in regulatory and economic decision-making. However, accurate results from these tools depends on the availability of materials parameters that are used as input. In recent years, atomistic and mesoscale simulation methods have emerged as a cost-effective, expedient means to obtain such input parameters. Phase-field simulations using INL’s Marmot application have been used to obtain microstructure-level parameters and to improve material models for fuel performance modeling using BISON. In this talk, recent examples of this process are given, including applications in UO2, U3Si2, and UZr fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1983

Numerical simulation of airfoil ice accretion

A mathematical model of glaze and rime ice accretion on a two-dimensional airfoil is presented. The model employs standard methods for calculating the flow field and cloud droplet trajectories. This data is used as input to a thermodynamic analysis of the ice accretion process which includes liquid runback. The structure of the model allows for dynamic updating of the droplet collection efficiency and heat and mass transfer processes as the accreted ice shape changes. This results in the improved shape prediction over previously described methods. Simulations made with the model are compared to experimental results obtained for a NACA 0012 airfoil in the NASA/Lewis Icing Research Tunnel. The agreement with experiment is found to be generally satisfactory for some cases, but a need for improvements is indicated by others. Suggestions for improvements are made.

Macarthur, C. D.↗

DOE OSTI2020

Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches

We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.

Tao, Xuecheng↗

DOE OSTI2016

Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT

A consistent “2D/1D” neutron transport method is derived from the 3D Boltzmann transport equation, to calculate fuel-pin-resolved neutron fluxes for realistic full-core Pressurized Water Reactor (PWR) problems. The 2D/1D method employs the Method of Characteristics to discretize the radial variables and a lower order transport solution to discretize the axial variable. This paper describes the theory of the 2D/1D method and its implementation in the MPACT code, which has become the whole-core deterministic neutron transport solver for the Consortium for Advanced Simulations of Light Water Reactors (CASL) core simulator VERA-CS. Several applications have been performed on both leadership-class and industry-class computing clusters. Results are presented for whole-core solutions of the Watts Bar Nuclear Power Station Unit 1 and compared to both continuous-energy Monte Carlo results and plant data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

SU-D-304-07: Application of Proton Boron Fusion Reaction to Radiation Therapy

Purpose: we present the introduction of a therapy method using the proton boron fusion reaction. The purpose of this study is to verify the theoretical validity of proton boron fusion therapy using Monte Carlo simulations. Methods: After boron is accumulated in the tumor region, the emitted from outside the body proton can react with the boron in the tumor region. An increase of the proton’s maximum dose level is caused by the boron and only the tumor cell is damaged more critically. In addition, a prompt gamma ray is emitted from the proton boron reaction point. Here we show that the effectiveness of the proton boron fusion therapy (PBFT) was verified using Monte Carlo simulations. Results: We found that a dramatic increase by more than half of the proton’s maximum dose level was induced by the boron in the tumor region. This increase occurred only when the proton’s maximum dose point was located within the boron uptake region (BUR). In addition, the 719 keV prompt gamma ray peak produced by the proton boron fusion reaction was positively detected. Conclusion: This therapy method features the advantages such as the application of Bragg-peak to the therapy, the accurate targeting of tumor, improved therapy effects, and the monitoring of the therapy region during treatment.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗

DOE OSTI2016

SU-F-J-123: CT-Based Determination of DIBH Variability and Its Dosimetric Impact On Post-Mastectomy Plus Regional Nodal Radiation Therapy

Purpose: Breast cancer radiotherapy delivered using voluntary deep inspiration breath-hold (DIBH) requires reproducible breath holds, particularly when matching supraclavicular fields to tangential fields. We studied the impact of variation in DIBHs on CTV and OAR dose metrics by comparing the dose distribution computed on two DIBH CT scans taken at the time of simulation. Methods: Ten patients receiving 50Gy in 25 fractions to the left chestwall and regional lymph nodes were studied. Two simulation CT scans were taken during separate DIBHs along with a free-breathing (FB) scan. The treatment was planned using one DIBH CT. The dose was recomputed on the other two scans using adaptive planning (Pinnacle 9.10) in which the scans are registered using a cross-correlation algorithm. The chestwall, lymph nodes and OARs were contoured on the scans following the RTOG consensus guidelines. The overall translational and rotational variation between the DIBH scans was used to estimate positional variation between breath-holds. Dose metrics between plans were compared using paired t-tests (p < 0.05) and means and standard deviations were reported. Results: The registration parameters were sub-millimeter and sub-degree. Although DIBH significantly reduced mean heart dose by 2.4Gy compared to FB (p < 0.01), no significant changes in dose were observed for targets or OARs between the two DIBH scans. Nodal coverage as assessed by V90% was 90%±8% and 89%±8% for supraclavicular and 99%±2% and 97%±22% for IM nodes. Though a significant decrease (10.5%±12.4%) in lung volume in the second DIBH CT was observed, the lung V20Gy was unchanged (14±2% and 14±3%) between the two DIBH scans. Conclusion: While the lung volume often varied between DIBHs, the CTV and OAR dose metrics were largely unchanged. This indicates that manual DIBH has the potential to provide consistent dose delivery to the chestwall and regional nodes targets when using matched fields.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗

DOE OSTI2014

Spin dynamics simulation of electron spin relaxation in Ni{sup 2+}(aq)

The ability to quantitatively predict and analyze the rate of electron spin relaxation of open-shell systems is important for electron paramagnetic resonance and paramagnetic nuclear magnetic resonance spectroscopies. We present a combined molecular dynamics (MD), quantum chemistry (QC), and spin dynamics simulation method for calculating such spin relaxation rates. The method is based on the sampling of a MD trajectory by QC calculations, to produce instantaneous parameters of the spin Hamiltonian used, in turn, to numerically solve the Liouville-von Neumann equation for the time evolution of the spin density matrix. We demonstrate the approach by simulating the relaxation of electron spin in an aqueous solution of Ni{sup 2+} ion. The spin-lattice (T{sub 1}) and spin-spin (T{sub 2}) relaxation rates are extracted directly from the simulations of the time dependence of the longitudinal and transverse magnetization, respectively. Good agreement with the available, indirectly obtained experimental data is obtained by our method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

SU-G-201-13: Investigation of Dose Variation Induced by HDR Ir-192 Source Global Shift Within the Varian Ring Applicator Using Monte Carlo Methods

Purpose: The Varian tandem and ring applicators are used to deliver HDR Ir-192 brachytherapy for cervical cancer. The source path within the ring is hard to predict due to the larger interior ring lumen. Some studies showed the source could be several millimeters different from planned positions, while other studies demonstrated minimal dosimetric impact. A global shift can be applied to limit the effect of positioning offsets. The purpose of this study was to assess the necessities of implementing a global source shift using Monte Carlo (MC) simulations. Methods: The MCNP5 radiation transport code was used for all MC simulations. To accommodate TG-186 guidelines and eliminate inter-source attenuation, a BrachyVision plan with 10 dwell positions (0.5cm step sizes) was simulated as the summation of 10 individual sources with equal dwell times for simplification. To simplify the study, the tandem was also excluded from the MC model. Global shifts of ±0.1, ±0.3, ±0.5 cm were then simulated as distal and proximal from the reference positions. Dose was scored in water for all MC simulations and was normalized to 100% at the normalization point 0.5 cm from the cap in the ring plane. For dose comparison, Point A was 2 cm caudal from the buildup cap and 2 cm lateral on either side of the ring axis. With seventy simulations, 108 photon histories gave a statistical uncertainties (k=1) <2% for (0.1 cm)3 voxels. Results: Compared to no global shift, average Point A doses were 0.0%, 0.4%, and 2.2% higher for distal global shifts, and 0.4%, 2.8%, and 5.1% higher for proximal global shifts, respectively. The MC Point A doses differed by < 1% when compared to BrachyVision. Conclusion: Dose variations were not substantial for ±0.3 cm global shifts, which is common in clinical practice.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Spatiotemporal pattern detection, generation, and computation with circuits

Abstract Implementations of neurons, delays, and synapse circuits are presented with simulations. These neural elements are used to create two small spiking neural networks, the Rate-Window and Order-Biased clusters, which are capable of detecting simple two-spike spatiotemporal patterns. A simple pattern detecting network (SPDN) is created by combining the Rate-Window and Order-Biased clusters, where clusters are small spiking neural networks, and its simple pattern detection ability is demonstrated in simulation. The SPDN is used to implement a complex pattern detecting network (CPDN) and its complex pattern detection ability is demonstrated in simulation. Methods for generating arbitrary spatiotemporal patterns are presented. The CPDN and spatiotemporal pattern generation methods are then used to implement a novel spatiotemporal computing paradigm based on detecting and responding to spatiotemporal symbols. A simulation of a spatiotemporal half adder is presented to demonstrate the computing paradigm.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2017

High-Performance First-Principles Molecular Dynamics for Predictive Theory and Modeling

This project focused on developing high-performance software tools for First-Principles Molecular Dynamics (FPMD) simulations, and applying them in investigations of materials relevant to energy conversion processes. FPMD is an atomistic simulation method that combines a quantum-mechanical description of electronic structure with the statistical description provided by molecular dynamics (MD) simulations. This reliance on fundamental principles allows FPMD simulations to provide a consistent description of structural, dynamical and electronic properties of a material. This is particularly useful in systems for which reliable empirical models are lacking. FPMD simulations are increasingly used as a predictive tool for applications such as batteries, solar energy conversion, light-emitting devices, electro-chemical energy conversion devices and other materials. During the course of the project, several new features were developed and added to the open-source Qbox FPMD code. The code was further optimized for scalable operation of large-scale, Leadership-Class DOE computers. When combined with Many-Body Perturbation Theory (MBPT) calculations, this infrastructure was used to investigate structural and electronic properties of liquid water, ice, aqueous solutions, nanoparticles and solid-liquid interfaces. Computing both ionic trajectories and electronic structure in a consistent manner enabled the simulation of several spectroscopic properties, such as Raman spectra, infrared spectra, and sum-frequency generation spectra. The accuracy of the approximations used allowed for direct comparisons of results with experimental data such as optical spectra, X-ray and neutron diffraction spectra. The software infrastructure developed in this project, as applied to various investigations of solids, liquids and interfaces, demonstrates that FPMD simulations can provide a detailed, atomic-scale picture of structural, vibrational and electronic properties of complex systems relevant to energy conversion devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗