Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulator”

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 433 records · Page 24

Comparison between simulated and observed LHC beam backgrounds in the ATLAS experiment at E beam =4 TeV

Results of dedicated Monte Carlo simulations of beam-induced background (BIB) in the ATLAS experiment at the Large Hadron Collider (LHC) are presented and compared with data recorded in 2012. During normal physics operation this background arises mainly from scattering of the 4 TeV protons on residual gas in the beam pipe. Methods of reconstructing the BIB signals in the ATLAS detector, developed and implemented in the simulation chain based on the \textsc{Fluka} Monte Carlo simulation package, are described. The interaction rates are determined from the residual gas pressure distribution in the LHC ring in order to set an absolute scale on the predicted rates of BIB so that they can be compared quantitatively with data. Through these comparisons the origins of the BIB leading to different observables in the ATLAS detectors are analysed. The level of agreement between simulation results and BIB measurements by ATLAS in 2012 demonstrates that a good understanding of the origin of BIB has been reached.

47 OTHER INSTRUMENTATION↗

NetSim Steer: A Runtime Steering Framework for Network Simulators

This paper presents NetSim-Steer, a runtime steering framework for network simulators. With this framework, users can specify execution constraints for the network protocols and network models. In addition to this, users can implement steering rules to be executed when a constraint is violated. These rules allows users to alter the parameters of the protocols/models during the simulation so that they stay within the boundaries of the constraints.

Network simulation↗

Advanced simulation capability for environmental management, integrated tool-sets and simulator to enhance public communication - 15186

The U.S. Department of Energy (US DOE) Office of Environmental Management (EM), Office of Soil and Groundwater Remediation, is supporting multi-laboratory development of the Advanced Simulation Capability for Environmental Management (ASCEM). ASCEM is an open source and modular computing framework that incorporates new advances and tools for predicting contaminant fate and transport in natural and engineered systems. ASCEM facilitates transparent, efficient and integrated approaches to modeling and site characterization, and provides robust and standardized assessments of performance and risk for EM cleanup and closure activities. Site subsurface simulation modeling and performance assessments are complex and involve sensitivity evaluations and uncertainty quantification to better understand the risks at a site. The ASCEM tool-set was designed in close collaboration with the sites and other potential end-users (both the public and regulators) to allow better communication and understanding of these complex modeling efforts. ASCEM's integrated tool-sets provide a user interface and tools for managing work-flows, including conceptual model development, management of data and metadata for model input, sensitivity analysis, model calibration and uncertainty analysis, model execution on diverse computational platforms, and processing of model output, including visualization. The hierarchical and modular design of ASCEM HPC simulator Amanzi easily lends itself to translating this conceptual model of a site into a numerical model that can easily be spatially and temporally discretized to run on the HPC simulator. Together this integrated work-flow allows for a seamless and integrated analysis tracking all provenance information and modeling uncertainties. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Algorithm for Building a Spectrum for NREL's One-Sun Multi-Source Simulator: Preprint

Historically, the tools used at NREL to compensate for the difference between a reference spectrum and a simulator spectrum have been well-matched reference cells and the application of a calculated spectral mismatch correction factor, M. This paper describes the algorithm for adjusting the spectrum of a 9-channel fiber-optic-based solar simulator with a uniform beam size of 9 cm square at 1-sun.The combination of this algorithm and the One-Sun Multi-Source Simulator (OSMSS) hardware reduces NREL's current vs. voltage measurement time for a typical three-junction device from man-days to man-minutes. These time savings may be significantly greater for devices with more junctions.

14 SOLAR ENERGY↗

A Hierarchical Evaluation of Regional Climate Simulations

Global climate models (GCMs) are the primary tools for predicting the evolution of the climate system. Through decades of development, GCMs have demonstrated useful skill in simulating climate at continental to global scales. However, large uncertainties remain in projecting climate change at regional scales, which limit our ability to inform decisions on climate change adaptation and mitigation. To bridge this gap, different modeling approaches including nested regional climate models (RCMs), global stretch-grid models, and global high-resolution atmospheric models have been used to provide regional climate simulations (Leung et al. 2003). In previous efforts to evaluate these approaches, isolating their relative merits was not possible because factors such as dynamical frameworks, physics parameterizations, and model resolutions were not systematically constrained. With advances in high performance computing, it is now feasible to run coupled atmosphere-ocean GCMs at horizontal resolution comparable to what RCMs use today. Global models with local refinement using unstructured grids have become available for modeling regional climate (e.g., Rauscher et al. 2012; Ringler et al. 2013). While they offer opportunities to improve climate simulations, significant efforts are needed to test their veracity for regional-scale climate simulations.

