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
NASA NTRS2020

InSight Aerothermal Environment Assessment

The InSight spacecraft was proposed to be a build-to-print copy of the Phoenix vehicle due to the knowledge that the lander payload would be similar and the trajectory would be similar. However, the InSight aerothermal analysts, based on tests performed in CO2 during the Mars Science Laboratory mission (MSL) and completion of Russian databases, considered radiative heat flux to the aftbody from the wake for the first time for a US Mars mission. The combined convective and radiative heat flux was used to determine if the as-flown Phoenix thermal protection system (TPS) design would be sufficient for InSight. All analyses showed that the design would be adequate. Once the InSight lander was successfully delivered to Mars on November 26, 2018, work began to reconstruct the atmosphere and trajectory in order to evaluate the aerothermal environments that were actually encountered by the spacecraft and to compare them to the design environments.The best estimated trajectory (BET) reconstructed for the InSight atmospheric entry fell between the two trajectories considered for the design, when looking at the velocity versus altitude values. The maximum heat rate design trajectory (MHR) flew at a higher velocity and the maximum heat load design trajectory (MHL) flew at a lower velocity than the BET. For TPS sizing, the MHL trajectory drove the design. Reconstruction has shown that the BET flew for a shorter time than either of the design environments, hence total heat load on the vehicle should have been less than used in design. Utilizing the BET, both DPLR and LAURA were first run to analyze the convective heating on the vehicle with no angle of attack. Both codes were run with axisymmetric, laminar flow in radiative equilibrium and vibrational non-equilibrium with a surface emissivity of 0.8. Eight species Mitcheltree chemistry was assumed with CO2, CO, N2, O2, NO, C, N, and O. Both codes agreed within 1% on the forebody and had the expected differences on the aftbody. The NEQAIR and HARA codes were used to analyze the radiative heating on the vehicle using full spherical ray-tracing. The codes agreed within 5% on most aftbody points of interest.The LAURA code was then used to evaluate the conditions at angle of attack at the peak heating and peak pressure times. Boundary layer properties were investigated to confirm that the flow over the forebody was laminar for the flight.Comparisons of the aerothermal heating determined for the reconstructed trajectory to the design trajectories showed that the as-flown conditions were less severe than design

Beck, Robin↗
NASA NTRS1991

Imaging the sun in hard X-rays - Spatial and rotating modulation collimators

Several approaches to imaging hard X-rays emitted from solar flares have been proposed or are planned for the nineties including the spatial modulation collimator (SMC) and the rotating modulation collimator (RMC). A survey of current solar flare theoretical literature indicates the desirability of spatial resolutions down to 1 arcsecond, field of views greater than the full solar disk (i.e., 32 arcminutes), and temporal resolutions down to 1 second. Although the sun typically provides relatively high flux levels, the requirement for 1 second temporal resolution raises the question as to the viability of Fourier telescopes subject to the aforementioned constraints. A basic photon counting, Monte Carlo 'end-to-end' model telescope was employed using the Astronomical Image Processing System (AIPS) for image reconstruction. The resulting solar flare hard X-ray images compared against typical observations indicated that both telescopes show promise for the future.

Campbell, Jonathan W.↗
NASA NTRS2001

Computed Tomography Support for Microgravity Materials Science Experiments

