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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2024

Estimates of global recycling coefficients for LTX- β discharges

We report the first observation of global recycling coefficient R near 0.5 in the Lithium Tokamak eXperiment-β (LTX-β), significantly below the minimum R previously reported in other devices. In a series of experiments with varied Li wall conditioning, estimates of the recycling coefficient have been made using a Lyman-α array and DEGAS2 modeling. A progressive reduction in Lyman-α emission with increased lithium and an increase in edge electron temperature are observed. It is also observed that with increasing Li coating thickness, the effective particle confinement time τp* is reduced and approaches TRANSP calculated energy confinement time (τE), with τp* near τE,TRANSP for the lowest recycling coefficients. Edge temperatures approaching core plasma temperatures, first reported in LTX, can now be directly connected to estimates of the recycling coefficient and qualitatively agree with previous UEDGE simulations. The particle flux to the limiting surfaces appears to be significantly reduced in comparison with fluid scrape-off layer (SOL) models, indicating that a large fraction of the SOL ions are mirror trapped. SOL collisionality drops more than an order of magnitude below the banana regime boundary, indicating the importance of kinetic effects. Full-f 1x2v gyrokinetic simulations of SOL field lines with the GKEYLL code indicate that the fraction of ions trapped along field lines increases as collisionality drops, as a result of increased lithium evaporation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Classification of River Catchments in the Contiguous United States: Code, Dataset, Similarity Patterns, and Resulting Classes

This dataset serves as supplementary information for the paper by Ciulla F. and Varadharajan C. A Network Approach for Multiscale Catchment Classification using Traits (see reference 1). It contains environmental and physical catchment traits, such as temperatures, precipitation, land use and human interference, from 9067 sites across the contiguous United States (CONUS). The purpose of this dataset is to provide information for a better trait-based categorization of river catchments in the CONUS using networks as an analytical tool. The traits variables match the ones present in the GAGES-II dataset and the preprocessing steps are described in the Methods section (processed_dataset.csv). Additionally we include the topologies (nodes, edges and clusters, also referred as classes) of the catchment network and traits network generated by said dataset (csv and json files). A series of tables support the information carried by the network providing more detailed descriptions of cluster components (SI1.pdf). A summary of all the plots of clusters of catchments with at least 50 nodes is provided (SI2.pdf). The characteristic traits for each cluster of catchments is presented as z-score (traits_categories_zscores_per_catchment_class.csv). The link to the hydrological behavior of clusters of catchments is displayed by boxplots, each describing a particular river discharge index (SI3.pdf). Both csv and json files can be read by common text editors but the data contained into them can be better handled using programming languages like python and database oriented libraries like pandas. Pdf files can be read by any pdf reader software.[02-23-2024] Update: The code and datasets necessary to reproduce the results of the study are available as a zipped repository (code_datasets_catchments_similarity.zip).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modelling of electrostatic ion-scale turbulence in divertor tokamaks with the gyrokinetic code $\mathrm{COGENT}$

Here, continuum gyrokinetic simulations of electrostatic ion scale turbulence are presented for the case of a diverted (single-null) tokamak geometry. The simulation model, implemented in the finite-volume code COGENT, solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to a vorticity equation for electrostatic potential variations, where a fluid model is used for an electron response. The model describes the ion scale ion temperature gradient (ITG) and resistive drift modes as well as neoclassical ion physics effects. Different turbulence regimes are observed depending on the plasma profiles, and the roles of a self-consistent background electric field and an X-point geometry are explored. In particular, increasing the pedestal density gradient and the corresponding radial electric field is demonstrated to suppress the ITG turbulence, whereas the same edge plasma background can still be destabilized by the resistive modes when the plasma resistivity is increased. The effects of X-point geometry are assessed by comparing cross-separatrix simulations with counterpart calculations performed for a toroidal annulus geometry. For the simulation parameters considered, similar global behaviour is observed in both cases, whereas strong local suppression of turbulence fluctuations is demonstrated near the X-point for the case of a single-null geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Versatile multi-energy hard x-ray camera to study confined and unconfined fast electron dynamics and anisotropies (Invited)

