Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “F 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

P2M Simulation Exercise on Past Fuel Melting Irradiation Experiments

This paper presents the results of the Power To Melt and Maneuverability (P2M) Simulation Exercise on past fuel melting irradiation experiments, organized within the Organisation for Economic Co-operation and Development/Nuclear Energy Agency Framework for IrraDiation ExperimentS (FIDES) framework by the Core Group (CEA, EDF, and SCK·CEN) and open to all FIDES members. The exercise consisted in simulating two past power ramps where fuel melting was detected: (1) the xM3 staircase power transient [ramp terminal level (RTL) 70 kW·m -1 , average burnup 27 GWd·tU -1 ], carried out in 2005 in the R2 reactor at Studsvik (Sweden), where the rodlet maintained its integrity, and (2) the HBC4 fast power transient (RTL 66 kW·m -1 , average burnup 48 GWd·tU -1 ), carried out in 1987 in the BR2 reactor at SCK·CEN (Belgium), where the cladding failed during the experiment. The exercise was joined by 13 organizations from 9 countries using 11 different fuel performance codes. In this paper, the main results of the Simulation Exercise are presented and compared to available postirradiation examinations (PIE) or on-line measurements during the power ramps (fuel and clad diameters, rod elongation, pellet-clad gap, and fission gas release). Since the focus of the Simulation Exercise is on fuel melting assessment, determination of the boundary between melted/nonmelted fuel and the consequent definition of a melting radius from PIE are first discussed. During the HBC4 ramp, fuel melting was predicted by most of the codes despite differences in the melting models. Higher discrepancies were observed for the xM3 rod that can be attributed partly to power uncertainty and partly to the limited capability of the models to describe partial melting of the fuel during this ramp. Finally, possible code developments to improve simulation results are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Sensitivity analysis for characterizing the impact of HNGD model on the prediction of hydrogen redistribution in Zircaloy cladding using BISON code

Hydrogen in zirconium cladding is able to precipitate into zirconium hydrides which impacts cladding integrity. The Hydride Nucleation-Growth-Dissolution (HNGD) model in the BISON code accounts for the precipitation and dissolution kinetics of hydride in Zircaloy material. This paper presents global sensitivity analyses of the HNGD model aiming to enhance our understanding of the hydride precipitation phenomena by quantifying the variance that key parameters have on the prediction of hydrogen behavior under various environmental conditions. Here, model predictions are compared to experimental data obtained under two different conditions: 1) with uniformly precharged specimens subjected to a linear thermal gradient, and 2) specimens precharged with a cathodically applied hydride rim at one end of the sample and subjected to an asymmetric thermal gradient. The Sobol sensitivity analysis identifies the key parameters in the HNGD model for both types of specimens. For linear temperature cases, the heat of transport dominates the accuracy of predictions when no precipitation occurs at the cold end, while Terminal Solid Solubility for Dissolution (TSSD) is the most important parameter when precipitation occurs. A large variation in the predicted hydrogen concentration profiles is found in the range of high TSSD due to the occurrence of precipitation. For asymmetric temperature cases, the solubility coefficient gives the largest impact on the predicted hydrogen distribution, as it determines the amount of solute hydrogen dissolved from the initially applied hydride rim. A large discrepancy in hydrogen distribution between simulations and experiments exists with the asymmetric specimens because BISON simulations fail to predict the precipitation of hydride at the cooler end. Comparative studies using former and updated models verifies the significant impact of the hydride growth mechanism on predicted hydrogen concentration profiles. In particular, when hydride initially exists, changes in TSSD generate a large variation in the predicted amount of precipitation by hydride growth, giving large uncertainty in predicting the hydrogen distribution over the sample length. The outputs characterize the significant impact of the hydride growth mechanism in the HNGD model on predicting hydrogen behavior, and improve the understanding of the precipitation of hydride in Zircaloy cladding within a range of expected environmental conditions. The analyses indicate work is still needed to improve the hydride solvus models in the BISON code to accurately predict experimentally observed hydride concentrations and distributions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

PARFUME/BISON Fission Product Release Predictions versus AGR-3/4 Heating Test Measurements

