2,277 research outputs found
IRIS: A Generic Three-Dimensional Radiative Transfer Code
We present IRIS, a new generic three-dimensional (3D) spectral radiative
transfer code that generates synthetic spectra, or images. It can be used as a
diagnostic tool for comparison with astrophysical observations or laboratory
astrophysics experiments. We have developed a 3D short-characteristic solver
that works with a 3D nonuniform Cartesian grid. We have implemented a piecewise
cubic, locally monotonic, interpolation technique that dramatically reduces the
numerical diffusion effect. The code takes into account the velocity gradient
effect resulting in gradual Doppler shifts of photon frequencies and subsequent
alterations of spectral line profiles. It can also handle periodic boundary
conditions. This first version of the code assumes Local Thermodynamic
Equilibrium (LTE) and no scattering. The opacities and source functions are
specified by the user. In the near future, the capabilities of IRIS will be
extended to allow for non-LTE and scattering modeling. IRIS has been validated
through a number of tests. We provide the results for the most relevant ones,
in particular a searchlight beam test, a comparison with a 1D plane-parallel
model, and a test of the velocity gradient effect. IRIS is a generic code to
address a wide variety of astrophysical issues applied to different objects or
structures, such as accretion shocks, jets in young stellar objects, stellar
atmospheres, exoplanet atmospheres, accretion disks, rotating stellar winds,
cosmological structures. It can also be applied to model laboratory
astrophysics experiments, such as radiative shocks produced with high power
lasers.Comment: accepted for publication in A&A; 17 pages, 9 figures, 2 table
Asymptotics for Hermite-Pade rational approximants for two analytic functions with separated pairs of branch points (case of genus 0)
We investigate the asymptotic behavior for type II Hermite-Pade approximation
to two functions, where each function has two branch points and the pairs of
branch points are separated. We give a classification of the cases such that
the limiting counting measures for the poles of the Hermite-Pade approximants
are described by an algebraic function of order 3 and genus 0. This situation
gives rise to a vector-potential equilibrium problem for three measures and the
poles of the common denominator are asymptotically distributed like one of
these measures. We also work out the strong asymptotics for the corresponding
Hermite-Pade approximants by using a 3x3 Riemann-Hilbert problem that
characterizes this Hermite-Pade approximation problem.Comment: 102 pages, 31 figure
From discretization to regularization of composite discontinuous functions
Discontinuities between distinct regions, described by different equation sets, cause difficulties for PDE/ODE solvers. We present a new algorithm that eliminates integrator discontinuities through regularizing discontinuities. First, the algorithm determines the optimum switch point between two functions spanning adjacent or overlapping domains. The optimum switch point is determined by searching for a “jump point” that minimizes a discontinuity between adjacent/overlapping functions. Then, discontinuity is resolved using an interpolating polynomial that joins the two discontinuous functions.
This approach eliminates the need for conventional integrators to either discretize and then link discontinuities through generating interpolating polynomials based on state variables or to reinitialize state variables when discontinuities are detected in an ODE/DAE system. In contrast to conventional approaches that handle discontinuities at the state variable level only, the new approach tackles discontinuity at both state variable and the constitutive equations level. Thus, this approach eliminates errors associated with interpolating polynomials generated at a state variable level for discontinuities occurring in the constitutive equations.
Computer memory space requirements for this approach exponentially increase with the dimension of the discontinuous function hence there will be limitations for functions with relatively high dimensions. Memory availability continues to increase with price decreasing so this is not expected to be a major limitation
- …