Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flux reconstruction”

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 181 records · Page 10

Tagging Cosmic Ray Background in the Icarus Detector

ICARUS T-600 will be the Far Detector of the upcoming SBN program, for which it will leverage a ∼ 1ns response and a high resolution for neutrino event reconstruction. ICARUS will respond to the challenge of the cosmic flux at the surface thanks to its Cosmic Ray Tagging (CRT) system. The work conducted during the Summer Students at Fermilab program contributed to the calibration activities that are ongoing in view of the second physics run in late 2022: the timing of CRT and PMT signals was analysed, highlighting the need for corrections to the PMTs. The timing delay for the global trigger signal for a section of the side CRT. The possibility of calibrating the light yield of the side CRTs with run data was finally explored.

47 OTHER INSTRUMENTATION↗

The Ultraviolet Radiation Environment around M Dwarf Exoplanet Host Stars

The spectral and temporal behavior of exoplanet host stars is a critical input to models of the chemistry and evolution of planetary atmospheres. Ultraviolet photons influence the atmospheric temperature profiles and production of potential biomarkers on Earth-like planets around these stars. At present, little observational or theoretical basis exists for understanding the ultraviolet spectra of M dwarfs, despite their critical importance to predicting and interpreting the spectra of potentially habitable planets as they are obtained in the coming decades. Using observations from the Hubble Space Telescope, we present a study of the UV radiation fields around nearby M dwarf planet hosts that covers both far-UV (FUV) and near-UV (NUV) wavelengths. The combined FUV+NUV spectra are publicly available in machine-readable format. We find that all six exoplanet host stars in our sample (GJ 581, GJ 876, GJ 436, GJ 832, GJ 667C, and GJ 1214) exhibit some level of chromospheric and transition region UV emission. No "UV-quiet" M dwarfs are observed. The bright stellar Lyman-alpha emission lines are reconstructed, and we find that the Lyman-alpha line fluxes comprise approximately 37%-75% of the total 1150-3100 A flux from most M dwarfs; approximately greater than 10(exp3) times the solar value. We develop an empirical scaling relation between Lyman-alpha and Mg II emission, to be used when interstellar H I attenuation precludes the direct observation of Lyman-alpha. The intrinsic unreddened flux ratio is F(Lyman-alpha)/F(Mg II) = 10(exp3). The F(FUV)/F(NUV) flux ratio, a driver for abiotic production of the suggested biomarkers O2 and O3, is shown to be approximately 0.5-3 for all M dwarfs in our sample, greather than 10(exp3) times the solar ratio. For the four stars with moderate signal-to-noise Cosmic Origins Spectrograph time-resolved spectra, we find UV emission line variability with amplitudes of 50%.500% on 10(exp2)-10(exp3) s timescales. This effect should be taken into account in future UV transiting planet studies, including searches for O3 on Earth-like planets. Finally, we observe relatively bright H2 fluorescent emission from four of the M dwarf exoplanetary systems (GJ 581, GJ 876, GJ 436, and GJ 832). Additional modeling work is needed to differentiate between a stellar photospheric or possible exoplanetary origin for the hot (T(H2) approximately equal to 2000-4000 K) molecular gas observed in these objects.

France, Kevin↗

Structural anisotropy in Sb thin films

Sb thin films have attracted wide interest due to their tunable band structure, topological phases, high electron mobility, and thermoelectric properties. We successfully grow epitaxial Sb thin films on a closely lattice-matched GaSb(001) surface by molecular beam epitaxy. We find a novel anisotropic directional dependence on their structural, morphological, and electronic properties. The origin of the anisotropic features is elucidated using first-principles density functional theory (DFT) calculations. The growth regime of crystalline and amorphous Sb thin films was determined by mapping the surface reconstruction phase diagram of the GaSb(001) surface under Sb 2 flux, with confirmation of structural characterizations. Crystalline Sb thin films show a rhombohedral crystal structure along the rhombohedral (211) surface orientation parallel to the cubic (001) surface orientation of the GaSb substrate. At this coherent interface, Sb atoms are aligned with the GaSb lattice along the [1̄10] crystallographic direction but are not aligned well along the [110] crystallographic direction, which results in anisotropic features in reflection of high-energy electron diffraction patterns, misfit dislocation formation, surface morphology, and transport properties. Our DFT calculations show that the preferential orientation of the rhombohedral Sb (211) plane may originate from the GaSb surface, where Sb atoms align with the Ga and Sb atoms on the reconstructed surface. The formation energy calculations confirm the stability of the experimentally observed structures. Our results provide optimal film growth conditions for further studies of novel properties of Bi 1-x Sb x thin films with similar lattice parameters and an identical crystal structure, as well as functional heterostructures of them with III–V semiconductor layers along the (001) surface orientation, supported by a theoretical understanding of the anisotropic film orientation.