climate↗

Efficient Simulation of Welding Distortion in Large Structures and Its Reduction by Jig Constraints

Two large construction machinery structures were welded separately with and without jig constraints. The welding distortion of the entire structure was measured by a 3D coordinate measuring system and simulated by elastic FEM using the inherent deformation method. To obtain an accurate inherent deformation, a thermal elastoplastic FE analysis of simple one-side fillet joints with and without jig constraints was performed. Efficient simulation of welding distortion in large structures was accomplished by applying inherent deformation in a localized region, and the effect of jig constraint on the reduction of welding distortion was clarified. The computation of inherent deformation, the weld zone definition and the conversion of inherent deformation into inherent strain were automated. Measured and computed welding distortions in large structures had a good correspondence with respect to both tendency and magnitude. Further investigation of jig configuration was performed for enhanced reduction of welding distortion. Alternative controlling techniques for common welded structures were also addressed.

36 MATERIALS SCIENCE↗

Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed

Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.

15 GEOTHERMAL ENERGY↗

Testing thermal gradient driving force for grain boundary migration using molecular dynamics simulations

Strong thermal gradients in low-thermal-conductivity ceramics may drive extended defects, such as grain boundaries and voids, to migrate in preferential directions. In this work, molecular dynamics simulations are conducted to study thermal gradient driven grain boundary migration and to verify a previously proposed thermal gradient driving force equation, using uranium dioxide as a model system. It is found that a thermal gradient drives grain boundaries to migrate up the gradient and the migration velocity increases under a constant gradient owing to the increase in mobility with temperature. Different grain boundaries migrate at very different rates due to their different intrinsic mobilities. The extracted mobilities from the thermal gradient driven simulations are compared with those calculated from two other well-established methods and good agreement between the three different methods is found, demonstrating that the theoretical equation of the thermal gradient driving force is valid, although a correction of one input parameter should be made. The discrepancy in the grain boundary mobilities between modeling and experiments is also discussed.

36 MATERIALS SCIENCE↗

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗

Multilevel Monte Carlo simulation of Coulomb collisions

We present a new, for plasma physics, highly efficient multilevel Monte Carlo numerical method for simulating Coulomb collisions. The method separates and optimally minimizes the finite-timestep and finite-sampling errors inherent in the Langevin representation of the Landau–Fokker–Planck equation. It does so by combining multiple solutions to the underlying equations with varying numbers of timesteps. For a desired level of accuracy ε, the computational cost of the method is O(ε{sup −2}) or O(ε{sup −2}(lnε){sup 2}), depending on the underlying discretization, Milstein or Euler–Maruyama respectively. This is to be contrasted with a cost of O(ε{sup −3}) for direct simulation Monte Carlo or binary collision methods. We successfully demonstrate the method with a classic beam diffusion test case in 2D, making use of the Lévy area approximation for the correlated Milstein cross terms, and generating a computational saving of a factor of 100 for ε=10{sup −5}. We discuss the importance of the method for problems in which collisions constitute the computational rate limiting step, and its limitations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Noiseless Vlasov–Poisson simulations with linearly transformed particles

We introduce a deterministic discrete-particle simulation approach, the Linearly-Transformed Particle-In-Cell (LTPIC) method, that employs linear deformations of the particles to reduce the noise traditionally associated with particle schemes. Formally, transforming the particles is justified by local first order expansions of the characteristic flow in phase space. In practice the method amounts of using deformation matrices within the particle shape functions; these matrices are updated via local evaluations of the forward numerical flow. Because it is necessary to periodically remap the particles on a regular grid to avoid excessively deforming their shapes, the method can be seen as a development of Denavit's Forward Semi-Lagrangian (FSL) scheme (Denavit, 1972 [8]). However, it has recently been established (Campos Pinto, 2012 [20]) that the underlying Linearly-Transformed Particle scheme converges for abstract transport problems, with no need to remap the particles; deforming the particles can thus be seen as a way to significantly lower the remapping frequency needed in the FSL schemes, and hence the associated numerical diffusion. To couple the method with electrostatic field solvers, two specific charge deposition schemes are examined, and their performance compared with that of the standard deposition method. Finally, numerical 1d1v simulations involving benchmark test cases and halo formation in an initially mismatched thermal sheet beam demonstrate some advantages of our LTPIC scheme over the classical PIC and FSL methods. Benchmarked test cases also indicate that, for numerical choices involving similar computational effort, the LTPIC method is capable of accuracy comparable to or exceeding that of state-of-the-art, high-resolution Vlasov schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational capacity in hydrodynamic real-time hybrid simulation applied to simulate the dynamic response of floating offshore wind turbines

