Testing weighted splitting schemes on a one-column transport-chemistry model

In many transport-chemistry models, a huge system of ODE’s of the advection-diffusion-reaction type has to be integrated in time. Typically, this is done with the help of operator splitting. Rosenbrock schemes combined with approximate matrix factorization (ROS-AMF) are an alternative to operator splitting which does not suffer from splitting errors. However, implementation of ROS-AMF schemes often requires serious changes in the code. In this paper we test another classical second order splitting introduced by Strang in 1963, which, unlike the popular Strang splitting, seemed to be forgotten and rediscovered recently (partially due to its intrinsic parallellism). This splitting, called symmetrically weighted sequential (SWS) splitting, is simple and straightforward to apply, independent of the order of the operators and has an operator-level parallelism. In the experiments, the SWS scheme compares favorably to the Strang splitting, but is less accurate than ROS-AMF

Hot-star wind models with magnetically split line blanketing

Fraction of hot stars posses strong magnetic fields that channel their radiatively driven outflows. We study the influence of line splitting in the magnetic field (Zeeman effect) on the wind properties. We use our own global wind code with radiative transfer in the comoving frame to understand the influence of the Zeeman splitting on the line force. We show that the Zeeman splitting has a negligible influence on the line force for magnetic fields that are weaker than about 100~kG. This means that the wind mass-loss rates and terminal velocities are not affected by the magnetic line splitting for magnetic fields as are typically found on the surface of nondegenerate stars. Neither have we found any strong flux variability that would be due to the magnetically split line blanketing.Comment: 4 pages, accepted for publication in Astronomy & Astrophysic

Geometrical relations between space time block code designs and complexity reduction

In this work, the geometric relation between space time block code design for the coherent channel and its non-coherent counterpart is exploited to get an analogue of the information theoretic inequality $I(X;S)\le I((X,H);S)$ in terms of diversity. It provides a lower bound on the performance of non-coherent codes when used in coherent scenarios. This leads in turn to a code design decomposition result splitting coherent code design into two complexity reduced sub tasks. Moreover a geometrical criterion for high performance space time code design is derived.Comment: final version, 11 pages, two-colum

Investigation of upwind, multigrid, multiblock numerical schemes for three dimensional flows. Volume 1: Runge-Kutta methods for a thin layer Navier-Stokes solver

A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, upwind numerical techniques, multigrid acceleration, and multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available: van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multi-grid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with multiple topologies. The results shown in this dissertation were chosen to validate the computer code and display its geometric flexibility, which is provided by the multi-block structure

GNSS signal acquisition in the presence of sign transitions

The next generation of Global Navigation Satellite Systems (GNSS), such as Galileo [1] and GPS modernization [2], will use signals with equal code and bit periods, which will introduce a potential sign transition in each segment of the signal processed in the acquisition block. If FFT is used to perform the correlations a sign transition occurring within the integration time may cause a splitting of the main peak of the Cross Ambiguity Function (CAF) into two smaller lobes along the Doppler shift axis [3]. In this paper a method to overcome the possible impairments due to the lobe splitting is proposed and validated by simulatio

A Multi-dimensional Code for Isothermal Magnetohydrodynamic Flows in Astrophysics

We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference method on an Eulerian grid, called the total variation diminishing (TVD) scheme. Recipes for building the one-dimensional IMHD code, including the normalized right and left eigenvectors of the IMHD Jacobian matrix, are presented. Then, we have extended the one-dimensional code to a multi-dimensional IMHD code through a Strang-type dimensional splitting. In the multi-dimensional code, an explicit cleaning step has been included to eliminate non-zero $\nabla\cdot B$ at every time step. To estimate the proformance of the code, one- and two-dimensional IMHD shock tube tests, and the decay test of a two-dimensional Alfv\'{e}n wave have been done. As an example of astrophysical applications, we have simulated the nonlinear evolution of the two-dimensional Parker instability under a uniform gravity.Comment: Accepted for publication in ApJ, using aaspp4.sty, 22 text pages with 10 figure

A Rate-Splitting Based Bound-Approaching Transmission Scheme for the Two-User Symmetric Gaussian Interference Channel with Common Messages

This paper is concerned with a rate-splitting based transmission strategy for the two-user symmetric Gaussian interference channel that contains common messages only. Each transmitter encodes its common message into multiple layers by multiple codebooks that drawn from one separate code book, and transmits the superposition of the messages corresponding to these layers; each receiver decodes the messages from all layers of the two users successively. Two schemes are proposed for decoding order and optimal power allocation among layers respectively. With the proposed decoding order scheme, the sum-rate can be increased by rate-splitting, especially at the optimal number of rate-splitting, using average power allocation in moderate and weak interference regime. With the two proposed schemes at the receiver and the transmitter respectively, the sum-rate achieves the inner bound of HK without time-sharing. Numerical results show that the proposed optimal power allocation scheme with the proposed decoding order can achieve significant improvement of the performance over equal power allocation, and achieve the sum-rate within two bits per channel use (bits/channel use) of the sum capacity

Evolution program for parton densities with perturbative heavy flavor boundary conditions

A new code for the scale evolution of modified-minimal-subtraction-scheme parton densities is described. Through next-to-leading order the program uses exact splitting functions. In next-to-next-to-leading order approximate splitting functions are used. For efficiency the program includes analytical results for the evaluation of the weights required for the integrations over the longitudinal momentum fractions of the partons. It also incorporates the operator matrix elements required for the matching conditions across heavy flavor thresholds in higher order perturbation theory. The more efficient handling of the weights implies that the code is faster than similar evolution codes in all modes of operation. The program is written in the C programming language.Comment: LaTeX with 4 figure

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 $\Lambda$ 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