The accurate measurement of density in both liquid and solid samples is of considerable interest to Principal Investigators with materials science experiments slated for the ISS. The work to be described is an innovative application of a conventional industrial nondestructive evaluation instrument. Traditional applications of industrial computed tomography (CT) rely on reconstructing cross sections of large structures to provide two-dimensional planar views which can identify defects such as porosity, or other material anomalies. This has been done on microgravity materials science experiments to check the integrity of ampoule-cartridge assemblies for safety purposes. With a substantially monoenergetic flux, as can be obtained with a radioactive cobalt source, there will be a direct correlation between absorption and density. Under such conditions it then becomes possible to make accurate measurements of density throughout a sample, and even when the sample itself is enclosed within a furnace and a safety required cartridge. Such a system has been installed at Kennedy Space Center (KSC) and is available to PIs to examine samples before and after flight. The CT system is being used to provide density information for two purposes. Firstly, the determination of density changes from liquid to solid is vital information to the PI for purposes of modeling the solidification behavior of his sample, and to engineers who have to design containment ampoules and must allow for shrinkage and other volume changes that may occur during processing. While such information can be obtained by pycnometric measurements, the possibility of using a furnace installed on the CT system enables one to examine potentially dangerous materials having high vapor pressures, while not needing visible access to the material. In addition, uniform temperature can readily be obtained, and the system can be controlled to ramp up, hold, and ramp down while collecting data over a wide range of parameters automatically. Results of initial tests on low melting point elements such as gallium, indium and tin will be presented, and the intent is to proceed to compounds such as InSb, HgCdTe and CdTe. Alloys such as Pb-Sb (PI - Poirier, U AZ) and Cu-Al (PI - Trivedi, Ames Lab.), which are the subjects of flight experiments, will also be examined. The second application is the conversion of measured density values directly to composition. This was successfully done with the mercury cadmium telluride alloys grown on the second and fourth United States Microgravity Payload (USMP-2 and USMP-4) missions by Lehoczky. CdTe values along the length of the boules were obtained at KSC prior to cutting the sample, and could have been obtained prior to its removal from the cartridge and ampoule. Examples of the data obtained will be shown. It is anticipated that several of the materials science PIs will avail themselves of the technique described, initially for determining densities prior to flight, and then to acquire early quantitative data on the compositional variation within their samples.

Gillies, Donald C.↗
NASA NTRS2016

Deep Morphological and Spectral Study of the SNR RCW 86 with Fermi-LAT

RCW 86 is a young supernova remnant (SNR) showing a shell-type structure at several wavelengths and is thought to be an efficient cosmic-ray (CR) accelerator. Earlier Fermi Large Area Telescope results reported the detection of gamma-ray emission coincident with the position of RCW 86 but its origin (leptonic or hadronic) remained unclear due to the poor statistics. Thanks to 6.5 years of data acquired by the Fermi-LAT and the new event reconstruction Pass 8, we report the significant detection of spatially extended emission coming from RCW 86. The spectrum is described by a power-law function with a very hard photon index (Gamma = 1.42 +/- 0.1(sub stat) +/- 0.06(sub syst)) in the 0.1-500 GeV range and an energy flux above 100 MeV of (2.91 +/- 0.8(sub stat) +/- 0.12(sub syst)) x 10(exp -11) erg/(sq cms). Gathering all the available multiwavelength (MWL) data, we perform a broadband modeling of the non-thermal emission of RCW 86 to constrain parameters of the nearby medium and bring new hints about the origin of the gamma-ray emission. For the whole SNR, the modeling favors a leptonic scenario in the framework of a two-zone model with an average magnetic field of 10.2 +/- 0.7 microG and a limit on the maximum energy injected into protons of 2 x 10(exp 49) erg for a density of 1 per cu cm. In addition, parameter values are derived for the north-east and south-west (SW) regions of RCW 86, providing the first indication of a higher magnetic field in the SW region.

acceleration of particles – cosmic rays – ISM:↗
NASA NTRSDate not supplied

Compressive Sensing Application for Transient Photometric Measurements

Compressive Sensing (CS) is a mathematical theory for simultaneous data acquisition and compression. Natural phenomena may be sparse in some physical or temporal domain. If we exploit this sparsity by applying the technique of CS to obtain information, how do our measurements change as a function of domain and measurement systematics? What are the specific implications for the science and for the sensing infrastructure? In this talk, we will discuss the generalized systematic effects consequent to the application of CS to time-series photometric measurements. We assess implications for observability, sparsification, and information loss in the detection, retrieval and reconstruction process. To study time-series photometry, we explore the field of gravitational microlensing. A source star, typically in the galactic bulge, gets microlensed when there is a precise alignment of a lensing star and its planetary system, with the source star. The microlensed source star changes in flux magnification as the lensing system crosses the precise path of alignment, resulting in a microlensing curve in time domain. A high-cadence, high-resolution system, which uses low power and bandwidth is essential to obtain valuable science measurements. Hence, we discuss application of CS to gravitational microlensing data sets, which in turn can be generalized to any time-ordered photometric measurements.