The fuel performance modeling codes PARFUME (PARticle FUel ModEl) and BISON were used to predict the release of fission products silver, cesium, and strontium from as-irradiated fuel compacts containing tristructural isotropic (TRISO) coated particles during heating tests post irradiation. The AGR-3/4 fuel compacts were irradiated as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) program. The heating tests were conducted at temperatures between 1200°C and 1700°C to simulate reactor accident conditions. The measured fission product release fractions from the heating tests were compared to modeling predictions calculated by PARFUME and BISON to evaluate how the codes compare to experimental results. Comparisons between the experimental measured fission product release fractions from silver, cesium and strontium indicate that both modeling codes overpredict the fission product release fractions demonstrating that the diffusivities used in the codes are overestimated. This results in a conservative estimate predicted by the codes when evaluating the fission product release relative to experimental data as it pertains to silver, cesium, and strontium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Shielding development for the spallation neutron source VENUS instrument

VENUS is a neutron imaging instrument that will use a broad range of neutron wavelengths, from epithermal to cold energies, and will include enhanced contrast mechanisms. It will offer novel energy-selective imaging techniques directly connecting complex engineering materials and systems’ structures, properties, and functions to reveal practical and fundamental answers about their real-world performance. Here, the instrument will be built at SNS beamline 10, facing the decoupled poisoned hydrogen moderator. The driving cost for the instrument is the beamline and instrument cave shielding. Final analyses were performed to evaluate the thickness and composition of shielding materials for the instrument cave and beamline to meet radiation safety criteria for the instrument to start up in 2024 after completing the SNS proton power upgrade.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

DESC: A stellarator equilibrium solver

In this paper, the new code DESC is presented to solve for fixed-boundary ideal magnetohydrodynamic equilibria in stellarators. The approach directly solves the equilibrium force balance as a system of nonlinear equations in the form f(x) = 0. The independent variables x represent nested magnetic flux surfaces expressed in the inverse representation with toroidal flux coordinates, and the equations f(x) quantify equilibrium force balance errors at discrete points in real space. Discretizing with global Fourier–Zernike basis functions properly treats the magnetic axis and minimizes the number of coefficients needed to describe the flux surfaces. The pseudospectral method provides great flexibility in where the errors are evaluated, and the system of equations is efficiently solved with a Newton–Raphson iteration. Equilibria are computed and compared against VMEC for both axisymmetric and non-axisymmetric examples. The results show fast convergence rates and solutions with low errors throughout the plasma volume.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Widespread misregulation of inter-species hybrid transcriptomes due to sex-specific and sex-chromosome regulatory evolution

When gene regulatory networks diverge between species, their dysfunctional expression in inter-species hybrid individuals can create genetic incompatibilities that generate the developmental defects responsible for intrinsic post-zygotic reproductive isolation. Both cis- and trans-acting regulatory divergence can be hastened by directional selection through adaptation, sexual selection, and inter-sexual conflict, in addition to cryptic evolution under stabilizing selection. Dysfunctional sex-biased gene expression, in particular, may provide an important source of sexually-dimorphic genetic incompatibilities. Here, we characterize and compare male and female/hermaphrodite transcriptome profiles for sibling nematode species Caenorhabditis briggsae and C. nigoni, along with allele-specific expression in their F 1 hybrids, to deconvolve features of expression divergence and regulatory dysfunction. Despite evidence of widespread stabilizing selection on gene expression, misexpression of sex-biased genes pervades F 1 hybrids of both sexes. This finding implicates greater fragility of male genetic networks to produce dysfunctional organismal phenotypes. Spermatogenesis genes are especially prone to high divergence in both expression and coding sequences, consistent with a “faster male” model for Haldane’s rule and elevated sterility of hybrid males. Moreover, underdominant expression pervades male-biased genes compared to female-biased and sex-neutral genes and an excess of cis-trans compensatory regulatory divergence for X-linked genes underscores a “large-X effect” for hybrid male expression dysfunction. Extensive regulatory divergence in sex determination pathway genes likely contributes to demasculinization of XX hybrids. The evolution of genetic incompatibilities due to regulatory versus coding sequence divergence, however, are expected to arise in an uncorrelated fashion. This study identifies important differences between the sexes in how regulatory networks diverge to contribute to sex-biases in how genetic incompatibilities manifest during the speciation process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Analytical formula for pellet fuel source density in toroidal plasma configurations based on an areal deposition model

