Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer codes”

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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Diffraction of radio frequency waves by spatially modulated interfaces in the plasma edge in tokamaks

We report the use of radio frequency (RF) waves in fusion plasmas for heating, for non-inductive current generation, for profile control and for diagnostics has been well established. The RF waves, excited by antenna structures placed near the wall of a fusion device, have to propagate through density fluctuations at the plasma edge. These fluctuations can modify the properties of the RF waves that propagate towards the core of the plasma. A full-wave electromagnetic computational code ScaRF based on the finite difference frequency domain (FDFD) method has been developed to study the effect of density turbulence on RF waves. The anisotropic plasma permittivity used in the scattering studies is that for a magnetized, cold plasma. The code is used to study the propagation of an RF plane wave through a modulated, spatially periodic density interface. Such an interface could arise in the edge region due to magnetohydrodynamic instability or drift waves. The frequency of the plane wave is taken to be in the range of the electron cyclotron frequency. The scattering analysis is applicable to ITER-like plasmas, as well as to plasmas in medium sized tokamaks such as TCV, ASDEX-U and DIII-D. The effect of different density contrasts across the interface and of different spatial modulations are discussed. While ScaRF is used to study a periodic density fluctuation, the code is general enough to include different varieties of density fluctuations in the edge region – such as blobs and filaments, and spatially random fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Microstability of β ~ 1 tokamak equilibria

High-power-density tokamaks offer a potential solution to design cost-effective fusion devices. One way to achieve high power density is to operate at a high β value (the ratio of thermal to magnetic pressure), i.e. β~1 . However, a β~1 state may be unstable to various pressure- and current-driven instabilities or have unfavourable microstability properties. To explore these possibilities, we generate β~1 equilibria and investigate their stability. First, we demonstrate the generation of high-β equilibria with the computer code VMEC. We then analyse these equilibria to determine their stability against the infinite-n ideal-ballooning mode. We follow that by engaging in a detailed microstability study using the GS2 code, beginning with assessments of electrostatic ion-temperature-gradient and trapped election mode instabilities. We observe interesting behaviour for the high-β equilibria – stabilization of these modes through two distinct mechanisms – large negative local shear and reversal of electron precession drift. Finally, we perform electromagnetic gyrokinetic simulations and observe enhanced stability in the outer core of high-β equilibria and absence of kinetic ballooning modes in the negative-triangularity, high-β equilibria. The enhanced outer-core stability of high-β equilibria is different from their lower-β counterparts and offers an alternative, potentially favourable regime of tokamak operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct Air Capture Using Aqueous Amino Acid Solvents in a Crossflow Absorber

Carbon dioxide (CO 2 ) is the most abundant of all greenhouse gases (GHGs). CO 2 levels in the atmosphere are 50% higher than in the preindustrial era, trapping heat. CO 2 removal from the atmosphere by direct air capture (DAC) is needed to achieve the internationally agreed global temperature goals. The most common CO 2 capture technology is absorption by aminebased solvents in packed columns. Amino acid solutions have recently gained attention due to their advantages over traditional amine solvents. To be implemented effectively, DAC industrial processes need to handle large airflow rates in separation absorbers. The large air and solvent flow rates preclude the use of countercurrent columns due to high-pressure drops and the occurrence of flooding. Crossflow air−liquid absorbers are used to handle large air and liquid volumes due to their lower-pressure drop. The objective of this work is to study the influence of crossflow absorber geometric parameters and operating conditions on product formation and process efficiency. An already derived theoretical model for countercurrent absorbers has been modified to simulate the operation of a crossflow DAC absorber. The predictive model was implemented into a computer code that was used to study the efficiency of the processes as geometric equipment dimensions and operating parameters vary. Practical suggestions are made to design more efficient DAC processes.

Absorption↗
DOE OSTI2020

High-Resolution Comonomer Sequencing of Blocky Brominated Syndiotactic Polystyrene Copolymers Using 13C NMR Spectroscopy and Computer Simulations

