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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Verification of global gyrokinetic simulation of low frequency mode excited by thermal plasma in spherical tokamak

A low frequency Alfven eigenmode in the ST40 spherical tokamak (ST) is found to be resonantly excited by thermal ions from global gyrokinetic GTC simulations. The analyses of phase-space resonances and wave-particle energy exchanges show that the Beta-induced Alfvén-Acoustic Eigenmode (BAAE) is primarily driven unstable by thermal ions, but is stabilized by energetic particles (EPs). When all kinetic effects are artificially suppressed, GTC simulations in the ideal magnetohydrodynamic (MHD) limit confirm the existence of the MHD eigenmode with a frequency in the BAAEs gap, as previously found by the ideal MHD code NOVA. The kinetic effects of thermal ions induce a frequency downshift compared to the stable BAAE in the MHD simulations. These results support the verification for global gyrokinetic simulations of the BAAE in the ST. However, the simulations do not explain the experimental observations of the BAAE in the ST40, possibly owing to the limitation of the collisionless gyrokinetic Monte Carlo model for the EPs.

Alfvén eigenmode↗
DOE OSTI2024

An international benchmark for wind plant wakes from the American WAKE ExperimeNt (AWAKEN)

This article introduces the first benchmark study within the International Energy Agency Wind Task 57 framework, focusing on wind plant wakes. Leveraging data from the American WAKE ExperimeNt (AWAKEN), the benchmark aims to assess the accuracy of simulation tools in modeling wind plant wakes and their impact on the downstream flow under diverse inflow conditions. The AWAKEN field campaign, conducted in Oklahoma from 2022 to 2024, provides unprecedented observations of wind plant-atmosphere interactions, thus offering a large dataset to validate numerical models of different complexity. The benchmark will include three phases—code calibration, blind comparison, and iteration—allowing participants to refine their numerical models based on the feedback from the benchmark team. This article describes the benchmark case study selected from observations providing details on atmospheric conditions, wake evidence, and wind turbine operation. The benchmark’s structure and timeline, along with the expected publication of results, are discussed as well. This collaborative effort aims to enhance the accuracy of wind plant wake simulations, thus contributing to the improvement of wind energy production estimates.

17 WIND ENERGY↗
NASA NTRS1999

Analysis of X-Ray Microradiographs of Al-Au Interface Quench Profile using Modeling of Solidification Including Double-Diffusion and Convection in the Melt

Experimental data on Al-0.8Au horizontal solidification of a 1 mm thick specimen in a BN crucible shows the effect of growth rate on the solidification interface shape. For translation rates below 0.5 micron/s the interface maintains a plain and flat shape. When the translation rate is 3 to 5 micron/s or more, the interface appearance changes to two planar zones, with the zone closer to the bottom having higher inclination. The interface shapes were measured by first quenching in place during growth. X-ray microscopy shows the interface shape within the quenched sample by viewing through the side of the specimen. In order to provide theoretical explanation of the phenomena, numerical modeling was undertaken using finite element code FIDAP. Double diffusion convection in Al-0.8Au melt and crystal-melt interface curvature during directional solidification was analyzed numerically. Actual thermophysical properties of Al-0.8Au including the binary Al-Au phase diagram were used. Although convection in the sample is weak, for the slower translation rate convection and diffusion is sufficient for the redistribution of initial compositional stratification caused by gravity. When translation rate is raised, neither convection nor diffusion can provide proper mixing so that solidification temperatures differ significantly near the bottom within the bulk of the sample. As a result, the solid-liquid interface appears to have two planar zones with different inclination.

Bune, Andris V.↗
NASA NTRS2001

Design of Critical Components

Critical component design is based on minimizing product failures that results in loss of life. Potential catastrophic failures are reduced to secondary failures where components removed for cause or operating time in the system. Issues of liability and cost of component removal become of paramount importance. Deterministic design with factors of safety and probabilistic design address but lack the essential characteristics for the design of critical components. In deterministic design and fabrication there are heuristic rules and safety factors developed over time for large sets of structural/material components. These factors did not come without cost. Many designs failed and many rules (codes) have standing committees to oversee their proper usage and enforcement. In probabilistic design, not only are failures a given, the failures are calculated; an element of risk is assumed based on empirical failure data for large classes of component operations. Failure of a class of components can be predicted, yet one can not predict when a specific component will fail. The analogy is to the life insurance industry where very careful statistics are book-kept on classes of individuals. For a specific class, life span can be predicted within statistical limits, yet life-span of a specific element of that class can not be predicted.