Asmita Abhay Korde↗
NASA NTRS2005

Two Successive Martian Years on the Orbit: Similarities and Differences of CO2 Seasonal Cycle from HEND/ODYSSEY Data

The three years of Mars Odyssey successful work on the martian orbit provide a lot of new information about peculiarities of long term variations of CO2 seasonal cycle. To start such analysis we have used observations of neutron albedo of Mars obtained by High Energy Neutron detector (HEND) mounted onboard Mars Odyssey spacecraft. The high latitude northern and southern regions of Mars are affected by global redistribution of atmospheric CO2 which resulted in 25% of atmospheric mass condensed on martian surface of these regions during winter period of time. The seasonal deposit is formed starting from 60N/60S latitudes and achieve its maximal thickness about 1 m at latitudes close to martian poles. Changes of CO2 deposit thickness is the reason for significant variations of neutron flux above martian poles from summer to winter seasons because CO2 frost effectively hides upper water rich surface layers from the orbit observations in neutrons and gamma-rays. This effect was used to estimate column density of CO2 deposit at different latitudes on North and South of Mars and reconstruct multidimensional model of CO2 deposit showing how snow depth varies as function of latitude, longitude and time. In this presentation we tried to make a next step in our study of martian seasonal CO2 cycle and look for similarities and differences between two successive martian years.

Litvak, M. L.↗
NASA NTRS2021

The Dark Side of Pluto

During its departure from Pluto, New Horizons used its LORRI camera to image a portion of Pluto's southern hemisphere that was in a decades-long seasonal winter darkness, but still very faintly illuminated by sunlight reflected by Charon. Recovery of this faint signal was technically challenging. The bright ring of sunlight forward-scattered by haze in the Plutonian atmosphere encircling the nightside hemisphere was severely overexposed, defeating the standard smeared-charge removal required for LORRI images. Reconstruction of the overexposed portions of the raw images, however, allowed adequate corrections to be accomplished. The small solar elongation of Pluto during the departure phase also generated a complex scattered-sunlight background in the images that was three orders of magnitude stronger than the estimated Charon-light flux (the Charon-light flux is similar to the flux of moonlight on Earth a few days before first quarter). A model background image was constructed for each Pluto image based on principal component analysis applied to an ensemble of scattered-sunlight images taken at identical Sun−spacecraft geometry to the Pluto images. The recovered Charon-light image revealed a high-albedo region in the southern hemisphere. We argue that this may be a regional deposit of N2 or CH4 ice. The Charon-light image also shows that the south polar region currently has markedly lower albedo than the north polar region of Pluto, which may reflect the sublimation of N2 ice or the deposition of haze particulates during the recent southern summer.

Pluto↗
NASA NTRS2010

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and complexity are studied for four nominally second-order accurate schemes: a node-centered scheme and three cell-centered schemes - a node-averaging scheme and two schemes with nearest-neighbor and adaptive compact stencils for least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Tests from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The tests of the second class are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes may degenerate on mixed grids, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to that of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping based on a distance function commonly available in practical schemes or modifying the scheme stencil to reflect the direction of strong coupling. The major conclusion is that accuracies of the node centered and the best cell-centered schemes are comparable at equivalent number of degrees of freedom.

Diskin, Boris↗
NASA NTRS1995

Solar Particle Induced Upsets in the TDRS-1 Attitude Control System RAM During the October 1989 Solar Particle Events

