Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable density”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO S2S Forecast System

Sea surface salinity (SSS) observations from space allow us to investigate if improved estimates of near-surface density stratification and associated mixing can positively impact seasonal to interannual variability of tropical Pacific Ocean dynamics as well as dynamical ENSO forecasts. For the first part of the presentation, we utilize our intermediate-complexity coupled model. Baseline experiments assimilate satellite sea level (multi-satellite gridded AVISO, 2013), SST (Reynolds et al., 2004), and in situ subsurface temperature and salinity observations (GTSPP NODC, 2006). These baseline experiments are then compared with experiments that additionally assimilate Aquarius (V5.0 Lilly and Lagerloef, 2008) and SMAP (V4.0 Fore et al., 2016) SSS. Twelve-month forecasts are initialized for each month from September 2011 to September 2017. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. This pattern enhances air/sea interaction and amplifies the equatorial Kelvin wave signal. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation over most forecast lead times. We next assess how different satellite SSS products impact the validation of ENSO forecasts. SMAP V4 reduces the salty bias in the western Pacific and so is an improvement upon SMAP V2 and SMOS V2 (Boutin et al., 2017) has similar validation characteristics as a combination of Aquarius and SMAP V4. Next we shift to present results from the NASA GMAO Sub-seasonal to seasonal (S2S_v2.1) production coupled model (i.e. the same model that contributes ENSO forecasts to the North American Multi-Model Ensemble Experiment). From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of satellite sea surface salinity (SSS) observations on dynamical ENSO forecasts for the big 2015 El Nino event. We will present distinct experiments for the overlap period that include 1) freely evolving SSS (i.e. no satellite SSS as the production system), 2) Aquarius, and 3) SMAP initialization. Our results show that using Aquarius slightly improves validation of the reanalysis (including sea level and temperature statistics). Our production system without SSS assimilation generated too warm forecasts for the 2015 El Nino from March initial conditions. Incorporating Aquarius into initialization of the coupled system leads to a deeper, more realistic MLD that acts to damp the downwelling Kelvin signal and slightly cool NINO3.4 SST. With Aquarius the forecasts better match the observed amplitude of the 2015 event. On the other hand, SMAP V2 relaxation generally degrades validation statistics. At forecast initialization, SMAP is much too salty within 10o of the equator, leading to deeper MLD east of 165W. This deeper MLD leads to over-damping of the downwelling signal (i.e. relative upwelling), in turn leading to relatively too cool ENSO forecasts.

Hackert, E.↗
DOE OSTI2018

The Kepler Light Curves of AGN: A Detailed Analysis

Here, we present a comprehensive analysis of 21 light curves of Type 1 active galactic nuclei (AGN) from the Kepler spacecraft. First, we describe the necessity and development of a customized pipeline for treating Kepler data of stochastically variable sources like AGN. We then present the light curves, power spectral density functions (PSDs), and flux histograms. The light curves display an astonishing variety of behaviors, many of which would not be detected in ground-based studies, including switching between distinct flux levels. Six objects exhibit PSD flattening at characteristic timescales that roughly correlate with black hole mass. These timescales are consistent with orbital timescales or free-fall accretion timescales. We check for correlations of variability and high-frequency PSD slope with accretion rate, black hole mass, redshift, and luminosity. We find that bolometric luminosity is anticorrelated with both variability and steepness of the PSD slope. We do not find evidence of the linear rms–flux relationships or lognormal flux distributions found in X-ray AGN light curves, indicating that reprocessing is not a significant contributor to optical variability at the 0.1%–10% level.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Calcium abundance variations in diffuse interstellar clouds

Observations are discussed which show that the Na I/H2 column-density ratio is relatively constant among different clouds where H2 is optically thin, but the Na I/Ca II column-density ratio varies. This implies that the gas-phase abundance of calcium is variable. It is proposed that the wide variation in the gas-phase calcium abundance can be explained by small variations in the calcium content of grains. The observed enhancement of the Ca II/Na I column-density ratio with cloud velocity is attributed to grain destruction in grain-grain collisions within the same cloud. The structure of shock waves propagating through the interstellar gas is computed, and the grain-destruction model is compared with observational data.

