Results for “Flux reconstruction” — Page 5

NASA NTRSDate not supplied

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

NASA NTRS1995

Adaptively Refined Euler and Navier-Stokes Solutions with a Cartesian-Cell Based Scheme

A Cartesian-cell based scheme with adaptive mesh refinement for solving the Euler and Navier-Stokes equations in two dimensions has been developed and tested. Grids about geometrically complicated bodies were generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells were created using polygon-clipping algorithms. The grid was stored in a binary-tree data structure which provided a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations were solved on the resulting grids using an upwind, finite-volume formulation. The inviscid fluxes were found in an upwinded manner using a linear reconstruction of the cell primitives, providing the input states to an approximate Riemann solver. The viscous fluxes were formed using a Green-Gauss type of reconstruction upon a co-volume surrounding the cell interface. Data at the vertices of this co-volume were found in a linearly K-exact manner, which ensured linear K-exactness of the gradients. Adaptively-refined solutions for the inviscid flow about a four-element airfoil (test case 3) were compared to theory. Laminar, adaptively-refined solutions were compared to accepted computational, experimental and theoretical results.

Coirier, William J.↗

NASA NTRSDate not supplied

Lessons Learned from Aerothermal Flight Data on NASA Mars Entry Vehicles

The NASA Mars Science Laboratory (MSL) and Mars 2020 entry vehicles included temperature measurements inside and direct heat flux on the surface of the thermal protection system (TPS) material. The instrumentation was included to begin addressing predictive capabilities for aerodynamic heating and ablative TPS material response. The MSL aeroshell included temperature measurements at seven heatshield locations, from which total aerodynamic heating was reconstructed. The Mars 2020 instrumentation suite included eleven heatshield measurement locations and new backshell measurements at nine locations: six with temperature sensors, two with sensors to directly measure total heat flux, and one to measure radiative heat flux. Both instrumentation suites returned full sets of data, with all temperature sensors nearest to the surface surviving entry, indicating minimal material recession. Heatshield boundary layer transition (BLT) was observed for both vehicles prior to the time of peak heating, with transition starting near the shoulder and progressing towards the stagnation area. MSL’s higher speed on approach to Mars resulted in higher total heating compared to Mars 2020, with the maximum occurring near the heatshield’s leeside shoulder. The Mars 2020 backshell heating data indicate significant contributions from shock layer radiation. Heating from Navier-Stokes flowfield calculations are in family with reconstructed heating. The paper covers the instrumentation, BLT and reconstructed heating on both heatshields, backshell heating for Mars 2020, TPS performance, and predicted versus reconstructed heat fluxes. Opportunities for improved design approaches for future Mars missions are discussed, especially for the upcoming Mars Sample Return mission.

Mars Science Laboratory↗

NASA NTRSDate not supplied

Lessons Learned from Aerothermal Flight Data on NASA Mars Entry Vehicles

The NASA Mars Science Laboratory (MSL) and Mars 2020 entry vehicles included temperature measurements inside and direct heat flux on the surface of the thermal protection system (TPS) material. The instrumentation was included to begin addressing predictive capabilities for aerodynamic heating and ablative TPS material response. The MSL aeroshell included temperature measurements at seven heatshield locations, from which total aerodynamic heating was reconstructed. The Mars 2020 instrumentation suite included eleven heatshield measurement locations and new backshell measurements at nine locations: six with temperature sensors, two with sensors to directly measure total heat flux, and one to measure radiative heat flux. Both instrumentation suites returned full sets of data, with all temperature sensors nearest to the surface surviving entry, indicating minimal material recession. Heatshield boundary layer transition (BLT) was observed for both vehicles prior to the time of peak heating, with transition starting near the shoulder and progressing towards the stagnation area. MSL’s higher speed on approach to Mars resulted in higher total heating compared to Mars 2020, with the maximum occurring near the heatshield’s leeside shoulder. The Mars 2020 backshell heating data indicate significant contributions from shock layer radiation. Heating from Navier-Stokes flowfield calculations are in family with reconstructed heating. The paper covers the instrumentation, BLT and reconstructed heating on both heatshields, backshell heating for Mars 2020, TPS performance, and predicted versus reconstructed heat fluxes. Opportunities for improved design approaches for future Mars missions are discussed, especially for the upcoming Mars Sample Return mission.

Mars Science Laboratory↗

NASA NTRS1992

