Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable Eddington factor”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

ODD (Version 0.X)

The Orthogonal Diffusion Discretizations (ODD) library is a collection of low order (Diffusion, P1, Ray Tracing, and Variable Eddington Factor (VEF)) thermal radiation solvers for orthogonal structured grids driven by an application programing interface (API).

Cleveland, Mathew↗
NASA NTRS1980

Supercritical time-dependent accretion onto compact objects. I - Neutron stars

The time-dependent flow resulting from spherically symmetric, supercritical accretion onto a nonmagnetic neutron star is computed for accretion rates between 3 and 30. In addition, the fully time-dependent self-consistent set of coupled radiation-hydrodynamic equations governing such flow has been solved. Effects that have not previously been treated are taken into account, including separate energy equations for ions and electrons allowing for nonequilibrium, variable Eddington factors in the radiation transfer to account for anisotropy of the radiation field due to highly spherically extended envelopes, and convection of ions, electrons and photons in the accretion envelope. It is found that supercritical flow, for material falling from a radius of less than 10 to the 10th cm, is characterized by stable accretion with no evidence of stellar wind outflow at any phase of the evolution.

Klein, R. I.↗
NASA NTRS1988

Cosmic-ray hydrodynamics at relativistic shocks

A test particle description of first-order Fermi acceleration of cosmic rays at plane-parallel relativistic shocks is presented in which the pitch angle distribution is obtained explicitly and used to determine the pressure anisotropy or variable Eddington factor as a function of position. The Eddington factor is in turn used to find the hydrodynamical cosmic-ray sound speed, taking the anisotropy into account. For upstream flows speeds above c/sq rt 3, it is shown that the cosmic-ray sound speed exceeds the shock speed in a region just upstream of the shock and decreases montonically with distance further upstream, reaching the value c/sq rt 3 far from the shock. Thus, the shock succeeds in maintaining an anisotropy in the cosmic-ray distribution which is sufficient to keep the flow subsonic with respect to the cosmic-ray sound speed both downstream and in a finite region upstream. The astrophysical implications of these findings are discussed.

Kirk, J. G.↗
NASA NTRS1991

Near-critical spherical accretion by neutron stars - General relativistic treatment

Time-independent spherical accretion by a neutron star is studied using general relativistic radiation hydrodynamics. Numerical integrations of the flow equations are presented. These show that when the luminosity is sufficiently close to (but below) the Eddington limit, the flow velocity increases with decreasing radius far from the neutron star, reaches a maximum at an intermediate radius, and decreases at small radii. A large fraction of the binding energy of the flow is transferred to the radiation through scattering before the flow strikes the surface of the neutron star. Following Miller's treatment of accretion at luminosites near the Eddington limit (which neglected general relativistic effects), analytic approximations for the decelerating phase of the flow's velocity profile are derived. The dependence of the solutions on the variable Eddington factor prescription chosen to close the radiation moment equations is also examined.

Park, Myeong-Gu↗
NASA NTRS1991

The shock breakout in SN1987A modelled with the time-dependent radiative transfer

A fully implicit high-order scheme was developed for the independent multi-group radiative transfer coupled with implicit hydrodynamics. The application of this scheme to the SN1987A explosion shows that shortly after the shock breakout a dense shell forms. A gasdynamic code describing the time dependent radiation transport in the multi-group approximation with variable Eddington factors was developed. Results and other aspects of the project are presented.

Blinnikov, S. I.↗
DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗
NASA NTRS1994

Neutrino bursts from gamma-ray bursts

If gamma-ray bursts originate at cosmological distances, as strongly indicated by the results from Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory (CGRO), then ultrarelativistic ejecta are the likely consequence of the highly super-Eddington luminosity of the sources. If the energy injection rate varies with time, then the Lorentz factor of the wind also varies, and the shells of ejected matter collide with each other. The collisions between baryons produce pions which decay into high-energy photons, electrons, electron positron pairs, and neutrino pairs. The bulk Lorentz factor of approximately 300 is required if our model is to be compatible with the observed millisecond variability. The strongest gamma-ray bursts are observed to deliver approximately 10(exp -4) ergs/sq cm in 100-200 keV photons. In our scenario more energy may be delivered in a neutrino burst. Typical neutrinos may be approximately 30 GeV if the protons have a Maxwellian energy distribution, and up to approximately TeV if the protons have a power-law distribution. Such neutrino bursts are close to the detection limit of the DUMAND II experiment.

Paczynski, Bohdan↗
DOE OSTI2020

Spectroscopic classification, variability, and SED of the Fermi -detected CSS 3C 286: the radio-loudest NLS1 galaxy?

ABSTRACT 3C 286 is a well-known calibrator source in radio astronomy. It is also one of very few compact steep-spectrum (CSS) sources detected in γ-rays. Here, we perform a detailed spectroscopic and variability analysis and present the first quasi-simultaneous optical to X-ray spectral energy distribution (SED) in order to reveal physical mechanisms that dominate its emission at different wavelengths, and arrive at a reliable optical source classification. The first main result of our study reveals several pitfalls when applying simple broad- or narrow-line Seyfert 1 (BLS1 or NLS1) classification criteria that only look at the [O iii]–Hβ complex. [O iii] and Hβ can be dominated by the same outflow components, in which case full width at half-maximum (Hβ) is no reliable classification criterion, and extinction by intrinsic or intervening material can make the highest velocity Hβ component undetectable. After careful combination of all information from ultraviolet (UV)–optical spectra along with multiwavelength data, we confirm that 3C 286 can be classified as NLS1 galaxy, with line properties and supermassive black hole mass (of order 108 M⊙ and accreting near the Eddington limit) close to the BLS1 regime, making it an important borderline object. The quasi-simultaneous SED taken with Swift shows a sharp rise in the optical–UV, implying the presence of a strong accretion-disc component with extreme UV excess, consistent with emission-line diagnostics. Finally, we report the discovery of X-ray variability of 3C 286, plausibly dominated by jet emission, and variable by at least a factor of ∼4. This result suggests to exercise caution when using 3C 286 as radio calibrator in high-resolution radio very long baseline interferometry observations.