This work demonstrates the first high-resolution comonomer sequencing of Blocky brominated syndiotactic polystyrene (sPS-co-sPS-Br) copolymers based on pentad assign- ments of the quaternary carbon region of the nuclear magnetic resonance spectrum. Copolymers containing p-bromostyrene (Br-Sty) units were prepared in matched sets using postpolymerization bromination methods carried out in the heterogeneous gel state (Blocky) and homogeneous solution state (Random). Quantitative information from the quaternary carbon spectra, heteronuclear multiple bond correlation spectroscopy, electronic structure calculations, and simulated statistically random copolymers was correlated to confirm the carbon resonance assignments for all 20 possible pentad comonomer sequences. Using the experimental pentad sequence prevalences, a computer code was developed to simulate chains with microstructures typical of each sample as a means to visually represent the copolymer blockiness with quantitative precision. Based on the microstructure and distribution of run lengths in these chains, the simulations revealed that the Blocky copolymers contain a high degree of blockiness. By comparing the run lengths in the simulated chains to the average number of styrene units in a crystalline segment of sPS (found by small-angle X-ray scattering), copolymer crystallizability was predicted. For the simulated Blocky B-21% (21 mol % Br-Sty) chain, the probability of randomly selecting a styrene unit in a crystallizable block was 25.8%, while that in the simulated Random R-18% was zero, in excellent agreement with the experimental crystallization behavior measured by differential scanning calorimetry. Additionally, these predictions confirmed that the simulated chains accurately represent the ensemble of chains in their respective copolymer samples. Furthermore, each simulated Blocky chain contained one or more long sPS blocks that paralleled the measured 38-40 styrene units spanning a crystalline segment within the sPS/CCl4 gel. This finding affirmed that the long sPS segments originated from the precise lamellar structure within the heterogeneous gel morphology (i.e., block length is correlated with lamellar thickness). Overall, the ability to tailor the copolymer microstructure through control of the semicrystalline gel morphology opens the door to synthesizing ordered copolymers by postpolymerization functionalization processes with unprecedented levels of compositional control.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Identifying microturbulence regimes in a TCV discharge making use of physical constraints on particle and heat fluxes

Reducing the uncertainty on physical input parameters derived from experimental measurements is essential towards improving the reliability of gyrokinetic turbulence simulations. This can be achieved by introducing physical constraints. Amongst them, the zero particle flux condition is considered here. A first attempt is also made to match as well the experimental ion/electron heat flux ratio. This procedure is applied to the analysis of a particular Tokamak à Configuration Variable discharge. A detailed reconstruction of the zero particle flux hyper-surface in the multi-dimensional physical parameter space at fixed time of the discharge is presented, including the effect of carbon as the main impurity. Both collisionless and collisional regimes are considered. Hyper-surface points within the experimental error bars are found. In conclusion, the analysis is done performing gyrokinetic simulations with the local version of the GENE code, computing the fluxes with a Quasi-Linear (QL) model and validating the QL results with non-linear simulations in a subset of cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Spall response of single-crystal copper

We performed a series of systematic spall experiments on single-crystal copper in an effort to determine and isolate the effects of crystal orientation, peak stress, and unloading strain rate on the tensile spall strength. Strain rates ranging from 0.62 to 2.2 × 10 6 s -1 and peak shock stresses in the 5–14 GPa range, with one additional experiment near 50 GPa, were explored as part of this work. Gun-driven impactors, called flyer plates, generated flat top shocks followed by spall. This work highlights the effect of crystal anisotropy on the spall strength by showing that the spall strength decreases in the following order: [100], [110], and [111]. Over the range of stresses and strain rates explored, the spall strength of [100] copper depends strongly on both the strain rate and shock stress. Except at the very highest shock stress, the results for the [100] orientation show linear relationships between the spall strength and both the applied compressive stress and the strain rate. In addition, hydrodynamic computer code simulations of the spall experiments were performed to calculate the relationship between the strain rate near the spall plane in the target and the rate of free surface velocity release during the pullback. As expected, strain rates at the spall plane are much higher than the strain rates estimated from the free surface velocity release rate. We have begun soft recovery experiments and molecular dynamics calculations to understand the unusual recompression observed in the spall signature for [100] crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Fourth-order vibrational perturbation theory with the Watson Hamiltonian: Report of working equations and preliminary results

