Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial refinement”

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 163 records · Page 9

Modeling Contamination Migration on the Chandra X-ray Observatory - II

During its first 14 years of operation, the cold (about -60C) optical blocking filter of the Advanced CCD Imaging Spectrometer (ACIS), aboard the Chandra X-ray Observatory, has accumulated a growing layer of molecular contamination that attenuates low-energy x rays. Over the past few years, the accumulation rate, spatial distribution, and composition have changed. This evolution has motivated further analysis of contamination migration within and near the ACIS cavity. To this end, the current study employs a higher-fidelity geometric model of the ACIS cavity, detailed thermal modeling based upon temperature data, and a refined model of the molecular transport.

O'Dell, Stephen L.↗

Simultaneous solution for core magnetic field and fluid flow beneath an electrically conducting mantle

The effects of laterally homogeneous mantle electrical conductivity were included in steady, frozen-flux core surface flow estimation along with refinements in method and weighting. The refined method allows simultaneous solution for both the initial radial geomagnetic field component at the core-mantle boundary (CMB) and the sub-adjacent fluid motion; it also features Gauss' method for solving the non-linear inverse problem associated with steady motional induction. The tradeoff between spatial complexity of the derived flows and misfit to the weighted Definitive Geomagnetic Reference Field models (DGRF's) is studied for various mantle conductivity profiles. For simple flow and a fixed initial geomagnetic condition, a fairly high deep-mantle conductivity performs better than either insulating or weakly conducting profiles; however, a thin, very high conductivity layer at the base of the mantle performs almost as well. Simultaneous solution for both initial geomagnetic field and flow reduces the misfit per degree of freedom even more than does changing the mantle conductivity profile. Moreover, when both core field and flow are estimated, the performance of the solutions and the derived flows become insensitive to the conductivity profile.

Voorhies, Coerte V.↗

Fractional Snowcover Estimates from Earth Observing System (EOS) Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS)

The Moderate Resolution Imaging Spectroradiometer (MODIS) on the NASA Earth Observing System (EOS) Terra and Aqua missions has shown considerable capability for mapping snowcover. The typical approach that has used, along with other criteria, the Normalized Snow Difference Index (NDSI) that takes the difference between 500 meter observations at 1.64 micrometers (MODIS band 6) and 0.555 micrometers (MODIS band 4) over the sum of these observations to determine whether MODIS pixels are snowcovered or not in mapping the extent of snowcover. For many hydrological and climate studies using remote sensing of snowcover, it is desirable to assess if the MODIS snowcover observations could not be enhanced by providing the fraction of snowcover in each MODIS observation (pixel). Pursuant to this objective studies have been conducted to assess whether there is sufficient "signal%o in the NDSI parameter to provide useful estimates of fractional snowcover in each MODIS 500 meter pixel. To accomplish this objective high spatial resolution (30 meter) Landsat snowcover observations were used and co-registered with MODIS 500 meter pixels. The NDSI approach was used to assess whether a Landsat pixel was or was not snowcovered. Then the number of snowcovered Landsat pixels within a MODIS pixel was used to determine the fraction of snowcover within each MODIS pixel. The e results were then used to develop statistical relationships between the NDSI value for each 500 meter MODIS pixel and the fraction of snowcover in the MODIS pixel. Such studies were conducted for three widely different areas covered by Landsat scenes in Alaska, Russia, and the Quebec Province in Canada. The statistical relationships indicate that a 10 percent accuracy can be attained. The variability in the statistical relationship for the three areas was found to be remarkably similar (-0.02 for mean error and less than 0.01 for mean absolute error and standard deviation). Independent tests of the relationships were accomplished by taking the relationship of fractional snow-cover to NDSI from one area (e.g., Alaska) and testing it on the other two areas (e.g. Russia and Quebec). Again the results showed that fractional snow-cover can be estimated to 10 percent. The results have been shown to have advantages over other published fractional snowcover algorithms applied to MODIS data. Most recently the fractional snow-cover algorithm has been applied using 500-meter observations over the state of Colorado for a period spanning 25 days. The results exhibit good behavior in mapping the spatial and temporal variability in snowcover over that 25-day period. Overall these studies indicate that robust estimates of fractional snow-cover can be attained using the NDSI parameter over areas extending in size from watersheds relatively large compared to MODIS pixels to global land cover. Other refinements to this approach as well as different approaches are being examined for mapping fractional snow-cover using MODIS observations.

Salomonson, Vincent V.↗

Multiscale Studies of the Three-Dimensional Dayside X-Line