Hendricks, Robert C.↗
DOE OSTI2023

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Model data for Flood Frequency Analysis using Stochastic Storm Transposition and an Integrated Surface-Subsurface Hydrological Model

This archived provides scripts and input files used for the implementation of a novel approach to conduct process-based Flood Frequency Analysis using a Stochastic Storm Transposition (SST) and an Integrated Surface-Subsurface Hydrological Model (ISSHM). As a proof-of-concept, this study uses the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct flood frequency analysis by simulating the flood response to 5,000 annual synthetic storm events in a ~2000 km2 Southeast Texas watershed.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include TXT, CSV, DAT, SBATCH, SHP, TIF, NetCDF, and HDF5 files, which can be read through Python scripts. The input files for the ATS model and RainyDay model have .XML and .SST extensions, respectively, and can be edited in any commonly used text editors.This archive contains:* Scripts and data files essential for generating the ATS model input. It uses the Watershed Workflow package to produce both mesh and ATS input files. * Jupyter notebooks designated for the ATS model evaluation, covering both long-term simulations and 40 rainfall-runoff events.* Input files required to simulate SST storm events using RainyDay.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Numerical computation of space shuttle heating and surface streamlines

Exact inviscid flow-field codes are used together with a quasi-three-dimensional boundary-layer analysis to provide estimates of the windward surface heating and streamline patterns of the shuttle orbiter vehicle. The accuracy and limitations of the methods are established by comparison with available wind-tunnel experiments and with more exact numerical solutions for simple flows. Flight predictions are presented showing the effects of finite-rate (nonequilibrium) chemical reactions, and the effects of varying boundary-layer edge conditions due to the growth of the boundary-layer into the inviscid flow (entropy layer swallowing). Differences between flow-field predictions at wind-tunnel and nominal flight conditions are discussed.

Rakich, J. V.↗
NASA NTRS1989

Prediction of the aerodynamic environment and heat transfer for rotor-stator configurations

A numerical study of the aerodynamic and thermal environment associated with axial turbine stages is presented. Computations were performed using a modification of the unsteady viscous code, ROTORI, and an improved version of the steady inviscid cascade system, MERIDL-TSONIC, coupled with boundary layer codes, BLAYER and STAN5. Two different turbine stages were analyzed: the first stage of the United Technologies Research Center Large Scale Rotating Rig (LSRR) and the first stage of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump turbine. The time-averaged airfoil midspan pressure and heat transfer profiles were predicted for numerous thermal boundary conditions including adiabatic wall, prescribed surface temperature, and prescribed heat flux. Computed solutions are compared with each other and with experimental data in the case of the LSRR calculations. Modified ROTORI predictions of unsteady pressure envelopes and instantaneous contour plots are also presented. Relative merits of the two computational approaches are discussed.

Griffin, L. W.↗
DOE OSTI2025

Characterization of W production during ICRF operations: experiments and modeling

