Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized 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 1,189 records · Page 66

VULCAN: An Open-source, Validated Chemical Kinetics Python Code for Exoplanetary Atmospheres

We present an open-source and validated chemical kinetics code for studying hot exoplanetary atmospheres, which we name VULCAN. It is constructed for gaseous chemistry from 500 to 2500 K, using a reduced C–H–O chemical network with about 300 reactions. It uses eddy diffusion to mimic atmospheric dynamics and excludes photochemistry. We have provided a full description of the rate coefficients and thermodynamic data used. We validate VULCAN by reproducing chemical equilibrium and by comparing its output versus the disequilibrium-chemistry calculations of Moses et al. and Rimmer and Helling. It reproduces the models of HD 189733b and HD 209458b by Moses et al., which employ a network with nearly 1600 reactions. We also use VULCAN to examine the theoretical trends produced when the temperature–pressure profile and carbon-to-oxygen ratio are varied. Assisted by a sensitivity test designed to identify the key reactions responsible for producing a specific molecule, we revisit the quenching approximation and find that it is accurate for methane but breaks down for acetylene, because the disequilibrium abundance of acetylene is not directly determined by transport-induced quenching, but is rather indirectly controlled by the disequilibrium abundance of methane. Therefore we suggest that the quenching approximation should be used with caution and must always be checked against a chemical kinetics calculation. A one-dimensional model atmosphere with 100 layers, computed using VULCAN, typically takes several minutes to complete. VULCAN is part of the Exoclimes Simulation Platform (ESP; exoclime.net) and publicly available at https://github.com/exoclime/VULCAN.

79 ASTRONOMY AND ASTROPHYSICS↗

Coordinated control design and verification of the HTR-PM plant - 18485

Through multi-modular scheme, i.e. the superheated steam flows produced by multiple MHTGR-based nuclear supplying system (NSSS) modules are combined together to drive a common thermal load, the inherent safety feature of MHTGR can be applicable to the large-scale nuclear plants at any desired power ratings. Since the plant power control technique of classical single-modular plants cannot be directly applied to multi-modular plants, it is necessary to develop the power control method of the multi-modular plants. In this paper, the plant power control strategy of the two-modular HTR-PM plant, i.e. the first under-constructed multi-modular high temperature gas-cooled nuclear plant in the world is given, which includes both the feedback loops and control algorithms. Then, the control strategy is verified through numerical simulation, and the simulation results in the cases of plant power step and ramp show the satisfactory performance of this control strategy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A simulation of a steam generator tube rupture in AP1000 nuclear power plant based on RELAP5-MV

In the operation of nuclear power plants, it is very important to evaluate different accident scenarios in actual plant conditions with different codes. In the field of nuclear safety, steam generator tube rupture (SGTR) is one of the main accidents. In this research, transient analysis of the primary system variation of thermo-hydraulics parameters in primary loop under SGTR accident in AP1000 nuclear power plant (NPP) is analyzed by Relap5-MV thermo-hydraulics code. The transient test scenario is modeled using the RELAP5-MV for the AP1000 NPP. The calculation was performed up to a transient time of 3,000.0 s. The results obtained from RELAP5 and RELAP5-MV are in agreement as they are based on the same models though in comparison with RELAP5, RELAP5-MV makes simulation of nuclear power systems easier and convenient for users most especially beginners. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Selective Shielding Bone Marrow: An Approach to Create Personal Protective Equipment for Gamma Radiation

