Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulations”

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 1,567 records · Page 87

Evaluation of UTLS Carbon Monoxide Simulations in GMI and GEOS-Chem Chemical Transport Models using Aura MLS Observations

This study evaluates the distribution and variation of carbon monoxide (CO) in the upper troposphere and lower stratosphere (UTLS) during 2004-2012 as simulated by two chemical transport models, using the latest version of Aura Microwave Limb Sounder (MLS) observations. The simulated spatial distributions, temporal variations and vertical transport of CO in the UTLS region are compared with those observed by MLS. We also investigate the impact of surface emissions and deep convection on CO concentrations in the UTLS over different regions, using both model simulations and MLS observations. Global Modeling Initiative (GMI) and GEOS-Chem simulations of UTLS CO both show similar spatial distributions to observations. The global mean CO values simulated by both models agree with MLS observations at 215 and 147 hPa, but are significantly underestimated by more than 40% at 100 hPa. In addition, the models underestimate the peak CO values by up to 70% at 100 hPa, 60% at 147 hPa and 40% at 215 hPa, with GEOS-Chem generally simulating more CO at 100 hPa and less CO at 215 hPa than GMI. The seasonal distributions of CO simulated by both models are in better agreement with MLS in the Southern Hemisphere (SH) than in the Northern Hemisphere (NH), with disagreements between model and observations over enhanced CO regions such as southern Africa. The simulated vertical transport of CO shows better agreement with MLS in the tropics and the SH subtropics than the NH subtropics. We also examine regional variations in the relationships among surface CO emission, convection and UTLS CO concentrations. The two models exhibit emission-convection- CO relationships similar to those observed by MLS over the tropics and some regions with enhanced UTLS CO.

spatial distribution↗

Uncertainty of Wheat Water Use: Simulated Patterns and Sensitivity to Temperature and CO2

Projected global warming and population growth will reduce future water availability for agriculture. Thus, it is essential to increase the efficiency in using water to ensure crop productivity. Quantifying crop water use (WU; i.e. actual evapotranspiration) is a critical step towards this goal. Here, sixteen wheat simulation models were used to quantify sources of model uncertainty and to estimate the relative changes and variability between models for simulated WU, water use efficiency (WUE, WU per unit of grain dry mass produced), transpiration efficiency (Teff, transpiration per kg of unit of grain yield dry mass produced), grain yield, crop transpiration and soil evaporation at increased temperatures and elevated atmospheric carbon dioxide concentrations ([CO2]). The greatest uncertainty in simulating water use, potential evapotranspiration, crop transpiration and soil evaporation was due to differences in how crop transpiration was modelled and accounted for 50 of the total variability among models. The simulation results for the sensitivity to temperature indicated that crop WU will decline with increasing temperature due to reduced growing seasons. The uncertainties in simulated crop WU, and in particularly due to uncertainties in simulating crop transpiration, were greater under conditions of increased temperatures and with high temperatures in combination with elevated atmospheric [CO2] concentrations. Hence the simulation of crop WU, and in particularly crop transpiration under higher temperature, needs to be improved and evaluated with field measurements before models can be used to simulate climate change impacts on future crop water demand.

Sensitivity↗

The Use of Microgravity Simulators for Space Research

The spaceflight environment is known to influence biological processes ranging from stimulation of cellular metabolism to possible impacts on cellular damage repair, suppression of immune functions, and bone loss in astronauts. Microgravity is one of the most significant stress factors experienced by living organisms during spaceflight, and therefore, understanding cellular responses to altered gravity at the physiological and molecular level is critical for expanding our knowledge of life in space. Since opportunities to conduct experiments in space are scarce, various microgravity simulators and analogues have been widely used in space biology ground studies. Even though simulated microgravity conditions have produced some, but not all of the biological effects observed in the true microgravity environment, they provide test beds that are effective, affordable, and readily available to facilitate microgravity research. Kennedy Space Center (KSC) provides ground microgravity simulator support to offer a variety of microgravity simulators and platforms for Space Biology investigators. Assistance will be provided by both KSC and external experts in molecular biology, microgravity simulation, and engineering. Comparisons between the physical differences in microgravity simulators, examples of experiments using the simulators, and scientific questions regarding the use of microgravity simulators will be discussed.

Simulation↗

