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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spatially calibrating polycyclic aromatic hydrocarbons (PAHs) as proxies of area burned by vegetation fires: Insights from comparisons of historical data and sedimentary PAH fluxes

Many regions worldwide have experienced increasing wildfire activity in recent years and climate changes are predicted to result in more frequent and severe fires. Reconstruction of past fire activity offers paleoenvironmental context for modern and future burning. Pyrogenic polycyclic aromatic hydrocarbons (PAHs) have been increasingly used as a molecular biomarker for fire occurrence in the paleorecord and offer opportunity for nuanced reconstructions of fire characteristics. A suite of PAHs are produced during combustion, and the emission amount and assemblage is influenced by many variables including fuel type, fire temperature, and oxygen availability. Despite recent advances in understanding the controls and taphonomy of these biomass burning markers, the spatial scale of this proxy is unknown. In this paper, measurements of PAH fluxes preserved in a lake sediment archive from the Sierra Nevada, California were compared with a historical geographic information system dataset of area burned up to 150 km distance from the lake to determine the spatial scales for which these biomarkers are reliable proxies of burning. Comparisons of PAH fluxes with charcoal accumulation rates in the same sediments suggest that pyrogenic particulate transport modulates low to mid-molecular weight PAHs via adsorption. Overall, the results indicate that PAH records integrate a combination of spatial signals of area burned and measurement of individual PAHs may enable cross-scale paleofire reconstructions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2023

Small-Scale Magnetic Flux Ropes in Stream Interaction Regions From Parker Solar Probe and Wind Spacecraft Observations

Using in situ measurements from the Parker Solar Probe and Wind spacecraft, we investigate the small-scale magnetic flux ropes (SFRs) and their properties inside stream interaction regions (SIRs). Within SIRs from ∼0.15 to 1 au, SFRs are found to exist in a wide range of solar wind speeds with more frequent occurrences after the stream interface, and the Alfvénicity of these structures decreases significantly with increasing heliocentric distances. Furthermore, we examine the variation of five corresponding SIRs from the same solar sources. The enhancements of suprathermal electrons within these SIRs persist at 1 au and are observed multiple times. An SFR appears to occur repeatedly with the recurring SIRs and is traversed by the Wind spacecraft at least twice. This set of SFRs has similarities in variations of the magnetic field components, plasma bulk properties, density ratio of solar wind alpha and proton particles, and unidirectional suprathermal electrons. We also show, through the detailed time-series plots and Grad–Shafranov reconstruction results, that they possess the same chirality and carry comparable amounts of magnetic flux. Lastly, we discuss the possibility for these recurring SFRs to be formed via interchange reconnection, maintain the connection with the Sun, and survive up to 1 au.

Solar wind↗
NASA NTRS2001

A High-Order Finite Spectral Volume Method for Conservation Laws on Unstructured Grids

A time accurate, high-order, conservative, yet efficient method named Finite Spectral Volume (FSV) is developed for conservation laws on unstructured grids. The concept of a 'spectral volume' is introduced to achieve high-order accuracy in an efficient manner similar to spectral element and multi-domain spectral methods. In addition, each spectral volume is further sub-divided into control volumes (CVs), and cell-averaged data from these control volumes is used to reconstruct a high-order approximation in the spectral volume. Riemann solvers are used to compute the fluxes at spectral volume boundaries. Then cell-averaged state variables in the control volumes are updated independently. Furthermore, TVD (Total Variation Diminishing) and TVB (Total Variation Bounded) limiters are introduced in the FSV method to remove/reduce spurious oscillations near discontinuities. A very desirable feature of the FSV method is that the reconstruction is carried out only once, and analytically, and is the same for all cells of the same type, and that the reconstruction stencil is always non-singular, in contrast to the memory and CPU-intensive reconstruction in a high-order finite volume (FV) method. Discussions are made concerning why the FSV method is significantly more efficient than high-order finite volume and the Discontinuous Galerkin (DG) methods. Fundamental properties of the FSV method are studied and high-order accuracy is demonstrated for several model problems with and without discontinuities.

Wang, Z. J.↗
NASA NTRS1987