Jura, M.↗
NASA NTRS2000

RXTE Observations of the Seyfert 2 Galaxy MrK 348

We present RXTE monitoring observations of the Seyfert 2 galaxy Mrk 348 spanning a 6 month period. The time-averaged spectrum in the 3-20 keV band shows many features characteristic of a Compton-thin Seyfert 2 galaxy, namely a hard underlying power-law continuum (Gamma approximately equal 1.8) with heavy soft X-ray absorption (N(sub H) approximately 10(exp 23)/sq cm) plus measurable iron K.alpha emission (equivalent width approximately 100 eV) and, at high energy, evidence for a reflection component (R approximately < 1). During the first half of the monitoring period the X-ray continuum flux from Mrk 348 remained relatively steady. However this was followed by a significant brightening of the source (by roughly a factor of 4) with the fastest change corresponding to a doubling of its X-ray flux on a timescale of about 20 days. The flux increase was accompanied by a marked softening of X-ray spectrum most likely attributable to a factor approximately 3 decline in the intrinsic line-of-sight column density. In contrast the iron K.alpha line and the reflection components showed no evidence of variability. These observations suggest a scenario in which the central X-ray source is surrounded by a patchy distribution of absorbing material located within about a light-week of the nucleus of Mrk 348. The random movement of individual clouds within the absorbing screen, across our line of sight, produces substantial temporal variations in the measured column density on timescales of weeks to months and gives rise to the observed X-ray spectral variability. However, as viewed from the nucleus the global coverage and typical thickness of the cloud layer remains relatively constant.

Smith, David A.↗
NASA NTRS2007

Hydrodynamical Simulations of the Jet in the Symbiotic Star MWC 560: Application to X-ray Jets in Symbiotic Stars - 3

In Papers I and II in this series, we presented hydrodynamical simulations of jet models with parameters representative of the symbiotic system MWC 560. These were simulations of a pulsed, initially underdense jet in a high-density ambient medium. Since the pulsed emission of the jet creates internal shocks and since the jet velocity is very high, the jet bow shock and the internal shocks are heated to high temperatures and should therefore emit X-ray radiation. In this paper, we investigate in detail the X-ray properties of the jets in our models. We have focused our study on the total X-ray luminosity and its temporal variability, the resulting spectra, and the spatial distribution of the emission. Temperature and density maps from our hydrodynamical simulations with radiative cooling presented in the second paper are used, together with emissivities calculated with the atomic database ATOMDB. The jets in our models show extended and variable X-ray emission, which can be characterized as a sum of hot and warm components with temperatures that are consistent with observations of CH Cyg and R Aqr. The X-ray spectra of our model jets show emission-line features that correspond to observed features in the spectra of CH Cyg. The innermost parts of our pulsed jets show iron line emission in the 6.4-6.7 keV range, which may explain such emission from the central source in R Aqr. We conclude that MWC 560 should be detectable with Chandra or XMM-Newton, and such X-ray observations will prove crucial for understanding jets in symbiotic stars.

hydrodynamics↗
DOE OSTI2024

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

High resolution density pedestal measurements during edge localized modes by short-pulse reflectometry in the TCV tokamak

