Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer code”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
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

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 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 OSTI2020

Correlation matrix renormalization theory in multi-band lattice systems

An appropriate treatment of electronic correlation effects plays an important role in accurate descriptions of physical and chemical properties of real materials. The recently proposed correlation matrix renormalization theory with sum rule correction (CMR) for studying correlated-electron materials has shown good performance in molecular systems and a periodic hydrogen chain in comparison with various quantum chemistry and quantum Monte Carlo calculations. Additionally, this work gives a detailed formulation and computational code implementation of CMR in multi-band periodic lattice systems. This lattice CMR ab initio theory is highly efficient, has no material specific adjustable parameters, and has no double counting issues faced by the hybrid approaches like LDA + U, DFT + DMFT and DFT + GA type theories. Benchmark studies on materials with s and p orbitals in this study show that CMR in its current implementation consistently performs well for these systems as the electron correlation increases from the bonding region to the bond breaking region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Calculation of RF sheath properties from surface wave-fields: a post-processing method

In ion cyclotron range of frequency (ICRF) experiments in fusion research devices, radio frequency (RF) sheaths form where plasma, strong RF wave fields and material surfaces coexist. These RF sheaths affect plasma material interactions such as sputtering and localized power deposition, as well as the global RF wave fields themselves. RF sheaths may be modeled by employing a sheath boundary condition (BC) in place of the more customary conducting wall BC; however, there are still many ICRF computer codes that do not implement the sheath BC. In this paper we present a method for post-processing results obtained with the conducting wall BC. The post-processing method produces results that are equivalent to those that would have been obtained with the RF sheath BC, under certain assumptions. We found the post-processing method is also useful for verification of sheath BC implementations and as a guide to interpretation and understanding of the role of RF sheaths and their interactions with the waves that drive them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Practical overview of the ICRP Publication 145 phantoms for use with MCNP

Abstract This paper provides a practical overview of the mesh-type reference phantoms (MRCPs) published by the ICRP by providing supplementary guidance on the proper implementation of using the MRCPs in computational codes, such as Monte Carlo N-Particle. Various validation scenarios for photon and neutron exposures were also examined where published values were reproduced using the MRCPs for both external and internal exposures. Generally, the MRCPs performed well in the validation scenarios with the values reproduced being within 10% of the published values.

Hetrick, Lucas (ORCID:0000000178978007)↗
DOE OSTI2025

Competing ionization and dissociation: Extension of the energy-dependent frame transformation to the gerade symmetry of H 2

This article solves two major tasks that frequently arise in the theory of electron collisions with a target molecular cation. First, it extends the energy-dependent frame transformation (EDFT) treatment, which is needed to map fixed-nuclei electron-molecule scattering matrices into an energy-dependent laboratory-frame scattering matrix with vibrational channel indices. The EDFT mapping can now be carried out even when the target molecule possesses multiple low-energy potential curves, significantly transcending previous applications. Second, it implements a method to extract the rest of the full laboratory-frame scattering matrix, i.e., the columns and rows describing input and/or output dissociation channels. The treatment is benchmarked in this article against the essentially exact solution of a refined two-dimensional model of the singlet gerade Σ symmetry of H 2 . Our tests demonstrate that the theory accurately maps fixed-nuclei scattering information, of the type provided by existing electron-molecule computer codes, into a laboratory-frame scattering matrix that includes both ionization and dissociation. Furthermore, this treatment can provide a general framework applicable to a broad class of electron collision processes involving diatomic target ions, suitable for an accurate description of challenging processes such as dissociative recombination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Neutron-neutron scattering length from the 6 He(p,pα)nn reaction

We propose a novel method to measure the neutron-neutron scattering length using the 6 He(p,pα)nn reaction in inverse kinematics at high energies. The method is based on the final-state interaction (FSI) between the neutrons after the sudden knockout of the α particle. We show that the details of the neutron-neutron relative-energy distribution allow for a precise extraction of the s-wave scattering length. Here, we present the state of the art in regard to the theory of this distribution. The distribution is calculated in two steps. First, we calculate the ground-state wave function of 6 He as a αnn three-body system. For this purpose we use Halo effective field theory, which also provides uncertainty estimates for the results. We compare our results at this stage to model calculations done with the computer code face. In a second step we determine the effects of the nn FSI using the nn t-matrix. We compare these FSI results to approximate FSI approaches based on standard FSI enhancement factors. While the final distribution is sensitive to the nn scattering length, it depends only weakly on the effective range. Throughout we emphasize the impact of theoretical uncertainties on the neutron-neutron relative-energy distribution, and discuss the extent to which those uncertainties limit the extraction of the neutron-neutron scattering length from the reaction 6 He(p,pα)nn.

6 ≤ A ≤ 19↗
DOE OSTI2024

Supersymmetry with scalar sequestering

Supersymmetric models with a strongly interacting superconformal hidden sector (HS) may drive soft supersymmetry (SUSY) breaking scalar masses, bilinear soft term B μ and Higgs combinations m H u , d 2 + μ 2 to small values at some intermediate scale, leading to unique sparticle mass spectra along with possibly diminished fine-tuning in spite of a large superpotential μ parameter. We set up a computer code to calculate such spectra, which are then susceptible to a variety of constraints: (1) possible charge-or-color breaking (CCB) minima in the scalar potential, (2) unbounded from below (UFB) scalar potential, (3) improper electroweak symmetry breaking, (4) a charged or sneutrino lightest SUSY particle (LSP), (5) generating m h ∼ 125 GeV , (6) consistency with LHC sparticle mass limits, and (7) naturalness. We find this bevy of constraints leaves little or no viable parameter space for the case where hidden sector dynamics dominates minimal supersymmetric standard model (MSSM) running, even for the case of nonuniversal gaugino masses. For the case with moderate HS running with comparable MSSM running, and with universal gaugino masses, then the fine-tuning is ameliorated, but nonetheless remains high. Viable spectra with moderate HS running and with low fine-tuning and large μ can be found for nonuniversal gaugino masses. Published by the American Physical Society 2024

Astronomy & Astrophysics↗