Credible practice of modeling and simulation in healthcare: ten rules from a multidisciplinary perspective

The complexities of modern biomedicine are rapidly increasing. Thus, modeling and simulation have become increasingly important as a strategy to understand and predict the trajectory of pathophysiology, disease genesis, and disease spread in support of clinical and policy decisions. In such cases, inappropriate or ill-placed trust in the model and simulation outcomes may result in negative outcomes, and hence illustrate the need to formalize the execution and communication of modeling and simulation practices. Although verification and validation have been generally accepted as significant components of a model’s credibility, they cannot be assumed to equate to a holistic credible practice, which includes activities that can impact comprehension and in-depth examination inherent in the devel-opment and reuse of the models. For the past several years, the Committee on Credible Practice of Modeling and Simulation in Healthcare, an interdisciplinary group seeded from a U.S. interagency initiative, has worked to codify best practices. Here, we provide Ten Rules for credible practice of modeling and simulation in healthcare developed from a comparative analysis by the Committee’s multidisciplinary membership, followed by a large stakeholder com-munity survey. These rules establish a unified conceptual framework for modeling and simulation design, implementation, evaluation, dissemination and usage across the modeling and simulation life-cycle. While biomedical science and clinical care domains have somewhat different requirements and expectations for credible practice, our study converged on rules that would be useful across a broad swath of model types. In brief, the rules are: (1) Define context clearly. (2) Use contextually appropriate data. (3) Evaluate within context. (4) List limitations explicitly. (5) Use version control. (6) Document appropriately. (7) Disseminate broadly. (8) Get independent reviews. (9) Test competing imple-mentations. (10) Conform to standards. Although some of these are common sense guidelines, we have found that many are often missed or misconstrued, even by seasoned practitioners. Computational models are already widely used in basic science to generate new biomedical knowledge. As they penetrate clinical care and healthcare policy, contributing to personalized and precision medicine, clinical safety will require established guidelines for the credible practice of modeling and simulation in healthcare.

Credibility↗

Numerical Simulation of the High-Speed Leg of the National Transonic Facility

Numerical simulations for the flow inside the high-speed leg of the National Transonic Facility was conducted. The NASA Tetrahedral Unstructured Software System (TetrUSS) with its USM3D_ME solver was used to perform the numerical simulations. USM3D_ME is developed and maintained by NASA Langley Research Center. Simulations were conducted for three configurations: empty tunnel, body of revolution installed, and NASA Common Research Model installed in the test section. Simulations were performed for a test section Mach number of 0.7 and 0.85 and a corresponding Reynolds number of 8 million per ft. A controller was developed that automated dynamic outflow boundary and streamlined the process of running multiple simulations. The use of a dynamic outflow boundary was the key parameter to drive simulation to the desired tunnel conditions. The numerical simulations captured expected flow features in NTF test section and in the tunnel plenum. The simulations depicted that separation of the flow inside the diffuser appears to be asymmetric and more extensive toward the top and bottom walls. Analysis of the flow field were conducted, and the computed drag coefficient compared to wind tunnel data.

NTF↗

Verification and Validation of Elastodynamic Simulation Software for Aerospace Research

Physics-based simulation of nondestructive evaluation (NDE) inspection can help to advance the inspectability and reliability of mechanical systems. However, NDE simulations applicable to non-idealized mechanical components often require large compute domains and long run times. This has prompted development of custom NDE simulation software tailored to high performance computing (HPC) hardware. Verification and validation (V&V) is an integral part of developing this software to ensure implementations are robust and applicable to inspection problems, producing tools and simulations suitable for computational NDE research. This presentation addresses factors common to V&V of several elastodynamic simulation codes applicable to ultrasonic NDE. Examples are drawn from in-house simulation software at NASA Langley Research Center, ranging from ensuring reliability in a 1D heterogeneous media wave equation solver to the V&V needs of 3D cluster-parallel elastodynamic software. Factors specific to a research environment are addressed, where individual simulation results can be as relevant as the software product itself. Distinct facets of V&V are discussed including testing to establish software reliability, employing systematic approaches for consistency with fundamental conservation laws, establishing the numerical stability of algorithms, and demonstrating concurrence with empirical data. This talk also addresses V&V practices for small groups of researchers. This includes establishing resources (e.g. time and personnel) for V&V during project planning to mitigate and control the risk of setbacks. Similarly, we identify ways for individual researchers to use V&V during simulation software development itself to both speed up the development process and reduce incurred technical debt.