Many casualties of the Hiroshima and Nagasaki atomic bombs, and many of the firefighters who first responded to the Chernobyl nuclear power plant accident, became ill with acute radiation syndrome (ARS). Most of the people who do not recover from ARS will die within a few weeks to a few months after exposure, with the primary cause of death being the destruction of the body's bone marrow (BM). This is exemplified by the Chernobyl firefighters, the majority of whom received whole body doses between 0.8 and 10 Gy, and the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) final report on the disaster concluded that the underlying cause of death of all 28 firefighters who succumbed to ARS was from bone marrow failure. BM is comprised of hematopoietic stem cells (HSCs) which are responsible for the constant renewal of blood cells. Due to their high rate of proliferation, HSCs are especially vulnerable to ionizing radiation, but are endowed with remarkable regenerative potential. Owing to their central role in blood production, lethal irradiation of HSCs leads to death from severe anemia, infection and internal bleeding. This relationship between high doses of radiation and HSC apoptosis has led to the use of HSC transplantation as a life-saving intervention in cases of acute exposure. This is demonstrable in medical practice, where thousands of individuals have undergone supra-lethal Total Body Irradiation (TBI) for purposes of cancer therapy, and were rescued by well-matched BM transplantation. Indeed, life-threatening damage may be reversed by BM transplantation in individuals receiving doses of radiation of up to 10 Gy. At doses > 10 Gy, damage to gastrointestinal (GI) tissue may become a limiting factor in survival. Remarkably, in the common procedure of BM transplantation, the number of HSCs extracted from a single active BM site containing less than 5% of the donor's BM tissue is sufficient to support the complete reconstitution of the HSC compartment in a lethally irradiated recipient. This capacity of BM to expand and replenish is due to the high regenerative potential of the HSCs it harbors. In contrast to mature blood cells, which are dispersed throughout the body, HSCs are confined to the bones, allowing for effective targeted shielding. This has been confirmed in several animal models, where sparing the HSC-rich area of a subject receiving otherwise full-body irradiation is sufficient to support hematopoietic functions and allow survival. Approximately 48% of the human body's active BM is contained within the lumbar vertebrae and the pelvic region. This high concentration of BM makes the pelvic region an ideal area to shield for preventing the serious hematopoietic effects of radiation exposure. Stem-cell rich tissues like active BM are also more sensitive to the stochastic effects (i.e. cancer) of radiation exposure. The high amount of BM in the pelvic region, combined with the presence of the sensitive gastrointestinal system and female gonads in this same area of the body allows for the possibility of significantly reducing stochastic effects of radiation by limiting radiation exposure to this vital area. Thus, shielding this region holds great promise for both acute and protracted exposures. Thus, a specialized radiation shielding device for the protection of active BM concentrations is presented here. This belt-like selective shielding device focuses on the protection of BM that is present in the pelvis and hip bones, protecting the medullary volume from which BM is commonly extracted for transplantation (i.e. the iliac bones) while allowing relatively unhindered movement of the wearer. In order to optimize the use of shielding materials towards the protection of active BM, the shielding uniquely brings into account the natural shielding properties of human tissue, with thickness being inversely related to the thickness and radiation attenuation of the underlying tissue at each point surrounding the area being protected. The findings resulting from both experimental and simulation testing of this selective shielding device are presented here.

61 RADIATION PROTECTION AND DOSIMETRY↗

Modeling of the SPERT Transients Using SERPENT2 with Time Dependent Capabilities

Over the last decade, multiphysics coupled calculations have become one of the most important research topics in the area of development and enhancement of reactor analysis and simulation tools. One of the main applications of these coupled simulations is transient analysis, which may involve modeling of experiments, normal plant operation or accident scenarios. Traditionally, the reactor physics side of coupled simulations has been represented by deterministic or hybrid Monte Carlo - deterministic methods. Monte Carlo techniques are limited by their high computational cost. However increased computing power makes their application to multiphysics modeling of transients feasible. Active research is underway to develop algorithms for time dependent calculations, implement them in Monte Carlo codes, and evaluate their performance. In this paper time dependent simulations of the SPERT experiments using the Serpent2 code are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transient analysis of a micro-reactor using the DireWolf code suite

Transient analyses of heat pipe micro-reactors are necessary to ensure that hypothetical accident scenarios do not comprise reactor safety. Due to its small size and reliance on heat-pipes for cooling, the micro-reactor design introduced in this paper is a tightly coupled system which requires multi-physics tools to accurately model transient events. Idaho National Laboratory's DireWolf code suite based on the MOOSE framework is tailor-built to model heat-pipe reactors. This paper demonstrates DireWolf's ability to simulate the coupled thermal-neutronics transient behavior of a heat-pipe micro-reactor. The transient events presented here include an inadvertent rotation of all control drums simultaneously and a sudden complete rotation of a single control drum. A detailed description of each event is provided along with simulation results, including time dependent power and temperature distributions, and discussion. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Mapping the current–current correlation function near a quantum critical point

The current–current correlation function is a useful concept in the theory of electron transport in homogeneous solids. The finite-temperature conductivity tensor as well as Anderson’s localization length can be computed entirely from this correlation function. Based on the critical behavior of these two physical quantities near the plateau–insulator or plateau–plateau transitions in the integer quantum Hall effect, we derive an asymptotic formula for the current–current correlation function, which enables us to make several theoretical predictions about its generic behavior. For the disordered Hofstadter model, we employ numerical simulations to map the current–current correlation function, obtain its asymptotic form near a critical point and confirm the theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A consistent and conservative scheme for MHD flows with complex boundaries on an unstructured Cartesian adaptive system