A powerful and flexible hard x-ray (HXR) camera has been recently installed and tested on the WEST tokamak (CEA, France) in collaboration with the Princeton Plasma Physics Laboratory. The diagnostic is a pinhole camera fielded with a 2D pixel detector equipped with a 1 mm thick CdTe sensor. The novelty of this diagnostic technique is the detector’s capability of adjusting the threshold energy at the pixel level. This innovation provides great flexibility in the energy configuration, allowing simultaneous space, energy, and time resolved x-ray measurements. The novel camera has been used to measure the core radiation from non-Maxwellian (fast) electrons accelerated by Lower Hybrid (LH) waves and also the beam–target emission of tungsten in the divertor region produced by fast electron losses interacting with the target. In addition, anisotropic hard x-ray emission has been detected for the first time at the WEST core and edge plasma, with opposite toroidal intensity trends. Finally, experimental vertical and toroidal HXR profiles have been successfully reproduced with the LH code LUKE.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Manufactured Home Equipment Sizing and Selection Tool

This tool is for use for manufactured homes built to the HUD Code (thermal provisions of the Manufactured Housing Construction and Safety Standards, Subpart F, Section §3280), including manufactured homes certified as ENERGY STAR® and DOE Efficient New Homes. It provides a simple look-up procedure to assist equipment specifiers, HVAC contractors, home installers, retailers, manufacturers, and electric utility staff. For new manufactured homes, it guides in selection of single speed heat pumps and air conditioners. For existing homes, it determines design conditions and nominal equipment size. The guidelines provide a simple but accurate procedure that reduces oversizing and its associated problems.

Source record↗
DOE OSTI2022

1D Heat Loss Models Validation Experiment

Contains data from the model validation in the 1D Heat Loss Models to Predict the Aquifer Temperature Profile during Hot/Cold Water Injection Project. The data include two COMSOL models (2D axisymmetric benchmark model and 2D Vinsome model), one python code (1D Vinsome based FEM numerical simulation), one matlab main code (1D Newton analytical solution and all results comparison visualization), and output files generated from the above models.

1D↗
DOE OSTI2024

Software stewardship and advancement of a high-performance computing scientific application: QMCPACK