Here, this publication presents high spatio-temporal resolution (mm µs -1 ) density profile measurements of the pedestal top during type I, III, and small edge localized mode (ELM) H-mode plasmas in the tokamak à configuration variable (TCV). These measurements were performed using a novel short-pulse reflectometer. Average inter-ELM density profiles are obtained via conditional averaging using the D α trace as ELM indicator. Changes to the pedestal density profile gradients prior to type-III ELMs reveal unique pedestal dynamics leading to the ELM crash which can provide important experimental data for validation of non-linear MHD ELM simulations. The small-ELM scenario is found to feature a ~25–35 kHz quasi-coherent density fluctuation near the separatrix ρ ψ 0.993–1.05 not observed during a similar type-I ELM discharge. This oscillation is also found in low-field-side magnetic pick-up probes displaying a ballooning character and n = +1 toroidal mode number. This oscillation could help explain the markedly different pedestal dynamics observed in the small-ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High resolution density pedestal measurements during edge localized modes by short-pulse reflectometry in the TCV tokamak

This publication presents high spatio-temporal resolution (mm/μs) density profile measurements of the pedestal top during type I, III, and small edge localized mode (ELM) H-mode plasmas in the Tokamak à Configuration Variable (TCV). These measurements were performed using a novel short-pulse reflectometer. Average inter-ELM density profiles are obtained via conditional averaging using the Dα trace as ELM indicator. Changes to the pedestal density profile gradients prior to type-III ELMs reveal unique pedestal dynamics leading to the ELM crash which can provide important experimental data for validation of non-linear MHD ELM simulations. The small-ELM scenario is found to feature a ∼25-35 kHz quasi-coherent density fluctuation near the separatrix rho_psi ∼0.993-1.05 not observed during a similar type-I ELM discharge. This oscillation is also found in low-field-side magnetic pick- up probes displaying a ballooning character and n=+1 toroidal mode number. This oscillation could help explain the markedly different pedestal dynamics observed in the small-ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Statistical Analyses of Post Irradiation Experiments of the BWR Forsmark

Computer codes which are to be used in nuclear safety research need to be validated prior to their application. For depletion codes, this validation can be done using data from post irradiation experiments (PIE). These experiments consist of a radiochemical analysis of small parts of a fuel rod (FR), which has been irradiated in a nuclear power plant as part of a fuel assembly (FA). This analysis allows the deduction of the isotopic composition of nuclides important for the desired application case. In addition the operation history at the sample location needs to be known and reported with the analysis results. In our work, we perform uncertainty and sensitivity analyses of a model based on a sample from the boiling water reactor (BWR) Forsmark-3 in order to deduce the influence of uncertain system parameters on the nuclide composition of the depleted fuel. In a second step we compare the resulting uncertainties of the nuclide densities with the deviation between their calculated and experimental values. The goal is to obtain quantitative uncertainty data which can be accounted for in a validation procedure. We use the Monte Carlo sampling approach, i.e. repeating the model calculation using statistically varied sets of the uncertain input parameters and analyzing the distributions of the result values statistically. BWR depletion calculations are especially demanding due to the increased complexity and heterogeneity compared to pressurized water reactors (PWR): extensively variable enrichments, complex FA geometries, partial length FRs, strong axially varying moderator density and the common use of gadolinium as burnable absorber in some of the FRs. However, the precise and correct calculation of BWR-depletion is extremely important in view of a possible application in burnup credit for BWR FAs. Therein credit is taken for the reduced reactivity of spent fuel due to the depletion of {sup 235}U for transport and intermediate or final storage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2017

Incorporating Planetary-Scale Waves Into the VTGCM: Understanding the Waves Impact on the Upper Atmosphere of Venus.

Venus has proven to have a very dynamic upper atmosphere. The upper atmosphere of Venus has been observed for many decades by multiple means of observation (e.g. ground-based, orbiters, probes, fly-by missions going to other planets). As of late, the European Space Agency Venus Express (VEX) orbiter has been a main observer of the Venusian atmosphere. Specifically, observations of Venus' O2 IR nightglow emission have been presented to show its variability. Nightglow emission is directly connected to Venus' circulation and is utilized as a tracer for the atmospheric global wind system. More recent observations are adding and augmenting temperature and density (e.g. CO, CO2, SO2) datasets. These additional datasets provide a means to begin analyzing the variability and study the potential drivers of the variability. A commonly discussed driver of variability is wave deposition. Evidence of waves has been observed, but these waves have not been completely analyzed to understand how and where they are important. A way to interpret the observations and test potential drivers is by utilizing numerical models.