A derivation of fourth-order vibrational perturbation theory (VPT4) based on the Watson Hamiltonian in dimensionless rectilinear normal coordinates is presented. Terms that are linear and cubic in the (n k + 1/2), with n k being the zeroth-order harmonic oscillator quantum numbers, appear in fourth order and extend the much simpler second-order vibrational perturbation theory model. The rather involved expressions for the fourth-order terms are derived with Rayleigh-Schrödinger perturbation theory, the process of verifying their correctness is described, and a computer code to generate the VPT4 constants from the potential energy surface derivatives is provided. Here, the paper concludes with numerical examples featuring the H 2 O, Si 2 C, and cyclic-C 3 H 2 molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

The role of rotation-vibration coupling in symmetric and asymmetric isotopomers of ozone

A theoretical framework and a computer code (SpectrumSDT) are developed for accurate calculations of coupled rotational–vibrational states in triatomic molecules using hyper-spherical coordinates and taking into account the Coriolis coupling effect. Concise final formulas are derived for the construction of the Hamiltonian matrix using an efficient combination of the variational basis representation and discrete variable representation methods with locally optimized basis sets and grids. First, the new code is tested by comparing its results with those of the APH3D program of Kendrick et al. [Kendrick, Pack, Walker, and Hayes, J. Chem. Phys. 110, 6673 (1999)]. Then, accurate calculations of the rovibrational spectra are carried out for doubly substituted symmetric ( 18 O 16 O 18 O) and asymmetric ( 18 O 18 O 16 O) ozone isotopomers for the total angular momentum up to J = 5. Together with similar data recently reported for the singly substituted symmetric ( 16 O 18 O 16 O) and asymmetric ( 16 O 16 O 18 O) ozone isotopomers, these calculations quantify the role of the Coriolis coupling effect in the large mass-independent isotopic enrichment of ozone, observed in both laboratory experiments and the atmosphere of the Earth. It is found that the Coriolis effect in ozone is relatively small, as evidenced by deviations of its rotational constants from the symmetric-top-rotor behavior, magnitudes of parity splittings (Λ-doubling), and ratios of rovibrational partition functions for asymmetric vs symmetric ozone molecules. It is concluded that all of these characteristics are influenced by the isotopic masses as much as they are influenced by the overall symmetry of the molecule. It is therefore unlikely that the Coriolis coupling effect could be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes in ozone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Yield degradation mechanisms for two-shock capsules evaluated through simulations

An investigation of twenty two-shock campaign indirectly driven capsules on the National Ignition Facility was conducted using the xRAGE computer code. The two-shock platform was developed to look at the sensitivity of fuel–ablator mix with shock timing, asymmetry, surface roughness, and convergence on roughly ignition size scale capsules. This platform used CH/CD (plastic/deuterated plastic) shell capsules that were about 685-μm outer radius and filled with D2 or hydrogen-tritium (HT) gas. The experimental radius and velocity vs time, neutron yield, burn averaged ion temperature (Tion), burn width, and self-emission image size were compared to one-dimensional (1D) and two-dimensional (2D) simulations. Our 2D simulations suggest that the mixing of glass from the fill tube was the dominant source of impurity in the gas region of the capsule during burn, along with fuel–ablator mix. The mass of glass mixed in is about 5–10 ng. Our 2D simulations capture most of the yield trends from different degradation mechanisms, and they match the observed burn width and Tion measurements. Our 2D models match all the available data to within 2.5 times the normalized experimental error for 19 of 20 capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The long road to ignition: An eyewitness account