We review recent experience from the Cluster, Double Star, and THEMIS missions for lessons that apply to the upcoming Magnetospheric Multiscale Mission (MMS) being developed for launch in 2014. On global scales, simulation and statistical studies lead to mean congurations of dayside reconnection, implying specific relative alignments of the inflow magnetic fields and X-line, with implications for MMS operations designed to maximize the number of close encounters with the diffusion region. At intermediate MHD-to-ion scales, reconstruction of features created by one or two X-lines have developed to the point where data from a cluster of spacecraft can determine their temporal trends and the approximate three-dimensional X-line structure. Recent petascale particle-in-cell (PIC) simulations of reconnection encompass three spatial dimensions with excellent resolution, and make striking predictions of electron scale physics that creates complex interacting flux ropes under component reconnection. High time resolution measurements from MMS will determine the detailed electron scale kinetics embedded within the global and MHDion scale contexts. These developments will lead to the refinement of our three-dimensional multiscale picture of reconnection, yielding improved understanding of the global, MHD, and local physics controlling the onset or quenching, variability, and mean rate of reconnection. This in turn will enable improved predictability of the structural features created by transient reconnection, and their space weather consequences.

Reconnection↗

AWP: NASA's Aerosol Wind Profiler Coherent Doppler Wind Lidar

A new airborne coherent Doppler Wind Lidar (DWL) instrument has been implemented leveraging NASA's development of the 2-micron Wind-Space Pathfinder (Wind-SP) lidar transceiver. A technology development project, Wind-SP refined and demonstrated numerous components required for a coherent-detection space wind lidar instrument. Aerosol Wind Profiler (AWP) transitioned this transceiver into an airborne wind lidar capable of providing full 3-D wind vector retrievals. AWP operates with laser pulse energy and repetition rate combination required for high spatial and vertical resolution wind profiling from space. Operation from aircraft platforms will yield very strong signal return required for detailed process studies at <2km spatial and <100m vertical resolution under most conditions. AWP will serve as NASA’s wind calibration and validation instrument for future space-based wind observations over the coming decades, while also providing data supporting space wind lidar simulation studies. The AWP instrument is introduced and preliminary results from recent AWP demonstration flights are presented.

AWP↗

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗

Clouds and the earth's radiant energy system (CERES) - Instrument design and development

Measurements of the earth's reflected shortwave and emitted longwave energy and of the effect of clouds on these quantities are planned using a refined version of the Earth Radiation Budget Experiment (ERBE) scanning instrument. The CERES instruments are being designed to accumulate earth radiance measurements with a repeatability of better than 0.5 percent over their five year life. Beginning in 1996, flights are planned on both polar and low earth orbit satellites to obtain the required temporal and spatial coverage. The design and development of CERES are discussed.

Kopia, Leonard P.↗

The X-Ray Lightcurve of Eta Carinae: Refinement of the Orbit and Evidence for Phase Dependent Mass Loss

We solve the RXTE X-ray lightcurve of the extremely luminous and massive star eta Carinae with a colliding wind emission model to refine the ground-based orbital elements. The sharp decline to X-ray minimum at the end of 1997 fixes the date of the last periastron passage at 1997.95 +/- 0.05, not 1998.13 as derived from ground-based radial velocities. This helps resolve a discrepancy between the ground-based radial velocities and spatially-resolved velocity measures obtained by STIS. The X-ray data are consistent with a mass function f(M) approx. = 1.5, lower than the value f(M) approx. = 7.5 previously reported, so that the masses of eta Carinae and the companion are M(sub eta) greater than or = 80 solar mass and M(sub c) approx. 30 solar mass respectively. In addition the X-ray data suggest that the mass loss rate from eta Carinae is generally less than 3 x 10(exp -4) solar mass/yr, about a factor of 5 lower than that derived from some observations in other wavebands. We could not match the duration of the X-ray minimum with any standard colliding wind model in which the wind is spherically symmetric and the mass loss rate is constant. However we show that we can match the variations around X-ray minimum if we include an increase of a factor of approx. 20 in the mass loss rate from eta Carinae for approximately 80 days following periastron. If real, this excess in M would be the first evidence of enhanced mass flow off the primary when the two stars are close (presumably driven by tidal interactions). Our interpretation of the X-ray data suggest that the ASCA and RXTE X-ray spectra near the X-ray minimum are significantly contaminated by unresolved hard emission (E greater than or = 2 keV) from sonic other nearby source, probably associated with scattering of tile colliding wind emission by circumstellar dust. Based on the X-ray fluxes the distance to n Carinae is 2300 pc with formal uncertainties of only approx. 10%.

Corcoran, M. F.↗

Quantitative Approach for the Accurate CFD Simulation of Hover in Turbulent Flow