Here, we provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions own runners, and NVIDIA and AMD GPUs used in pre-exascale systems, (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Docker containers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the stewardship and advancement of community HPC codes to enable scientific discovery at scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Deuterium Extraction from Helium with a Vanadium Vacuum Permeator: Commissioning of the Tritium Extraction eXperiment (TEX)

Here, the Tritium Extraction eXperiment (TEX) is a forced-convection lead-lithium (PbLi) loop in the Safety and Tritium Applied Research (STAR) facility at Idaho National Laboratory (INL) with the purpose of providing validation data for the vacuum permeator tritium extraction concept. A vanadium tube of 1000 mm length, 12.7 mm outside diameter, and 0.50 mm wall thickness is installed in the test section of TEX. The installed vanadium tube is characterized to quantify impurity concentrations, surface chemistry, and microstructure to elucidate permeation phenomena observed in experimentation. Herein, the permeation properties of the vanadium tube are characterized by measuring deuterium permeation at 300 °C, 325 °C, and 350 °C at 100 kPa, 125 kPa, and 150 kPa total pressures with 5000 ppm deuterium in helium gas mixture in a once-through flow configuration. The hydrogen isotope permeation through the vanadium tube in the test section is measured with quadrupole mass spectrometers and the hydrogen isotope concentration in the feed and retentate gas stream is measured with gas chromatography. The transient permeation results are modeled with MELCOR-TMAP, a thermal-hydraulic tritium transport code, and compared well with literature data.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evaluation of Low Environmental Impact Distributed Scroll Booster Technology for Supermarket Refrigeration

A distributed scroll booster refrigeration system is proposed as an alternative for current refrigeration systems used in supermarket applications. The distributed scroll booster system is configured to enable the use of low-pressure refrigerants such as R-1234yf, which has a global warming potential (GWP) of less than 1, in a cycle configuration that employs a two-stage compression process. This scroll booster system is comprised of medium-temperature compressors (high stage) coupled with low-temperature compressors (low stage). This two-stage system architecture allows both low-temperature and medium-temperature refrigerated display cases to be integrated in the same system. By design, this system reduces the refrigerant charge, produces higher efficiencies (two-stage compression), uses existing components minimizing cost impact, and will be able to use ultra-low GWP class A2L (low toxicity and low flammability) refrigerants such as R-1234yf when standards and building codes are updated (expected by 2025).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

AMPX and SCALE Nuclear Data Libraries for Depletion [Slides]

This presentation touches on AMPX and SCALE nuclear data libraries for depletion. AMPX is a nuclear data processing code suite that generates continuous-energy and multigroup neutron and gamma libraries used in SCALE and VERA, covariance data used in S/U calculations, and depletion, activation and decay origin data for Oak Ridge Isotope Generation code (ORIGEN).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Calculated covariance matrices for fission product yields using BeoH

Fission product yields (FPY) are important for a variety of applications (reactor neutronics, spent fuel, dosimetry, radiochemistry, etc.) and are currently included in many of the evaluated libraries around the world. The FPYs in the current US evaluation, ENDF/B-VIII.0, are mainly based on the 1994 evaluation of England and Rider and have only had slight updates—such as the inclusion of a 2 MeV point for 239 Pu—since their development. Additionally, only mean values and uncertainties are included in the evaluation, not full correlations. Los Alamos National Laboratory, in collaboration with several other institutes, has been working on an updated evaluation for the FPYs of 239 Pu(n,f), 235 U(n,f), 238 U(n,f), and 252 Cf(sf) using the deterministic, Hauser-Feshbach, fission fragment decay code, BeoH. BeoH calculates the FPYs consistently with many other prompt and delayed fission observables, explicitly taking into account multi-chance fission and ensuring consistency between observables. In addition to providing updated means and uncertainties for the FPYs on a pointwise energy grid from thermal to 20 MeV, we calculate correlations between all FPYs at each incident energy and across incident energies. Here, we discuss the development of these covariance matrices, differences in the correlations between FPYs based on the parameters that are included in the model optimization, and correlations across incident energies for neutron-induced fission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electric field and turbulence in global Braginskii simulations across the ASDEX Upgrade edge and scrape-off layer

Turbulence simulations in diverted geometry across the edge and scrape-off layer (SOL) of ASDEX Upgrade are performed with the GRILLIX code (Stegmeir et al 2019 Phys. Plasmas 26 052517). The underlying global (full-f) drift-reduced Braginskii model allows to concurrently study the self-consistent dynamics of the turbulence and the background as well as the evolution of toroidal and zonal flows. Different contributions to the radial electric field are identified. The dominant contribution on closed flux surfaces comes from the ion pressure gradient, due to the diamagnetic drift in the curved magnetic field. Large deviations can be induced, in particular, by the polarization particle flux, leading to zonal flows. The latter are driven by small-scale eddies, but do not exhibit much impact on the overall transport which is driven by ballooning modes at larger scales. Ion viscosity is found to be important in damping poloidal rotation through adjusting of the parallel velocity profile, but not via direct vorticity damping. The zonal flow drive peaks at the separatrix, where a strong shear layer forms due to the sheath-induced counter-propagating SOL flow, allowing for the formation of a transport barrier. The temperature profile across the separatrix is determined by the competition between cross-field transport and outflow in the SOL, the latter being largely controlled by the parallel heat conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗
DOE OSTI2023

The DESC stellarator code suite. Part 2. Perturbation and continuation methods

A new perturbation and continuation method is presented for computing and analysing stellarator equilibria. The method is formally derived from a series expansion about the equilibrium condition $\boldsymbol {F} \equiv \boldsymbol {J}\times \boldsymbol {B} - \boldsymbol {\nabla } p = 0$ , and an efficient algorithm for computing solutions to second- and third-order perturbations is developed. The method has been implemented in the DESC stellarator equilibrium code, using automatic differentiation to compute the required derivatives. Examples are shown demonstrating its use for computing complicated equilibria, perturbing a tokamak into a stellarator and performing parameter scans in pressure, rotational transform and boundary shape in a fraction of the time required for a full solution.

Conlin, Rory (ORCID:0000000183662111)↗
DOE OSTI2024

Aeroelastic code comparison using the IEA 22MW reference turbine

Reference wind turbine designs and the associated aeroelastic models are widely used in both research and industry. Reference models representing future concepts are of particular interest. Current state of the art aeroelastic tools are relied upon to design the next generation of large wind turbines. However, modelling assumptions may be invalidated by upcoming very large turbines, and different aeroelastic tools may give inconsistent results. A 22MW turbine model has been defined as part of International Energy Agency (IEA) Wind Task 55 on Reference Wind Turbines and Farms to represent future turbines to be deployed in the 2030s. In this study, an aeroelastic model of this turbine has been created in four tools; Bladed, HAWC2, OpenFAST, and QBlade. Code comparisons are presented for steady state operation, linear stability analysis, and time domain power production simulations in steady and turbulent wind. Generally, the codes show a good agreement, but with some differences present in the linear stability analysis, periodic azimuthal variation, and time domain simulations. The models are a good basis for further study with the IEA 22MW turbine, and further code comparison exercises.

17 WIND ENERGY↗
DOE OSTI2025

Describing hadronization via histories and observables for Monte-Carlo event reweighting

We introduce a novel method for extracting a fragmentation model directly from experimental data without requiring an explicit parametric form, called Histories and Observables for Monte-Carlo Event Reweighting (HOMER), consisting of three steps: the training of a classifier between simulation and data, the inference of single fragmentation weights, and the calculation of the weight for the full hadronization chain. We illustrate the use of HOMER on a simplified hadronization problem, a q\bar{q} q q ‾ string fragmenting into pions, and extract a modified Lund string fragmentation function f(z) f ( z ) . We then demonstrate the use of HOMER on three types of experimental data: (i) binned distributions of high-level observables, (ii) unbinned event-by-event distributions of these observables, and (iii) full particle cloud information. After demonstrating that f(z) f ( z ) can be extracted from data (the inverse of hadronization), we also show that, at least in this limited setup, the fidelity of the extracted f(z) f ( z ) suffers only limited loss when moving from (i) to (ii) to (iii). Public code is available at https://gitlab.com/uchep/mlhad.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Report on PTT Imaging of Defects in AM Metallic Materials-Part 2

Metal Additive Manufacturing (AM) is a promising method for cost-efficient fabrication of complex shape structures for applications in harsh environment, such as in a nuclear reactor. However, internal defects (pores) occur in high-strength AM alloys, which are manufactured with Laser Powder Bed Fusion (LPBF) AM method. Pulsed Infrared Thermography (PIT) is an efficient nondestructive evaluation (NDE) method to examine actual structures, because this method offers one-sided non-contact measurements, and fast processing of large sample areas. However, imaging of material defects, particularly defects with sizes at microscopic level, is challenging. In this report, we benchmark the performance of several Unsupervised Learning (UL) algorithms designed to enhance imaging of microscopic defects in metals with PIT. UL aims to learn the latent principal patterns (dictionaries) in PIT data to detect defects with minimal human supervision. Performance of Independent Component Analysis (ICA), Sparse Coding (SC), Principal Component Analysis (PCA) and Exploratory Factor Analysis (EFA) was compared using F-score, UL model training time and defects reconstruction time. We obtained the average F-score of 0.75, and a highest F-score of 0.89 for the EFA algorithm. Overall, EFA outperforms other UL algorithms considered in this study.

36 MATERIALS SCIENCE↗