Local Two-Dimensional Splitting Schemes for 3D Suspended Matter Transport Problem Parallel Solution

A 3D model of suspended matter transport in coastal marine systems is considered, which takes into account many factors, including the hydraulic size or the rate of particle deposition, the propagation of suspended matter, sedimentation, the intensity of distribution of suspended matter sources, etc. The difference operators of diffusion transport in the horizontal and vertical directions for this problem have significantly different characteristic spatiotemporal scales of processes, as well as spectra. With typical sampling, applied to shallow-water systems in the South of Russia (the Sea of Azov, the Tsimlyansk reservoir), the steps in horizontal directions are 200-1000 meters, the coefficients of turbulent exchange (turbulent diffusion) are (103-104) m2/sec; in the vertical direction - - - steps of 0.1 m-1 m, and the coefficients of microturbulent exchange in the vertical — (0.1-1) m2/sec. If we focus on the use of explicit locally twodimensional - - - locally one-dimensional splitting schemes, then the permissible values of the time step for a two-dimensional problem will be about 10-100 seconds, and for a one-dimensional problem in the vertical direction - - - 0.1 – 1 sec. This motivates us to construct an additive locally-two-dimensional-locallyonedimensional splitting scheme in geometric directions. The paper describes a parallel algorithm that uses both explicit and implicit schemes to approximate the two-dimensional diffusion-convection problem in horizontal directions and the one-dimensional diffusion-convection problem in the vertical direction. The two-dimensional implicit diffusion-convection problem in horizontal directions is numerically solved by the adaptive alternating-triangular method. The numerical implementation of the one-dimensional diffusion-convection problem in the vertical direction is carried out by a sequential run-through method for a series of independent one-dimensional three-point problems in the vertical direction on a given layer. To increase the efficiency of parallel calculations, the decomposition of the calculated spatial grid and all grid data in one or two spatial directions - in horizontal directions-is also performed. The obtained algorithms are compared taking into account the permissible values of time steps and the actual time spent on performing calculations and exchanging information on each time layer.

A prescription for projectors to compute helicity amplitudes in D dimensions

This article discusses a prescription to compute polarized dimensionally regularized amplitudes, providing a recipe for constructing simple and general polarized amplitude projectors in D dimensions that avoids conventional Lorentz tensor decomposition and avoids also dimensional splitting. Because of the latter, commutation between Lorentz index contraction and loop integration is preserved within this prescription, which entails certain technical advantages. The usage of these D-dimensional polarized amplitude projectors results in helicity amplitudes that can be expressed solely in terms of external momenta, but different from those defined in the existing dimensional regularization schemes. Furthermore, we argue that despite being different from the conventional dimensional regularization scheme (CDR), owing to the amplitude-level factorization of ultraviolet and infrared singularities, our prescription can be used, within an infrared subtraction framework, in a hybrid way without re-calculating the (process-independent) integrated subtraction coefficients, many of which are available in CDR. This hybrid CDR-compatible prescription is shown to be unitary. We include two examples to demonstrate this explicitly and also to illustrate its usage in practice.

Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition

Research on seepage flow in the vadose zone has largely been driven by engineering and environmental problems affecting many fields of geotechnics, hydrology, and agricultural science. Mathematical modeling of the subsurface flow under unsaturated conditions is an essential part of water resource management and planning. In order to determine such subsurface flow, the two-dimensional (2D) Richards equation can be used. However, the computation process is often hampered by a high spatial resolution and long simulation period as well as the non-linearity of the equation. A new highly efficient and accurate method for solving the 2D Richards equation has been proposed in the paper. The developed algorithm is based on dimensional splitting, the result of which means that 1D equations can be solved more efficiently than as is the case with unsplit 2D algorithms. Moreover, such a splitting approach allows any algorithm to be used for space as well as time approximation, which in turn increases the accuracy of the numerical solution. The robustness and advantages of the proposed algorithms have been proven by two numerical tests representing typical engineering problems and performed for typical properties of soil.

TWO-DIMENSIONAL SPLITTING SCHEMES FOR HYPERBOLIC EQUATIONS

The article considers splitting schemes in geometric directions that approximate the initial-boundary value problem for p-dimensional hyperbolic equation by chain of two-dimensional-one-dimensional problems. Two ways of constructing splitting schemes are considered with an operator factorized on the upper layer, algebraically equivalent to the alternating direction scheme, and additive schemes of total approximation. For the first scheme, the restrictions on the shape of the region G at p=3 can be weakened in comparison with schemes of alternating directions, which are a chain of three-point problems on the upper time layer, the region can be a connected union of cylindrical regions with generators parallel to the axis OX3. In the second case, for a three-dimensional equation of hyperbolic type, an additive scheme is constructed, which is a chain «two-dimensional problem – one-dimensional problem» and approximates the original problem in a summary sense (at integer time steps). For the numerical implementation of the constructed schemes – the numerical solution of two-dimensional elliptic problems – one can use fast direct methods based on the Fourier algorithm, cyclic reduction methods for three-point vector equations, combinations of these methods, and other methods. The proposed two-dimensional splitting schemes in a number of cases turn out to be more economical in terms of total time expenditures, including the time for performing computations and exchanges of information between processors, compared to traditional splitting schemes based on the use of three-point difference problems for multiprocessor computing systems, with different structures of connections between processors type «ruler», «matrix», «cube», with universal switching.