Scattering Polarization and Hanle Effect in Stellar Atmospheres with Horizontal Inhomogeneities

Scattering of light from an anisotropic source produces linear polarization in spectral lines and the continuum. In the outer layers of a stellar atmosphere the anisotropy of the radiation field is typically dominated by the radiation escaping away, but local horizontal fluctuations of the physical conditions may also contribute, distorting the illumination and hence, the polarization pattern. Additionally, a magnetic field may perturb and modify the line scattering polarization signals through the Hanle effect. Here, we study such symmetry-breaking effects. We develop a method to solve the transfer of polarized radiation in a scattering atmosphere with weak horizontal fluctuations of the opacity and source functions. It comprises linearization (small opacity fluctuations are assumed), reduction to a quasi-planeparallel problem through harmonic analysis, and numerical solution by generalized standard techniques. We apply this method to study scattering polarization in atmospheres with horizontal fluctuations in the Planck function and opacity. We derive several very general results and constraints from considerations on the symmetries and dimensionality of the problem, and we give explicit solutions of a few illustrative problems of especial interest. For example, we show (a) how the amplitudes of the fractional linear polarization signals change when considering increasingly smaller horizontal atmospheric inhomogeneities, (b) that in the presence of such inhomogeneities even a vertical magnetic field may modify the scattering line polarization, and (c) that forward scattering polarization may be produced without the need of an inclined magnetic field. These results are important to understand the physics of the problem and as benchmarks for multidimensional radiative transfer codes.Comment: 27 pages, 13 figures, to appear in Ap

Optimized Schwarz waveform relaxation for Primitive Equations of the ocean

In this article we are interested in the derivation of efficient domain decomposition methods for the viscous primitive equations of the ocean. We consider the rotating 3d incompressible hydrostatic Navier-Stokes equations with free surface. Performing an asymptotic analysis of the system with respect to the Rossby number, we compute an approximated Dirichlet to Neumann operator and build an optimized Schwarz waveform relaxation algorithm. We establish the well-posedness of this algorithm and present some numerical results to illustrate the method

High performance computing of explicit schemes for electrofusion jointing process based on message-passing paradigm

The research focused on heterogeneous cluster workstations comprising of a number of CPUs in single and shared architecture platform. The problem statements under consideration involved one dimensional parabolic equations. The thermal process of electrofusion jointing was also discussed. Numerical schemes of explicit type such as AGE, Brian, and Charlies Methods were employed. The parallelization of these methods were based on the domain decomposition technique. Some parallel performance measurement for these methods were also addressed. Temperature profile of the one dimensional radial model of the electrofusion process were also given

Algorithms and data structures for adaptive multigrid elliptic solvers

Adaptive refinement and the complicated data structures required to support it are discussed. These data structures must be carefully tuned, especially in three dimensions where the time and storage requirements of algorithms are crucial. Another major issue is grid generation. The options available seem to be curvilinear fitted grids, constructed on iterative graphics systems, and unfitted Cartesian grids, which can be constructed automatically. On several grounds, including storage requirements, the second option seems preferrable for the well behaved scalar elliptic problems considered here. A variety of techniques for treatment of boundary conditions on such grids are reviewed. A new approach, which may overcome some of the difficulties encountered with previous approaches, is also presented