An analytical formula for the fuel source density η F (Ψ) over a magnetic flux surface Ψ following pellet injection is derived. During pellet motion the ablated and ionized fuel will be initially deposited nearby the pellet, prior to subsequent homogenization within the magnetic flux surfaces. Here in our research, it is assumed that the deposition area projected normal to the magnetic field is on the order of the cross-sectional area of the ionized ablation column whose axis is aligned with the magnetic field. Previous pellet deposition codes [1] employed a point source deposition approximation, which overpredicts the surface-averaged density, particularly when the obliquity between the pellet trajectory and the flux surface normal approaches 90 degrees. A key finding of this study is that as the obliquity approaches 90 degrees, grazing trajectory, the areal model predicts a well behaved η F and removes the unphysical divergence stemming from the point model. To illustrate the difference exhibited between areal and point models for near grazing trajectories, we considered vertical pellet injection into a CFETR-like tokamak with a negative triangularity plasma cross section using a pellet ablation/deposition module.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

The cis-regulatory codes of response to combined heat and drought stress in Arabidopsis thaliana

Abstract Plants respond to their environment by dynamically modulating gene expression. A powerful approach for understanding how these responses are regulated is to integrate information about cis-regulatory elements (CREs) into models called cis-regulatory codes. Transcriptional response to combined stress is typically not the sum of the responses to the individual stresses. However, cis-regulatory codes underlying combined stress response have not been established. Here we modeled transcriptional response to single and combined heat and drought stress in Arabidopsis thaliana. We grouped genes by their pattern of response (independent, antagonistic and synergistic) and trained machine learning models to predict their response using putative CREs (pCREs) as features (median F-measure = 0.64). We then developed a deep learning approach to integrate additional omics information (sequence conservation, chromatin accessibility and histone modification) into our models, improving performance by 6.2%. While pCREs important for predicting independent and antagonistic responses tended to resemble binding motifs of transcription factors associated with heat and/or drought stress, important synergistic pCREs resembled binding motifs of transcription factors not known to be associated with stress. These findings demonstrate how in silico approaches can improve our understanding of the complex codes regulating response to combined stress and help us identify prime targets for future characterization.

Azodi, Christina B.↗
DOE OSTI2020

Electron transport in gaseous detectors with a Python-based Monte Carlo simulation code

Understanding electron drift and diffusion in gases and gas mixtures is a topic of central importance for the development of modern particle detection instrumentation. The industry-standard MagBoltz code has become an invaluable tool during its 20 years of development, providing capability to solve for electron transport (‘swarm’) properties based on a growing encyclopedia of built-in collision cross sections. We have made a refactorization of this code from FORTRAN into Cython, and studied a range of gas mixtures of interest in high energy and nuclear physics. The results from the new open source PyBoltz package match the outputs from the original MagBoltz code, with comparable simulation speed. An extension to the capabilities of the original code is demonstrated, in implementation of a new Modified Effective Range Theory interface. We hope that the versatility afforded by the new Python code-base will encourage continued use and development of the MagBoltz tools by the particle physics community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Charge breeders: Development of diagnostic tools to probe the underlying physics

Charge breeders were developed more than 20 years ago in the context of radioactive ion beam (RIB) production. The main goal is to boost the charge state of a singly charged RIB to a higher value matching the A/Q of a post-accelerator. In that way, the RIB produced at some tens of keV can efficiently be accelerated to energies in a range of few MeV/u up to several tens of MeV/u, which is of interest for nuclear structure and nuclear astrophysics experiments. Two families of charge breeders have been developed: one based on Electron Cyclotron Resonance Ion Sources (ECRISs) suitable for Continuous Wave (CW) operation and another built around Electron Beam Ion Sources (EBISs) used for pulsed post-accelerators. For many years, experimental studies have been carried out to enhance the charge breeding process and improve the extracted beam quality, i.e., purity, intensity, emittance, and time structure. The evolution of the charge breeders is also closely related to the emergence of new facilities delivering even more exotic beams. Diagnostic tools, such as Faraday cups, beam profile monitors, and emittance scanners, are mandatory to tune the machines and evaluate their performances. However, to go beyond that, a better understanding of plasma physics (in ECRIS) and electron-ion interactions (in EBIS), as well as the 1+ beam characteristics, for instance, energy spread transverse emittance, is needed. Therefore, simulation codes modeling the physics phenomena in both type of charge breeders were developed to support the understanding of experimental outputs and to gain insights into non-observable parameters. Nevertheless, more sophisticated diagnostic tools are required to improve the charge breeding performance and to extend the potential applications toward new RIB facilities. This Review will present the two kinds of charge breeders and their technical evolutions. It will review the diagnostic tools and simulation codes employed for operation and evaluation of the internal physics processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Modeling of advanced accelerator concepts