NDE↗

Regolith Simulant Preparation and Geotechnical Characterization for Plume Surface Interaction Testing

Descent engine plumes interact with the lunar surface and accelerate regolith particles to potentially high velocities. These ejecta create risks to surface assets that have yet to be fully assessed. To better understand these risks, plume surface interactions can be simulated on the ground by firing a test engine plume into a bin of lunar regolith simulant under vacuum conditions. The dynamics of the resultant ejecta can then be recorded. In this technical memorandum we discuss the processes used in preparing a 150 kg bin of lunar regolith simulant for plume surface interaction ground tests under vacuum conditions for the NASA STMD Plume Surface Interaction project. We present our approach to mitigating regolith simulant eruptions during pump-down, the methods used to fill and reset the regolith simulant bin for each test, and the techniques used to characterize the consistency of regolith simulant geotechnical properties before each new firing. The challenges of preparing a regolith simulant test bin below an ambient pressure of one atmosphere, particularly on the large scale, could largely be overcome with a system that could fill the test bin with simulant inside the chamber and under vacuum conditions.

lunar regolith↗

Initial Testing the Exploration Extravehicular Mobility Unit (xEMU) in Lunar Environment Simulation at the Neutral Buoyancy Lab (NBL) in 2022

Understanding how to effectively train for Extravehicular Activities (EVAs) for Artemis missions is critical. Tasked by the new Extravehicular Activity and Human Surface Mobility Program(EHP), a team of stakeholders from engineering and flight operations worked together to test the Exploration Extravehicular Mobility (xEMU) in the Neutral Buoyancy Lab(NBL) simulating the lunar environment. The NBL is a pool in which astronauts are typically trained to complete simulated International Space Station (ISS) EVA tasks in a neutral weigh-out. The overall focus for this test series was evaluating how well the NBL could be used as an effective simulated environment for Artemis lunar EVAs. New NBL support hardware was manufactured to support this test series and update the interfaces to be more representative of xEMU hardware such as a new NBL Portable Life Support System(PLSS)mock-up and umbilical, lights and cameras, and display and control unit (DCU).To simulate the 1/6th-gravity (1/6-g) environment, extra weights were added around the suit to create an accurate center of gravity representation. The Partial Gravity Weigh-out Stand (PGWS) predicted where and how much weight would be required for the 1/6th-G simulation and predicted the center of gravity placement. After the weigh-out attempt was completed, the suited participant would complete a set of directed tasks such as kneeling, object pick-up, and walking to evaluate the weigh-out. To create the lunar environment, sand, boulders, and ramps were added to the pool floor. If time allowed after weigh-out was determined acceptable, the flight operations team would direct the suited participants to complete simulated Artemis EVA tasks using flight-like tools and tool cart. The team also evaluated turning out the overhead lights in the NBL for a portion of some of the test runs and using a solar simulator to create amore realistic lunar south pole lighting environment. The test series completed one unmanned test and nine manned tests. With four of the tests, dual xEMU suited operations were completed with subjects in both the small and large xEMU sizes. Overall, the NBL was able to create an acceptable lunar landscape and ability to have two subjects work side by side in pressurized xEMUs which currently cannot be achieved in other test environments. Achieving a stable 1/6th-G lunar weigh-out and clear, consistent dual suit communications was a challenge for this series. Future work is planned to continue to improve this simulation environment.

spacesuit↗

Coding Structures for Seated Row Simulation of an Active Controlled Vibration Isolation and Stabilization System for Astronaut’s Exercise Platform

Simulation for seated aerobic row exercise was a continued task to assist NASA in analyzing a one-dimensional vibration isolation and stabilization system for astronaut’s exercise platform. Feedback delay and signal noise were added to the simulation model. Simulation runs for this study were conducted in two software simulation tools, Trick and MBDyn, software simulation environments developed at the NASA Johnson Space Center. The exciter force in the simulation was calculated from motion capture of an exerciser during a seated aerobic row exercise. The simulation runs include passive control, active control using a Proportional, Integral, Derivative (PID) controller, and active control using a Piecewise Linear Integral Derivative (PWLID) controller. Output parameters include displacements of the exercise platform, the exerciser, and the counterweight; transmitted force to the wall of spacecraft; and actuator force to the platform. The simulation results showed excellent force reduction in the active controlled system compared to the passive controlled system, which resulted in less force reduction.

