Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D Simulation”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Electron acceleration in three-dimensional magnetic reconnection with a guide field

Kinetic simulations of 3D collisionless magnetic reconnection with a guide field show a dramatic enhancement of energetic electron production when compared with 2D systems. In the 2D systems, electrons are trapped in magnetic islands that limit their energy gain, whereas in the 3D systems the filamentation of the current layer leads to a stochastic magnetic field that enables the electrons to access volume-filling acceleration regions. The dominant accelerator of the most energetic electrons is a Fermi-like mechanism associated with reflection of charged particles from contracting field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling, analysis, and code/data validation of DIII-D tokamak divertor experiments on ELM and non-ELM plasma tungsten sputtering erosion

Abstract We analyzed recent DIII-D tokamak tungsten divertor probe experiments using advanced, coupled, sputter erosion/redeposition, plasma, and surface response code packages. Modeling is done for ELMing H-mode, and L-mode plasmas, impinging on various size tungsten deposits on Divertor Material Evaluation System (DiMES) carbon probes. The simulations compute 3D, full kinetic, sub-gyromotion, impurity sputtering and transport, including changes in tungsten surface composition and response due to mixed deuterium and carbon ions irradiation. Per our analysis, ELM (edge localized mode) plasma sputtering in DIII-D mostly involves free-streaming high energy (∼500–1000 eV) D + and C +6 ions, with high near-surface plasma density. L-Mode sputtering is due to impurity sputtering (C, W) only, with lower density. All cases show complete redeposition of tungsten on the divertor, with significant redeposition on the tungsten spots themselves, and low self-sputtering. Comparison of ELM plasma gross tungsten erosion simulation results with in-situ spectroscopic data is good, as are code/data comparisons of net erosion using post-exposure Rutherford backscattering (RBS) data for the L-mode probes. The analysis, extrapolated to a full tungsten divertor, implies low net erosion and negligible plasma contamination from sputtering. These results support the use of high-Z plasma facing surfaces in ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Slipping magnetic reconnection during an X-class solar flare observed by SDO/AIA

We present SDO/AIA observations of an eruptive X-class flare of 2012 July 12, and compare its evolution with the predictions of a three-dimensional (3D) numerical simulation. We focus on the dynamics of flare loops that are seen to undergo slipping reconnection during the flare. In the Atmospheric Imaging Assembly (AIA) 131 Å observations, lower parts of 10 MK flare loops exhibit an apparent motion with velocities of several tens of km s{sup –1} along the developing flare ribbons. In the early stages of the flare, flare ribbons consist of compact, localized bright transition-region emission from the footpoints of the flare loops. A differential emission measure analysis shows that the flare loops have temperatures up to the formation of Fe XXIV. A series of very long, S-shaped loops erupt, leading to a coronal mass ejection observed by STEREO. The observed dynamics are compared with the evolution of magnetic structures in the 'standard solar flare model in 3D.' This model matches the observations well, reproducing the apparently slipping flare loops, S-shaped erupting loops, and the evolution of flare ribbons. All of these processes are explained via 3D reconnection mechanisms resulting from the expansion of a torus-unstable flux rope. The AIA observations and the numerical model are complemented by radio observations showing a noise storm in the metric range. Dm-drifting pulsation structures occurring during the eruption indicate plasmoid ejection and enhancement of the reconnection rate. The bursty nature of radio emission shows that the slipping reconnection is still intermittent, although it is observed to persist for more than an hour.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A new method for analyzing and visualizing plasma simulations using a phase-space tessellation

We apply a novel phase-space interpolation technique referred to as the simplex-in-cell (SIC) method to analyze two- and three-dimensional particle-in-cell (PIC) simulations of electromagnetic plasmas. SIC relies on a discretization of the initial phase-space distribution function into simplices, which allows an approximation to the full, continuously defined distribution function to be constructed at any later time in the simulation. This allows densities, currents, and even full momentum distribution functions to be measured at any point in the simulation domain without averaging over control volumes. The SIC approach applies to any PIC simulation for which a tessellation of the initial particle distribution can be constructed. In this study, we use outputs from standard PIC simulations of the Weibel instability and compare physical quantities such as charge and current densities calculated in postprocessing using SIC and standard particle deposits. Using 2D simulations with 1–65 536 particles-per-cell, we find that SIC eliminates discrete particle noise and in some cases can reach a given noise level using ~1000 times fewer simulation particles than with standard particle deposition schemes. In regions of low density, such as between current filaments, SIC is able to capture small amplitude features even with fewer particles than gridpoints due to the deformable nature of the SIC volume elements. Here, by calculating momentum distributions, we show how SIC can capture low density tails in the spectrum using far fewer particles than are necessary for standard particle deposits. We calculate the charge density on spatial grids of increasing resolution to demonstrate the ability of SIC to reveal fine-scale details that are not accessible with standard particle deposits. Finally, we show how SIC can be extended to 3D and give an example of its use to calculate the charge density from 3D PIC simulations of the Weibel instability. These results motivate the future implementation of SIC directly in the simulation force calculation for a novel low-noise electromagnetic plasma simulation method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2021

