Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “V 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 649 records · Page 36

Hybrid simulations of thin current sheets

A one-dimensional, hybrid simulation code is used to study current sheets with a nonzero normal magnetic field B(sub z) and a dawn-to-dusk electric field E(sub y). Such configurations are dependent upon only two parameters: we use the normalized normal magnetic field B-normalized (sub z) = B(sub z)/(4(pi)(n(sub b)) (v(exp 2 sub T))(exp 1/2) and normalized electric field V-normalized (sub D) = (1/V(sub T)(cE(sub y)/B(sub z)), where V(sub T) is the thermal velocity of ions prior to their interaction with the current sheet and n(sub b) is the number density outside the current sheet (at the simulation boundary). A third parameter that is relevant to the motion of particles in current sheets is kappa(sub A), the value of kappa = (R(sub min)/rho(sub max))(exp 1/2) for particles of average energy. We find that if either B-normalized (sub z) is close to or greater than 1, or if kappa(sub A) is close to 1, a rotational mode develops in which the z = 0 current rotates with the ion sense about the normal magnetic field, while for small values of both B-normalized (sub z) or kappa(sub A), the configuration is quasi-steady. To achieve values of kappa(sub A) of the order of or larger than 1, we decrease the value of V-normalized (sub D) uniformly. We find that the magnetic field fluctuations and particle distribution functions are similar in many respects to what was observed in the day 240, 1986, Active Magnetospheric Particle Tracer Explorer (AMPTE)/CCE current disruption event, an event that appears to be located at the site of initiation of current disruption and related particle energization.

Burkhart, G. R.↗

Reacting Multi-Species Gas Capability for USM3D Flow Solver

The USM3D Navier-Stokes flow solver contributed heavily to the NASA Constellation Project (CxP) as a highly productive computational tool for generating the aerodynamic databases for the Ares I and V launch vehicles and Orion launch abort vehicle (LAV). USM3D is currently limited to ideal-gas flows, which are not adequate for modeling the chemistry or temperature effects of hot-gas jet flows. This task was initiated to create an efficient implementation of multi-species gas and equilibrium chemistry into the USM3D code to improve its predictive capabilities for hot jet impingement effects. The goal of this NASA Engineering and Safety Center (NESC) assessment was to implement and validate a simulation capability to handle real-gas effects in the USM3D code. This document contains the outcome of the NESC assessment.

Frink, Neal T.↗

Modernizing Fortran 77 legacy codes

Legacy software has great value since it is generally well debugged, produces results that are trusted and is actively meeting end-user goals. The amount of hidden expert knowledge embedded in such software can be significant, making its preservation important.

Fortran scientific programming↗

Discovering invariant spatial features in electron energy loss spectroscopy images on the mesoscopic and atomic levels

Over the last two decades, Electron Energy Loss Spectroscopy (EELS) imaging with a scanning transmission electron microscope has emerged as a technique of choice for visualizing complex chemical, electronic, plasmonic, and phononic phenomena in complex materials and structures. The availability of the EELS data necessitates the development of methods to analyze multidimensional data sets with complex spatial and energy structures. Traditionally, the analysis of these data sets has been based on analysis of individual spectra, one at a time, whereas the spatial structure and correlations between individual spatial pixels containing the relevant information of the physics of underpinning processes have generally been ignored and analyzed only via the visualization as 2D maps. Here, we develop a machine learning-based approach and workflows for the analysis of spatial structures in 3D EELS data sets using a combination of dimensionality reduction and multichannel rotationally invariant variational autoencoders. This approach is illustrated for the analysis of both the plasmonic phenomena in a system of nanowires and in the core excitations in functional oxides using low loss and core-loss EELS, respectively. The code developed in this manuscript is open sourced and freely available and provided as a Jupyter notebook for the interested reader.

36 MATERIALS SCIENCE↗

Early results in the search for extreme coronal line emitters with the Dark Energy Spectroscopic Instrument

Here, we present the results of our search through the early data release (EDR) of the Dark Energy Spectroscopic Instrument (DESI) for extreme coronal line emitters (ECLEs) – a rare classification of galaxies displaying strong, high-ionization iron coronal emission lines within their spectra. With the requirement of a strong X-ray continuum to generate the coronal emission, ECLEs have been linked to both active galactic nuclei (AGNs) and tidal disruption events (TDEs). We focus our search on identifying TDE-linked ECLEs. We identify three such objects within the EDR sample, highlighting DESI’s effectiveness for discovering new nuclear transients, and determine a galaxy-normalized TDE-linked ECLE rate of RG=5 −3+5×10−6 galaxy−1 yr−1 at a median redshift of z=0.2 – broadly consistent with previous work. Additionally, we identify more than 200 AGNs displaying coronal emission lines, which serve as the primary astrophysical contaminants in searches for TDE-related events. We also include an outline of the custom python code developed for this search.

Clark, Peter [Southampton U.; Portsmouth U., ICG] ↗

Toward a 2D Local Implementation of Quantum Low-Density Parity-Check Codes

Geometric locality is an important theoretical and practical factor for quantum low-density parity-check (qLDPC) codes that affects code performance and ease of physical realization. For device architectures restricted to two-dimensional (2D) local gates, naively implementing the high-rate codes suitable for low-overhead fault-tolerant quantum computing incurs prohibitive overhead. In this work, we present an error-correction protocol built on a bilayer architecture that aims to reduce operational overheads when restricted to 2D local gates by measuring some generators less frequently than others. We investigate the family of bivariate-bicycle qLDPC codes and show that they are well suited for a parallel syndrome-measurement scheme using fast routing with local operations and classical communication (LOCC). Through circuit-level simulations, we find that in some parameter regimes, bivariate-bicycle codes implemented with this protocol have logical error rates comparable to the surface code while using fewer physical qubits. Published by the American Physical Society 2025

Berthusen, Noah (ORCID:0000000275862786)↗

Celeritas: GPU-accelerated particle transport for detector simulation in High Energy Physics experiments

Within the next decade, experimental High Energy Physics (HEP) will enter a new era of scientific discovery through a set of targeted programs recommended by the Particle Physics Project Prioritization Panel (P5), including the upcoming High Luminosity Large Hadron Collider (LHC) HL-LHC upgrade and the Deep Underground Neutrino Experiment (DUNE). These efforts in the Energy and Intensity Frontiers will require an unprecedented amount of computational capacity on many fronts including Monte Carlo (MC) detector simulation. In order to alleviate this impending computational bottleneck, the Celeritas MC particle transport code is designed to leverage the new generation of heterogeneous computer architectures, including the exascale computing power of U.S. Department of Energy (DOE) Leadership Computing Facilities (LCFs), to model targeted HEP detector problems at the full fidelity of Geant4. This paper presents the planned roadmap for Celeritas, including its proposed code architecture, physics capabilities, and strategies for integrating it with existing and future experimental HEP computing workflows.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

GDSA Framework Development and Process Model Integration FY2022

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is disposal system modeling (Sassani et al. 2021). The SFWST Geologic Disposal Safety Assessment (GDSA) work package is charged with developing a disposal system modeling and analysis capability for evaluating generic disposal system performance for nuclear waste in geologic media. This report describes fiscal year (FY) 2022 advances of the Geologic Disposal Safety Assessment (GDSA) performance assessment (PA) development groups of the SFWST Campaign. The common mission of these groups is to develop a geologic disposal system modeling capability for nuclear waste that can be used to assess probabilistically the performance of generic disposal options and generic sites. The modeling capability under development is called GDSA Framework (pa.sandia.gov). GDSA Framework is a coordinated set of codes and databases designed for probabilistically simulating the release and transport of disposed radionuclides from a repository to the biosphere for post-closure performance assessment. Primary components of GDSA Framework include PFLOTRAN to simulate the major features, events, and processes (FEPs) over time, Dakota to propagate uncertainty and analyze sensitivities, meshing codes to define the domain, and various other software for rendering properties, processing data, and visualizing results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Type 1 Diabetes Genetic Risk in 109,954 Veterans With Adult-Onset Diabetes: The Million Veteran Program (MVP)

OBJECTIVE To characterize high type 1 diabetes (T1D) genetic risk in a population where type 2 diabetes (T2D) predominates. RESEARCH DESIGN AND METHODS Characteristics typically associated with T1D were assessed in 109,594 Million Veteran Program participants with adult-onset diabetes, 2011–2021, who had T1D genetic risk scores (GRS) defined as low (0 to <45%), medium (45 to <90%), high (90 to <95%), or highest (≥95%). RESULTS T1D characteristics increased progressively with higher genetic risk (P < 0.001 for trend). A GRS ≥90% was more common with diabetes diagnoses before age 40 years, but 95% of those participants were diagnosed at age ≥40 years, and their characteristics resembled those of individuals with T2D in mean age (64.3 years) and BMI (32.3 kg/m2). Compared with the low-risk group, the highest-risk group was more likely to have diabetic ketoacidosis (low GRS 0.9% vs. highest GRS 3.7%), hypoglycemia prompting emergency visits (3.7% vs. 5.8%), outpatient plasma glucose <50 mg/dL (7.5% vs. 13.4%), a shorter median time to start insulin (3.5 vs. 1.4 years), use of a T1D diagnostic code (16.3% vs. 28.1%), low C-peptide levels if tested (1.8% vs. 32.4%), and glutamic acid decarboxylase antibodies (6.9% vs. 45.2%), all P < 0.001. CONCLUSIONS Characteristics associated with T1D were increased with higher genetic risk, and especially with the top 10% of risk. However, the age and BMI of those participants resemble those of people with T2D, and a substantial proportion did not have diagnostic testing or use of T1D diagnostic codes. T1D genetic screening could be used to aid identification of adult-onset T1D in settings in which T2D predominates.

Yang, Peter K. (ORCID:0000000193796981)↗

Analysis of the Tanana River Basin using LANDSAT data

Digital image classification techniques were used to classify land cover/resource information in the Tanana River Basin of Alaska. Portions of four scenes of LANDSAT digital data were analyzed using computer systems at Ames Research Center in an unsupervised approach to derive cluster statistics. The spectral classes were identified using the IDIMS display and color infrared photography. Classification errors were corrected using stratification procedures. The classification scheme resulted in the following eleven categories; sedimented/shallow water, clear/deep water, coniferous forest, mixed forest, deciduous forest, shrub and grass, bog, alpine tundra, barrens, snow and ice, and cultural features. Color coded maps and acreage summaries of the major land cover categories were generated for selected USGS quadrangles (1:250,000) which lie within the drainage basin. The project was completed within six months.

Morrissey, L. A.↗

On Flowfield Periodicity in the NASA Transonic Flutter Cascade: Numerical Study - Part 2

The transonic flutter cascade facility at NASA Glenn Research Center was redesigned based on a combined program of experimental measurements and numerical analyses. The objectives of the redesign were to improve the periodicity of the cascade in steady operation, and to better quantify the inlet and exit flow conditions needed for CFD predictions. Part I of this paper describes the experimental measurements, which included static pressure measurements on the blade and endwalls made using both static taps and pressure sensitive paints, cobra probe measurements of the endwall boundary layers and blade wakes, and shadowgraphs of the wave structure. Part II of this paper describes three CFD codes used to analyze the facility, including a multibody panel code, a quasi-three-dimensional viscous code, and a fully three-dimensional viscous code. The measurements and analyses both showed that the operation of the cascade was heavily dependent on the configuration of the sidewalls. Four configurations of the sidewalls were studied and the results are described. For the final configuration, the quasi-three-dimensional viscous code was used to predict the location of mid-passage streamlines for a perfectly periodic cascade. By arranging the tunnel sidewalls to approximate these streamlines, sidewall interference was minimized and excellent periodicity was obtained.

Chima, Rodrick V.↗

A review on charged-particle transport modeling for laser direct-drive fusion

Inertial confinement fusion (ICF) with the laser-indirect-drive scheme has recently made a tremendous breakthrough recently after decades of intensive research effort. Taking this success to the next step, the ICF community is coming to a general consensus that laser direct-drive (LDD) fusion might be the viable way for enabling inertial fusion energy (IFE) and high-gain targets for other applications. Designing and understanding LDD fusion targets heavily rely on radiation-hydrodynamic code simulations, in which charged-particle transport plays an essential role in modeling laser-target energy coupling and bootstrap heating of fusion-produced α-particles. To better simulate charged-particle transport in LDD targets, over the past four decades the plasma physics community has advanced transport calculations from simple plasma physics models to sophisticated computations based on first-principles methods. In this review, we give an overview of the current status of charged-particle transport modeling for LDD fusion, including what challenges we still face and the possible paths moving forward to advance transport modeling for ICF simulations. We hope this review will provide a summary of exciting challenges to stimulate young minds to enter the field, facilitate further progress in understanding warm-dense matter physics, and ultimately bridge toward the success of reliable LDD fusion designs for IFE and other high-gain ICF applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Geometric Discontinuities on Creep Behavior of Alloy 617

Due to the excellent long-term high-temperature performance of Alloy 617, the Alloy has been identified as the primary structural material for the intermediate heat exchanger of the Very High Temperature Reactor (VHTR). Both the base and welded alloy has been qualified to be used up to 950 °C up to 100,000 hours by the American Society of Mechanical Engineer (ASME) Boiler and Pressure Vessel Code (BPVC). The qualification properties of the material were exclusively determined by uniaxial creep tests. The components may, however, experience multiaxial stress conditions in service due to the reactor geometry. This study used notched creep tests to investigate the effect of geometric discontinuities on creep behavior of both the base and welded Alloy 617. The V-notch tests reveal that the geometric discontinuities do not reduce the creep rupture live under high and intermediate stress conditions. The U-notch tests show that a larger radius notch behave in a similar way as a straight gauge section of a uniaxial test. The creep rupture live of the welded alloy was slightly longer than the base metal.

36 - MATERIALS SCIENCE↗

A Non-perturbative Approach to Computing Seismic Normal Modes in Rotating Planets

In this work, a continuous Galerkin method based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial filtering eigensolver. The relevant elastic-gravitational system of equations, including the Coriolis force, is subjected to a mixed finite-element method, while self-gravitation is accounted for with the fast multipole method. Our discretization utilizes fully unstructured tetrahedral meshes for both solid and fluid regions. The relevant eigenvalue problem is solved by a combination of several highly parallel and computationally efficient methods. We validate our three-dimensional results in the non-rotating case using analytical results for constant elastic balls, as well as numerical results for an isotropic Earth model from standard “radial” algorithms. We also validate the computations in the rotating case, but only in the slowly-rotating regime where perturbation theory applies, because no other independent algorithms are available in the general case. The algorithm and code are used to compute the point spectra of eigenfrequencies in several Earth and Mars models studying the effects of heterogeneity on a large range of scales.

58 GEOSCIENCES↗

Transient magnetic field diffusion considerations relevant to magnetically assisted indirect drive inertial confinement fusion

We report application of a magnetic field to an indirect drive inertial confinement fusion target requires diffusion of the field through the high-Z and electrically conducting Hohlraum. The onset of the external field generates eddy currents in the Hohlraum wall that result in (1) a reduction of the peak field at the capsule, (2) heating of the Hohlraum wall through Ohmic dissipation, and (3) wall movement due to the inward force from the eddy current interacting with the field. Heating of the wall causes an increase in blackbody radiation which can preheat the capsule and frozen deuterium–tritium fuel, while wall motion leads to potential misalignment of the lasers at the Hohlraum wall. Limiting these detrimental effects sets requirements on the tolerable magnitude of each effect. We present a nonlinear model for B-field diffusion through an infinitely long thin-walled cylinder with a temperature dependent resistivity, to show that a 15 μm thick wall of pure gold fails to meet these requirements. A new Hohlraum material made from an alloy of Au and Ta has a measured resistivity of ≥60 times that of Au and is shown with the nonlinear model to meet the requirements for magnetization. We compare the nonlinear model to simulations of the actual Hohlraum target using a finite element code which includes temperature-dependent Hohlraum resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite element implementation of viscoplastic models

A brief description of the implementation in MARK, the general purpose finite element structural analysis code, of two viscoplastic models developed by Robinson is given. One model is for isotropic materials and the other is for metal matrix composites. Also presented are analytical results obtained for hot section components using these models.

Arya, V. K.↗

Simultaneous Estimation of Model State Variables and Observation and Forecast Biases Using a Two-Stage Hybrid Kalman Filter

In this paper, we present a two-stage hybrid Kalman filter to estimate both observation and forecast bias in hydrologic models, in addition to state variables. The biases are estimated using the discrete Kalman filter, and the state variables using the ensemble Kalman filter. A key issue in this multi-component assimilation scheme is the exact partitioning of the difference between observation and forecasts into state, forecast bias and observation bias updates. Here, the error covariances of the forecast bias and the unbiased states are calculated as constant fractions of the biased state error covariance, and the observation bias error covariance is a function of the observation prediction error covariance. In a series of synthetic experiments, focusing on the assimilation of discharge into a rainfall-runoff model, it is shown that both static and dynamic observation and forecast biases can be successfully estimated. The results indicate a strong improvement in the estimation of the state variables and resulting discharge as opposed to the use of a bias-unaware ensemble Kalman filter. Furthermore, minimal code modification in existing data assimilation software is needed to implement the method. The results suggest that a better performance of data assimilation methods should be possible if both forecast and observation biases are taken into account.

Data Assimilation↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D

Full seven-equation Reynolds stress turbulence models are a relatively new and promising tool for todays aerospace technology challenges. This paper uses two stress-omega full Reynolds stress models to evaluate challenging flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSG/LRR full second-moment Reynolds stress models have been implemented into the FUN3D (Fully Unstructured Navier-Stokes Three Dimensional) unstructured Navier-Stokes code and are evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST-V) two-equation turbulence models.

reynolds averaging↗