The three large solar particle events, beginning on October 19, 1989 and lasting approximately six days, were characterized by high fluences of solar protons and heavy ions at 1 AU. During these events, an abnormally large number of upsets (243) were observed in the random access memory of the attitude control system (ACS) control processing electronics (CPE) on-board the geosynchronous TDRS-1 (Telemetry and Data Relay Satellite). The RAM unit affected was composed of eight Fairchild 93L422 memory chips. The Galileo spacecraft, launched on October 18, 1989 (one day prior to the solar particle events) observed the fluxes of heavy ions experienced by TDRS-1. Two solid-state detector telescopes on-board Galileo, designed to measure heavy ion species and energy, were turned on during time periods within each of the three separate events. The heavy ion data have been modeled and the time history of the events reconstructed to estimate heavy ion fluences. These fluences were converted to effective LET spectra after transport through the estimated shielding distribution around the TDRS-1 ACS system. The number of single event upsets (SEU) expected was calculated by integrating the measured cross section for the Fairchild 93L422 memory chip with average effective LET spectrum. The expected number of heavy ion induced SEU's calculated was 176. GOES-7 proton data, observed during the solar particle events, were used to estimate the number of proton-induced SEU's by integrating the proton fluence spectrum incident on the memory chips, with the two-parameter Bendel cross section for proton SEU'S. The proton fluence spectrum at the device level was gotten by transporting the protons through the estimated shielding distribution. The number of calculated proton-induced SEU's was 72, yielding a total of 248 predicted SEU'S, very dose to the 243 observed SEU'S. These calculations uniquely demonstrate the roles that solar heavy ions and protons played in the production of SEU's during the October 1989 solar particle events.

Croley, D. R.↗
NASA NTRS2009

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes - Part 1

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling.

Diskin, Boris↗
NASA NTRS1995

Coupled Ablation, Heat Conduction, Pyrolysis, Shape Change and Spallation of the Galileo Probe

The Galileo probe enters the atmosphere of Jupiter in December 1995. This paper presents numerical methodology and detailed results of our final pre-impact calculations for the heat shield response. The calculations are performed using a highly modified version of a viscous shock layer code with massive radiation coupled with a surface thermochemical ablation and spallation model and with the transient in-depth thermal response of the charring and ablating heat shield. The flowfield is quasi-steady along the trajectory, but the heat shield thermal response is dynamic. Each surface node of the VSL grid is coupled with a one-dimensional thermal response calculation. The thermal solver includes heat conduction, pyrolysis, and grid movement owing to surface recession. Initial conditions for the heat shield temperature and density were obtained from the high altitude rarefied-flow calculations of Haas and Milos. Galileo probe surface temperature, shape, mass flux, and element flux are all determined as functions of time along the trajectory with spallation varied parametrically. The calculations also estimate the in-depth density and temperature profiles for the heat shield. All this information is required to determine the time-dependent vehicle mass and drag coefficient which are necessary inputs for the atmospheric reconstruction experiment on board the probe.

Milos, Frank S.↗
NASA NTRS2018

The Zwicky Transient Facility: Data Processing, Products, and Archive

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ∼3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5σ in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The real-time pipeline employs a novel image-differencing algorithm, optimized for the detection of point-source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. We present some preliminary results of the calibration performance delivered by the real-time pipeline. The reconstructed astrometric accuracy per science image with respect to Gaia DR1 is typically 45 to 85 milliarcsec. This is the RMS per-axis on the sky for sources extracted with photometric S/N ≥10 and hence corresponds to the typical astrometric uncertainty down to this limit. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. The high end of these ranges corresponds to an airmass approaching ∼2—the limit of the public survey. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae; rare flux transients; variable stars; eclipsing binaries; variability from active galactic nuclei; counterparts to gravitational wave sources; a more complete census of Type Ia supernovae; and solar-system objects.

Frank J. Masci↗
NASA NTRS2022

The Influence of Recent and Future Climate Change on Spring Arctic Cyclones