36 MATERIALS SCIENCE↗

Application of the metabolic modeling pipeline in KBase to categorize reactions, predict essential genes, and predict pathways in an isolate genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

DOE knowledgebase↗

Solution-Adaptive Cartesian Cell Approach for Viscous and Inviscid Flows

A Cartesian cell-based approach for adaptively refined solutions of the Euler and Navier-Stokes equations in two dimensions is presented. Grids about geometrically complicated bodies are generated automatically, by the recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal cut cells are created using modified polygon-clipping algorithms. The grid is stored in a binary tree data structure that provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite volume formulation. The convective terms are upwinded: A linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The results of a study comparing the accuracy and positivity of two classes of cell-centered, viscous gradient reconstruction procedures is briefly summarized. Adaptively refined solutions of the Navier-Stokes equations are shown using the more robust of these gradient reconstruction procedures, where the results computed by the Cartesian approach are compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

LOFTID Surface Heating Reconstruction

On November 10, 2022, the Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) reentry vehicle launched to low-Earth orbit aboard a United Launch Alliance Atlas V rocket out of Vandenberg Air Force Base. The aeroshell, the largest Hypersonic Inflatable Aerodynamic Decelerator (HIAD) ever flown, was inflated to its full 6-meter diameter before successful re-entry into the atmosphere. The aeroshell was heavily instrumented in order to understand its behavior during entry. There were 82 thermocouples (TCs) distributed across the aeroshell, with 22 integrated into the flexible thermal protection system (FTPS) on the rigid nose, 36 in the FTPS on the deployable structure, and 24 on the inflatable structure. TCs were placed at different depths throughout the FTPS. Those nearest to the surface were located just beneath the two SiC outer fabric layers. The near-surface TCs on the rigid nose were Type R with flame spray alumina insulation, while those on the flank were Type N with mica/ceramic insulation. Additionally, a radiometer was placed at the center of the nose surrounded by four total heat flux gauges in a cruciform configuration at a radius of 0.41 m. The nose instrumentation is shown in Fig. 1 and a cross-section of the aeroshell with all TC locations is shown in Fig. 2. The objective of this work was to use the temperatures measured by the TCs during flight to estimate the surface heat rate across the aeroshell throughout the period of re-entry by inverse analysis methodology. The results were used to evaluate the fidelity of measurements from the total heat flux gauges on the nose, determine the surface heat flux at aeroshell locations where gauges were not present, and compare to pre-flight CFD-based heating predictions. Inversely estimated surface heat flux continues to be used to correlate FTPS thermal models to reconstruct in-flight thermal response.

H S Alpert↗

LOFTID Surface Heating Reconstruction

On November 10, 2022, the Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) reentry vehicle launched to low-Earth orbit aboard a United Launch Alliance Atlas V rocket out of Vandenberg Air Force Base. The aeroshell, the largest Hypersonic Inflatable Aerodynamic Decelerator (HIAD) ever flown, was inflated to its full 6-meter diameter before successful re-entry into the atmosphere. The aeroshell was heavily instrumented in order to understand its behavior during entry. There were 82 thermocouples (TCs) distributed across the aeroshell, with 22 integrated into the flexible thermal protection system (FTPS) on the rigid nose, 36 in the FTPS on the deployable structure, and 24 on the inflatable structure. TCs were placed at different depths throughout the FTPS. Those nearest to the surface were located just beneath the two SiC outer fabric layers. The near-surface TCs on the rigid nose were Type R with flame spray alumina insulation, while those on the flank were Type N with mica/ceramic insulation. Additionally, a radiometer was placed at the center of the nose surrounded by four total heat flux gauges in a cruciform configuration at a radius of 0.41 m. The nose instrumentation is shown in Fig. 1 and a cross-section of the aeroshell with all TC locations is shown in Fig. 2. The objective of this work was to use the temperatures measured by the TCs during flight to estimate the surface heat rate across the aeroshell throughout the period of re-entry by inverse analysis methodology. The results were used to evaluate the fidelity of measurements from the total heat flux gauges on the nose, determine the surface heat flux at aeroshell locations where gauges were not present, and compare to pre-flight CFD-based heating predictions. Inversely estimated surface heat flux continues to be used to correlate FTPS thermal models to reconstruct in-flight thermal response.

