Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Immersive Simulations and Engineering Environment (iSEE): Increasing the Realism of Simulations

The Immersive Simulations and Engineering Environment (iSEE) lab at Kennedy Space Center has been operating with its current team since 2017. The lab provides National Aeronautics and Space Administration (NASA) with simulations that show the stress and strain that a strenuous task can cause on a body. Previously the lab was using software that was efficient in ergonomic analysis of a human body, but limited in creating a realistic simulation. To produce a more realistic simulation, the iSEE invested in new software that could potentially create an environment that could show the viewer a simulation that was closer to reality.

Hodge, Brianna↗
NASA NTRSDate not supplied

Comparison of Physical Workload Across Eva-Simulation Analog Environments: Hybrid Space Suit Simulator and Pressurized Suit Testing

NASA conducts research, testing, and training across a variety of analog environments to support characterization of human performance during Extravehicular Activities (EVAs). Time utilizing pressurized suits in an offloaded environment is both limited and expensive, making it challenging to carry out extensive research with these suits. The Human Physiology, Performance, Protection, and Operations (H-3PO) Laboratory at NASA Johnson Space Center designed a Hybrid Space Suit Simulator (HS3) as a low-cost, workload approximator and easy access research tool to provide relevant physical and cognitive workloads during simulated EVAs. A pilot study was conducted in a 1g analog environment where six healthy subjects (3 male, 3 female) underwent simulated 5-hour EVAs in the HS3. This study used a COSMED K5 portable metabolic analyzer and a Polar H10 heart rate monitor to evaluate physical workload during the EVAs. For direct comparison, EVA tasks and timelines were modeled after a similar study conducted in pressurized suits (Mark-III spacesuit, n=3 male; small xPGS spacesuit, n=3 female) at NASA’s Active Response Gravity Offload System (ARGOS) offloaded to Lunar gravity (1/6 g). The 1g HS3 data demonstrated increased metabolic rate when compared to pressurized, Lunar-offloaded suited simulated EVAs for certain tasks including the kneeling scoop (HS3: 1194±199 BTU/hr, ARGOS: 860±170 BTU/hr, p < 0.05) and 20 lb object relocation (HS3: 1651±136 BTU/hr, ARGOS: 1048±252.8 BTU/hr, p < 0.05), but it did not demonstrate any significant differences in heart rate (p > 0.05). Other simulated EVA tasks, such as a 500 m traverse (HS3: 1431±177 BTU/hr, ARGOS: 1212±294 BTU/hr, p > 0.05), showed similar workload profiles with no significant workload differences between HS3 and ARGOS. The HS3 is a useful research tool for increasing workloads in EVA research without pressurized suits, though it may overestimate workload during certain EVA tasks when used in 1g environments.

Zachary Wusk↗
DOE OSTI2018

Towards a theory of economic value for modeling and simulation: incremental cost of parallel simulation (wip)

A useful perspective on the economic value of simulation in engineering comes from considering the role simulation plays in constructing a system. We propose a conceptual model that combines an embryonic concept of the value of an information bit with cost estimation tools used in software engineering. The model provides an economic justification for the common conception that models should be built to address specific questions rather than to comprehensively represent a referent system. Extending these ideas, we offer a model for the incremental value of parallel execution to explain why simulation applications that take advantage of sophisticated parallel discrete event simulation tools form a relatively narrow niche despite the ubiquity of parallel computing hardware. The proposed model enables us to bring value of information criteria into model development and to make predictions about the future of parallel simulation applications.

Nutaro, James J↗
DOE OSTI2021

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗
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

Numerical simulations of the two-dimensional multimode Richtmyer-Meshkov instability

We show the two-dimensional Richtmyer-Meshkov instability occurs as shock waves pass through a perturbed material interface, triggering transition to an inhomogeneous turbulence variable density flow. This paper presents a series of large-eddy-simulations of the two dimensional turbulent RM instability and compares the results to the fully three dimensional simulations. There are two aims for this paper, the first is to explore what numerical resolution is required for a statistically converged solution for a two dimensional inhomogeneous flow field. The second aim is to elucidate the key differences in flow physics between the two dimensional and three dimensional Richtmyer-Meshkov instabilities, particularly their asymptotic self-similar regime. Convergence is achieved using 64 independent realisations and grid resolutions up to 4096 2 in the plane. It is shown that for narrowband cases the growth rate θ = 0.48 which is substantially higher than the three-dimensional equivalent. Mix measures are consistently lower compared to three-dimensional, and the kinetic energy distribution is homogeneous at late time. The broadband case has a similar initial growth rate as the three-dimensional case, with a marginally lower θ = 0.63. Mix is similar in magnitude, but is reducing at late time. The spectra in both cases exhibit the dual-cascade expected from two-dimensional turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Combining experiment and optical simulation in coherent X-ray nanobeam characterization of Si/SiGe semiconductor heterostructures

The highly coherent and tightly focused x-ray beams produced by hard x-ray light sources enable the nanoscale characterization of the structure of electronic materials but are accompanied by significant challenges in the interpretation of diffraction and scattering patterns. X-ray nanobeams exhibit optical coherence combined with a large angular divergence introduced by the x-ray focusing optics. The scattering of nanofocused x-ray beams from intricate semiconductor heterostructures produces a complex distribution of scattered intensity. We report here an extension of coherent x-ray optical simulations of convergent x-ray beam diffraction patterns to arbitrary x-ray incident angles to allow the nanobeam diffraction patterns of complex heterostructures to be simulated faithfully. These methods are used to extract the misorientation of lattice planes and the strain of individual layers from synchrotron x-ray nanobeam diffraction patterns of Si/SiGe heterostructures relevant to applications in quantum electronic devices. The systematic interpretation of nanobeam diffraction patterns from semiconductor heterostructures presents a new opportunity in characterizing and ultimately designing electronic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

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↗
DOE OSTI2015

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↗
DOE OSTI2017

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↗
DOE OSTI2025

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↗
DOE OSTI2015

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↗
DOE OSTI2026

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↗
DOE OSTI2014

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↗
DOE OSTI2014

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↗
DOE OSTI2026

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↗
DOE OSTI2013

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↗