In recent decades, the Arctic has experienced rapid atmospheric warming and sea ice loss, with an ice-free Arctic projected by the end of this century. Cyclones are synoptic weather events that transport heat and moisture into the Arctic, and have complex impacts on sea ice, and the local and global climate. However, the effect of a changing climate on Arctic cyclone behavior remains poorly understood. This study uses high resolution (4 km), regional modeling techniques and downscaled global climate reconstructions and projections to examine how recent and future climatic changes alter cyclone behavior. Results suggest that recent climate change has not yet had an appreciable effect on Arctic cyclone characteristics. However, future sea ice loss and increasing surface temperatures drive large increases in the near-surface temperature gradient, sensible and latent heat fluxes, and convection during cyclones. The future climate can alter cyclone trajectories and increase and prolong intensity with greatly augmented wind speeds, temperatures, and precipitation. Such changes in cyclone characteristics could exacerbate sea ice loss and Arctic warming through positive feedbacks. The increasing extreme nature of these weather events has implications for local ecosystems, communities, and socio-economic activities.

Chelsea L Parker↗
NASA NTRS1998

Polar Hydra Data Analysis

The science activities are: 1) Hydra is still operating successfully on orbit. 2) A large amount of analysis and discovery has occurred with the Hydra ground data processing this past year. 3) Full interdetector calibration has been implemented and documented. This intercalibration was necessitated by the incorrect installation of bias resistors in the pre-acceleration stage to the electron channeltrons. This had the effect of making the counting efficiency for electrons energy dependent as well as channeltron specific. The nature of the error had no impact on the ion detection efficiency since they have a different bias arrangement. This intercalibration is so effective, that the electron and ion moment densities are routinely produced with a level of agreement better than 20%. 4) The data processing routinely removes glint in the sensors and produces public energy time spectrograms on the web overnight. 6) Routine, but more intensive computer processing codes are operational that determine for electrons and ions, the density, the flow vector, the pressure tensor and the heat flux by numerical integration. These codes use the magnetic field to sustain the quality of their output. To gain access to this high quality magnetic field within our data stream we have monitored Russell's web page for zero levels and timing files (since his data acquisition is not telemetry synchronous) and have a local reconstruction of B for our use. We have also detected a routine anomaly in the magnetometer data stream that we have documented to Chris Russell and developed an editing algorithm to intercept these "hits" and remove them from the geophysical analysis.

Scudder, J. D.↗
NASA NTRS1993

Solar flux forecasting using mutual information with an optimal delay

Solar flux F(sub 10.7) directly affects the atmospheric density, thereby changing the lifetime and prediction of satellite orbits. For this reason, accurate forecasting of F(sub 10.7) is crucial for orbit determination of spacecraft. Our attempts to model and forecast F(sub 10.7) uncovered highly entangled dynamics. We concluded that the general lack of predictability in solar activity arises from its nonlinear nature. Nonlinear dynamics allow us to predict F(sub 10.7) more accurately than is possible using stochastic methods for time scales shorter than a characteristic horizon, and with about the same accuracy as using stochastic techniques when the forecasted data exceed this horizon. The forecast horizon is a function of two dynamical invariants: the attractor dimension and the Lyapunov exponent. In recent years, estimation of the attractor dimension reconstructed from a time series has become an important tool in data analysis. In calculating the invariants of the system, the first necessary step is the reconstruction of the attractor for the system from the time-delayed values of the time series. The choice of the time delay is critical for this reconstruction. For an infinite amount of noise-free data, the time delay can, in principle, be chosen almost arbitrarily. However, the quality of the phase portraits produced using the time-delay technique is determined by the value chosen for the delay time. Fraser and Swinney have shown that a good choice for this time delay is the one suggested by Shaw, which uses the first local minimum of the mutual information rather than the autocorrelation function to determine the time delay. This paper presents a refinement of this criterion and applies the refined technique to solar flux data to produce a forecast of the solar activity.

Ashrafi, S.↗
NASA NTRS2020

Reconstructing Magma Storage Depths for the 2018 Kilauean Eruption from melt inclusion CO2 Contents: The importance of Vapor Bubbles

