对流扩散问题的Crank-Nicolson差分-流线扩散法

1　引　言Streamline- Diffusion method (SD方法 )是近年来 Hughes和 Brooks提出的一种求解定常的对流占优和对流扩散问题的人工粘性有限元方法[1 ] ,[2 ] ,它具有标准有限元方法的高阶精度特点和人工粘性 Galerkin方法的稳定性特点 ,因此越来越受到人们的重视 .现在 ,SD方法已被推广到 Euler方程和 Navier- Stokes方程等发展型对流扩散问题[3 ] [4] ,但是常常采用时空有限元 [3 ] [5] ,这样能把时间和空间的精度很好地统一起来 ,却增大了数值计算的复杂性 ,基于此 [6 ]对非线性的对流占优扩散问题提出一种 Finite Difference- Strea…     

调水工程线路方案优选的灰色系统模型及其应用

调水工程为解决区域间水资源分布不均衡问题提供了一种有效途径 .调水工程线路方案的优选问题涉及社会、经济、环境、技术等多方面的因素 ,是一个复杂的多目标决策问题 .本文建立了其评价指标体系 ,应用 AHP法和信息熵方法相结合确定评价指标权重 ,建立了调水工程线路方案优选的灰色系统评价模型 .算例显示方法的有效性 .     

基于光频调节的干涉型光纤水听器相位补偿检测方法

提出了一种基于光源频率调节的干涉型光纤水听器主动相位补偿的信号检测方法。详细介绍了该方法的基本检测原理,对信号解调误差进行了理论分析与仿真。为了验证该方法的可行性,搭建了实验系统,并编写了实时的信号采集、处理程序,对某一干涉型光纤水听器的声压灵敏度进行了测试。在频带20 Hz-1.3 kHz上,平均声压灵敏度为-162.2 dB(0 dB=1 rad/μPa),波动小于±0.8 dB,与采用相位载波调制解调方法测量的结果基本吻合。实验结果证明该方法是有效的。由于传感部分不含有源器件,便于实现全光纤化,且解调算法简单、检测频带宽,该方法能被广泛应用到各种干涉型光纤传感器的信号检测当中。     

磁单极与磁洛伦兹力

概述了磁单极概念的历史发展,从洛伦兹变换出发,利用电磁场张量和四维力的协变性以及电荷相对论不变,直接证明了运动磁单极受磁洛伦兹力,建议了一个磁洛伦兹力的验证方案,并用磁洛伦兹力公式导出狄拉克电荷量子化条件．证明了磁洛伦兹力公式具有与库仑定律相同的精确度．     

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.