Venus↗
NASA NTRS2009

Lunar Neutral Exposphere Properties from Pickup Ion Analysis

Composition and structure of neutral constituents in the lunar exosphere can be determined through measurements of phase space distributions of pickup ions borne from the exosphere [1]. An essential point made in an early study [ 1 ] and inferred by recent pickup ion measurements [2, 3] is that much lower neutral exosphere densities can be derived from ion mass spectrometer measurements of pickup ions than can be determined by conventional neutral mass spectrometers or remote sensing instruments. One approach for deriving properties of neutral exospheric source gasses is to first compare observed ion spectra with pickup ion model phase space distributions. Neutral exosphere properties are then inferred by adjusting exosphere model parameters to obtain the best fit between the resulting model pickup ion distributions and the observed ion spectra. Adopting this path, we obtain ion distributions from a new general pickup ion model, an extension of a simpler analytic description obtained from the Vlasov equation with an ion source [4]. In turn, the ion source is formed from a three-dimensional exospheric density distribution, which can range from the classical Chamberlain type distribution to one with variable exobase temperatures and nonthermal constituents as well as those empirically derived. The initial stage of this approach uses the Moon's known neutral He and Na exospheres to deriv e He+ and Na+ pickup ion exospheres, including their phase space distributions, densities and fluxes. The neutral exospheres used are those based on existing models and remote sensing studies. As mentioned, future ion measurements can be used to constrain the pickup ion model and subsequently improve the neutral exosphere descriptions. The pickup ion model is also used to estimate the exosphere sources of recently observed pickup ions on KAGUYA [3]. Future missions carrying ion spectrometers (e.g., ARTEMIS) will be able to study the lunar neutral exosphere with great sensitivity, yielding the necessary ion velocity spectra needed to further analysis of parent neutral exosphere properties.

Hartle, R. E.↗
NASA NTRS1999

A Cutoff in the X-Ray Fluctuation Power Density Spectrum of the Seyfert 1 Galaxy NGC 3516

During 1997 March-July, RXTE observed the bright, strongly variable Seyfert 1 galaxy NGC 3516 once every approx. 12.8 hr for 4.5 months and nearly continuously (with interruptions due to SAA passage but not Earth occultation) for a 4.2 day period in the middle. These were followed by ongoing monitoring once every approx. 4.3 days. These data are used to construct the first well-determined X-ray fluctuation power density spectrum (PDS) of an active galaxy to span more than 4 decades of usable temporal frequency. The PDS shows no signs of any strict or quasi-periodicity, but does show a progressive flattening of the power-low slope from -1.74 at short time scales to -0.73 at longer time scales. This is the clearest observation to date of the long-predicted cutoff in the PDS. The characteristic variability time scale corresponding to this cutoff temporal frequency is approx. 1 month. Although it is unclear how this time scale may be interpreted in terms of a physical size or process, there are several promising candidate models. The PDS appears similar to those seen for Galactic black hole candidates such as Cyg X-1, suggesting that these two classes of objects with very different luminosities and putative black hole masses (differing by more than a factor of 10(exp 5)) may have similar X-ray generation processes and structures.

Edelson, Rick↗
NASA NTRS2016

The Nature of the Torus in the Heavily Obscured AGN Markarian 3: an X-Ray Study