Simulation↗

Initial Testing of the Exploration Extravehicular Mobility Unit (xEMU) in Lunar Environment Simulation at the Neutral Buoyancy Lab (NBL) in 2022

Understanding how to effectively train for Extravehicular Activities (EVAs) for Artemis missions is critical. Tasked by the new Extravehicular Activity and Human Surface Mobility Program(EHP), a team of stakeholders from engineering and flight operations worked together to test the Exploration Extravehicular Mobility (xEMU) in the Neutral Buoyancy Lab(NBL) simulating the lunar environment. The NBL is a pool in which astronauts are typically trained to complete simulated International Space Station (ISS) EVA tasks in a neutral weigh-out. The overall focus for this test series was evaluating how well the NBL could be used as an effective simulated environment for Artemis lunar EVAs. New NBL support hardware was manufactured to support this test series and update the interfaces to be more representative of xEMU hardware such as a new NBL Portable Life Support System(PLSS)mock-up and umbilical, lights and cameras, and display and control unit (DCU).To simulate the 1/6th-gravity (1/6-g) environment, extra weights were added around the suit to create an accurate center of gravity representation. The Partial Gravity Weigh-out Stand (PGWS) predicted where and how much weight would be required for the 1/6th-G simulation and predicted the center of gravity placement. After the weigh-out attempt was completed, the suited participant would complete a set of directed tasks such as kneeling, object pick-up, and walking to evaluate the weigh-out. To create the lunar environment, sand, boulders, and ramps were added to the pool floor. If time allowed after weigh-out was determined acceptable, the flight operations team would direct the suited participants to complete simulated Artemis EVA tasks using flight-like tools and tool cart. The team also evaluated turning out the overhead lights in the NBL for a portion of some of the test runs and using a solar simulator to create amore realistic lunar south pole lighting environment. The test series completed one unmanned test and nine manned tests. With four of the tests, dual xEMU suited operations were completed with subjects in both the small and large xEMU sizes. Overall, the NBL was able to create an acceptable lunar landscape and ability to have two subjects work side by side in pressurized xEMUs which currently cannot be achieved in other test environments. Achieving a stable 1/6th-G lunar weigh-out and clear, consistent dual suit communications was a challenge for this series. Future work is planned to continue to improve this simulation environment.

spacesuit↗

Verification of the Generalized Aerospace Simulation in Simulink (R)

NASA uses six-degrees-of-freedom (6-DOF) simulations tools to design, test, develop Guidance Navigation and Control (GN&C) software, and certify vehicle performance prior to flight. Therefore, it is critical that the 6-DOF tools used for vehicle design and certification are validated. The focus of this work is the vali-dation of the NASA Marshall Space Flight Center 6-DOF “GeneraLized Aero-space Simulation in Simulink” (GLASS) framework tool. The GLASS tool framework is currently used to support NASA GN&C insight for the Human Landing System (HLS) project, simulating vehicle dynamics during lunar descent and as-cent. The GLASS framework utilizes the off-the-shelf Mathworks ® Simscape Multibody® toolbox to model vehicle multi-body dynamics. NASA’s Engineering and Safety Center (NESC) provides a set of 6-DOF simulation verification “check cases” that are available to any user needing to verify 6-DOF tools. The check cases contain seventeen atmospheric and twenty-six orbital test scenarios are provided to validate equations of motion, environmental models (e.g., atmosphere, gravitation, and geodesy) and tool propagators. This paper compares GLASS 6-DOF simulation results against the NESC check cases’ results via simulation-to-simulation comparisons. The comparisons demonstrate that GLASS simulation results are “in family” with the outputs of the applicable NASA NESC check cases and verifies the GLASS core framework dynamics and the correct implementation of the check case scenario models.

6-Dof↗

Verification of the Generalized Aerospace Simulation in Simulink