The numerical simulation of Magnetohydrodynamics (MHD) flows with complex boundaries has been a topic of great interest in the development of a fusion reactor blanket for the difficulty to accurately simulate the Hartmann layers and side layers along arbitrary geometries. An adaptive version of a consistent and conservative scheme has been developed for simulating the MHD flows. Besides, the present study forms the first attempt to apply the cut-cell approach for irregular wall-bounded MHD flows, which is more flexible and conveniently implemented under adaptive mesh refinement (AMR) technique. It employs a Volume-of-Fluid (VOF) approach to represent the fluid–conducting wall interface that makes it possible to solve the fluid–solid coupling magnetic problems, emphasizing at how electric field solver is implemented when conductivity is discontinuous in cut-cell. For the irregular cut-cells, the conservative interpolation technique is applied to calculate the Lorentz force at cell-center. On the other hand, it will be shown how consistent and conservative scheme is implemented on fine/coarse mesh boundaries when using AMR technique. Then, the applied numerical schemes are validated by five test simulations and excellent agreement was obtained for all the cases considered, simultaneously showed good consistency and conservative properties.

42 ENGINEERING↗

A spectral unaveraged algorithm for free electron laser simulations

We propose and discuss a numerical method to model electromagnetic emission from the oscillating relativistic charged particles and its coherent amplification. The developed technique is well suited for free electron laser simulations, but it may also be useful for a wider range of physical problems involving resonant field–particles interactions. The algorithm integrates the unaveraged coupled equations for the particles and the electromagnetic fields in a discrete spectral domain. Using this algorithm, it is possible to perform full three-dimensional or axisymmetric simulations of short-wavelength amplification. In this paper we describe the method, its implementation, and we present examples of free electron laser simulations comparing the results with the ones provided by commonly known free electron laser codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An improved random walk algorithm for the implicit Monte Carlo method

In this work, we introduce a modified Implicit Monte Carlo (IMC) Random Walk (RW) algorithm, which increases simulation efficiency for multigroup radiative transfer problems with strongly frequency-dependent opacities. To date, the RW method has only been implemented in “fully-gray” form; that is, the multigroup IMC opacities are group-collapsed over the full frequency domain of the problem to obtain a gray diffusion problem for RW. This formulation works well for problems with large spatial cells and/or opacities that are weakly dependent on frequency; however, the efficiency of the RW method degrades when the spatial cells are thin or the opacities are a strong function of frequency. To address this inefficiency, we introduce a RW frequency group cutoff in each spatial cell, which divides the frequency domain into optically thick and optically thin components. In the modified algorithm, opacities for the RW diffusion problem are obtained by group-collapsing IMC opacities below the frequency group cutoff. Particles with frequencies above the cutoff are transported via standard IMC, while particles below the cutoff are eligible for RW. This greatly increases the total number of RW steps taken per IMC time-step, which in turn improves the efficiency of the simulation. We refer to this new method as Partially-Gray Random Walk (PGRW). We present numerical results for several multigroup radiative transfer problems, which show that the PGRW method is significantly more efficient than standard RW for several problems of interest. In general, PGRW decreases runtimes by a factor of ∼2–4 compared to standard RW, and a factor of ∼3–6 compared to standard IMC. While PGRW is slower than frequency-dependent Discrete Diffusion Monte Carlo (DDMC), it is also easier to adapt to unstructured meshes and can be used in spatial cells where DDMC is not applicable. This suggests that it may be optimal to employ both DDMC and PGRW in a single simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dense plasma heating and Gbar shock formation by a high intensity flux of energetic electrons

Process of shock ignition in inertial confinement fusion implies creation of a high pressure shock with a laser spike having intensity of the order of a few PW/cm{sup 2}. However, the collisional (Bremsstrahlung) absorption at these intensities is inefficient and a significant part of laser energy is converted in a stream of energetic electrons. The process of shock formation in a dense plasma by an intense electron beam is studied in this paper in a planar geometry. The energy deposition takes place in a fixed mass target layer with the areal density determined by the electron range. A self-similar isothermal rarefaction wave of a fixed mass describes the expanding plasma. Formation of a shock wave in the target under the pressure of expanding plasma is described. The efficiency of electron beam energy conversion into the shock wave energy depends on the fast electron energy and the pulse duration. The model is applied to the laser produced fast electrons. The fast electron energy transport could be the dominant mechanism of ablation pressure creation under the conditions of shock ignition. The shock wave pressure exceeding 1 Gbar during 200–300 ps can be generated with the electron pulse intensity in the range of 5–10 PW/cm{sup 2}. The conclusions of theoretical model are confirmed in numerical simulations with a radiation hydrodynamic code coupled with a fast electron transport module.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulation of ion charge breeding in electron beam ion source