Computer modeling is essential to research on Advanced Accelerator Concepts (AAC), as well as to their design and operation. This paper summarizes the current status and future needs of AAC systems and reports on several key aspects of (i) high-performance computing (including performance, portability, scalability, advanced algorithms, scalable I/Os and In-Situ analysis), (ii) the benefits of ecosystems with integrated workflows based on standardized input and output and with integrated frameworks developed as a community, and (iii) sustainability and reliability (including code robustness and usability).

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Opacities of singly and doubly ionized neodymium and uranium for kilonova emission modeling

ABSTRACT Even though the electromagnetic counterpart AT2017gfo to the binary neutron star merger GW170817 is powered by the radioactive decay of r-process nuclei, only few tentative identifications of light r-process elements have been made so far. One of the major limitations for the identification of heavy nuclei is incomplete or missing atomic data. While substantial progress has been made on lanthanide atomic data over the last few years, for actinides there has been less emphasis, with the first complete set of opacity data only recently published. We perform atomic structure calculations of neodymium (Z = 60) as well as the corresponding actinide uranium (Z = 92). Using two different codes [flexible atomic code (fac) and hartree–fock-relativistic (hfr)] for the calculation of the atomic data, we investigate the accuracy of the calculated data (energy levels and electric dipole transitions) and their effect on kilonova opacities. For the fac calculations, we optimize the local central potential and the number of included configurations and use a dedicated calibration technique to improve the agreement between theoretical and available experimental atomic energy levels (AELs). For ions with vast amounts of experimental data available, the presented opacities agree quite well with previous estimations. On the other hand, the optimization and calibration method cannot be used for ions with only few available AELs. For these cases, where no experimental nor benchmarked calculations are available, a large spread in the opacities estimated from the atomic data obtained with the various atomic structure codes is observed. We find that the opacity of uranium is almost double the neodymium opacity.

Flörs, A. (ORCID:0000000320242819)↗
DOE OSTI2022

PETSc/TAO Users Manual: Revision 3.18

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.20)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

PETSc/TAO Users Manual Revision 3.22

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations (PDEs) and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for implementing large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Propagation of radio frequency waves through turbulent plasmas

The practical and economic viability of tokamak fusion reactors depends, in a significant way, on the efficiency of radio frequency (RF) waves to deliver energy and momentum to the plasma in the core of the reactor. The RF electromagnetic waves, excited by antenna structures placed near the wall of a tokamak, have to propagate through the turbulent edge plasma along their path to the core of the fusion device. In present day experiments, the radial width of the edge region and scrape-off layer is of the order of a few centimeters. In ITER, and in future fusion reactors, this width will be of the order of tens of centimeters. Any effects on RF waves due to plasma turbulence have to be properly understood in order to optimize the delivery of RF energy and momentum into the core. This paper is on a multi-pronged, theoretical and computational, approach that is being pursued to quantify the effect of edge plasma turbulence on the propagation of RF waves. The theoretical and analytical models are based on solutions of the Faraday-Ampere equation in a magnetized plasma and on the Kirchhoff tangent plane approximation. An effective medium approach has been developed so as to approximate the permittivity of a turbulent plasma by analytical expressions. The computations are being carried out with a newly developed code ScaRF that is based on the finite difference finite domain technique for solving Maxwell's equations. The plasma permittivity can be assigned as desired. The code is being used to validate the analytical and theoretical models and to evaluate their limitations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