Yao, Su↗
NASA NTRS1985

Solar radiance models for determination of ERBE scanner filter factor

Shortwave spectral radiance models for use in the spectral correction algorithms for the ERBE Scanner Instrument are provided. The required data base was delivered to the ERBe Data Reduction Group in October 1984. It consisted of two sets of data files: (1) the spectral bidirectional angular models and (2) the spectral flux modes. The bidirectional models employ the angular characteristics of reflection by the Earth-atmosphere system and were derived from detailed radiance calculations using a finite difference model of the radiative transfer process. The spectral flux models were created through the use of a delta-Eddington model to economically simulate the effects of atmospheric variability. By combining these data sets, a wide range of radiances may be approximated for a number of scene types.

Arduini, R. F.↗
NASA NTRS2003

A Chandra X-Ray Study of NGC 1068 IL the Luminous X-Ray Source Population

We present an analysis of the compact X-ray source population in the Seyfert 2 galaxy NGC 1068, imaged with a approx. 50 ks Chandra observation. We find a total of 84 compact sources on the S3 chip, of which 66 are located within the 25.0 B-mag/arcsec isophote of the galactic disk of NGC 1068. Spectra have been obtained for the 21 sources with at least 50 counts and modeled with both multicolor disk blackbody and power-law models. The power-law model provides the better description of the spectrum for 18 of these sources. For fainter sources, the spectral index has been estimated from the hardness ratio. Five sources have 0.4 - 8 keV intrinsic luminosities greater than 10(exp 39)ergs/ s, assuming that their emission is isotropic and that they are associated with NGC 1068. We refer to these sources as intermediate-luminosity X-ray objects (ISOs). If these five sources are X-ray binaries accreting with luminosities that are both sub-Eddington and isotropic, then the implied source masses are approx greater than 7 solar mass, and so they are inferred to be black holes. Most of the spectrally modeled sources have spectral shapes similar to Galactic black hole candidates. However, the brightest compact source in NGC 1068 has a spectrum that is much harder than that found in Galactic black hole candidates and other ISOs. The brightest source also shows large amplitude variability on both short-term and long-term timescales, with the count rate possibly decreasing by a factor of 2 in approx. 2 ks during our Chundra observation, and the source flux decreasing by a factor of 5 between our observation and the grating observations taken just over 9 months later. The ratio of the number of sources with luminosities greater than 2.1 x 10(exp 38) ergs/s in the 0.4 - 8 keV band to the rate of massive (greater than 5 solar mass) star formation is the same, to within a factor of 2, for NGC 1068, the Antennae, NGC 5194 (the main galaxy in M51), and the Circinus galaxy. This suggests that the rate of production of X-ray binaries per massive star is approximately the same for galaxies with currently active star formation, including "starbursts."

Smith, David A.↗
NASA NTRS1991

Radius-luminosity and mass-luminosity relationships for active galactic nuclei

Broad-line region (BLR) sizes derived from spectral variability and BLR line widths are used to directly derive the mass (M) of the central objects of ten active galactic nuclei (AGNs) in a uniform manner. It is shown that the luminosity-weighted C IV 1549-emitting BLR radius (R) correlates with the bolometric luminosity L(Bol) and is consistent with R about sq rt L(Bol). The measurements also permit a verification of the Dibai mass-luminosity (M-L) relationship (previously derived indirectly). It is found that L(Bol) is proportional to M exp (1.1 + or - 0.3). It is found that the efficiency factor epsilon, defined as the ratio of L(Bol) to the Eddington luminosity increases from 0.03 in the low-luminosity Seyferts up to 0.06 in the most luminous objects in the sample.

Koratkar, Anuradha P.↗
NASA NTRS2020

Discovery of Thermonuclear (Type 1) X-ray Bursts in the X-ray Binary Swift J1858.6-0814 Observed with NICER and NuSTAR

Swift J1858.6–0814 is a recently discovered X-ray binary notable for extremely strong variability (by factors of >100 in soft X-rays) in its discovery state. We present the detection of five thermonuclear (Type I) X-ray bursts from Swift J1858.6–0814, implying that the compact object in the system is a neutron star (NS). Some of the bursts show photospheric radius expansion, so their peak flux can be used to estimate the distance to the system. The peak luminosity, and hence distance, can depend on several system parameters; for the most likely values, a high inclination and a helium atmosphere, D = 12.8+0.8−0.6 kpc, although systematic effects allow a conservative range of 9–18 kpc. Before one burst, we detect a QPO at 9.6 ± 0.5 mHz with a fractional rms amplitude of 2.2 ± 0.2 per cent (0.5–10 keV), likely due to marginally stable burning of helium; similar oscillations may be present before the other bursts but the light curves are not long enough to allow their detection. We also search for burst oscillations but do not detect any, with an upper limit in the best case of 15 per cent fractional amplitude (over 1–8 keV). Finally, we discuss the implications of the NS accretor and this distance on other inferences which have been made about the system. In particular, we find that Swift J1858.6–0814 was observed at super-Eddington luminosities at least during bright flares during the variable stage of its outburst.

D J K Buisson↗