解二维扩散方程的高精度多重网格方法

本文提出了数值求解二维扩散方程的一种高精度加权平均隐式差分格式,理论分析结果表明其为无条件稳定的。为了克服传统迭代法在求解隐格式方面的困难,采用了多重网格算法,大大加快了迭代收敛速度,提高了求解效率。数值模拟了二维方腔内溶质的浓度扩散问题,数值实验结果验证了方法的精确性和可靠性。     

用代数多重网格法求解一维分裂格式的Euler方程

提出了代数多重网格法(AMG)的一种新算法。新算法改进了插值公式和粗网格方程,并把它应用到求解一维的分裂格式Euler方程。数值结果表明,对于具有高CFL条件数的Euler方程,代数多重网格法可以求解;对于Gaus-Seidel方法求解不能收敛的代数方程组,代数多重网格法求解可以收敛。新算法改进了代数多重网格法的收敛性和扩展了它的应用范围,数值结果表明了它的有效性和强壮性。     

隐-显积分因子间断Galerkin方法求解二维辐射扩散方程

提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是一种非常有效的数值算法.     

基于二阶格式的有限体积保正格式

从二阶线性格式出发,通过对法向通量进行重构,得到非线性两点通量,获得四面体网格上的单元中心型有限体积保正格式.该格式适用于求解间断和各向异性扩散系数问题.无需假设辅助未知量非负,避免了辅助未知量计算出负时"遇负置零"的人为处理方式;并且证明该格式在每个非线性Picard迭代步具有强保正性,即当源项和边界条件非负时,线性化格式的非平凡解是严格大于零的.数值算例验证该格式具有二阶收敛性且是保正的.     

杂质扩散方程的数值解法

研究杂质在等离子体中的扩散(非定态)时,提出了求解二阶非线性抛物型偏微分方程组的问题。对于这类方程组的数值求解,不少人进行过研究,但是在理论上至今还不够完善。仅就常用的差分方法而言,对于具体问题仍有采用哪种差分格式、如何线性化以及如何迭代等问题。这些问题的解决带有一定的经验性质。     

非定常对流扩散方程的高精度多重网格方法

由已有的求解定常对流扩散方程的高阶紧致差分格式出发,直接推导出了数值求解非定常对流扩散方程的一种高阶隐式紧致差分格式,其时间为二阶精度,空间为四阶精度,并且是无条件稳定的。为了加快传统迭代法在求解隐格式时在每一个时间步上的迭代收敛速度,采用了多重网格加速技术。数值实验结果验证了本文方法的高阶精度、高效性及高稳定性。     

带阻尼项定常Navier-Stokes方程的并行两水平有限元算法

基于两重网格离散和区域分解技术，提出数值求解带阻尼项定常Navier-Stokes方程的三种并行两水平有限元算法。其基本思想是首先在粗网格上求解完全的非线性问题，以获得粗网格解，然后在重叠的局部细网格子区域上并行求解Stokes、 Oseen和Newton线性化的残差问题，最后在非重叠的局部细网格子区域上校正近似解。数值算例验证了算法的有效性。     

修正的变分迭代法在四阶Cahn-Hilliard方程和BBM-Burgers方程中的应用

变分迭代法是一种基于变分原理,具有高数值精度的数值格式,目前已广泛应用于各类强非线性孤立波方程的数值求解中.本文利用修正的变分迭代法对两类非线性方程进行研究.该格式是对原数值方法的一种改进,即在变分项前引入了参数h.通过定义误差函数的离散二范数并在定义域内绘出h-曲线,从而确定出使误差达到最小的h,再返回原迭代过程进行求解.同时,参数的引入也扩大了原数值解的收敛域,在迭代次数一定的情况下达到了数值最优.在数值实验中,将上述结果应用于四阶的Cahn-Hilliard方程和BenjaminBona-Mahoney-Burgers方程.对于四阶的Cahn-Hilliard方程,普通的变分迭代法绝对误差在10~(-1)左右,经过修正后,绝对误差降为10~(-4),而且修正后的方法扩大了原数值解的收敛域.对于Benjamin-Bona-MahonyBurgers方程,利用带有辅助参数的变分迭代法将数值解的精度提高到10~(-3),对真解的逼近效果优于原始的变分迭代法.此数值方法也为其他强非线性孤立波微分方程的数值求解提供了方法和参考.     