Time-dependent Navier-Stokes simulations have been carried out for a V22 rotor in hover using an improved, low-diffusion, HLLE++ upwind algorithm in the OVERFLOW code. Emphasis is placed on lessons learned over the past decade regarding the effects of high-order spatial accuracy, grid resolution, and the use of detached eddy simulation in predicting the rotor performance, i.e., figure-of-merit (FM). A quick-start procedure is described together with a statistical measure of FM convergence that reduces hover computations by 5 times, similar in computational work for forward flight. Furthermore, Cartesian adaptive mesh refinement in the rotor wake revealed a complex turbulent flow with worm-like structures of various scales. These turbulent worms, found computationally more than a decade ago, have recently been observed in a separate German Aerospace Center (DLR) experiment. Moreover, adaptive mesh refinement has been used to resolve the tip-vortex to its correct physical size.

Neal M. Chaderjian↗

Dynamic measurement of bulk modulus of dielectric materials using a microwave phase shift technique

A microwave Doppler shift technique was developed for measuring the dynamic bulk modulus of dielectric materials such as solid propellants. The system has a demonstrated time resolution on the order of milliseconds and a theoretical spatial resolution of a few microns. Accuracy of the technique is dependent on an accurate knowledge of the wavelength of the microwave in the sample being tested. Such measurement techniques are discussed. Preliminary tests with two solid propellants, one non-aluminized and one containing 16% aluminum, yielded reasonable, reproducible results. It was concluded that with refinements the technique holds promise as a practical means for obtaining accurate dynamic bulk modulus data over a variety of transient conditions.

Barker, B. J.↗

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

Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.

Casper, Jay↗

Multiplexed Holographic Data Storage in Bacteriorhodopsin

High density optical data storage, driven by the information revolution, remains at the forefront of current research areas. Much of the current research has focused on photorefractive materials (SBN and LiNbO3) and polymers, despite various problems with expense, durability, response time and retention periods. Photon echo techniques, though promising, are questionable due to the need for cryogenic conditions. Bacteriorhodopsin (BR) films are an attractive alternative recording medium. Great strides have been made in refining BR, and materials with storage lifetimes as long as 100 days have recently become available. The ability to deposit this robust polycrystalline material as high quality optical films suggests the use of BR as a recording medium for commercial optical disks. Our own recent research has demonstrated the suitability of BR films for real time spatial filtering and holography. We propose to fully investigate the feasibility of performing holographic mass data storage in BR. Important aspects of the problem to be investigated include various data multiplexing techniques (e.g. angle- amplitude- and phase-encoded multiplexing, and in particular shift-multiplexing), multilayer recording techniques, SLM selection and data readout using crossed polarizers for noise rejection. Systems evaluations of storage parameters, including access times, memory refresh constraints, erasure, signal-to-noise ratios and bit error rates, will be included in our investigations.

Mehrl, David J.↗

Determination of an Infrared Reference Frame

The mission of the European Space Agency's Hipparcos satellite was to provide high precision astrometric and photometric data for over 100,000 stars and somewhat lower precision data on over a million additional stars. These observations were made in the optical band (340 - 850 nm) and resulted in median precisions of better than 1 milli-arcsec (mas) in position and parallax and I mas/yr in proper motion. The Hipparcos frame is inertial within an estimated uncertainty of +/- 0.25 mas/yr. The other important wavelength band for astrometry has been the radio band. A reference frame has been defined based on absolute Very Long Baseline Interferometry (VLBI) positions of several hundreds of radio sources. There have been a number of projects to refine the link between the optical and radio reference frames. In recent years there has been increasing interest in the astrometry of infrared sources. For example, stars with circumstellar shells can be observed in both the optical and radio (e.g. SiO maser emission). The infrared emission from such sources is dominated by thermal dust emission. The spatial distribution of radiation in the three bands is expected to be quite different. In 1982 we proposed a number of infrared-selected sources (S(sub 10 micron) greater than or approx. equal to 100 J(sub y)) for inclusion in the Hipparcos Input Catalog. Most of those sources were observed, and in Sutton (1997) we report the optical astrometric results (position, proper motion, and parallax) for 87 such sources. These sources were selected to have reasonable agreement between infrared and optical positions, and so may be used as primary standards for future infrared astrometry. They are well distributed across the sky, but exhibit some bias towards the galactic plane and the northern hemisphere. We have also obtained a set of 1480 secondary standards (unpublished) A comparison between optical positions and 86 GHz SiO maser positions for 10 sources indicates coincidence at the 0.15 arcsec level, consistent with current uncertainties in the SiO positions. Improved measurements of the SiO maser positions are needed. VLBI data show that SiO maser spots can be distributed in broken rings of radii 10 - 30 mas, presumably centered on the stellar position. However the SiO flux may be distributed asymmetrically, in which case the flux-weighted center of the SiO emission can be different from the stellar position. We have also published a related paper on phase fluctuations in millimeter interferometry (Sutton and Hueckstaedt 1996). Radiometric monitoring of emission from atmospheric water vapor may be a method for correcting for such fluctuations. This is relevant to positional determinations at centimeter and millimeter wavelengths. The integration of visible, infrared, and radio astrometry is one of our principal goals.