LOFTID↗

Comprehensive Measurement of the Reactor Antineutrino Spectrum and Flux at Daya Bay

This Letter reports the precise measurement of the reactor antineutrino spectrum and flux based on the full dataset of 4.7 ×10 6 inverse-beta-decay candidates collected at Daya Bay near detectors. Expressed in terms of the inverse-beta-decay yield per fission, the antineutrino spectra from all reactor fissile isotopes and the specific 235 U and 239 Pu isotopes are measured with 1.3%, 3%, and 8% uncertainties, respectively, near the 3 MeV spectrum peak in reconstructed energy, reaching the best precision in the world. The total antineutrino flux and isotopic 235 U and 239 Pu fluxes are precisely measured to be 5.84 ± 0.07, 6.16 ± 0.12, and 4.16 ± 0.21 in units of 10 −43 cm 2 /fission. These measurements are compared with the Huber-Mueller model, the reevaluated conversion model based on the Kurchatov Institute measurement, and the latest summation model (SM2023). Furthermore, the Daya Bay flux shows good consistency with the Kurchatov Institute and SM2023 models but disagrees with the Huber-Mueller model. The Daya Bay spectrum, however, disagrees with all model predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Calibration and Timing Performance of the Light Detection System in the ICARUS Detector

ICARUS is the largest Liquid Argon Time Projection Chamber (LArTPC) in operation and serves as the Far Detector of the Short Baseline Neutrino (SBN) program at Fermilab. It aims to investigate the possible existence of sterile neutrinos with $\Delta m^2 \approx \SI{1}{eV^2}$ using the Booster Neutrino Beam (BNB) and explore physics beyond the Standard Model with the Neutrinos at the Main Injector (NuMI) beam. The ICARUS light detection system, comprising 360 TPB-coated large-area Photo-Multiplier Tubes (PMTs), is crucial for triggering and event reconstruction. Due to its shallow installation, the detector is exposed to a high flux of cosmic rays, necessitating precise timing to reject background events and align neutrino interactions with the beam time profile. This talk will detail the timing inter-calibration procedures for the ICARUS light detection system, which achieve sub-nanosecond resolution. Additionally, the performance of the system in reconstructing the timing of neutrino interactions from the BNB and NuMI beams will be discussed. The results highlight the effectiveness of the ICARUS light detection system in enhancing the detector's capability for precise and reliable neutrino selection.

43 PARTICLE ACCELERATORS↗

Calibration and Timing Performance of the Light Detection System in the ICARUS Detector

ICARUS is the largest Liquid Argon Time Projection Chamber (LArTPC) in operation and serves as the Far Detector of the Short Baseline Neutrino (SBN) program at Fermilab. It aims to investigate the possible existence of sterile neutrinos with $\Delta m^2 \approx 1 \,\mathrm{eV^2}$ using the Booster Neutrino Beam (BNB) and explore physics beyond the Standard Model with the Neutrinos at the Main Injector (NuMI) beam. The ICARUS light detection system, comprising 360 TPB-coated large-area Photo-Multiplier Tubes (PMTs), is crucial for triggering and event reconstruction. Due to its shallow installation, the detector is exposed to a high flux of cosmic rays, necessitating precise timing to reject background events and align neutrino interactions with the beam time profile. This talk will detail the timing inter-calibration procedures for the ICARUS light detection system, which achieve sub-nanosecond resolution. Additionally, the performance of the system in reconstructing the timing of neutrino interactions from the BNB and NuMI beams will be discussed. The results highlight the effectiveness of the ICARUS light detection system in enhancing the detector's capability for precise and reliable neutrino selection.

43 PARTICLE ACCELERATORS↗

Reconciling LCROSS and Orbital Neutron Water Abundance Estimates in Cabeus Crater