For successfully heating plasma with waves in the ion cyclotron range of frequencies (ICRFs), mitigating impurity production is just as crucial as maximizing power coupling, especially in high-Z environments. ICRF can effectively deposit energy on ions, modify turbulence-driven transport, and enhance fusion reaction efficiency, but only when its power coupling has minimal impact on impurity production. To do so, one must rely on a toroidal array of at least three active elements excited with appropriate phasing and power ratio to reduce the currents induced on the antenna frame below levels critical for physical sputtering. In contrast to classic two-strap antennas, which are optimized for dipole phasing with equal power on both straps, three-strap antennas in ASDEX Upgrade (AUG)—but also four-strap antennas in JET, Alcator C-Mod, SPARC and ITER—offer the possibility to act also on the power ratio between the central and outer straps. With optimal settings, impurity production can be reduced substantially, making the ICRF compatible with the high-Z wall. This paper explores the characteristics of the AUG three-strap antennas in terms of impurity production, as well as the key role of plasma composition in this process. Numerical simulations were performed using SSWICH and Petra-M (finite element codes) to quantify impurity production and compare with experimental results. Energies of ions falling on antenna limiters (measured with probes) are well predicted by both codes. These tools are then used to further describe the source of the impurity, namely the gross erosion of tungsten from an ICRF antenna, for different plasma mixtures. Results are also compared to spectroscopy data. Ultimately, we show that deleterious effects of the ICRF on plasma surface interactions will be weaker in plasmas containing larger fractions of highly ionized heavier low-Z impurity, which is typically relevant for experiments relying on impurity seeding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Speeding up particle track reconstruction using a parallel Kalman filter algorithm

One of the most computationally difficult problems expected for the High-Luminosity Large Hadron Collider (HL-LHC) is determining the trajectory of charged particles during event reconstruction. Algorithms used at the LHC today rely on Kalman filtering, which builds physical trajectories incrementally while incorporating material effects and error estimation. Recognizing the need for faster computational throughput, we have adapted Kalman-filter-based methods for highly parallel, many-core SIMD architectures that are now prevalent in high-performance hardware. In this paper, we discuss the design and performance of the improved tracking algorithm, referred to as mkFit. A key piece of the algorithm is the Matriplex library, containing dedicated code to optimally vectorize operations on small matrices. The physics performance of the mkFit algorithm is comparable to the nominal CMS tracking algorithm when reconstructing tracks from simulated proton-proton collisions within the CMS detector. We study the scaling of the algorithm as a function of the parallel resources utilized and find large speedups both from vectorization and multi-threading. mkFit achieves a speedup of a factor of 6 compared to the nominal algorithm when run in a single-threaded application within the CMS software framework.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Vertex-finding and reconstruction of contained two-track neutrino events in the MicroBooNE detector

In this work, we describe algorithms developed to isolate and accurately reconstruct two-track events that are contained within the MicroBooNE detector. This method is optimized to reconstruct two tracks of lengths longer than 5cm. This code has applications to searches for neutrino oscillations and measurements of cross sections using quasi-elastic-like charged current events. The algorithms we discuss will be applicable to all detectors running in Fermilab's Short Baseline Neutrino program (SBN), and to any future liquid argon time projection chamber (LArTPC) experiment with beam energies ~ 1 GeV. The algorithms are publicly available on a GITHUB repository. This reconstruction offers a complementary and independent alternative to the Pandora reconstruction package currently in use in LArTPC experiments, and provides similar reconstruction performance for two-track events.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗
NASA NTRS1999

Development, Analysis and Testing of the High Speed Research Flexible Semispan Model

This report presents the work performed by Lockheed Martin Engineering and Sciences (LMES) in support of the High Speed Research (HSR) Flexible Semispan Model (FSM) wind-tunnel test. The test was conducted in order to assess the aerodynamic and aeroelastic character of a flexible high speed civil transport wing. Data was acquired for the purpose of code validation and trend evaluation for this type of wing. The report describes a number of activities in preparing for and conducting the wind-tunnel test. These included coordination of the design and fabrication, development of analytical models, analysis/hardware correlation, performance of laboratory tests, monitoring of model safety issues, and wind-tunnel data acquisition and reduction. Descriptions and relevant evaluations associated with the pretest data are given in sections 1 through 6, followed by pre- and post-test flutter analysis in section 7, and the results of the aerodynamics/loads test in section 8. Finally, section 9 provides some recommendations based on lessons learned throughout the FSM program.

Schuster, David M.↗
NASA NTRS2013

The Effectiveness of Magnetic Shielding in High-Isp Hall Thrusters