The Electron Beam Ion Source particle-in-cell code (EBIS-PIC) tracks ions in an EBIS electron beam while updating electric potential self-consistently and atomic processes by the Monte Carlo method. Recent improvements to the code are reported in this paper. The ionization module has been improved by using experimental ionization energies and shell effects. The acceptance of injected ions and the emittance of extracted ion beam are calculated by extending EBIS-PIC to the beam line transport region. An EBIS-PIC simulation is performed for a Cs charge-breeding experiment at BNL. The charge state distribution agrees well with experiments, and additional simulation results of radial profiles and velocity space distributions of the trapped ions are presented.

43 PARTICLE ACCELERATORS↗

Full wave simulation of waves in ECRIS plasmas based on the finite element method

This paper describes the modeling and the full wave numerical simulation of electromagnetic waves propagation and absorption in an anisotropic magnetized plasma filling the resonant cavity of an electron cyclotron resonance ion source (ECRIS). The model assumes inhomogeneous, dispersive and tensorial constitutive relations. Maxwell's equations are solved by the finite element method (FEM), using the COMSOL Multiphysics{sup ®} suite. All the relevant details have been considered in the model, including the non uniform external magnetostatic field used for plasma confinement, the local electron density profile resulting in the full-3D non uniform magnetized plasma complex dielectric tensor. The more accurate plasma simulations clearly show the importance of cavity effect on wave propagation and the effects of a resonant surface. These studies are the pillars for an improved ECRIS plasma modeling, that is mandatory to optimize the ion source output (beam intensity distribution and charge state, especially). Any new project concerning the advanced ECRIS design will take benefit by an adequate modeling of self-consistent wave absorption simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear simplified model to study localization of kinetic Alfvén wave

We have presented the numerical simulation of the coupled equations governing the dynamics of kinetic Alfvén wave (KAW) and ion acoustic wave in the intermediate β plasma, where β is the ratio of thermal pressure to the background magnetic pressure. We have also developed a simplified model for this nonlinear interaction using the results obtained from the simulation to understand the physics of nonlinear evolution of KAW. Localization of magnetic field intensity of KAW has been studied by means of the simplified model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interactions of mobile helium clusters with surfaces and grain boundaries of plasma-exposed tungsten

We report results of atomistic computations for the interactions of small mobile helium clusters (He{sub n}) with free surfaces and grain boundaries (GBs) in tungsten toward development of continuum drift-diffusion-reaction models for the dynamics of mobile helium clusters in plasma-exposed tungsten. Molecular-statics (MS) simulations based on reliable many-body interatomic potentials are carried out for He{sub n} (1 ≤ n ≤ 7) clusters near sinks to obtain the potential energy profiles of the He{sub n} clusters as a function of the clusters' center-of-mass distance from a sink. Sinks investigated include surfaces, GBs, and regions in the vicinity of junctions where GBs intersect free surfaces. Elastic interaction potentials based on elastic inclusion theory provide an excellent description of the MS results for the cluster-sink interactions. The key parameter in the elastic models is the sink segregation strength, which is found to increase with increasing cluster size. Such cluster-sink interactions are responsible for the migration of small helium clusters by drift and for helium segregation on surfaces and grain boundaries in tungsten. Such helium segregation on sinks is observed in large-scale molecular-dynamics simulations of helium aggregation in model polycrystalline tungsten at 933 K upon helium implantation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Properties of a soft-core model of methanol: An integral equation theory and computer simulation study

Thermodynamic and structural properties of a coarse-grained model of methanol are examined by Monte Carlo simulations and reference interaction site model (RISM) integral equation theory. Methanol particles are described as dimers formed from an apolar Lennard-Jones sphere, mimicking the methyl group, and a sphere with a core-softened potential as the hydroxyl group. Different closure approximations of the RISM theory are compared and discussed. The liquid structure of methanol is investigated by calculating site-site radial distribution functions and static structure factors for a wide range of temperatures and densities. Results obtained show a good agreement between RISM and Monte Carlo simulations. The phase behavior of methanol is investigated by employing different thermodynamic routes for the calculation of the RISM free energy, drawing gas-liquid coexistence curves that match the simulation data. Preliminary indications for a putative second critical point between two different liquid phases of methanol are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy resolved actinometry for simultaneous measurement of atomic oxygen densities and local mean electron energies in radio-frequency driven plasmas

A diagnostic method for the simultaneous determination of atomic oxygen densities and mean electron energies is demonstrated for an atmospheric pressure radio-frequency plasma jet. The proposed method is based on phase resolved optical emission measurements of the direct and dissociative electron-impact excitation dynamics of three distinct emission lines, namely, Ar 750.4 nm, O 777.4 nm, and O 844.6 nm. The energy dependence of these lines serves as basis for analysis by taking into account two line ratios. In this frame, the method is highly adaptable with regard to pressure and gas composition. Results are benchmarked against independent numerical simulations and two-photon absorption laser-induced fluorescence experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