NASA uses six-degrees-of-freedom (6-DOF) simulations tools to design, test, develop Guidance Navigation and Control (GN&C) software, and certify vehicle performance prior to flight. Therefore, it is critical that the 6-DOF tools used for vehicle design and certification are validated. The focus of this work is the validation of the NASA Marshall Space Flight Center 6-DOF “GeneraLized Aerospace Simulation in Simulink” (GLASS) framework tool. The GLASS tool framework is currently used to support NASA GN&C insight for the Human Landing System (HLS) project, simulating vehicle dynamics during lunar descent and ascent. The GLASS framework utilizes the off-the-shelf Mathworks (R) Simscape (TM) Multibody (TM) toolbox to model vehicle multi-body dynamics. NASA’s Engineering and Safety Center (NESC) provides a set of 6-DOF simulation verification “check cases” that are available to any user needing to verify 6-DOF tools. The check cases contain seventeen atmospheric and twenty-six orbital test scenarios are provided to validate equations of motion, environmental models (e.g., atmosphere, gravitation, and geodesy) and tool propagators. This paper compares GLASS 6-DOF simulation results against the NESC check cases’ results via simulation-to-simulation comparisons. The comparison demonstrates that GLASS simulation results are “in family” with the outputs of the applicable NASA NESC check-cases and verify the GLASS core framework dynamics and the implementation of the check case scenario models.

6-Dof↗

Immersive Technologies for Human-in-the-Loop Lunar Surface Simulations

NASA, the National Aeronautics and Space Administration, continually seeks innovative solutions to enhance its operations, particularly in the realms of testing, evaluation, and training for future missions. Immersive technologies, such as virtual, augmented, and mixed reality have proven to be powerful tools for immersing users in realistic, interactive, and engaging environments. This paper explores how the Simulation and Graphics Branch at NASA’s Johnson Space Center (JSC) leverages immersive technology, modern commercial rendering engines, and physics-based systems simulations to develop human-in-the-loop systems for humanity’s return to the Moon through the Artemis program. When NASA returns to the Moon, astronauts will travel to the Moon’s South Pole where lighting conditions will cause a more complex operational environment. Human-in-the-loop simulations play a crucial role in NASA’s mission planning, spacecraft and space systems development, and evaluation of operational scenarios. The development of immersive environments such as a lunar rover mockup at a video wall enables engineers and astronauts to simulate and experience mission scenarios, integrated spacecraft systems, and operational procedures in a relevant environment before deployment. By integrating realistic virtual environments, immersive technology allows for the visualization and interaction with virtual spacecraft models, mission landscapes, and complex operational tasks. This approach helps identify potential design flaws, operational challenges, and safety considerations. It also provides valuable insights for risk reduction and helps improve mission efficiency and effectiveness. With advanced motion tracking systems and custom virtual environments, data can be gathered and evaluated to help NASA refine training protocols, develop specialized training procedures, and optimize human-robotic interactions for future space missions. Furthermore, immersive technology offers opportunities for future training initiatives at NASA. The Virtual Reality Laboratory at JSC has pioneered training with Virtual Reality (VR) since the Hubble Space Telescope repair missions in the early 1990’s. Extended Reality (XR) simulations enable astronauts to rehearse complex spacewalks, spacecraft maneuvers, and extravehicular activities in a safe and controlled environment. By replicating the physical and cognitive challenges of space missions, immersive training experiences enhance astronauts’ situational awareness, decision-making abilities, and adaptability to unexpected scenarios. Additionally, immersive technology facilitates collaborative training, allowing geographically dispersed crew and mission control personnel to engage in synchronized simulations, fostering teamwork and effective communication. The adoption of immersive technology in NASA’s testing, evaluation, and future training programs has yielded significant benefits. By incorporating human-in-the-loop simulations for studies involving Extra Vehicular Activities (EVA), surface mobility and landing systems, NASA can identify and mitigate risks, optimize operational procedures, and enhance mission success. Ultimately, immersive simulation experiences can empower astronauts to better navigate the complexities of space missions, ensuring their safety, productivity, and success in the dynamic and challenging environments they will experience at the Lunar South Pole.

HITL↗

High Temperature Dielectric Properties and Differential Scanning Calorimetry of Lunar Simulants

