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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Physics-informed neural networks (PINNs) for fluid mechanics: a review

Despite the significant progress over the last 50 years in simulating flow problems using numerical discretization of the Navier–Stokes equations (NSE), we still cannot incorporate seamlessly noisy data into existing algorithms, mesh-generation is complex, and we cannot tackle high-dimensional problems governed by parametrized NSE. Moreover, solving inverse flow problems is often prohibitively expensive and requires complex and expensive formulations and new computer codes. Here, we review flow physics-informed learning, integrating seamlessly data and mathematical models, and implement them using physics-informed neural networks (PINNs). Here, we demonstrate the effectiveness of PINNs for inverse problems related to three-dimensional wake flows, supersonic flows, and biomedical flows.

42 ENGINEERING↗
DOE OSTI2018

Metropolis–Hastings Monte Carlo Method for Neutron Emissivity Tomographic Inversion in Tokamak Plasma

The article presents a new approach to plasma neutron emissivity reconstruction based on the Metropolis–Hastings Monte Carlo algorithm. The algorithm is based on a biased random walk. A dedicated computer code generates pseudo-random samples within the domain that contains the solution. The properly chosen objective function ensures convergence to the desired solution. The method allows to overcome the problem of trapping in local minima by accepting, with some finite probability, changes in the direction opposite to that indicated by the objective function. The developed method was tested using a synthetic data set. The results show that the method provides reliable reconstructions and thus can be used as an alternative or complementary method to the commonly applied ones. The method capabilities were also tested against noise present in the line-integrated projections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Energy balance in the DSSAT-CSM-CROPGRO model