High resolution upwind schemes for the three-dimensional incompressible Navier-Stokes equations

Based on flux-difference splitting, implicit high resolution schemes are constructed for efficient computations of steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. These schemes use first-order accurate Euler backward-time differencing and second-order central differencing for the viscous shear fluxes. Up to third-order accurate upwind differencing is achieved through a reconstruction of the solution from its cell averages. The reconstruction is accomplished by linear interpolation, where the node stencils are selected such that in regions of smooth solution the flow is highly resolved while spurious oscillations in regions of rapid changes in gradient are still suppressed. Fairly rapid convergence to steady-state solutions is attained with a completely vectorizable hybrid time-marching method. Flows around a sharp-edged delta wing are computed with the maximum accuracy of the upwind-differencing restricted to first-, second-, and third-order, to illustrate the effect of accuracy on the global and on the local vortical flow fields. The results are validated with experimental data.

Hartwich, PETER-M.↗
NASA NTRS1990

Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction

High order accurate finite-volume schemes for solving the Euler equations of gasdynamics are developed. Central to the development of these methods are the construction of a k-exact reconstruction operator given cell-averaged quantities and the use of high order flux quadrature formulas. General polygonal control volumes (with curved boundary edges) are considered. The formulations presented make no explicit assumption as to complexity or convexity of control volumes. Numerical examples are presented for Ringleb flow to validate the methodology.

Barth, Timothy J.↗
NASA NTRS1991

A fast upwind solver for the Euler equations on three-dimensional unstructured meshes

An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a novel cell reconstruction process which results in computational times per cell comparable to those of structured codes. The approach yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicit limiting. Solutions are advanced in time by a 3-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Solutions are presented for a range of configurations in the transonic speed regime to demonstrate code accuracy, speed, and robustness. The results include an assessment of grid sensitivity and convergence acceleration by mesh sequencing.

Frink, Neal T.↗
NASA NTRS1992

Upwind scheme for solving the Euler equations on unstructured tetrahedral meshes

An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a multidimensional linear reconstruction process. The solution gradients required for the higher-order differenes are computed by a novel approach that yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicitly applying a limiter. Solutions are advanced in time by a three-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Transonic solutions are presented for two meshes around the ONERA M6 wing and demonstrate substantial accuracy and insensitivity to mesh size.

Frink, Neal T.↗
NASA NTRS2004

A Fast Upwind Solver for the Euler Equations on Three-Dimensional Unstructured Meshes

An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a novel cell reconstruction process which results in computational times per cell comparable to those of structured codes. The approach yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicit limiting. Solutions are advanced in time by a 3-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Solutions are presented for a range of configurations in the transonic speed regime to demonstrate code accuracy, speed, and robustness. The results include an assessment of grid sensitivity and convergence acceleration by mesh sequencing.

Frink, Neal T.↗
NASA NTRS2008

Modelling a Simultaneous Radio/X-Ray Flare from Cyg X-1

The long-term monitoring campaign of Cyg X-1 has provided the detection of the first simultaneous radio/X-ray flare seen from that source. We investigate the physical characteristics of the event in terms of emission from a homogeneous, expanding blob of pair-plasma, superimposed on a baseline flux in both bands. We find that while the radio flare can be reconstructed under various configurations of a cooling blob, continuous (re)acceleration of particles inside the jet is necessary to sustain X-ray emission at the levels implied by the data, for the observed duration. We present major results of the modelling and discuss implications for the role of microquasar jets.

Leventis, Konstantinos↗
NASA NTRS2021

Overview of Dragonfly Entry Aerosciences Measurements (DrEAM)