Application of Radon Transform to Multi-Angle Measurements Made by the Research Scanning Polarimeter: A New Approach to Cloud Tomography. Part I: Theory and Tests on Simulated Data.

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.

clouds↗
DOE OSTI2023

Transition in particle transport under resonant magnetic perturbations in a tokamak

Abstract Nonlinear 3D MHD simulations and validations reveal that the hybrid particle-MHD transport is a key process for driving the pump-out in the presence of Resonant Magnetic Perturbations (RMPs) in the KSTAR tokamak. Particle transport and the resulting density pump-out by RMPs are shown to be composed of not only the classical flow convection near magnetic islands due to polarization but also the neoclassical ion diffusion across perturbed magnetic surfaces. The latter is known as the Neoclassical Toroidal Viscosity (NTV) and is integrated into nonlinear MHD simulations here for the first time, revealing that the two-stage pump-outs observed in KSTAR experiments are reproduced only with such integrated nonlinear MHD and transport evolution. Near-resonant responses, which have received less attention than the resonant response, play distinct roles in the pump-out along with the island formation. In addition, this modeling is used to investigate the pump-outs in double-null-like plasmas and numerically capture the effect of the double-null shape on the pump-outs, which may explain the difficulty of Edge Localized Mode (ELM) suppression access in double-like plasmas. This reveals new aspects of the impact toroidal geometry and mode coupling have on 3D physics and reveals the importance of near-resonant components in suppressing ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On the physical and chemical details of alumina atomic layer deposition: A combined experimental and numerical approach

Alumina thin film is typically studied as a model atomic layer deposition (ALD) process due to its high dielectric constant, high thermal stability, and good adhesion on various wafer surfaces. Despite extensive applications of alumina ALD in microelectronics industries, details on the physical and chemical processes are not yet well understood. ALD experiments are not able to shed adequate light on the detailed information regarding the transient ALD process. Most of current numerical approaches lack detailed surface reaction mechanisms, and their results are not well correlated with experimental observations. In this paper, the authors present a combined experimental and numerical study on the details of flow and surface reactions in alumina ALD using trimethylaluminum and water as precursors. Results obtained from experiments and simulations are compared and correlated. By experiments, growth rate on five samples under different deposition conditions is characterized. The deposition rate from numerical simulation agrees well with the experimental results. Details of precursor distributions in a full cycle of ALD are studied numerically to bridge between experimental observations and simulations. The 3D transient numerical model adopts surface reaction kinetics and mechanisms based on atomic-level studies to investigate the surface deposition process. Surface deposition is shown as a strictly self-limited process in our numerical studies. ALD is a complex strong-coupled fluid, thermal and chemical process, which is not only heavily dependent on the chemical kinetics and surface conditions but also on the flow and material distributions.

42 ENGINEERING↗
DOE OSTI2014

Characterizing the convective velocity fields in massive stars

We apply the mathematical formalism of vector spherical harmonics decomposition to convective stellar velocity fields from multidimensional hydrodynamics simulations and show that the resulting power spectra furnish a robust and stable statistical description of stellar convective turbulence. Analysis of the power spectra helps identify key physical parameters of the convective process such as the dominant scale of the turbulent motions that influence the structure of massive evolved pre-supernova stars. We introduce the numerical method that can be used to calculate vector spherical harmonics power spectra from two-dimensional (2D) and three-dimensional (3D) convective shell simulation data. Using this method we study the properties of oxygen shell burning and convection for a 15 M {sub ☉} star simulated by the hydrodynamics code FLASH in 2D and 3D. We discuss the importance of realistic initial conditions to achieving successful core-collapse supernova explosions in multidimensional simulations. We show that the calculated power spectra can be used to generate realizations of the velocity fields of presupernova convective shells. We find that the slope of the solenoidal mode power spectrum remains mostly constant throughout the evolution of convection in the oxygen shell in both 2D and 3D simulations. We also find that the characteristic radial scales of the convective elements are smaller in 3D than in 2D, while the angular scales are larger in 3D.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Thermodynamic balance of three-dimensional stratospheric winds derived from a data assimilation procedure