In this paper, we report the results of an X-ray monitoring campaign on the heavily obscured Seyfert galaxy, Markarian 3, carried out between the fall of 2014 and the spring of 2015 with NuSTAR, Suzaku and XMMNewton. The hard X-ray spectrum of Markarian 3 is variable on all the time-scales probed by our campaign, down to a few days. The observed continuum variability is due to an intrinsically variable primary continuum seen in transmission through a large, but still Compton-thin column density (N(sub H) approx. 0.8-1.1 x 10(exp 24)/sq cm). If arranged in a spherical-toroidal geometry, the Compton scattering matter has an opening angle approx. 66deg, and is seen at a grazing angle through its upper rim (inclination angle approx. 70deg). We report a possible occultation event during the 2014 campaign. If the torus is constituted by a system of clouds sharing the same column density, this event allows us to constrain their number (17 +/- 5) and individual column density, [approx. (4.9 +/- 1.5) x 10(exp 22)/ sq cm]. The comparison of IR and X-ray spectroscopic results with state-of-the art torus models suggests that at least two-thirds of the X-ray obscuring gas volume might be located within the dust sublimation radius. We report also the discovery of an ionized absorber, characterized by variable resonant absorption lines due to He- and H-like iron. This discovery lends support to the idea that moderate column density absorbers could be due to clouds evaporated at the outer surface of the torus, possibly accelerated by the radiation pressure due to the central AGN emission leaking through the patchy absorber.

X-rays: individual: Markarian 3↗
NASA NTRS2004

An Expanded RXTE Survey of Long-Term X-ray Variability in Seyfert 1 Galaxies

The first seven years of RXTE monitoring of Seyfert 1 active galactic nuclei have been systematically analyzed to yield five homogenous samples of 2-12 keV light curves, probing hard X-ray variability on successively longer durations from approx. 1 day to approx. 3.5 years. 2-10 keV variability on time scales of approx. 1 day, as probed by ASCA, are included. All sources exhibit stronger X-ray variability towards longer time scales, with variability amplitudes saturating at the longest time scales, but the increase is greater for relatively higher luminosity sources. The well-documented anticorrelation between variability amplitude and luminosity is confirmed on all time scales. However, anticorrelations between variability amplitude and black hole mass estimate are evident on only the shortest time scales probed. The data are consistent with the models of power spectral density (PSD) movement described in Markowitz et al. (2003) and McHardy et al. (2004), whereby Seyfert 1 galaxies variability can be described by a single, universal PSD shape whose cutoff frequency scales with black hole mass. The best-fitting scaling relations between variability time scale, black hole mass and X-ray luminosity support an average accretion rate of 2% of the Eddington limit for the sample. Nearly all sources exhibit stronger variability in the relatively soft 2-4 keV band compared to the 7-12 keV band on all time scales. Color-flux diagrams support also Seyfert 1s' softening as they brighten. There are indications that relatively less luminous or less massive sources exhibit a greater degree of spectral variability for a given increase in overall flux.

Markowitz, A.↗
DOE OSTI2017

Cusps in the center of galaxies: a real conflict with observations or a numerical artefact of cosmological simulations?

Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of observations favor variable profiles with a core in the center. In this paper, we investigate the convergency of standard N-body simulations, especially in the cusp region, following the approach proposed by [1]. We simulate the well known Hernquist model using the SPH code Gadget-3 and consider the full array of dynamical parameters of the particles. We find that, although the cuspy profile is stable, all integrals of motion characterizing individual particles suffer strong unphysical variations along the whole halo, revealing an effective interaction between the test bodies. This result casts doubts on the reliability of the velocity distribution function obtained in the simulations. Moreover, we find unphysical Fokker-Planck streams of particles in the cusp region. The same streams should appear in cosmological N-body simulations, being strong enough to change the shape of the cusp or even to create it. Our analysis, based on the Hernquist model and the standard SPH code, strongly suggests that the UDPs generally found by the cosmological N-body simulations may be a consequence of numerical effects. A much better understanding of the N-body simulation convergency is necessary before a 'core-cusp problem' can properly be used to question the validity of the CDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Lagrangian and Dirac constraints for the ideal incompressible fluid and magnetohydrodynamics