NASA Ames Research Center (ARC) leads the Dragonfly Entry Aerosciences Measurements (DrEAM) project, which is an aeroshell instrumentation suite on the Dragonfly mission that fulfills the Engineering Science Investigation requirement for the New Frontiers mission. NASA ARC is partnering with NASA Langley Research Center (LaRC) and the German Aerospace Center (DLR) to provide a comprehensive sensor suite, including a DLR-provided Data Acquisition System (DAS). DrEAM will provide key aerothermodynamic data and performance analysis for Dragonfly’s forebody and backshell Thermal Protection System (TPS). Titan’s atmosphere predominantly consists of ni-trogen (~98% by mole) with small amounts of me-thane (~2% by mole) and other trace gases. CN is a strong radiator, and is found in nonequilibrium con-centrations for Titan entry. The accurate modeling of nonequilibrium CN radiation has proven to be a diffi-cult task. Prompted by the Huygens mission, many experimental campaigns and analyses were performed to better understand the aerothermal environments experienced by the probe during Titan entry [1]. How-ever, the Huygens probe carried no heatshield instru-mentation. Therefore, the DrEAM sensor suite will sig-nificantly advance the state-of-the-art not only by documenting the environment and performance of Dragonfly’s entry system but also by making key in situ measurements in Titan’s atmosphere for the first time. Aerothermal environments and TPS response will be measured using sensors whose flight heritage is tak-en from the Mars Entry, Descent, and Landing In-strumentation 2 (MEDLI2) thermocouple plugs and the COMbined Aerothermal and Radiometer Sensor (COMARS) suite [2], with the latter supplied by DLR. The MEDLI2 project used embedded thermocouples to directly measure the in-depth TPS temperature-time history at several locations on the heat shield and backshell of the Mars 2020 entry vehicle. These tem-perature measurements, in turn, can be used to infer surface environments via an inverse analysis proce-dure analogous to that used for MEDLI. For DrEAM, the thermocouple plug subsystem will be known as Dragonfly Sensors for Aero-Thermal Reconstruction (DragSTR). On Schiaparelli, the COMARS suite in-cluded three total surface-mounted heat flux sensors, three pressure sensors, and one radiometer. For DrEAM, the COMARS package will be known as the COmbined Sensor System for Titan Atmosphere (COSSTA). Since the methane concentration in the Titan atmos-phere is directly proportional to the radiative heat flux, the COSSTA measurements will be used to reduce the current uncertainty in the methane volume fraction. Atmospheric density measurements and capsule aero-dynamic data will be obtained through the onboard Inertial Measurement Unit (IMU), supplemented by pressure transducers similar to those used by the MEDLI and MEDLI2 projects. The DrEAM pressure sensors will be known as the Dragonfly Atmospheric Flight Transducers. (DrAFT). The pressure measure-ments, when combined with data from the on-board IMU, will allow for reconstruction of such quantities as vehicle Mach number, freestream density, and atmos-pheric winds. DrAFT measurements will enhance Dragonfly trajectory reconstruction and enable a sepa-ration of the aerodynamics from the atmosphere, as was done for MEDLI [3] and is currently in process for the MEDLI2 flight data set.

EDL↗
NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗
NASA NTRSDate not supplied

Mars Entry, Descent, and Landing Instrumentation 2 Trajectory, Aerodynamics, and Atmosphere Reconstruction

On February 18th, 2021, the Mars 2020 entry system successfully delivered the Perseverance rover to the surface of Mars at Jezero Crater. The entry capsule carried a set of instrumentation installed on the heat shield and backshell, named the Mars Entry, Descent, and Landing Instrumentation 2. The instruments include pressure transducers, thermocouples, heat flux gauges, and radiometers to measure the aerodynamic and aerothermodynamic performance of the entry vehicle. This paper describes the trajectory and atmosphere reconstruction results based on the pressure sensor measurements. The process uses a Kalman filter approach to estimate the freestream atmospheric properties from the pressure measurements combined with a model of the pressure distribution of the heatshield and other sensor inputs, including an inertial measurement unit and other on-board navigation sensors, and several external atmospheric observations. The results indicate upper altitude density was up to 150% higher than nominal, which is consistent with the observed early entry guidance start time. The density below 40 km was within 12% the pre-flight predictions. The reconstructed axial force coefficient was approximately 2% lower than the pre-flight prediction across the flight range.

Christopher D Karlgaard↗
DOE OSTI2023

A High-Quality Genome-Scale Model for Rhodococcus opacus Metabolism