The NASA/Goddard three-dimensional chemistry and transport model is driven by winds from a stratospheric data assimilation system. Synoptic- and planetary-scale patterns, apparent in satellite observations of trace constituents, are successfully reproduced for seasonal integrations. As model integrations proceed, however, the quality of simulations decreases, and systematic differences between calculation and measurement appear. The differences are explained by examining the zonal-mean residual circulation. The vertical velocity w-bar (sup star) is calculated two ways: (1) from the diabatic heating rates and temperature tendency and (2) from the Eulerian vertical velocity and the horizontal eddy heat flux convergence. The results from these calculations differ substantially. Periodic insertion of observational data during the assimilation process continually shocks the general circulation model and produces these differences, which leads to an overestimate of the mean vertical heat and cconstituent transport. Such differences are expected to be general to all data assimilation products. This interpretation is corroborated by two-dimensional (2D) model calculations. When w-bar(sup star) is calculated from (2), the 2D ozone evolution is unrealistic and qualitatively similar to the 3D model simulation. The 2D ozone evolution is reasonable when w-bar (sup star) from (1).

Weaver, Clark J.↗
DOE OSTI2016

Modeling Heat Loss through Piston and Effects of Thermal Boundary Coatings in Diesel Engine Simulations using Conjugate Heat Transfer models

Heat loss through wall boundaries play a dominant role in the overall performance and efficiency of internal combustion engines. Typical engine simulations use constant temperature wall boundary conditions. These boundary conditions cannot be estimated accurately from experiments due to the complexities involved with engine combustion. As a result they introduce a large uncertainty in engine simulations and serve as a tuning parameter. Modeling the process of heat transfer through the solid walls in an unsteady engine computational fluid dynamics (CFD) simulation can lead to the development of higher fidelity engine calculations. These models can be used to study the impact of heat loss on engine efficiency and explore new design methodologies that can reduce heat losses. In this work, a single cylinder diesel engine is modeled along with the solid piston coupled to the fluid domain. Conjugate heat transfer (CHT) modeling techniques were implemented to model heat losses for a full cycle of a Navistar diesel engine. This CFD model is then validated against experimental data available from thermocouples embedded inside the piston surface. The overall predictions from the model match closely with the experimental observations. The validated model is further used to explore the benefits of thermal barrier coatings (TBC) on piston bowls. The effect of TBC coatings were modeled as a thermal resistance in the heat transfer models. Full cycle 3D engine simulations provide quantitative insights into heat loss and thus calculate the efficiency gain by the use of TBC coatings. The work establishes a validated modeling framework for CHT modeling in reciprocating engine simulations.

Kundu, Prithwish↗
DOE OSTI2016

Scaled laboratory experiments explain the kink behaviour of the Crab Nebula jet

X-ray images from the Chandra X-ray Observatory show that the South-East jet in the Crab nebula changes direction every few years. This remarkable phenomenon is also observed in jets associated with pulsar wind nebulae and other astrophysical objects, and therefore is a fundamental feature of astrophysical jet evolution that needs to be understood. Theoretical modeling and numerical simulations have suggested that this phenomenon may be a consequence of magnetic fields (B) and current-driven magnetohydrodynamic (MHD) instabilities taking place in the jet, but until now there has been no verification of this process in a controlled laboratory environment. Here we report the first such experiments, using scaled laboratory plasma jets generated by high-power lasers to model the Crab jet and monoenergetic-proton radiography to provide direct visualization and measurement of magnetic fields and their behavior. The toroidal magnetic field embedded in the supersonic jet triggered plasma instabilities and resulted in considerable deflections throughout the jet propagation, mimicking the kinks in the Crab jet. We also demonstrated that these kinks are stabilized by high jet velocity, consistent with the observation that instabilities alter the jet orientation but do not disrupt the overall jet structure. We successfully modeled these laboratory experiments with a validated three-dimensional (3D) numerical simulation, which in conjunction with the experiments provide compelling evidence that we have an accurate model of the most important physics of magnetic fields and MHD instabilities in the observed, kinked jet in the Crab nebula. The experiments initiate a novel approach in the laboratory for visualizing fields and instabilities associated with jets observed in various astrophysical objects, ranging from stellar to extragalactic systems. We expect that future work along this line will have important impact on the study and understanding of such fundamental astrophysical phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

