Search NASA⌕ Search

SEARCH · Search NASA

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 127 records · Page 7

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↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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↗

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↗

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↗

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↗

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.↗

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↗

Artificial Neural Networks for Determining Magnetospheric Conditions

This chapter presents a neural-network-based technique that allows for the reconstruction of the global, time-varying distribution of some physical quantity Q, that has been sparsely sampled at various locations within the magnetosphere, and at different times. We begin with a general introduction to the problem of prediction and specification, and why it is important and difficult to achieve with existing methods. We then provide a basic introduction to neural networks, and describe our technique using the specific example of reconstructing the electron plasma density in the Earth's inner magnetosphere on the equatorial plane. We then show more advanced uses of the technique, including 3D reconstruction of the plasma density, specification of chorus and hiss waves, and energetic particle fluxes. We summarize and conclude with a general discussion of how machine learning techniques might be used to advance the state-of-the-art in space weather prediction, and insight discovery.

Bortnik, Jacob↗

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present the first fifth order, semi-discrete central upwind method for approximating solutions of multi-dimensional Hamilton-Jacobi equations. Unlike most of the commonly used high order upwind schemes, our scheme is formulated as a Godunov-type scheme. The scheme is based on the fluxes of Kurganov-Tadmor and Kurganov-Tadmor-Petrova, and is derived for an arbitrary number of space dimensions. A theorem establishing the monotonicity of these fluxes is provided. The spacial discretization is based on a weighted essentially non-oscillatory reconstruction of the derivative. The accuracy and stability properties of our scheme are demonstrated in a variety of examples. A comparison between our method and other fifth-order schemes for Hamilton-Jacobi equations shows that our method exhibits smaller errors without any increase in the complexity of the computations.

Bryson, Steve↗

Timing Analysis with INTEGRAL: Comparing Different Reconstruction Algorithms

INTEGRAL is one of the few instruments capable of detecting X-rays above 20keV. It is therefore in principle well suited for studying X-ray variability in this regime. Because INTEGRAL uses coded mask instruments for imaging, the reconstruction of light curves of X-ray sources is highly non-trivial. We present results from the comparison of two commonly employed algorithms, which primarily measure flux from mask deconvolution (ii-lc-extract) and from calculating the pixel illuminated fraction (ii-light). Both methods agree well for timescales above about 10 s, the highest time resolution for which image reconstruction is possible. For higher time resolution, ii-light produces meaningful results, although the overall variance of the lightcurves is not preserved.

Grinberg, V.↗

ANTARES Constrains a Blazar Origin of Two IceCube PeV Neutrino Events

Context. The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Objects like these are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. Aims. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope. Methods. The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons - and hence their neutrino progenitors - from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for di erent possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin. Results. Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653-329 and 1714-336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC 14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than -2.4.

general↗

Time-Resolved Ultraviolet Spectroscopy of The M-Dwarf GJ 876 Exoplanetary System

Extrasolar planets orbiting M-stars may represent our best chance to discover habitable worlds in the coming decade. The ultraviolet spectrum incident upon both Earth-like and Jovian planets is critically important for proper modeling of their atmospheric heating and chemistry. In order to provide more realistic inputs for atmospheric models of planets orbiting low-mass stars, we present new near- and far-ultraviolet (NUV and FUV) spectroscopy of the M-dwarf exoplanet host GJ 876 (M4V). Using the COS and STIS spectrographs on board the Hubble Space Telescope, we have measured the 1150-3140 A spectrum of GJ 876. We have reconstructed the stellar H1 Ly alpha emission line profile, and find that the integrated Ly alpha flux is roughly equal to the rest of the integrated flux (1150-1210 A + 1220-3140 A) in the entire ultraviolet bandpass (F(Ly alpha)/F(FUV+NUV) equals approximately 0.7). This ratio is approximately 2500x greater than the solar value. We describe the ultraviolet line spectrum and report surprisingly strong fluorescent emission from hot H2 (T(H2) greater than 2000 K). We show the light curve of a chromospheric + transition region flare observed in several far-UV emission lines, with flare/quiescent flux ratios greater than or equal to 10. The strong FUV radiation field of an M-star (and specifically Ly alpha) is important for determining the abundance of O2--and the formation of biomarkers-in the lower atmospheres of Earth-like planets in the habitable zones of low-mass stars.

France, Kevin↗