Rhodococcus opacus is a bacterium that has a high tolerance to aromatic compounds and can produce significant amounts of triacylglycerol (TAG). Here, we present iGR1773, the first genome-scale model (GSM) of R. opacus PD630 metabolism based on its genomic sequence and associated data. The model includes 1773 genes, 3025 reactions, and 1956 metabolites, was developed in a reproducible manner using CarveMe, and was evaluated through Metabolic Model tests (MEMOTE). We combine the model with two Constraint-Based Reconstruction and Analysis (COBRA) methods that use transcriptomics data to predict growth rates and fluxes: E-Flux2 and SPOT (Simplified Pearson Correlation with Transcriptomic data). Growth rates are best predicted by E-Flux2. Flux profiles are more accurately predicted by E-Flux2 than flux balance analysis (FBA) and parsimonious FBA (pFBA), when compared to 44 central carbon fluxes measured by 13C-Metabolic Flux Analysis (13C-MFA). Under glucose-fed conditions, E-Flux2 presents an R2 value of 0.54, while predictions based on pFBA had an inferior R2 of 0.28. We attribute this improved performance to the extra activity information provided by the transcriptomics data. For phenol-fed metabolism, in which the substrate first enters the TCA cycle, E-Flux2’s flux predictions display a high R2 of 0.96 while pFBA showed an R2 of 0.93. We also show that glucose metabolism and phenol metabolism function with similar relative ATP maintenance costs. These findings demonstrate that iGR1773 can help the metabolic engineering community predict aromatic substrate utilization patterns and perform computational strain design.

Roell, Garrett W.↗
DOE OSTI2021

Use of Neutrino Scattering Events with Low Hadronic Recoil to Inform Neutrino Flux and Detector Energy Scale

Charged-current neutrino interactions with low hadronic recoil ("low-nu") have a cross-section that is approximately constant versus neutrino energy. These interactions have been used to measure the shape of neutrino fluxes as a function of neutrino energy at accelerator-based neutrino experiments such as CCFR, NuTeV, MINOS and MINERvA. In this paper, we demonstrate that low-nu events can be used to measure parameters of neutrino flux and detector models and that utilization of event distributions over the upstream detector face can discriminate among parameters that affect the neutrino flux model. From fitting a large sample of low-nu events obtained by exposing MINERvA to the NuMI medium-energy beam, we find that the best-fit flux parameters are within their a priori uncertainties, but the energy scale of muons reconstructed in the MINOS detector is shifted by 3.6% (or 1.8 times the a priori uncertainty on that parameter). These fit results are now used in all MINERvA cross-section measurements, and this technique can be applied by other experiments operating at MINERvA energies, such as DUNE.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Examining a hadronic γ-ray scenario for the radiative shell and molecular clouds of the old GeV supernova remnant G298.6−0.0

Abstract Based on the 13.7 yr Fermi-LAT data, Yeung, Bamba, and Sano (2023, PASJ, 75, 384) claimed detection of two γ-ray sources (namely Src-NE and Src-NW) associated with the supernova remnant (SNR) G298.6−0.0, and interpreted it as an old GeV SNR interacting with molecular clouds (MCs). In this follow-up study, we refine the flux measurements below 2 GeV with Fermi-LAT event types of better angular reconstruction. Then, we report our γ-ray spectral fittings and cosmic-ray phenomenology in a hadronic scenario, considering both the shell and MC regions of SNR G298.6−0.0. We confirm that the γ-ray spectra of both Src-NE and Src-NW exhibit spectral breaks, at $1.50_{-0.50}^{+0.60}$ and $0.68_{-0.11}^{+0.32}\:$GeV, respectively. Src-NW has a harder broad-band photon index than Src-NE, suggesting an appreciable difference between the physical separations of their respective emission sites from SNR G298.6−0.0. The cosmic-ray spectrum responsible for Src-NE starts with a minimum energy $E_\mathrm{CR,min}=1.38_{-0.16}^{+0.47}\:$GeV, and has a proton index $\Gamma _\mathrm{CR}=2.57_{-0.21}^{+0.18}$ below the exponential cutoff energy $E_\mathrm{CR,max}=240_{-150}^{+240}\:$GeV. Accordingly, we argue that Src-NE is dominated by the SNR shell, while only a minor portion of lower-energy emission is contributed by the MCs interacting with the SNR. The cosmic-ray population for Src-NW starts at a higher energy such that the ECR, min ratio of Src-NW to Src-NE is ≳2. The high ECR, min, as well as the high cosmic-ray energy density required (∼26 eV cm−3), supports the interpretation that Src-NW is predominantly the γ-ray emission from the farther MCs being bombarded by protons that had earlier escaped from SNR G298.6−0.0. By comparing the high-energy features of G298.6−0.0 with those of analogical SNRs, especially SNR W 28 and SNR W 44, we further constrain the age of SNR G298.6−0.0 to be 10–30 kyr, under the assumption of a purely hadronic scenario.