Real-time hybrid simulation (RTHS) mitigates similitude distortions in model-scale tests of floating offshore wind turbines (FOWTs) by coupling physical experiments with numerical models in real time. The coupling requires faster-than-real-time numerical computations to satisfy temporal similitude with the physical experiment, presenting a bottleneck for using more complex numerical models in RTHS. This paper presents a hydrodynamic-RTHS (hydro-RTHS) framework for FOWTs that simulates the hydrodynamics physically and the aerodynamics numerically with sensor feedback from the physical testing. The framework adapts the three-loop hardware architecture to leverage greater computational resources and mitigate strict temporal requirements, enabling more computationally demanding numerical analyses in hydro-RTHS. The three-loop hardware architecture integrates multiple machines, each dedicated to either numerical analysis or RTHS controls, with a rate-transition algorithm to synchronize the tasks executed across the different machine processors. Virtual and physical tests verified and validated the hydro-RTHS framework, respectively. The ”virtual” tests, which approximates the physical domain numerically, verified the RTHS framework with respect to a numerical full-scale complete FOWT model simulated in the open-source software, OpenFAST. The virtual tests were able to maintain comparable control signals while enabling greater computational resources for the numerical calculations. Real-world physical tests demonstrated that the hydro-RTHS framework computes aerodynamic forces similar to the complete OpenFAST model, validating the hydro-RTHS framework using the three-loop hardware architecture. Findings show that the hydro-RTHS framework with the three-loop hardware architecture is computationally efficient, with reserve capacity to simulate more complex problems due to the customized software, hardware, and rate-transition algorithm.

17 WIND ENERGY↗

Direct simulation of electron transfer using ring polymer molecular dynamics: Comparison with semiclassical instanton theory and exact quantum methods

The use of ring polymer molecular dynamics (RPMD) for the direct simulation of electron transfer (ET) reaction dynamics is analyzed in the context of Marcus theory, semiclassical instanton theory, and exact quantum dynamics approaches. For both fully atomistic and system-bath representations of condensed-phase ET, we demonstrate that RPMD accurately predicts both ET reaction rates and mechanisms throughout the normal and activationless regimes of the thermodynamic driving force. Analysis of the ensemble of reactive RPMD trajectories reveals the solvent reorganization mechanism for ET that is anticipated in the Marcus rate theory, and the accuracy of the RPMD rate calculation is understood in terms of its exact description of statistical fluctuations and its formal connection to semiclassical instanton theory for deep-tunneling processes. In the inverted regime of the thermodynamic driving force, neither RPMD nor a related formulation of semiclassical instanton theory capture the characteristic turnover in the reaction rate; comparison with exact quantum dynamics simulations reveals that these methods provide inadequate quantization of the real-time electronic-state dynamics in the inverted regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct simulation of proton-coupled electron transfer across multiple regimes