The incompressibility constraint for fluid flow was imposed by Lagrange in the so-called Lagrangian variable description using his method of multipliers in the Lagrangian (variational) formulation. An alternative is the imposition of incompressibility in the Eulerian variable description by a generalization of Dirac’s constraint method using noncanonical Poisson brackets. Here it is shown how to impose the incompressibility constraint using Dirac’s method in terms of both the canonical Poisson brackets in the Lagrangian variable description and the noncanonical Poisson brackets in the Eulerian description, allowing for the advection of density. Both cases give the dynamics of infinite-dimensional geodesic flow on the group of volume preserving diffeomorphisms and explicit expressions for this dynamics in terms of the constraints and original variables is given. Since Lagrangian and Eulerian conservation laws are not identical, comparison of the various methods is made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Venus lower thermosphere studies

Studies undertaken in this project have sought to understand lower thermospheric structure and dynamics (less than or equal to 145 km), particularly the processes responsible. This is a region just below the reach of in-situ instruments onboard PVO (Pioneer Venus Orbiter) during the first few diurnal cycles. PVO remote airglow observations (nitric oxide, O2, visible, O 1304A) have been coupled with ground-based observations (CO densities, winds, temperatures, O2 IR nightglow) to address the behavior of lower thermospheric winds and chemistry over 95 to 150 km. This interpretation of PVO and related data is accomplished by using the NCAR Venus thermospheric general circulation model (VTGCM) (Bougher et al., 1988; 1990). This model has been modified over the last two years to improve its ability to calculate O, CO, and O2 densities, temperatures, nightglow, and subsolar-to-antisolar and zonal winds over 95 to 150 km (Bougher and Borucki, 1992). Our VTGCM studies show that: (1) O2 visible and IR nightglow distributions can be used to trace lower thermosphere / upper mesosphere winds over 100-130 km. Typically, weak zonal winds (less than or equal to 25 m/sec) and nightglow maximum patches near 0100 LT prevail. Occasionally, strong zonal winds (30-60 m/sec) and airglow patches peaking near 0300 LT characterize the Venus lower thermosphere. (2) It is clear that the dynamics of the Venus 90-130 km region is highly variable on time scales as short as an hour. This is most likely due to the time variable nature of upward propagating gravity waves, which grow in amplitude and eventually break. The resulting turbulence gives rise to local time variable eddy diffusion and momentum drag, both of which strongly impact global density and nightglow distributions. (3) The oxygen chemistry (O, O2, etc.) over 90-120 km is strongly dependent on HO(x) and CLO(x) tracer species that must be properly included in any coupled chemical dynamical model. (4) The density profiles of light species (O, CO, N, He) are strongly affected by large-scale transport by the winds. Strong eddy diffusion is not a suitable model parameterization for approximating these light species, especially for extrapolation into the region below 140 km where PVO in-situ data is lacking. Instead, the fully coupled chemical dynamical VTGCM model should be used to improve estimates of densities within the VTS3 empirical model (Hedin et al., 1983) below 140 km.

Bougher, Stephen W.↗
NASA NTRS1992

The relationship between pressure-based and density-based algorithms

The PISO, pressure-based algorithm, is compared with implicit time-marching systems to ascertain their similarities and differences. Both methods are expressed in vector form for comparison purposes. The vector form of the PISO method is triangular, allowing an uncoupled solution procedure, while the Euler implicit method requires the simultaneous solution of all equations. Upwind differencing is performed according to the direction of the eigenvalues in both systems, but this is the particle velocity in the PISO method and the acoustic velocity in the Euler method. Vector stability calculations show that the PISO method is conditionally stable depending on the time step, but that it provides adequate damping at small time steps to enable good convergence. Unconditional stability can be provided by retaining the energy coupling terms in the time derivative. Transforming the Euler implicit equations to the PISO variables along with a modification of the time derivatives gives the density-based method the same amplification factors at low speeds as the pressure-based method without affecting their behavior at high speeds.

Merkle, Charles L.↗