Astronomy & Astrophysics↗
DOE OSTI2025

GRB 221009A: The B.O.A.T. Burst that Shines in Gamma Rays

We present a complete analysis of Fermi Large Area Telescope (LAT) data of GRB 221009A, the brightest gamma-ray burst (GRB) ever detected. The burst emission above 30 MeV detected by the LAT preceded, by 1 s, the low-energy (<10 MeV) pulse that triggered the Fermi Gamma-Ray Burst Monitor (GBM), as has been observed in other GRBs. The prompt phase of GRB 221009A lasted a few hundred seconds. It was so bright that we identify a bad time interval of 64 s caused by the extremely high flux of hard X-rays and soft gamma rays, during which the event reconstruction efficiency was poor and the dead time fraction quite high. The late-time emission decayed as a power law, but the extrapolation of the late-time emission during the first 450 s suggests that the afterglow started during the prompt emission. We also found that high-energy events observed by the LAT are incompatible with synchrotron origin, and, during the prompt emission, are more likely related to an extra component identified as synchrotron self-Compton (SSC). A remarkable 400 GeV photon, detected by the LAT 33 ks after the GBM trigger and directionally consistent with the location of GRB 221009A, is hard to explain as a product of SSC or TeV electromagnetic cascades, and the process responsible for its origin is uncertain. Because of its proximity and energetic nature, GRB 221009A is an extremely rare event.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Grid-independent upwind scheme for multidimensional flow

Recent advances in the development of a grid-independent finite volume scheme for the Euler equations of gas dynamics are described. In the proposed method, flowfield gradient data are reconstructed locally (on a triangle) using five elementary planar waves, and an upwind numerical flux function for grid-oblique waves is developed to model the effect of the passage of these waves on the data in a cell. Numerical examples for several two-dimensional test problems are included. These results show high wave resolution and nearly monotone strong-wave transitions.

Parpia, Ijaz H.↗
NASA NTRS2013

Surface Ocean pCO2 Seasonality and Sea-Air CO2 Flux Estimates for the North American East Coast

Underway and in situ observations of surface ocean pCO2, combined with satellite data, were used to develop pCO2 regional algorithms to analyze the seasonal and interannual variability of surface ocean pCO2 and sea-air CO2 flux for five physically and biologically distinct regions of the eastern North American continental shelf: the South Atlantic Bight (SAB), the Mid-Atlantic Bight (MAB), the Gulf of Maine (GoM), Nantucket Shoals and Georges Bank (NS+GB), and the Scotian Shelf (SS). Temperature and dissolved inorganic carbon variability are the most influential factors driving the seasonality of pCO2. Estimates of the sea-air CO2 flux were derived from the available pCO2 data, as well as from the pCO2 reconstructed by the algorithm. Two different gas exchange parameterizations were used. The SS, GB+NS, MAB, and SAB regions are net sinks of atmospheric CO2 while the GoM is a weak source. The estimates vary depending on the use of surface ocean pCO2 from the data or algorithm, as well as with the use of the two different gas exchange parameterizations. Most of the regional estimates are in general agreement with previous studies when the range of uncertainty and interannual variability are taken into account. According to the algorithm, the average annual uptake of atmospheric CO2 by eastern North American continental shelf waters is found to be between 3.4 and 5.4 Tg C/yr (areal average of 0.7 to 1.0 mol CO2 /sq m/yr) over the period 2003-2010.

interannual variability↗