Development of a High Resolution Liquid Xenon Imaging Telescope for Medium Energy Gamma Ray Astrophysics

In the third year of the research project, we have (1) tested a 3.5 liter prototype of the Liquid Xenon Time Projection Chamber, (2) used a prototype having a 4.4 cm drift gap to study the charge and energy resolution response of the 3.5 liter chamber, (3) obtained an energy resolution as good as that previously measured by us using chambers with drift gaps of the order of millimeters, (4) observed the induction signals produced by MeV gamma rays, (4) used the 20 hybrid charge sensitive preamplifiers for a nondestructive readout of the electron image on the induction wires, (5) performed extensive Monte Carlo simulations to obtain results on efficiency, background rejection capability, and source flux sensitivity, and (6) developed a reconstruction algorithm for events with multiple interaction points.

Aprile, Elena↗

NASA NTRS2003

Finite Volume Methods: Foundation and Analysis

Finite volume methods are a class of discretization schemes that have proven highly successful in approximating the solution of a wide variety of conservation law systems. They are extensively used in fluid mechanics, porous media flow, meteorology, electromagnetics, models of biological processes, semi-conductor device simulation and many other engineering areas governed by conservative systems that can be written in integral control volume form. This article reviews elements of the foundation and analysis of modern finite volume methods. The primary advantages of these methods are numerical robustness through the obtention of discrete maximum (minimum) principles, applicability on very general unstructured meshes, and the intrinsic local conservation properties of the resulting schemes. Throughout this article, specific attention is given to scalar nonlinear hyperbolic conservation laws and the development of high order accurate schemes for discretizing them. A key tool in the design and analysis of finite volume schemes suitable for non-oscillatory discontinuity capturing is discrete maximum principle analysis. A number of building blocks used in the development of numerical schemes possessing local discrete maximum principles are reviewed in one and several space dimensions, e.g. monotone fluxes, E-fluxes, TVD discretization, non-oscillatory reconstruction, slope limiters, positive coefficient schemes, etc. When available, theoretical results concerning a priori and a posteriori error estimates are given. Further advanced topics are then considered such as high order time integration, discretization of diffusion terms and the extension to systems of nonlinear conservation laws.

Barth, Timothy↗

NASA NTRS2018

Increased Snowfall over the Antarctic Ice Sheet Mitigated Twentieth-Century Sea-Level Rise

Changes in accumulated snowfall over the Antarctic Ice Sheet have an immediate and time-delayed impact on global mean sea level. The immediate impact is due to the instantaneous change in freshwater storage over the ice sheet, whereas the time-delayed impact acts in opposition through enhanced ice dynamic flux into the ocean. Here, we reconstruct 200 years of Antarctic-wide snow accumulation by synthesizing a newly compiled database of ice core records using reanalysis derived spatial coherence patterns. The results reveal that increased snow accumulation mitigated twentieth-century sea-level rise by approximately 10 millimeters since 1901, with rates increasing from 1.1 mm per decade between 1901 and 2000 to 2.5 mm per decade after 1979. Reconstructed accumulation trends are highly variable in both sign and magnitude at the regional scale, and linked to the trend towards a positive Southern Annular Mode since 1957. Because the observed Southern Annular Mode trend is accompanied by a decrease in Antarctic Ice Sheet accumulation, changes in the strength and location of the circumpolar westerlies cannot explain the reconstructed increase, which may instead be related to stratospheric ozone depletion. However, our results indicate that a warming atmosphere cannot be excluded as a dominant force in the underlying increase.

Antarctic↗

NASA NTRS2018

Increased snowfall over the Antarctic Ice Sheet mitigated twentieth-century sea-level rise

Changes in accumulated snowfall over the Antarctic Ice Sheet have an immediate and time-delayed impact on global mean sea level. The immediate impact is due to the instantaneous change in freshwater storage over the ice sheet, whereas the time-delayed impact acts in opposition through enhanced ice-dynamic flux into the ocean1. Here, we reconstruct 200 years of Antarctic-wide snow accumulation by synthesizing a newly compiled database of ice core records2 using reanalysis-derived spatial coherence patterns. The results reveal that increased snow accumulation mitigated twentieth-century sea-level rise by ~10 mm since 1901, with rates increasing from 1.1 mm decade−1 between 1901 and 2000 to 2.5 mm decade−1 after 1979. Reconstructed accumulation trends are highly variable in both sign and magnitude at the regional scale, and linked to the trend towards a positive Southern Annular Mode since 19573. Because the observed Southern Annular Mode trend is accompanied by a decrease in Antarctic Ice Sheet accumulation, changes in the strength and location of the circumpolar westerlies cannot explain the reconstructed increase, which may instead be related to stratospheric ozone depletion4. However, our results indicate that a warming atmosphere cannot be excluded as a dominant force in the underlying increase.