Sutton, Edmund C.↗

CFD Predictions of Fluctuating Pressure Environments on NASA’s SLS Ascent Unsteady Aerodynamics Wind Tunnel Test

This paper details work that is being conducted on predicting fluctuating pressure environments on NASA launch vehicles using computational fluid dynamics (CFD). The accurate characterization of these aeroacoustics environments is necessary in the analysis of vehicle structural health and aerodynamic performance and in the vehicle design process. High-fidelity simulations of the unsteady flow over NASA’s SLS Ascent Unsteady Aerodynamics wind tunnel Test (AUAT) using the hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulations (LES) methodologies in Mississippi States’ Loci/CHEM solver are presented. The two available methods, Nelson-Nichols and dynamic hybrid RANS-LES (DHRL), are compared head-to-head on identical grids at two Mach numbers. At the high-subsonic Mach number, locally supersonic flow expanding over a shoulder induces a separation-reattachment system that is predicted poorly by the Nelson-Nichols method due to an under-prediction in unsteady content. The DHRL solution predicts the spatial extent and frequency domain response of this system well but under-predicts the peak noise levels due to a delay in the RANS-to-LES transition. At the low-supersonic Mach number, compression corner dynamics dominate the unsteady flow. While the Nelson Nichols solution captures this phenomenon well, another solution on a refined grid indicates that some model stress depletion is occurring. The DHRL solution again performs well at this Mach number and matches the wind tunnel data fairly closely, but perhaps most impressively, is fairly agnostic to increasing grid resolution. The findings in this paper are expected to be generally applicable and will be used to guide the prediction of aeroacoustics environments of other NASA launch vehicles.

CFD↗

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗

Parallel Adaptive Mesh Refinement Library

Parallel Adaptive Mesh Refinement Library (PARAMESH) is a package of Fortran 90 subroutines designed to provide a computer programmer with an easy route to extension of (1) a previously written serial code that uses a logically Cartesian structured mesh into (2) a parallel code with adaptive mesh refinement (AMR). Alternatively, in its simplest use, and with minimal effort, PARAMESH can operate as a domain-decomposition tool for users who want to parallelize their serial codes but who do not wish to utilize adaptivity. The package builds a hierarchy of sub-grids to cover the computational domain of a given application program, with spatial resolution varying to satisfy the demands of the application. The sub-grid blocks form the nodes of a tree data structure (a quad-tree in two or an oct-tree in three dimensions). Each grid block has a logically Cartesian mesh. The package supports one-, two- and three-dimensional models.

Mac-Neice, Peter↗

Evaluation of spatial, radiometric and spectral thematic mapper performance for coastal studies

An area along the southeastern shore of the Chesapeake Bay was subsetted from TM imagery. The subsetted image was then enhanced and classified using an ERDAS 400 system. Results obtained were compared with a chart showing the distribution of both Zolsters marina and Rupplia martime in the Vaucluse Shores and which supports a large community of SAV. Radiative transfer models describing the irradiance reflectance of a water column containing SAV are being refined. Radiative transfer theory was used to model upwelling radiance for an orbiting sensor viewing an estuarine environment. Upwelling radiance was calculated for a clear maritime atmosphere, an optically shallow estuary of either clear or turbid water, and one of three bottom types: vegetation, sand, or mud using TM bands 1, 2, and 3 and MSS bands 4 and 5. A spectral quality index was defined similar to the equation for apparent contrast and used to evaluate the relative effectiveness of TM and MSS bands in detecting submerged vegetation.

Klemas, V.↗

An evaluation of three spatial differencing schemes for the discrete ordinates method in participating media

Three popular spatial differencing schemes for the discrete ordinates method are examined for two-dimensional Cartesian coordinates system. These are a positive, the step, and the diamond schemes. Contrary to the common belief that negative intensities will not occur when fine spatial discretization is used with the diamond scheme, under certain conditions, the diamond scheme will produce negative intensities irrespective of the number of control volums employed. The positive scheme can produce physically unrealistic trends. The diamond and positive schemes are also capable of producing physically unrealistic overshoots. In absorbing-emitting or absorbing-emitting-scattering media, grid refinement can result in negative intensities when the diamond or positive scheme is used.

Chai, John C.↗