随机Ginzburg-Landau方程的数值模拟

对随机Ginzburg-Landau方程进行数值研究,构造一个非线性差分格式和一个线性化差分格式.通过对确定性和随机Ginzburg-Landau方程的计算,表明所构造的格式具有较高的精度和较快的计算效率.对随机Ginzburg-Landau方程就噪声振幅的不同取值进行了数值模拟,并对由此引发的各种行为进行了描述.     

二维翼型跨音速绕流的数值模拟

本文结合高精度TVD格式的数值通量和时间进展多步法给出了一种求解定常流问题的数值方法。同时给出一些特殊处理来加快数值解的收敛速度。本文用以上方法计算了翼型跨音速绕流问题,结果表明此方法具有分辨率高,收敛速度较快之性质。     

一种最优控制数值计算方法

本文为解决非线性微分方程自治系统的周期轨道的数值计算,提出了一种最优控制的数值计算方法。其方法的实质是利用连续的打靶法,在非线性最小二乘公式的监护下,附加惩罚函数进行控制来实现的。并利用该法计算了典型的化学反应模型。     

Using the shooting method to solve boundary-value problems involving nonlinear coupled-wave equations

Using the shooting method, one of the numerical methods to solve two-point boundary-value problems, numerical solutions of the nonlinear coupled-wave equations in degenerate two-wave and four-wave mixing can be obtained. In this first part of the paper the general shooting method is described, and then applied to two-wave mixing in a reflection geometry. Computed results are presented in graphical form. Comparison between the shooting method and the direct numerical method [1] is made also. In the second part of the paper, numerical solutions for four-wave mixing in a reflection geometry will be given.     

Using the seventh-order numerical method to solve first-order nonlinear coupled-wave equations for degenerate two-wave and four-wave mixing

Using a new seventh-order numerical method [theO(h ⁷) method] for solving two-point boundary value problems, numerical solutions of the first-order nonlinear coupledwave equations for degenerate two-wave and four-wave mixing in a reflection geometry have been obtained. A computer program employing the Gauss-Jordan elimination technique has also been adopted to effectively solve the resultant large, sparse and unsymmetric matrix, obtained from theO(h ⁷) method and the Newton-Raphson iteration method. Numerical results from the computer calculations are presented graphically. A comparison between thisO(h ⁷) method and the shooting method, mainly from the viewpoint of computational efficiency, is also made.     

A fully second order implicit/explicit time integration technique for hydrodynamics plus nonlinear heat conduction problems

We present a fully second order implicit/explicit time integration technique for solving hydrodynamics coupled with nonlinear heat conduction problems. The idea is to hybridize an implicit and an explicit discretization in such a way to achieve second order time convergent calculations. In this scope, the hydrodynamics equations are discretized explicitly making use of the capability of well-understood explicit schemes. On the other hand, the nonlinear heat conduction is solved implicitly. Such methods are often referred to as IMEX methods [2], [1], [3]. The Jacobian-Free Newton Krylov (JFNK) method (e.g. [10], [9]) is applied to the problem in such a way as to render a nonlinearly iterated IMEX method. We solve three test problems in order to validate the numerical order of the scheme. For each test, we established second order time convergence. We support these numerical results with a modified equation analysis (MEA) [21], [20]. The set of equations studied here constitute a base model for radiation hydrodynamics.     

Effects of Thermal Radiation and Slip Mechanism on Mixed Convection Flow of Williamson Nanofluid Over an Inclined Stretching Cylinder