B Medley↗

NASA NTRS2012

Effects of Mesh Irregularities on Accuracy of Finite-Volume Discretization Schemes

The effects of mesh irregularities on accuracy of unstructured node-centered finite-volume discretizations are considered. The focus is on an edge-based approach that uses unweighted least-squares gradient reconstruction with a quadratic fit. For inviscid fluxes, the discretization is nominally third order accurate on general triangular meshes. For viscous fluxes, the scheme is an average-least-squares formulation that is nominally second order accurate and contrasted with a common Green-Gauss discretization scheme. Gradient errors, truncation errors, and discretization errors are separately studied according to a previously introduced comprehensive methodology. The methodology considers three classes of grids: isotropic grids in a rectangular geometry, anisotropic grids typical of adapted grids, and anisotropic grids over a curved surface typical of advancing layer grids. The meshes within the classes range from regular to extremely irregular including meshes with random perturbation of nodes. Recommendations are made concerning the discretization schemes that are expected to be least sensitive to mesh irregularities in applications to turbulent flows in complex geometries.

Diskin, Boris↗

NASA NTRSDate not supplied

Reconstruction of Thermal Protection System Aeroheating using a Green’s Function Approach

Inverse heat transfer (IHT) techniques are often used to reconstruct the surface heating conditions on spacecraft thermal protection systems (TPS) during atmospheric entry. Current IHT techniques for entry spacecraft applications, however, demand substantial computational resources, and are impractical for analyses such as uncertainty quantification and real-time health monitoring. In this paper, a Green’s function sensor fusion approach is used to reconstruct the TPS surface aeroheating conditions on experimental spaceflight and ground test systems from collocated temperature and heat flux sensors embedded in the TPS. The algorithm leverages Green’s functions to model the heat conduction within the spacecraft TPS and stabilizes the recovery of the surface heating condition using the direct heat flux sensor measurement. The algorithm is validated using arc-jet ground test data and applied to the reconstruction of the Mars 2020 backshell heating during Martian atmospheric entry. The performance of the algorithm is benchmarked against a current state-of-the-art IHT framework, FIAT_Opt. The Green’s function-based reconstruction algorithm recovers the net hot-wall heat flux absorbed by the TPS and the incident heat flux from the atmospheric entry environment in close agreement with FIAT_Opt. Notably, computation of the surface heating condition is completed in three orders of magnitude less time with the Green’s function sensor fusion approach using a consumer-grade PC, versus with FIAT_Opt running on a high performance computer cluster. The efficiency of the algorithm is leveraged to compute the uncertainty contributions of input parameters to the total uncertainty in reconstructed Mars 2020 backshell heating for the full atmospheric entry heat pulse. The sensitivity analysis uncovers that, at different times throughout the entry heat pulse, uncertainties in the TPS specific heat, thermal conductivity, and emissivity are all dominant drivers of the reconstruction uncertainty. These results demonstrate Green’s functions and sensor-fusion techniques as promising IHT approaches to reconstruct atmospheric entry environments from TPS-embedded measurements, and highlight how these techniques may give access to post-flight analyses previously hindered by the prohibitive cost of current methods.

Kenneth McAfee↗

NASA NTRSDate not supplied

Reconstruction of Thermal Protection System Aeroheating using a Green’s Function Approach