The coupled transfer of electrons and protons is a central feature of biological and molecular catalysis, yet fundamental aspects of these reactions remain poorly understood. In this study, we extend the ring polymer molecular dynamics (RPMD) method to enable direct simulation of proton-coupled electron transfer (PCET) reactions across a wide range of physically relevant regimes. In a system-bath model for symmetric, co-linear PCET in the condensed phase, RPMD trajectories reveal distinct kinetic pathways associated with sequential and concerted PCET reaction mechanisms, and it is demonstrated that concerted PCET proceeds by a solvent-gating mechanism in which the reorganization energy is mitigated by charge cancellation among the transferring particles. We further employ RPMD to study the kinetics and mechanistic features of concerted PCET reactions across multiple coupling regimes, including the fully non-adiabatic (both electronically and vibrationally non-adiabatic), partially adiabatic (electronically adiabatic, but vibrationally non-adiabatic), and fully adiabatic (both electronically and vibrationally adiabatic) limits. Comparison of RPMD with the results of PCET rate theories demonstrates the applicability of the direct simulation method over a broad range of conditions; it is particularly notable that RPMD accurately predicts the crossover in the thermal reaction rates between different coupling regimes while avoiding a priori assumptions about the PCET reaction mechanism. Finally, by utilizing the connections between RPMD rate theory and semiclassical instanton theory, we show that analysis of ring-polymer configurations in the RPMD transition path ensemble enables the a posteriori determination of the coupling regime for the PCET reaction. This analysis reveals an intriguing and distinct 'transient-proton-bridge' mechanism for concerted PCET that emerges in the transition between the proton-mediated electron superexchange mechanism for fully non-adiabatic PCET and the hydrogen atom transfer mechanism for partially adiabatic PCET. Taken together, these results provide a unifying picture of the mechanisms and physical driving forces that govern PCET across a wide range of physical regimes, and they raise the possibility for PCET mechanisms that have not been previously reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Verification of particle simulation of radio frequency waves in fusion plasmas

Radio frequency (RF) waves can offer heating, current and flow drive, as well as instability control for steady state operations of fusion experiments. A particle simulation model has been designed in this work to provide a first-principles tool for studying the RF nonlinear interactions with plasmas. In this model, ions are considered as fully kinetic particles using the Vlasov equation and electrons are treated as guiding centers using the drift kinetic equation. This model has been implemented in a global gyrokinetic toroidal code using real electron-to-ion mass ratio. To confirm the model, linear simulations of ion plasma oscillation, ion Bernstein wave, and lower hybrid wave are carried out in cylindrical geometry and found to agree well with analytic predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomistic simulation study of short pulse laser interactions with a metal target under conditions of spatial confinement by a transparent overlayer

The distinct characteristics of short pulse laser interactions with a metal target under conditions of spatial confinement by a solid transparent overlayer are investigated in a series of atomistic simulations. The simulations are performed with a computational model combining classical molecular dynamics (MD) technique with a continuum description of the laser excitation, electron-phonon equilibration, and electronic heat transfer based on two-temperature model (TTM). Two methods for incorporation of the description of a transparent overlayer into the TTM-MD model are designed and parameterized for Ag-silica system. The material response to the laser energy deposition is studied for a range of laser fluences that, in the absence of the transparent overlayer, covers the regimes of melting and resolidification, photomechanical spallation, and phase explosion of the overheated surface region. In contrast to the irradiation in vacuum, the spatial confinement by the overlayer facilitates generation of sustained high-temperature and high-pressure conditions near the metal-overlayer interface, suppresses the generation of unloading tensile wave, decreases the maximum depth of melting, and prevents the spallation and explosive disintegration of the surface region of the metal target. At high laser fluences, when the laser excitation brings the surface region of the metal target to supercritical conditions, the confinement prevents the expansion and phase decomposition characteristic for the vacuum conditions leading to a gradual cooling of the hot compressed supercritical fluid down to the liquid phase and eventual solidification. The target modification in this case is limited to the generation of crystal defects and the detachment of the metal target from the overlayer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A simulation assessment of the thermodynamics of dense ion-dipole mixtures with polarization

Molecular dynamics (MD) simulations are employed to ascertain the relative importance of various electrostatic interaction contributions, including induction interactions, to the thermodynamics of dense, hot ion-dipole mixtures. In the absence of polarization, we find that an MD-constrained free energy term accounting for the ion-dipole interactions, combined with well tested ionic and dipolar contributions, yields a simple, fairly accurate free energy form that may be a better option for describing the thermodynamics of such mixtures than the mean spherical approximation (MSA). Polarization contributions induced by the presence of permanent dipoles and ions are found to be additive to a good approximation, simplifying the thermodynamic modeling. We suggest simple free energy corrections that account for these two effects, based in part on standard perturbative treatments and partly on comparisons with MD simulation. Even though the proposed approximations likely need further study, they provide a first quantitative assessment of polarization contributions at high densities and temperatures and may serve as a guide for future modeling efforts.

74 ATOMIC AND MOLECULAR PHYSICS↗