Related papers: Radiative transfer in decomposed domains
A fast method is proposed for solving the high frequency Helmholtz equation. The building block of the new fast method is an overlapping source transfer domain decomposition method for layered medium, which is an extension of the source…
We simulate convection near the solar surface, where the continuum optical depth is of order unity. Hence, to determine the radiative heating and cooling in the energy conservation equation, we must solve the radiative transfer equation…
We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core…
Domain decomposition is a technique used to reduce memory overhead on large neutron transport problems. Currently, the optimal load-balanced processor allocation for these domains is typically determined through small-scale simulations of…
We present a numerical method and computer code to calculate the radiative transfer and excitation of molecular lines. Formulating the Monte Carlo method from the viewpoint of cells rather than photons allows us to separate local and…
We present describe a new computer code that solves the radiative transfer problem on multi-resolution grids. If the cloud model is from an MHD simulation on a regular cartesian grid, criteria based for example on local density or velocity…
With tens of petaflops supercomputers already in operation and exaflops machines expected to appear within the next 10 years, efficient parallel computational methods are required to take advantage of such extreme-scale machines. In this…
We present a general formalism for computing self-consistent, numerical solutions to the time-dependent radiative transfer equation in low velocity, multi-level ions undergoing radiative interactions. Recent studies of time-dependent…
A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick…
We shall derive and propose several efficient overlapping domain decomposition methods for solving some typical linear inverse problems, including the identiffication of the flux, the source strength and the initial temperature in second…
This work proposes a discretization of the acoustic wave equation with possibly oscillatory coefficients based on a superposition of discrete solutions to spatially localized subproblems computed with an implicit time discretization. Based…
We present a ray-tracing technique for radiative transfer modeling of complex three-dimensional (3D) structures which include dense regions of high optical depth like in dense molecular clouds, circumstellar disks, envelopes of evolved…
We develop innovative algorithms for solving the strong-constraint formulation of four-dimensional variational data assimilation in large-scale applications. We present a space-time decomposition approach that employs domain decomposition…
In this article we propose a numerical algorithm to compute the intensity and polarization of a polychromatic electromagnetic radiation crossing a medium with graded refractive index and modeled by the Vector Radiative Refractive Transfer…
A methodology to reduce the computational cost of time domain computations of eddy currents problems is proposed and implemented in a parallel computing environment. It is based on the modal decomposition of the current density and it is…
A domain decomposition method is proposed based on carefully chosen impedance transmission operators for a hybrid formulation of the eddy current problem. Preliminary analysis and numerical results are provided in the spherical case showing…
Solving the null geodesic equations for a ray of light is a difficult task even considering a stationary spacetime. The problem becomes even more difficult if the electromagnetic signal propagates through a flowing optical medium. Indeed,…
We develop a representation of reverse-time migration in terms of Fourier integral operators the canonical relations of which are graphs. Through the dyadic parabolic decomposition of phase space, we obtain the solution of the wave equation…
This paper presents a multiscale methodology for efficient unsteady conjugate heat transfer simulations. The solid domain is modelled by coupling a global representation of the temperature field, based on the eigenfunctions of the unsteady…
Context. The numerical modeling of the generation and transfer of polarized radiation is a key task in solar and stellar physics research and has led to a relevant class of discrete problems that can be reframed as linear systems. In order…