To guide development of microwave process technology that could be used during in situ construction on the Moon, we measured the high-temperature basic dielectric properties (εʹ and εʺ) of 17 lunar simulants and related materials. In order to confidently use these data one needs to understand the data’s strengths and weaknesses. Therefore, a goal of this publication is to provide insights into the comparative effects of sample composition, pre-treatments, experimental variables, high temperatures, and other factors on the measured response. The dielectric measurements were performed using the cavity perturbation method over a temperature range between room temperature to 1000 °C, or higher, and provided the real and imaginary components of permittivity at six frequencies. The utility of the original values was limited by the varying density of the pellets used in the measurement. Therefore, all of the εʹ and εʺ measurements at the frequency of 2466 MHz have been scaled to a constant density, 1.75 g/cm 3 . Here the data are presented as graphs chosen to aid analysis within and across simulant groups. To gain additional insight into the processes happening at the elevated temperatures in the dielectric measurements, heat capacity data was obtained using differential scanning calorimetry (DSC) on several of the simulant materials. Our data show that over the frequency range 397 MHz – 2985 MHz a material’s behavior does not greatly change, as compared to the scale of differences observed between lunar mare and highland simulants at high temperatures. For example at 1000 °C, the mare simulant JSC-1A absorbs 10 times more power than the highland simulant NUW-LHT-5M. We observe that as melting temperatures are reached both permittivity and dielectric loss rise non-linearly, helping to explain thermal runaway during microware heating. Our data show that even less than a few weight % of many non-lunar minerals, and the use of mixtures in simulants can affect the dielectric behavior at higher temperatures. A comparison of our results with published dielectric data for Apollo samples and with remote sensing of the Moon supports the conclusion the simulants and lunar material at room temperature have very similar dielectric values.

Differential Scanning Calorimetry↗

Evaluation of the Inertial Response of Variable-Speed Wind Turbines Using Advanced Simulation: Preprint

In this paper, we focus on the temporary frequency support effect provided by wind turbine generators (WTGs) through the inertial response. With the implemented inertial control methods, the WTG is capable of increasing its active power output by releasing parts of the stored kinetic energy when the frequency excursion occurs. The active power can be boosted temporarily above the maximum power points, but the rotor speed deceleration follows and an active power output deficiency occurs during the restoration of rotor kinetic energy. In this paper, we evaluate and compare the inertial response induced by two distinct inertial control methods using advanced simulation. In the first stage, the proposed inertial control methods are analyzed in offline simulation. Using an advanced wind turbine simulation program, FAST with TurbSim, the response of the researched wind turbine is comprehensively evaluated under turbulent wind conditions, and the impact on the turbine mechanical components are assessed. In the second stage, the inertial control is deployed on a real 600-kW wind turbine, the three-bladed Controls Advanced Research Turbine, which further verifies the inertial control through a hardware-in-the-loop simulation. Various inertial control methods can be effectively evaluated based on the proposed two-stage simulation platform, which combines the offline simulation and real-time hardware-in-the-loop simulation. The simulation results also provide insights in designing inertial control for WTGs.

17 WIND ENERGY↗

Modeling Integrated Electricity and Natural Gas Systems Using Co-Simulation

One approach for modeling natural gas and electricity networks as integrated systems is co-simulation, which involves separating a single model into subproblems that are solved iteratively, with coupling information passed between then. Although a fully integrated problem might be more direct, co-simulation offers advantages in that it can sometimes be more computationally tractable, and can also leverage existing simulation tools to represent the subproblems. Here we present a framework for co-simulation of natural gas and electricity networks using the HELICS tool. We provide validation of the co-simulation approach by comparing results with a fully integrated approach that uses the same simulation tool, finding that the two approaches provide identical results for a range of test cases. We then present results from natural gas and electricity network co-simulation case study of the Colorado Front Range, with an evaluation of the present-day system and also a future system with higher shares of renewable power. The presentation will close with a discussion on future avenues of work for co-simulation of gas and power networks.

co-simulation↗

Kinetic Monte Carlo Simulation of Hydrogen Diffusion in Tungsten

