Study of a millimeter-wave squint indirect holographic algorithm suitable for imaging with large field-of-view

In this paper a millimeter-wave （MMW） squint indirect holographic method is presented, which is suitable for imaging with a large field-of-view. The proposed system employs the squint operation mode to remove the background and twin- image interferences, which achieves a similar effect to off-axis holography but leaves out the large-aperture quasi-optical component. The translational scanning manner enables a large field of view and ensures the image uniformity, which is difficult to realize in off-axis holography. In addition, a corresponding imaging algorithm for the presented scheme is developed to reconstruct the image from the recorded hologram. Some imaging results on typical objects, obtained with electromagnetic simulation, demonstrate good performance of the imaging scheme and validate the effectiveness of the image reconstruction algorithm.     

A cell functional minimization scheme for parabolic problem

We are interested in a robust and accurate finite volume scheme for 2-D parabolic problems derived from the cell functional minimization approach. The scheme has a local stencil, is locally conservative, treats discontinuity rigorously and leads to a symmetric positive definite linear system. Since the scheme has both cell centered unknowns and cell edge unknowns, the computational cost is an issue and a parallel algorithm is then suggested based on nonoverlapping domain decomposition approach. The interface condition is of the Dirichlet–Robin type and has a parameter λ. By choosing this parameter properly, the convergence of the iteration process could be sped up. Numerical results for linear and nonlinear problems demonstrate the good performance of the cell functional minimization scheme and its parallel version on distorted meshes.     

A class of deformational flow test cases for linear transport problems on the sphere

A class of new benchmark deformational flow test cases for the two-dimensional horizontal linear transport problems on the sphere is proposed. The scalar field follows complex trajectories and undergoes severe deformation during the simulation; however, the flow reverses its course at half-time and the scalar field returns to its initial position and shape. This process makes the exact solution available at the end of the simulation, and facilitates assessment of the accuracy of the underlying transport scheme. A procedure to eliminate possible cancellations of errors when the flow reverses is proposed.     

A well‐balanced upstream flux‐splitting finite‐volume scheme for shallow‐water flow simulations with irregular bed topography

This study extends the upstream flux‐splitting finite‐volume (UFF) scheme to shallow water equations with source terms. Coupling the hydrostatic reconstruction method (HRM) with the UFF scheme achieves a resultant numerical scheme that adequately balances flux gradients and source terms. The proposed scheme is validated in three benchmark problems and applied to flood flows in the natural/irregular river with bridge pier obstructions. The results of the simulations are in satisfactory agreement with the available analytical solutions, experimental data and field measurements. Comparisons of the present results with those obtained by the surface gradient method (SGM) demonstrate the superior stability and higher accuracy of the HRM. The stability test results also show that the HRM requires less CPU time (up to 60%) than the SGM. The proposed well‐balanced UFF scheme is accurate, stable and efficient to solve flow problems involving irregular bed topography. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical simulation of core convection by a multi-layer semi-implicit spherical spectral method

A semi-implicit multi-layer spherical spectral method for simulating stellar core convection is described. The fully compressible three-dimensional hydrodynamic equations with rotation and energy generation are solved. Prognostic variables are expressed as finite sums of spherical harmonics in the horizontal directions and handled by the finite difference method in the radial direction. The stratified approximation is used to simplify the nonlinearity to quadratic. A multi-layer scheme is employed to overcome the time step problem arising from shrinking grid sizes in the physical space near the center of the star. Despite of the different spectral truncations in different layers, round-off conservation of the total mass and total angular momentum of the whole domain can be maintained, and were confirmed numerically. The code is parallelized; with 12 processors the speedup factor is about 9. The solutions of model core convection with and without rotation are discussed.     

APPROXIMATING COMMON FIXED POINTS OF NEARLY ASYMPTOTICALLY NONEXPANSIVE MAPPINGS

We use an iteration scheme to approximate common fixed points of nearly asymptotically nonexpansive mappings.We generalize corresponding theorems of [1] to the case of two nearly asymptotically nonexpansive mappings and those of [9] not only to a larger class of mappings but also with better rate of convergence.     

Normal halogen dependence of 13C NMR chemical shifts of halogenomethanes revisited at the four‐component relativistic level

The ‘Normal Halogen Dependence’ of ¹³C NMR chemical shifts in the series of halogenomethanes is revisited at the four‐component relativistic level. Calculations of ¹³C NMR chemical shifts of 70 halogenomethanes have been carried out at the density functional theory (DFT) and MP2 levels with taking into account relativistic effects using the four‐component relativistic theory of Dirac‐Coulomb within the different computational methods (4RPA, 4OPW91) and hybrid computational schemes (MP2 + 4RPA, MP2 + 4OPW91). The most efficient computational protocols are derived for practical purposes. Relativistic shielding effect reaches as much as several hundreds of ppm for heavy halogenomethanes, and to account for this effect in comparison with experiment at the qualitative level, relativistic Dyall's basis sets of triple‐zeta quality or higher are to be used within the framework of the four‐component relativistic theory taking into account solvent effects. Relativistic geometrical optimization (as compared with the non‐relativistic level) is essential for the molecules containing at least two iodines at one carbon atom. Copyright © 2016 John Wiley & Sons, Ltd.     

求解Schrdinger方程的高阶紧致差分格式

基于Richardson外推法提出了一种求解Schrdinger方程的高阶紧致差分方法.该方法首先利用二阶微商的四阶精度紧致差分逼近公式对原方程进行求解,然后利用Richardson外推技术外推一次,得到了Schrdinger方程具有O(r~4+h~4)精度的数值解.通过Fourier分析方法证明了该格式是无条件稳定的.数值实验验证了该方法的高阶精度及有效性.     

The parameter uniform numerical method for singularly perturbed parabolic reaction–diffusion problems on equidistributed grids

This article presents the study of singularly perturbed parabolic reaction–diffusion problems with boundary layers. To solve these problems, we use a modified backward Euler finite difference scheme on layer adapted nonuniform meshes at each time level. The nonuniform meshes are obtained by equidistribution of a positive monitor function, which involves the second-order spatial derivative of the singular component of the solution. The equidistributing monitor function at each time level allows us to use this technique to non-linear parabolic problems. The truncation error and the stability analysis are obtained. Parameter–uniform error estimates are derived for the numerical solution. To support the theoretical results, numerical experiments are carried out.     

A phase transition for the heights of a fragmentation tree

We provide information about the asymptotic regimes for a homogeneous fragmentation of a finite set. We establish a phase transition for the asymptotic behavior of the shattering times, defined as the first instants when all the blocks of the partition process have cardinality less than a fixed integer. Our results may be applied to the study of certain random split trees. © 2010 Wiley Periodicals, Inc. Random Struct. Alg., 39, 247‐274, 2011