This paper reviews the many twists and turns in the long journey that culminated in ignition in late 2022 using the laser heated indirect-drive approach to imploding DT filled targets at the National Ignition Facility (NIF), located at the Lawrence Livermore National Laboratory (LLNL). We describe the early origins of the Laser Program at LLNL and key developments such as the paradigm shifting birth of high energy density physics (HEDP) studies with lasers, changes in choice of laser wavelength, and the development of key diagnostics and computer codes. Fulfilling the requirements of the multi-faceted Nova Technical Contract was a necessary condition for the approval of the NIF, but more importantly, the end of the Cold War and the cessation of nuclear testing were key catalysts in that approval, along with the ready-and-waiting field of HEDP. The inherent flexibility of the field of laser driven inertial confinement fusion played a fundamental role in achieving success at the NIF. We describe how the ultimately successful ignition target design evolved from the original “point design” target, through the lessons of experiment. All key aspects of that original design changed: The capsule's materials and size were changed; the hohlraum's materials, size, laser entrance hole size, and gas fills were also all changed, as were the laser pulse shapes that go along with all those changes. The philosophy to globally optimize performance for stability (by raising the adiabat and thus lowering the implosion convergence) was also key, as was progress in target fabrication, and in increasing NIF's energy output. The persistence of the research staff and the steadfast backing of our supporters were also necessary elements in this success. We gratefully acknowledge seven decades of researcher endeavors and four decades of the dedicated efforts of many hundreds of personnel across the globe who have participated in NIF construction, operation, target fabrication, diagnostic, and theoretical advances that have culminated in ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

MELCOR Analysis of OSU Multi-Application Small Light Water Reactor (MASLWR) Experiment

A multi-application small light water reactor (MASLWR) conceptual design was developed by Oregon State University (OSU) with emphasis on passive safety systems. The passive containment safety system employs condensation and natural circulation to achieve the necessary heat removal from the containment in case of postulated accidents. Containment condensation experiments at the MASLWR test facility at OSU are modeled and analyzed with MELCOR, a system-level reactor accident analysis computer code. The analysis assesses its ability to predict condensation heat transfer in the presence of noncondensable gas for accidents where high-energy steam is released into the containment. This work demonstrates MELCOR’s ability to predict the pressure-temperature response of the scaled containment. Our analysis indicates that the heat removal rates are underestimated in the experiment due to the limited locations of the thermocouples and applies corrections to these measurements by conducting integral energy analyses along with CFD simulation for confirmation. Furthermore, the corrected heat removal rate measurements and the MELCOR predictions on the heat removal rate from the containment show good agreement with the experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Dose Rate Evaluation for the ES-3100 Package with HEU Content Using MCNP, ADVANTG, Monaco, and MAVRIC

Here, this paper presents a comparative study of dose rate calculations for the ES-3100 package with highly enriched uranium (HEU) content for different source configurations using computer codes: (1) MCNP; (2) Automated Variance Reduction Generator (ADVANTG)/MCNP; (3) Monaco; and (4) Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC). The Model ES-3100 package was developed at the Y-12 National Security Complex for domestic and international transportation of Type B fissile radioactive material. In this study, six different source configurations (i.e., solid cylinder, cylindrical hemi-shell, cylindrical shell, rectangular plate, cylindrical rod, or cylindrical segment form) having 36 kg of HEU metal inside the package containment vessel (based on configurations in the ES-3100/HEU safety analysis report for packaging) are evaluated. Dose rates at the package surfaces and 1 m from the package surfaces are calculated for these different source configurations. MCNP and Monaco cases are run without any biasing options to accelerate the convergence. Consistent Adjoint Driven Importance Sampling (CADIS) and Forward-Weighted CADIS (FW-CADIS) methods developed at the Oak Ridge National Laboratory are implemented in ADVANTG/MCNP and MAVRIC codes to accelerate the convergence. ADVANTG generates variance reduction parameters using Denovo code, and MCNP is used with the variance reduction parameters to accelerate the convergence. MAVRIC uses Denovo code to construct an importance map and a biased source distribution that are supplied to Monaco to accelerate the Monte Carlo simulation. The FW-CADIS option in ADVANTG and MAVRIC is used to accelerate the convergence in this study. The accelerated convergence cases (ADVANTG/MCNP and MAVRIC) are about 100 times faster with 100 times less particle simulation than those cases run without biasing options (analog MCNP and analog Monaco). Lastly, the MCNP, ADVANTG/MCNP, Monaco, and MAVRIC calculated dose rates at the package surfaces and at 1 m from the package surfaces for the different source configurations are compared and are found to be in general agreement.