Nuclear fusion is considered as the most promising sustainable energy source. If achieved, it could provide electricity to the mankind for millions of years. During the reaction, the fusion ash is continuously removed from the magnetic gap, the divertor region, and the divertor facing component is subjected to intense plasma bombardment. It is proposed to use tungsten as the divertor facing material, because it has low physical sputtering yield and superior thermomechanical properties. However, the high mobility of hydrogen in tungsten and the high bombardment flux bring the hydrogen isotopes deep into the tungsten tile. The presence of hydrogen in tungsten tends to weaken the gram boundaries and produce surface blistering. The ejection of tungsten particulates into the plasma core will cool the reaction and cease the fusion reactions. Therefore, it is important to have a detailed understanding of the hydrogen diffusion in tungsten bulk. The KMC method is a stochastic computer simulation to predict the time evolution of a certain system. Typically, at a given time, the system will evolve through a few possible paths with known probabilities. In each step, the KMC method analyzes the entire system and randomly determines the system evolving events based on all transition rates of all possible events and advance the physical time. This process is repeated till the simulation of the entire physical time is completed. In previous work, the authors developed a kinetic Monte Carlo (KMC) algorithm and program to simulate hydrogen diffusion process on tungsten reconstructed surface. In this article, this KMC method is extended to model the hydrogen diffusion in tungsten bulk. The diffusion coefficients are calculated and fitted to the Arrhenius equation. The result matches the experiment measurement very well. KMC diffusion simulation needs the locations of the hydrogen trapping sites, diffusion paths and the corresponding diffusion energy barriers. The hydrogen interstitial sites in tungsten crystal have been thoroughly studied by density function theory (DFT) and molecular dynamics (MD). Johnson and Carter used DFT to determine the preferred location of the hydrogen interstitial sites m tungsten crystal, which are the tetrahedral sites (T-sites) in body-centered-cubic (bcc) lattice surfaces. All T-sites are located on the surfaces of a lattice cell. There are four T-sites per surface and twelve T-sites per lattice cell. The hydrogen atoms diffuse by traveling from one T-site to one of its nearest four T-site neighbors. Thus, each T-site is surrounded by two neighbors on the same surface, and one neighbor on each of the two adjacent surfaces perpendicular to the surface holding the center T-site. The diffusion barrier between two neighboring T-sites was calculated to be 0.38 eV. Knowing the diffusion path, energy barrier, and the crystal structure, the KMC simulation is readily to be performed. Based on other researchers' DFT calculations, the hydrogen stable interstitial site in tungsten bulk is determined as the tetrahedral site on the surface of the bee lattice. Each T-site has 4 neighbors and the energy barrier between two adjacent T-sites are estimated to be 0.38 eV. Based on those facts and the KMC theory, a MATLAB program is developed to model the hydrogen diffusion in tungsten bulk. In the MATLAB program, substrate temperature, hydrogen concentration, physical time to be simulated, and the number of cycles are all customizable. To compare the KMC determined diffusion coefficients with experiment results in the literature, the simulations are performed for temperature 300 K - 2500 K and low hydrogen concentration (0.1%), mimicking the experiment condition. The calculated diffusion coefficients match the reference values measured by the experiment under low hydrogen concentration condition. The results are also fitted to the Arrhenius equation. The pre-exponential factor and the average activation energy in Arrhenius equation are obtained. The good agreement between simulation and experiment proves that the developed KMC algorithm truly reflect the well-established theory of interstitial hydrogens in tungsten. The developed program could easily be used to model the diffusion process and predict the diffusion coefficients under other concentrations and temperatures. This work is part of the ongoing development of hydrogen diffusion in tungsten integrating surface diffusion (previous work), bulk diffusion (this work), and adsorption/desorption on surface (future work)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of MAVRIC/MONACO and TRIPOLI-4{sup R} simulations of the LLNL pulsed spheres benchmark experiments - 25403

Simulations of the Lawrence Livermore Pulsed Spheres with MAVRIC/Monaco and TRIPOLI-4{sup R} are presented. Both codes use ENDF/B-VII.1 cross section data, but TRIPOLI-4{sup R} also uses cross section data based on JEFF-3.1.1. The MAVRIC/Monaco simulations were performed with the SCALE 6.2.2 release while the TRIPOLI-4{sup R} simulations used an internal development version of the code. For the most part, the simulated and measured neutron counts agree within 30%. Simulated gamma flux data is compared, as well. There are no measured gamma data. The agreement between the simulated gamma flux spectra is often within 10% or less, but there some notable differences, especially when comparing simulations using ENDF/B-VII.1 and JEFF 3.1.1 based cross section data. The comparison of photon simulations suggests that improvements are needed in the cross-section data, cross section processing codes, and/or transport codes and gamma measurement data is needed for experimental validation. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