THE STRUCTURE OF SPIRAL SHOCKS EXCITED BY PLANETARY-MASS COMPANIONS

Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three-dimensional (3D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the nonlinear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3D hydrodynamical models than images based on two-dimensional models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion’s position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion’s mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

On Differential Rotation and Overshooting in Solar-like Stars

We seek to characterize how the change of global rotation rate influences the overall dynamics and large-scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such a change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the last profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number R{sub of}= ω-tilde /2Ω{sub ∗}∼ v-tilde /2Ω{sub ∗}R{sub ∗}, which we calibrate based on our 3D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map (R {sub of} ≳ 1.5–2) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Felippa, C. A.↗
DOE OSTI2016

Study on air ingress phenomena during a depressurization accident of a VHTR - 18632

A depressurization accident is the one of the design-basis accidents of a Very-High-Temperature Reactor (VHTR). When a depressurization accident occurs, air is expected to enter into the reactor pressure vessel from the breach and oxidize in-core graphite structures. Therefore, it is important to understand the mixing processes of different kinds of gases in the stable and unstable stratified fluid layers. In particular, it is also important to examine the influence of localized natural convection and molecular diffusion on the mixing process from a safety viewpoint. Therefore, in order to predict or analyze the air ingress phenomena during a depressurization accident, it is important to develop a method for the prevention of air ingress during an accident. We carried out experiments and numerical analysis using three-dimensional (3D) CFD code to obtain the mixing process of two-component gases and the flow characteristics of localized natural convection. This study also investigated a control method for the natural circulation of air through the injection of helium gas. The numerical model consists of a storage tank and a reverse U-shaped vertical rectangular passage. They are separated by a horizontal partition plate. One sidewall of the high-temperature side vertical passage is heated and the other sidewall is cooled. The low-temperature vertical passage is cooled by ambient air. The experimental procedures and numerical analysis process are as follows. Firstly, the storage tank is filled with heavy gas and the reverse U-shaped vertical passage is filled with a light gas. In the vertical passage of the high-temperature side, localized natural convection is generated by the temperature difference between the vertical walls. These flow characteristics were obtained from the experiments and steady-state analysis. The experiment and analysis began after the partition plate was opened. The results obtained from the experiments were quantitatively simulated using 3D numerical analysis. The two component gases were mixed via molecular diffusion and natural convection. After some time elapsed, natural circulation occurred through the reverse U-shaped vertical passage. These flow characteristics are the same as those of phenomena generated in the passage between a permanent reflector and a pressure vessel wall of the VHTR. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1996

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗
DOE OSTI2022

The HUNTER Dynamic Human Reliability Analysis Tool: Coupling an External Plant Code

The Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) is a framework to support dynamic human reliability analysis (HRA) with the aim to develop standalone software to perform the dynamic HRA calculations. Within the HRA, human actions in nuclear power plants (NPPs) are predicated by plant states, and human actions influence the plant. In other words, plant operations are necessarily recursive, and it becomes challenging to model complex human-plant interactions. Consequently, we have linked two software simulations that complement those shortcomings. RELAP5-3D—the Reactor Excursion and Leak Analysis Program (RELAP; Aumiller, Tomlinson, and Bauer 2001) is the foundational thermal-hydraulic software used to model nuclear systems. Using RELAP5-3D, we have simulated the plant operations proceeding according to procedures developed to address emergent situations in NPPs. Plant operations include various actions such as the operator checking plant parameters, as well as actions that are continuously performed over time until a specific parameter reaches certain criteria. That means that HUNTER and RELAP5-3D exchange information with each other and should be carried out simultaneously over time. To simulate plant operations, which represent the actual operator checks of plant parameters and corresponding manual control actions, changes in plant status are identified through simulation and performed according to the criteria and order of the procedure. Thus, the goal of coupling HUNTER with RELAP5-3D is to facilitate synchronous coupling, where human and plant models provide iterative feedback loops that drive the course of actions. The advantage of coupling with RELAP5-3D to serve as the external environment module in HUNTER is the ability to customize the plant model and streamline for particular model applications. In this paper, we will address the key features of the coupling and the coupling structure built to perform the feedback loops.

99 GENERAL AND MISCELLANEOUS↗