The Lunar Prospector Neutron Spectrometer (LPNS) first revealed Cabeus crater (84.9 deg S, 35.5degW) as having the highest inferred hydrogen on the Moon. Because of the broad LPNS footprint (approximately 40 km FWHM), the apparent peak water-equivalent hydrogen (WEH) concentration is only approximately 0.25 wt%, but could be much higher in smaller areas than the spectrometer footprint. Earlier image reconstruction work suggested that areas within permanent shadow have abundances approximately 1 wt% WEH. However, the LCROSS impact yielded total water estimates, ice plus vapor, of between 3 and 10 wt%. The large disagreement between LCROSS and apparent orbital values imply that either the ice is buried, by perhaps as much as 50 to 100 cm; or the ice distribution within Cabeus is spatially inhomogeneous, or both. Modeling reveals that the areal extent of a "shallow permafrost zone" is far greater than the area of permanent shadow. Ice can be virtually stable for billions of years within a few tens of centimeters of the surface in these areas. However, the LCROSS impact took place in an area of permanent shadow. If stably-trapped volatiles can be found in locales that receive occasional, oblique sunlight, landed missions may target these sites and eventual resource exploitation may be done more easily. Are orbital neutron data consistent with areally-extensive, volatile-rich cold traps? Orbital epithermal neutron data over the northern half of Cabeus (near the LCROSS impact site) are consistent with 0.2 wt% WEH or less in the "permafrost zone" near the crater. On the other hand, pixon reconstructions that confine the hydrogen enhancements to permanent shadow result in higher abundance estimates -- around 1 wt% if homogeneously mixed. But if the PSR abundance is increased to 10 wt%, consistent with the sum of all H-bearing compounds seen by LCROSS, a much larger-than-observed reduction in neutron count rate would be seen from orbit. It is likely that volatiles are inhomogeneously distributed, due to both impact processes and emplacement history. Two possibilities may bring consistency to the orbital and LCROSS measurements. Inhomogeneous lateral distribution: Consider the extreme case of a bimodal distribution within the crater -- dry and wet. In this case the epithermal leakage flux seen from orbit is a mixture of two different values, weighted according to fractional areas. Two possible outcomes, depending on whether the inferred leakage flux for the PSR or "permafrost" areas are considered. In the first case, approximately 40% of the PSR may be "wet", the remainder dry (and LCROSS was slightly lucky). However, if the whole area of permafrost is considered, then as little as 20% of the area will be as "wet" as the LCROSS results (and LCROSS was quite lucky). Inhomogeneous depth distribution: The leakage flux of thermal and epithermal neutrons depends on depth of burial of an icy layer beneath dry ferroan anorthosite soil (FAn). For the Cabeus PSR, the pixon reconstruction values for the epithermal flux allows a range of abundance and burial depth, while that of the thermal+epi detector constrains this range. (Uncertainties in iron abundance in the FAn can have significant impact on thermal neutron leakage flux estimates.) Between 20% and 40% of the Cabeus floor may be "wet", or alternatively a 5-10 wt% "wet" layer exists between 50 and 100 cm beneath a layer of dry regolith within the PSR. But volatile abundances of 5 wt% or more, distributed uniformly and homogeneously throughout the Cabeus PSR do not agree with orbital measurements

Elphic, Richard↗

A Stochastic Approach to Reconstruct Gamma-Ray-burst Light Curves

Gamma-ray bursts (GRBs), as they are observed at high redshift ( z = 9.4), are vital to cosmological studies and investigating Population III stars. To tackle these studies, we need correlations among relevant GRB variables with the requirement of small uncertainties on their variables. Thus, we must have good coverage of GRB light curves (LCs). However, gaps in the LC hinder the precise determination of GRB properties and are often unavoidable. Therefore, extensive categorization of GRB LCs remains a hurdle. We address LC gaps using a stochastic reconstruction, wherein we fit two preexisting models (the Willingale model; W07; and a broken power law; BPL) to the observed LC, then use the distribution of flux residuals from the original data to generate data to fill in the temporal gaps. We also demonstrate a model-independent LC reconstruction via Gaussian processes. At 10% noise, the uncertainty of the end time of the plateau, its correspondent flux, and the temporal decay index after the plateau decreases by 33.3%, 35.03%, and 43.32% on average for the W07, and by 33.3%, 30.78%, 43.9% for the BPL, respectively. The uncertainty of the slope of the plateau decreases by 14.76% in the BPL. After using the Gaussian process technique, we see similar trends of a decrease in uncertainty for all model parameters for both the W07 and BPL models. These improvements are essential for the application of GRBs as standard candles in cosmology, for the investigation of theoretical models, and for inferring the redshift of GRBs with future machine-learning analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine learning for single-ended event reconstruction in PROSPECT experiment