The current investigation highlights the mixed convection slip flow and radiative heat transport of uniformly electrically conducting Williamson nanofluid yield by an inclined circular cylinder in the presence of Brownian motion and thermophoresis parameter.A Lorentzian magnetic body force model is employed and magnetic induction effects are neglected.The governing equations are reduced to a system of nonlinear ordinary differential equations with associated boundary conditions by applying scaling group transformations.The reduced nonlinear ordinary differential equations are then solved numerically by Runge-Kutta-Fehlberg fifth-order method with shooting technique.The effects of magnetic field,Prandtl number,mixed convection parameter,buoyancy ratio parameter,Brownian motion parameter,thermophoresis parameter,heat generation/absorption parameter,mass transfer parameter,radiation parameter and Schmidt number on the skin friction coefficient and local Nusselt are analyzed and discussed.It is found that the velocity of the fluid decreases with decrease in curvature parameter,whereas it increases with mixed convection parameter.Further,the local Nusselt number decreases with an increase in the radiation parameter.The numerical comparison is also presented with the existing published results and found that the present results are in excellent agreement which also confirms the validity of the present methodology.     

MHD Flow and Heat Transfer Characteristics in a Casson Liquid Film Towards an Unsteady Stretching Sheet with Temperature-Dependent Thermal Conductivity

Theoretical and numerical outcomes of the non-Newtonian Casson liquid thin film fluid flow owing to an unsteady stretching sheet which exposed to a magnetic field, Ohmic heating and slip velocity phenomena is reported here. The non-Newtonian thermal conductivity is imposed and treated as it vary with temperature. The nonlinear partial differential equations governing the non-Newtonian Casson thin film fluid are simplified into a group of highly nonlinear ordinary differential equations by using an adequate dimensionless transformations. With this in mind, the numerical solutions for the ordinary conservation equations are found using an accurate shooting iteration technique together with the Runge-Kutta algorithm. The lineaments of the thin film flow and the heat transfer characteristics for the pertinent parameters are discussed through graphs. The results obtained here detect many concern for the local Nusselt number and the local skin-friction coefficient in which they may be beneficial for the material processing industries. Furthermore, in some special conditions, the present problem has an excellent agreement with previously published work.     

An efficient shooting method for fiber amplifiers and lasers

A simple shooting method is proposed for the design of distributed multi-pumped Raman fiber amplifiers (RFAs) and Yb-doped double-clad fiber lasers (DCFLs). Using the proposed method a distributed RFA with 10 pump sources and 1700 mW total input pump power is simulated and high-power Yb-doped DCFL rate equations are solved numerically. The numerical simulation results show that the proposed method has good convergence in the condition of increasing numbers and power of input pump sources.     

Computing for second harmonic generation in one-dimensional nonlinear photonic crystals

A numerical study of second harmonic generation (SHG) in one-dimensional nonlinear photonic crystals based on full nonlinear system of equations, implemented by a combination of the method of finite elements and fixed-point iterations, is reported. This model is derived from a nonlinear system of Maxwell’s equations, which partly overcomes the known shortcoming of some existing models relied on the undepleted-pump approximation. We derive a general solution of SHG in one-dimensional nonlinear photonic crystals structures. The convergence of our method is fast. Numerical simulations also show the conversion efficiency of SHG can be significantly enhanced when the frequencies of the fundamental wave are located at the photonic band edges or are assigned to the designed defect states.     

A Robust and Efficient Method for Steady State Patterns in Reaction-Diffusion Systems

An inhomogeneous steady state pattern of nonlinear reaction-diffusion equations with no-flux boundary conditions is usually computed by solving the corresponding time-dependent reaction-diffusion equations using temporal schemes. Nonlinear solvers (e.g., Newton's method) take less CPU time in direct computation for the steady state; however, their convergence is sensitive to the initial guess, often leading to divergence or convergence to spatially homogeneous solution. Systematically numerical exploration of spatial patterns of reaction-diffusion equations under different parameter regimes requires that the numerical method be efficient and robust to initial condition or initial guess, with better likelihood of convergence to an inhomogeneous pattern. Here, a new approach that combines the advantages of temporal schemes in robustness and Newton's method in fast convergence in solving steady states of reaction-diffusion equations is proposed. In particular, an adaptive implicit Euler with inexact solver (AIIE) method is found to be much more efficient than temporal schemes and more robust in convergence than typical nonlinear solvers (e.g., Newton's method) in finding the inhomogeneous pattern. Application of this new approach to two reaction-diffusion equations in one, two, and three spatial dimensions, along with direct comparisons to several other existing methods, demonstrates that AIIE is a more desirable method for searching inhomogeneous spatial patterns of reaction-diffusion equations in a large parameter space.