42 ENGINEERING↗
DOE OSTI2023

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Neutronics and Thermo-Fluids Simulation of Generic Pebble-Bed Fluoride-Salt-Cooled High-Temperature Reactor

The fluoride-salt-cooled high-temperature reactor (FHR) is one of the advanced reactors that has been attracting considerable interest from both the research community and the nuclear industry. To help facilitate the nuclear community's familiarity with the FHR, Kairos Power has developed a generic FHR (gFHR) benchmark. In the research performed here, this benchmark was used to assess innovative modeling methods that combine stochastic and deterministic computer codes to perform the design and analysis of the gFHR. Further, the Monte Carlo code Serpent 2 was used to generate few-group cross sections that were then used in the neutron diffusion and thermal-fluids code AGREE to perform full-core neutronics and thermal-fluids steady-state and transient core analysis. The Argonne National Laboratory code SAM was then used to model the gFHR system and to simulate the load-follow operation of the gFHR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

SCALE Shielding Calculations for Advanced Reactor Accident Scenarios

The study results presented in this paper demonstrate the capabilities of the SCALE computer code for non–light water reactor (non-LWR) radiation source term and shielding calculations in support of the U.S. Nuclear Regulatory Commission non-LWR fuel cycle demonstration project. Representative non-LWR types, including the sodium-cooled fast reactor (SFR), the molten salt reactor (MSR), and the heat pipe microreactor (HPMR), were analyzed to evaluate dose rates associated with postulated accident scenarios (SFR and MSR), reactor operation (MSR and HPMR), and fuel self-protecting characteristics (HPMR). New features were implemented in SCALE depletion codes to better simulate MSR operation.

MSR↗
DOE OSTI2019

Simulation of carbon dioxide absorption by amino acids in two-phase batch and bubble column reactors

The absorption of carbon dioxide (CO 2 ) is a critical process in many practical applications. The use of amino acid solutions as absorption solvents has the potential to reduce the amount of energy required by the regeneration process. The goal of this research project is to develop dynamic models to simulate CO 2 absorption by using amino acid solutions as absorption solvents. A reaction scheme is introduced to represent the chemical reactions between the amino acid and CO 2 . Two reaction models, for a two-phase batch reactor and a bubble column, based upon transient mass and energy balances for the chemical species found in CO 2 gas-liquid absorption are presented. Computer codes have been written to implement the proposed models. Simulation results are presented and discussed. The proposed models can be used to optimize and control CO 2 absorption in practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Integrated Liquid Metal Flowing First Wall and Open-Surface Divertor for Fusion Nuclear Science Facility: Concept, Design, and Analysis

The Fusion Energy System Studies (FESS) Fusion Nuclear Science Facility (FNSF) project team in the United States is examining the use of liquid metals (LMs) for plasma-facing components (PFCs). Our approach has been to utilize an already established fusion design, FESS-FNSF, which is a tokamak-based machine with 518 MW fusion power, a 4.8-m major radius, a 1.2-m minor radius, and a machine average neutron wall loading of ~1 MW/m 2 . For this design, we propose a PFC concept that integrates a flowing LM first wall (FW) and an open-surface divertor. The flowing LM first removes the surface heat flux from the FW and then proceeds to the lower section of the vacuum chamber to form a large area LM surface for absorbing high peak surface heat flux in the divertor region. In pursuing the application of large open LM surfaces in the FNSF, two new computer codes have been developed and then applied to the analysis of free-surface magnetohydrodynamic flows and heat transfer, including fast thin flowing liquid layers over the solid FW (liquid wall), a tublike divertor, and a fast flow divertor. The analysis is aimed at optimization of the liquid wall design by matching certain proposed design criteria and also at evaluation of the maximum heat fluxes, using liquid lithium (Li) as a working fluid. It was demonstrated that the flowing Li FW (at ~2 cm and ~10 m/s) can tolerate a surface heat flux of ~1 MW/m 2 , while the open-surface Li divertor can remove a maximum high peak heat flux of 10 MW/m 2 . Furthermore, the paper also focuses on the underlying science. One such example is the evaluation and characterization of heat transfer mechanisms and heat transfer intensification in the tublike Li divertor.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