One potential way to improve crop growth models is for the models to predict energy balance and evapotranspiration (ET) from first principles, thus serving as a check on “engineered” ET methodology. In this paper, we present new implementations and the results of an energy balance model (EBL) developed by Jagtap and Jones (1989) and then implemented in DSSAT’s CROPGRO (CG-EBL) model by Pickering et al. (1995) as a linked energy balance-photosynthesis model that has not been field-tested until now. The energy balance code computes evapotranspiration and other energy balance components, as well as a canopy air temperature, based on three sources (sunlit leaves, shaded leaves, soil surface). Model performance was evaluated with measured biomass and energy fluxes from two sites in Nebraska, namely, the US-Ne2 irrigated maize-soybean rotation field and the US-Ne3 rainfed maize-soybean rotation field, which are part of the Ameriflux eddy covariance network (htt ps://ameriflux.lbl.gov/sites). After implementing new aerodynamic resistances and the stomatal conductance model of the Ball–Berry–Leuning, crop growth, evapotranspiration and soil temperature were simulated well by the EBL model. The EBL improved ET predictions slightly over the often-used FAO56 method [Penman–Monteith (Allen et al., 1998)] for 4 of the 5 years evaluated for both irrigated and rainfed conditions. Further, a significant improvement was achieved using EBL for the simulation of soil temperature at the various depths compared to STEMP, the original subroutine in DSSAT for simulating soil temperature. Compared to the other available DSSAT methods, the EBL explicitly simulates the impacts of crop morphology, physiology and management on the crop’s environment and energy and mass exchange, which in turn directly affect the water use and irrigation requirements, phenology, photosynthesis, growth, sterility, and yield of the crop.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Mass transport analysis for tritium removal in FHRs

Tritium control and mitigation is one of the most significant issues in Fluoride Salt-cooled Hightemperature Reactors (FHRs). To address this issue, a cross-flow tritium removal facility has been proposed to remove molecular tritium, T 2 , from the primary coolant. To model the performance of the design, a method is developed based on the logarithmic mean difference of the square root of the tritium partial pressures on the upstream side and the purging gas side. A computer code has been developed based on this global analysis method to study the effectiveness of tritium removal in such a cross-flow tritium removal facility. The code is benchmarked against experimental data from a hydrogen separation experiment in the literature and compared with the results from a code based on traditional finite volume method. Evaluation of the cross-flow tritium removal facility designed for FHR systems has also been performed. The results show how the key factors, i.e., the fluid flow rates and tritium inlet concentration affect the performance of the facility.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A fast multi-fidelity method with uncertainty quantification for complex data correlations: Application to vortex-induced vibrations of marine risers

Here we develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.

42 ENGINEERING↗
DOE OSTI2023

Predicting displacement damage for ion irradiation: Origin of the overestimation of vacancy production in SRIM full-cascade calculations

Ion irradiation and implantation have wide applications that demand accurate determination of displacement damage profile and distribution of implanted ion concentration. The prediction of vacancies is especially important to determine displacements per atom (dpa), the standard parameter of primary radiation damage in materials. However, significant discrepancies exist in estimations of vacancies between full-cascade (F-C) and quick calculation (Q-C) options in the popular computer code SRIM. This study inspected the SRIM code and a relatively new code called Iradina, which uses a similar methodology, to develop an understanding of the origin of vacancy overestimation in the F-C options for SRIM and Iradina. Here, we found that the default values of thresholds (namely final energy in SRIM and replacement energy in Iradina) in displacement production calculations results in excessively large number of calculated vacancies and very few replacements. After conducting multiple calculations using SRIM, Iradina, and MARLOWE (all based on the binary collision approximation), a comparison of the results indicates that there is a shortcoming in the SRIM and Iradina F-C methodology for treating near-threshold collisions. This issue is responsible for the deficiency of replacements and excess of calculated vacancies in the SRIM and Iradina F-C results. Drawing on the principles of collision physics, we propose recommendations for modifying the source codes to address these issues.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Algorithms and algebraic solutions of decay chain differential equations for stable and unstable nuclide fractionation

Radioactive decay processes, such as alpha decay, produce decay chains where the mass numbers of nuclides decrease as larger nuclides expel energetic particles to form smaller nuclides. Under these conditions, the coefficient matrix that describes the differential rate expressions for radioactive decay can be made lower triangular. With this special structure, formulating an algebraic solution to the decay chains can be done by first formulating the eigenvectors that make up the coefficient matrix, which can then be solved using forward substitution for a lower triangular matrix. This work details the derivation of algebraic solutions for decay chains of any number of stable and unstable nuclides with any number of branching based on this eigenvector analysis. A prototype computational code was developed to validate and compare this methodology against a number of other methods for solving similar systems. A two-phase sorting algorithm yielding the lower triangular matrix structure was established to apply the developed algebraic solutions for decay chains involving beta-emitting radionuclides transformed into daughter nuclides without change in their mass number. In conclusion, the methodologies produced in this work provide an efficient way to estimate nuclide fractions from natural decay processes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

2p x-ray absorption spectroscopy of 3d transition metal systems

This review provides an overview of the different methods and computer codes that are used to interpret 2p x-ray absorption spectra of 3d transition metal ions. We first introduce the basic parameters and give an overview of the methods used. We start with the semi-empirical multiplet codes and compare the different codes that are available. A special chapter is devoted to the user friendly interfaces that have been written on the basis of these codes. Next we discuss the first principle codes based on band structure, including a chapter on Density Functional theory based approaches. We also give an overview of the first-principle multiplet codes that start from a cluster calculation and we discuss the wavefunction based methods, including multi-reference methods. We end the review with a discussion of the link between theory and experiment and discuss the open issues in the spectral analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerical simulation of coupled THM behaviour of full-scale EBS in backfilled experimental gallery in the Horonobe URL

Bentonite-based engineered barrier system (EBS) is a key component of many repository designs for the geological disposal of high-level radioactive waste. Given the complexity and interaction of the phenomena affecting the barrier system, coupled thermo-hydro-mechanical (THM) numerical analyses are a potentially useful tool for a better understanding of their behaviour. In this context, a Task (the Horonobe EBS experiment) was undertaken to study, using numerical analyses, the thermo-hydro-mechanical (and thermo-hydro) interactions in bentonite based engineered barriers within the international cooperative project DECOVALEX 2023. One full-scale in-situ experiment and four laboratory experiments, largely complementary, were selected for modelling. The Horonobe EBS experiment is a temperature-controlled non-isothermal experiment combined with artificial groundwater injection. The Horonobe EBS experiment consists of the heating and cooling phases. Six research teams performed the THM or TH (depended on research team approach) numerical analyses using a variety of computer codes, formulations and constitutive laws. Finally, for each experiment, the basic features of the analyses are described and the comparison between calculations and laboratory experiments and field observations are presented and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A constitutive model incorporating grain refinement strengthening on metallic alloys

Surface nano-crystallization techniques have been recently developed as one of the most efficient ways to optimize materials’ structure, and therefore develop the local and global mechanical behavior as to increase strength without compromising ductility. Here, in this work, we present a constitutive model incorporating grain refinement hardening to simulate the nano-crystallization technique, specifically, surface mechanical attrition treatment. The computation is implemented using user-defined VUMAT subroutines. As an example of its application, a geometry model with full coverage of random impacts are employed. The results show that the model has rather precise predictability of grain size evolution during plastic deformation. The readily embedded with a computational code of material dynamics enables this novel model to be a promising tool to study the dynamic evolution of microstructures under plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nuclear Data Sheets for A=203

Evaluated nuclear structure and decay data for all nuclei with mass number A=203 ( 203 Os, 203 Ir, 203 Pt, 203 Au, 203 Hg, 203 Tl, 203 Pb, 203 Bi, 203 Po, 203 At, 203 Rn, 203 Fr, 203 Ra) are presented. All available experimental data are compiled and evaluated, and best values for level and γ-ray energies, quantum numbers, lifetimes, γ-ray intensities and transition probabilities, as well as other nuclear properties, are recommended. Inconsistencies and discrepancies that exist in the literature are discussed. A number of computer codes (https://www-nds.iaea.org/public/ensdf_pgm/index.htm) developed by members of the NSDD network were used during the evaluation process. For example, the reported absolute γ-ray emission probabilities and their uncertainties were determined using the GABS code. Here, the γ-ray transition probabilities were determined using the RULER code and the corresponding uncertainties were estimated by means of a Monte-Carlo approach. This work supersedes the earlier evaluation by F.G. Kondev (2005Ko20), published in Nuclear Data Sheets 105, 1 (2005).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Nuclear Data Sheets for A=201

Evaluated nuclear structure and decay data for all nuclei with mass number A=201 ( 201 Os, 201 Ir, 201 Pt, 201 Au, 201 Hg, 201 Tl, 201 Pb, 201 Bi, 201 Po, 201 At, 201 Rn, 201 Fr, 201 Ra) are presented. All available experimental data are compiled and evaluated, and best values for level and γ-ray energies, quantum numbers, lifetimes, γ-ray intensities and transition probabilities, as well as other nuclear properties, are recommended. Inconsistencies and discrepancies that exist in the literature are discussed. A number of computer codes (https://wwwnds.iaea.org/public/ensdf pgm/index.htm) developed by members of the NSDD network were used during the evaluation process. For example, the reported absolute γ-ray emission probabilities and their uncertainties in various decay data sets were determined using the GABS code. The γ-ray transition probabilities were determined using the RULER code and the corresponding uncertainties were determined using a Monte-Carlo approach. In conclusion, this work supersedes the earlier evaluation by F.G. Kondev (2007Ko06), published in Nuclear Data Sheets 108, 365 (2007).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Recommended Nuclear Structure and Decay Data for A = 200 Isobars

Evaluated nuclear structure and decay data for all nuclei with mass number A = 200 ( 200 Os, 200 Ir, 200 Pt, 200 Au, 200 Hg, 200 Tl, 200 Pb, 200 Bi, 200 Po, 200 At, 200 Rn, and 200 Fr), are presented. All available experimental data are compiled and evaluated, and best values for level and gamma-ray energies, quantum numbers, lifetimes, gamma-ray intensities and transition probabilities, as well as other nuclear properties, are recommended. Inconsistencies and discrepancies that exist in the literature are discussed. A number of computer codes (https://www-nds.iaea.org/public/ensdf pgm/) developed by members of the NSDD network were used during the evaluation process. In conclusion, this work supersedes the earlier evaluation by F.G. Kondev and S. Lalkovski (2007Ko42), published in Nuclear Data Sheets 108, 1471 (2007).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Recommended Nuclear Structure and Decay Data for A=206 Isobars

Here, evaluated nuclear structure and decay data for all nuclei with mass number A=206 ( 206 Pt, 206 Au, 206 Hg, 206 Tl, 206 Pb, 206 Bi, 206 Po, 206 At, 206 Rn, 206 Fr, 206 Ra and 206 Ac), are presented. All available experimental data are compiled and evaluated, and best values for level and γ-ray energies, quantum numbers, lifetimes, γ-ray intensities and transition probabilities, as well as other nuclear properties, are recommended. Inconsistencies and discrepancies that exist in the literature are discussed. A number of computer codes (https://www-nds.iaea.org/public/ensdf_pgm/) developed by members of the NSDD network were used during the evaluation process. This work supersedes the earlier evaluation by F.G. Kondev (2008Ko21), published in Nuclear Data Sheets 109, 1527 (2008).

Kondev, F. G. [Argonne National Laboratory (ANL), ↗
DOE OSTI2024

Scaling methodologies and similarity analysis for thermal hydraulics test facility development for water-cooled small modular reactor

Small modular reactors (SMRs) represent a promising option for providing clean and sustainable energy due to their potential for enhanced safety, reduced capital costs, and increased siting flexibility. However, new reactor systems require the development and operation of representative scaled-down test facilities to support the verification and validation of system computer codes and models. Here this study reviews the research on scaling methodologies and similarity principles pivotal in developing non-nuclear integral effects test and separate effects test facilities for water-cooled SMRs. The study focuses on a review of the scaling methods, similarity approaches, and possible challenges posed by the unique and compact design features of integral-pressurized water reactor-type SMRs, and their representative test facilities. This study also reviews previous research related to scaling and similarity methodologies and provides insights into design considerations for achieving prototypic conditions in test facilities. The findings and recommendations emphasize the broader impact of appropriate scaling and similarity principles to ensure meaningful and transferable results from non-nuclear test facilities to accelerate the safe and efficient deployment of next-generation water-cooled SMRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A seamless multiscale operator neural network for inferring bubble dynamics

Modelling multiscale systems from nanoscale to macroscale requires the use of atomistic and continuum methods and, correspondingly, different computer codes. Here, we develop a seamless method based on DeepONet, which is a composite deep neural network (a branch and a trunk network) for regressing operators. In particular, we consider bubble growth dynamics, and we model tiny bubbles of initial size from 100 nm to 10 $\mathrm {\mu }\textrm {m}$ , modelled by the Rayleigh–Plesset equation in the continuum regime above 1 $\mathrm {\mu }\textrm {m}$ and the dissipative particle dynamics method for bubbles below 1 $\mathrm {\mu }\textrm {m}$ in the atomistic regime. After an offline training based on data from both regimes, DeepONet can make accurate predictions of bubble growth on-the-fly (within a fraction of a second) across four orders of magnitude difference in spatial scales and two orders of magnitude in temporal scales. The framework of DeepONet is general and can be used for unifying physical models of different scales in diverse multiscale applications.

Mechanics↗
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↗