The Precision Reactor Oscillation and Spectrum Experiment, PROSPECT, was a segmented antineutrino detector that successfully operated at the High Flux Isotope Reactor in Oak Ridge, TN, during its 2018 run. Despite challenges with photomultiplier tube base failures affecting some segments, innovative machine learning approaches were employed to perform position and energy reconstruction, and particle classification. This work highlights the effectiveness of convolutional neural networks and graph convolutional networks in enhancing data analysis. By leveraging these techniques, a 3.3% increase in effective statistics was achieved compared to traditional methods, showcasing their potential to improve analysis performance. Furthermore, these machine learning methodologies offer promising applications for other segmented particle detectors, underscoring their versatility and impact.

47 OTHER INSTRUMENTATION↗

High Flux Isotope Reactor Neutron Spectrum Shape Estimation From Activation Experiment Data

Here, this article provides a comprehensive review of historical irradiation dosimetry available for different locations within the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). This article includes a summary of the available measured activation target data covering a span of over 15 years and 39 experimental campaigns, including 200 individual sample locations evaluated. Using this broad set of data, we reconstruct historic average neutron spectra shapes for HFIR at various locations, including within the flux trap region, beryllium reflectors, and hydraulic tube (HT) regions, at both the 100- and 85-MW operational power. Our findings indicate that the general axial flux distribution shows a relatively small change in transition from 100- to 85-MW operating power, with differences of -6% to +15% for the thermal energy range and around -16% to +8% for the fast range, indicating a sharper drop-off of the thermal neutron flux away from the axial center. Compared with historical dosimetry estimates of the HFIR flux shape, we generally find a moderately broader axial profile shape for thermal neutrons in the interior target regions for the 100-MW samples evaluated but relatively close agreement for the present 85-MW flux shape for both thermal and fast fluxes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Source shape estimation for neutron imaging systems using convolutional neural networks

Neutron imaging systems are important diagnostic tools for characterizing the physics of inertial confinement fusion reactions at the National Ignition Facility (NIF). In particular, neutron images give diagnostic information on the size, symmetry, and shape of the fusion hot spot and surrounding cold fuel. Images are formed via collection of neutron flux from the source using a system of aperture arrays and scintillator-based detectors. Currently, reconstruction of fusion source geometry from the collected neutron images is accomplished by solving a computationally intensive maximum likelihood estimation problem via expectation maximization. In contrast, it is often useful to have simple representations of the overall source geometry that can be computed quickly. In this work, we develop convolutional neural networks (CNNs) to reconstruct the outer contours of simple source geometries. We compare the performance of the CNN for penumbral and pinhole data and provide experimental demonstrations of our methods on both non-noisy and noisy data.

Machine learning, neutron imaging, source reconstr↗

Diverse electron carriers drive syntrophic interactions in an enriched anaerobic acetate-oxidizing consortium

Abstract In many anoxic environments, syntrophic acetate oxidation (SAO) is a key pathway mediating the conversion of acetate into methane through obligate cross-feeding interactions between SAO bacteria (SAOB) and methanogenic archaea. The SAO pathway is particularly important in engineered environments such as anaerobic digestion (AD) systems operating at thermophilic temperatures and/or with high ammonia. Despite the widespread importance of SAOB to the stability of the AD process, little is known about their in situ physiologies due to typically low biomass yields and resistance to isolation. Here, we performed a long-term (300-day) continuous enrichment of a thermophilic (55 °C) SAO community from a municipal AD system using acetate as the sole carbon source. Over 80% of the enriched bioreactor metagenome belonged to a three-member consortium, including an acetate-oxidizing bacterium affiliated with DTU068 encoding for carbon dioxide, hydrogen, and formate production, along with two methanogenic archaea affiliated with Methanothermobacter_A. Stable isotope probing was coupled with metaproteogenomics to quantify carbon flux into each community member during acetate conversion and inform metabolic reconstruction and genome-scale modeling. This effort revealed that the two Methanothermobacter_A species differed in their preferred electron donors, with one possessing the ability to grow on formate and the other only consuming hydrogen. A thermodynamic analysis suggested that the presence of the formate-consuming methanogen broadened the environmental conditions where ATP production from SAO was favorable. Collectively, these results highlight how flexibility in electron partitioning during SAO likely governs community structure and fitness through thermodynamic-driven mutualism, shedding valuable insights into the metabolic underpinnings of this key functional group within methanogenic ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗

Tetrahedral Finite-Volume Solutions to the Navier-Stokes Equations on Complex Configurations

A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.

Frink, Neal T.↗

Tetrahedral Finite-Volume Solutions to the Navier-Stokes Equations on Complex Configurations

A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.

Frink, N. T.↗