A series of numerical simulations and experiments have been performed to assess the effectiveness of magnetic shielding in a Hall thruster operating in the discharge voltage range of 300-700 V (Isp 2000-2700 s) at 6 kW, and 800 V (Isp 3000) at 9 kW. In this paper we report on the simulation results and their validation with experimental measurements. At 6 kW the magnetic field topology with which we recently demonstrated highly effective magnetic shielding at 300 V was retained for all other discharge voltages; only the magnitude of the field was changed to achieve optimum thruster performance. It is found that magnetic shielding remains highly effective for all discharge voltages studied. Maximum erosion rates that remain fairly constant across the range of 300-700 V are computed, with values not exceeding 10-2 mm/kh. Such rates are ~3 orders of magnitude less than those observed in the unshielded version of the same thruster at 300 V. At 9 kW and 800 V, saturation of the magnetic circuit did not permit us to attain precisely the same magnetic shielding topology as that employed during the 6-kW operation since this thruster was not designed to operate at this condition. Consequently, the maximum erosion rate at the inner wall is found to be ~1 order of magnitude higher (~10-1 mm/kh) than that at the 6-kW level. At the outer wall the ion energy is below the sputtering yield threshold so no measurable erosion is expected. The reasons behind the effectiveness of magnetic shielding at higher discharge voltages are discussed.

magnetic-field-aligned-mesh (MFAM) code↗
NASA NTRS2018

Three-Dimensional, Space-Dependent Mesoscale Diffusivity: Derivation and Implications

Recently, we presented a parameterization of an arbitrary tracer 3D mesoscale flux that describes both diabatic and adiabatic regimes without using arbitrary tapering functions. However, we did not parameterize the mesoscale diffusivity, which is the subject of this work. A key difference between the present and previous diffusivity parameterizations is that in the latter, the two main ingredients, mesoscale drift velocity and eddy kinetic energy, were not parameterized but determined using present data, which deprives the models of predictive power. Since winds, stratification, etc., are predicted to change in the future, use of these parameterizations to study future climate scenarios becomes questionable. In this work, we parameterize drift velocity and eddy kinetic energy (vertical–horizontal components), which we first assess with data [WOCE, TOPEX/Poseidon (T/P), and North Atlantic Tracer Release Experiment (NATRE)] and then use in a coarse-resolution stand-alone ocean code under Coordinated Ocean-Ice Reference Experiment I (CORE-I) forcing. We present results for the global ocean temperature and salinity, Atlantic overturning circulation, meridional heat transport, and Drake Passage transport, which we compare with several previous studies. The temperature drift is less than that of five of seven previous OGCMs, and the salinity drift is among the smallest in those studies. The predicted winter Antarctic Circumpolar Current mixed layer depths (MLDs) are in good agreement with the data. Predicting the correct MLD is important in climate studies since models that predict very deep mixed layers transfer more of the radiative perturbation to the deep ocean, reducing surface warming (and vice versa).

Eddies; Mesoscale processes; Mixing; Ocean models;↗
DOE OSTI2022

Morphology of shocked lateral outflows in colliding hydrodynamic flows

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. Here, we introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design and validation of a gamma-ray scanning system for measuring irradiated nuclear fuel

We report a prototype gamma-ray scanning system was designed to perform post-irradiation measurements of nuclear fuel at Idaho National Laboratory (INL). The system is composed of a coaxial high-purity germanium (HPGe) detector, a collimator, and mechanical positioning stages that translate the fuel sample across the front of the collimator. This paper documents the design process, modeling and the laboratory characterization performed to validate the system. The modeling of the system was done using the Monte Carlo N-Particle code with vendor-supplied design specifications, dimensional measurements, and x-ray radiographs of the HPGe detector. The characterization included measurements of calibrated radiation sources in different geometries, and calculations for absolute efficiency and dead time. Benchmark measurements were performed scanning an irradiated fuel rodlet containing ten pellets of 4.9% enriched UO2 in zirconium alloy. These results were compared to simulated spectra to further help characterize the detector model. One goal of the system is to determine the number of fissions per gram of UO2 in the fuel; comparison of the results from this prototype system's assay of the irradiated fuel rodlet are in agreement with previously estimated results.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

AutoSourceID-Classifier: Star-galaxy classification using a convolutional neural network with spatial information

Aims.Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification’s reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. Methods.The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results. Results.We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C’s direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.

Astronomy & Astrophysics↗