17,955 research outputs found
A 3D radiative transfer framework: IV. spherical & cylindrical coordinate systems
We extend our framework for 3D radiative transfer calculations with a
non-local operator splitting methods along (full) characteristics to spherical
and cylindrical coordinate systems. These coordinate systems are better suited
to a number of physical problems than Cartesian coordinates. The scattering
problem for line transfer is solved via means of an operator splitting (OS)
technique. The formal solution is based on a full characteristics method. The
approximate operator is constructed considering nearest neighbors
exactly. The code is parallelized over both wavelength and solid angle using
the MPI library. We present the results of several test cases with different
values of the thermalization parameter for the different coordinate systems.
The results are directly compared to 1D plane parallel tests. The 3D results
agree very well with the well-tested 1D calculations.Comment: A&A, in pres
Numerical Solution of the Expanding Stellar Atmosphere Problem
In this paper we discuss numerical methods and algorithms for the solution of
NLTE stellar atmosphere problems involving expanding atmospheres, e.g., found
in novae, supernovae and stellar winds. We show how a scheme of nested
iterations can be used to reduce the high dimension of the problem to a number
of problems with smaller dimensions. As examples of these sub-problems, we
discuss the numerical solution of the radiative transfer equation for
relativistically expanding media with spherical symmetry, the solution of the
multi-level non-LTE statistical equilibrium problem for extremely large model
atoms, and our temperature correction procedure. Although modern iteration
schemes are very efficient, parallel algorithms are essential in making large
scale calculations feasible, therefore we discuss some parallelization schemes
that we have developed.Comment: JCAM, in press. 28 pages, also available at
ftp://calvin.physast.uga.edu:/pub/preprints/CompAstro.ps.g
A 3D radiative transfer framework IX. Time dependence
Context. Time-dependent, 3D radiation transfer calculations are important for
the modeling of a variety of objects, from supernovae and novae to simulations
of stellar variability and activity. Furthermore, time-dependent calculations
can be used to obtain a 3D radiative equilibrium model structure via relaxation
in time. Aims. We extend our 3D radiative transfer framework to include direct
time dependence of the radiation field; i.e., the
terms are fully considered in the solution of radiative transfer problems.
Methods. We build on the framework that we have described in previous papers in
this series and develop a subvoxel method for the
terms. Results. We test the implementation by comparing the 3D results to our
well tested 1D time dependent radiative transfer code in spherical symmetry. A
simple 3D test model is also presented. Conclusions. The 3D time dependent
radiative transfer method is now included in our 3D RT framework and in
PHOENIX/3D.Comment: A&A in press, 7 pages, 14 figure
A 3D radiative transfer framework: VII. Arbitrary velocity fields in the Eulerian frame
A solution of the radiative-transfer problem in 3D with arbitrary velocity
fields in the Eulerian frame is presented. The method is implemented in our 3D
radiative transfer framework and used in the PHOENIX/3D code. It is tested by
comparison to our well- tested 1D co-moving frame radiative transfer code,
where the treatment of a monotonic velocity field is implemented in the
Lagrangian frame. The Eulerian formulation does not need much additional memory
and is useable on state-of-the-art computers, even large-scale applications
with 1000's of wavelength points are feasible
General Relativistic Radiative Transfer
We present a general method to calculate radiative transfer including
scattering in the continuum as well as in lines in spherically symmetric
systems that are influenced by the effects of general relativity (GR). We
utilize a comoving wavelength ansatz that allows to resolve spectral lines
throughout the atmosphere. The used numerical solution is an operator splitting
(OS) technique that uses a characteristic formal solution. The bending of
photon paths and the wavelength shifts due to the effects of GR are fully taken
into account, as is the treatment of image generation in a curved spacetime. We
describe the algorithm we use and demonstrate the effects of GR on the
radiative transport of a two level atom line in a neutron star like atmosphere
for various combinations of continuous and line scattering coefficients. In
addition, we present grey continuum models and discuss the effects of different
scattering albedos on the emergent spectra and the determination of effective
temperatures and radii of neutron star atmospheres
Preliminary Spectral Analysis of SN 1994I
We present optical spectra of the Type Ic supernova 1994I in M51 and
preliminary non-LTE analysis of the spectra. Our models are not inconsistent
with the explosions of C+O cores of massive stars. While we find no direct
evidence for helium in the optical spectra, our models cannot rule out small
amounts of helium. More than 0.1~\msol\ of helium seems unlikely.Comment: LaTeX, MN style, psfig, and natbib substyles, 7 pages, 4 figures, to
appear in MNRAS. Postscript file available from
http://www.nhn.uoknor.edu/~baro
- …