The 2018 Lower East Rift Zone (LERZ) eruption of Kīlauea Volcano and the accompanying collapse of the summit caldera marked the most destructive phase of activity on Hawai’i in the last 200 years. The integration of petrological data extracted from lava samples collected throughout the eruption with geodetic data examining the caldera collapse event, and estimates of the co-erupted flux of SO2 from the main eruptive fissure (Fissure 8), provides an exceptional opportunity to determine the reservoir geometry and magma transport paths supplying Kīlauea’s LERZ. The forsterite contents of erupted olivines and the degree of disequilibrium with their carrier melts indicate that two distinct olivine populations were erupted from Fissure 8. Melt inclusion entrapment pressures reveal that more evolved olivines (Fo<81.5) crystallized at ~2 km depth within the shallower Halema’uma’u reservoir, while more primitive olivines (Fo>81.5)crystallized within the deeper South Caldera reservoir at ~3–5 km depth. Crucially, primitive olivines experienced extensive post-entrapment crystallization, driving the growth of a vapor bubble. Raman spectroscopy reveals that this bubble contains up to 99% of the total inclusionCO2 budget (median=93%). Measurements of CO2 in only the glass phase would have underestimated entrapment depths by up to 60× (median=11×), and the importance of the SC reservoir as a source of magma to Fissure 8 would have been overlooked. Overall, we demonstrate that Raman measurements of bubbles, along with careful choice of suitably-calibrated H2O-CO2 solubility model, is vital to place accurate constraints on the depths of magma storage regions supplying volcanic eruptions.

SIMS↗
NASA NTRS2010

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes

Cell-centered and node-centered approaches have been compared for unstructured finite-volume discretization of inviscid fluxes. The grids range from regular grids to irregular grids, including mixed-element grids and grids with random perturbations of nodes. Accuracy, complexity, and convergence rates of defect-correction iterations are studied for eight nominally second-order accurate schemes: two node-centered schemes with weighted and unweighted least-squares (LSQ) methods for gradient reconstruction and six cell-centered schemes two node-averaging with and without clipping and four schemes that employ different stencils for LSQ gradient reconstruction. The cell-centered nearest-neighbor (CC-NN) scheme has the lowest complexity; a version of the scheme that involves smart augmentation of the LSQ stencil (CC-SA) has only marginal complexity increase. All other schemes have larger complexity; complexity of node-centered (NC) schemes are somewhat lower than complexity of cell-centered node-averaging (CC-NA) and full-augmentation (CC-FA) schemes. On highly anisotropic grids typical of those encountered in grid adaptation, discretization errors of five of the six cell-centered schemes converge with second order on all tested grids; the CC-NA scheme with clipping degrades solution accuracy to first order. The NC schemes converge with second order on regular and/or triangular grids and with first order on perturbed quadrilaterals and mixed-element grids. All schemes may produce large relative errors in gradient reconstruction on grids with perturbed nodes. Defect-correction iterations for schemes employing weighted least-square gradient reconstruction diverge on perturbed stretched grids. Overall, the CC-NN and CC-SA schemes offer the best options of the lowest complexity and secondorder discretization errors. On anisotropic grids over a curved body typical of turbulent flow simulations, the discretization errors converge with second order and are small for the CC-NN, CC-SA, and CC-FA schemes on all grids and for NC schemes on triangular grids; the discretization errors of the CC-NA scheme without clipping do not converge on irregular grids. Accurate gradient reconstruction can be achieved by introducing a local approximate mapping; without approximate mapping, only the NC scheme with weighted LSQ method provides accurate gradients. Defect correction iterations for the CC-NA scheme without clipping diverge; for the NC scheme with weighted LSQ method, the iterations either diverge or converge very slowly. The best option in curved geometries is the CC-SA scheme that offers low complexity, second-order discretization errors, and fast convergence.

Diskin, Boris↗
NASA NTRS1990

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