Inverse heat transfer (IHT) techniques are often used to reconstruct the surface heating conditions on spacecraft thermal protection systems (TPS) during atmospheric entry. Current IHT techniques for entry spacecraft applications, however, demand substantial computational resources, and are impractical for analyses such as uncertainty quantification and real-time health monitoring. In this paper, a Green’s function sensor fusion approach is used to reconstruct the TPS surface aeroheating conditions on experimental spaceflight and ground test systems from collocated temperature and heat flux sensors embedded in the TPS. The algorithm leverages Green’s functions to model the heat conduction within the spacecraft TPS and stabilizes the recovery of the surface heating condition using the direct heat flux sensor measurement. The algorithm is validated using arc-jet ground test data and applied to the reconstruction of the Mars 2020 backshell heating during Martian atmospheric entry. The performance of the algorithm is benchmarked against a current state-of-the-art IHT framework, FIAT_Opt. The Green’s function-based reconstruction algorithm recovers the net hot-wall heat flux absorbed by the TPS and the incident heat flux from the atmospheric entry environment in close agreement with FIAT_Opt. Notably, computation of the surface heating condition is completed in three orders of magnitude less time with the Green’s function sensor fusion approach using a consumer-grade PC, versus with FIAT_Opt running on a high performance computer cluster. The efficiency of the algorithm is leveraged to compute the uncertainty contributions of input parameters to the total uncertainty in reconstructed Mars 2020 backshell heating for the full atmospheric entry heat pulse. The sensitivity analysis uncovers that, at different times throughout the entry heat pulse, uncertainties in the TPS specific heat, thermal conductivity, and emissivity are all dominant drivers of the reconstruction uncertainty. These results demonstrate Green’s functions and sensor-fusion techniques as promising IHT approaches to reconstruct atmospheric entry environments from TPS-embedded measurements, and highlight how these techniques may give access to post-flight analyses previously hindered by the prohibitive cost of current methods.

Kenneth McAfee↗

NASA NTRS2023

Understanding the Deflection of the “Cartwheel CME”: Data Analysis and Modeling

We study the low corona evolution of the “Cartwheel” coronal mass ejection (CME; 2008 April 9) by- reconstructing its three-dimensional path and modeling it with magnetohydrodynamic simulations. This event exhibited a double deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The Cartwheel CME traveled toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full-trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the preeruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight into the motion exhibited by the 2008 April 9 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR toward the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.

Solar prominences↗

NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/sq cm for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/sq cm, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8-species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/cm2 for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/cm2, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8- species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗

NASA NTRS2022

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes OPTIM3D and RADMC3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N-band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.

stars: atmospheres↗

NASA NTRS1999

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗

NASA NTRSDate not supplied

Direct Reconstruction of Ablative Thermal Protection System Aeroheating Using A Green's Function Approach

A Green’s function inverse heat transfer (IHT) approach is used to reconstruct the surface heating conditions on ablative thermal protection system (TPS) materials from embedded heat flux sensor and temperature probe measurements. The approach models the temperature time-history at the measurement location as a discrete linear system, allowing for the heat flux boundary condition to be recovered directly without the need for time-marching schemes. The effects of material decomposition and pyrolysis gas transport are modeled using an energy source/sink analogue. The performance of the reconstruction approach is analyzed on a 1D test case representative of an atmospheric entry heating scenario. The approach can recover the TPS surface heat flux to within 4% of the input heating condition with a computation time of 2-3 seconds (>3 orders of magnitude faster than current time-marching IHT methods). As a byproduct of the surface heating reconstruction, the algorithm also captures the surface pyrolysis gas mass flux and solid decomposition at multiple through-thickness locations within the TPS.

Kenneth McAfee↗

NASA NTRS2018

Elliptic-Cylindrical Analytical Flux Rope Model for Magnetic Clouds

In this paper, we present the elliptic-cylindrical analytical flux rope model, which constitutes the first level of complexity above that of a circular-cylindrical geometry. The framework of this series of models was established by Nieves-Chinchilla et al. with the circular-cylindrical analytical flux rope model. The model describes the magnetic flux rope topology with distorted cross section as a possible consequence of the flux rope interaction with the solar wind. In this model, for the first time, a flux rope is completely described by a nonorthogonal geometry. The Maxwell equations can be consistently solved using tensorial analysis, and relevant physical quantities can be derived, such as magnetic fluxes, number of turns, or Lorentz force distribution. The model is generalized in termsof the radial dependence of the poloidal and axial current density components. The circular-cylindrical reconstruction technique has been adapted to the new geometry for a specific case of the model and tested againstan interplanetary coronal mass ejection observed by the Wind spacecraft on 2005 June 12. In this specific case,from the comparative analysis between the circular-cylindrical and elliptic-cylindrical models, the inclusion ofthe cross-section distortion in the 3D reconstruction results in significant changes in the derived axis orientation,size, central magnetic field, magnetic fluxes, and force-freeness. The case studied in this paper exemplifies the useof the model and reconstruction technique developed. Furthermore, the novel mathematical formulation to modelflux ropes in heliophysics paves the way to the inclusion of more complex magnetic field configurations.

CMEs↗