三维非线性多层渗流方程耦合系统的差分方法

对三维多层非线性渗流方程耦合系统提出适合并行计算的二阶迎风分数步差分格式, 利用变分形式、能量方法、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶的l²误差估计. 该方法已成功地应用到油资源运移聚集数值模拟的生产实际中.     

三维非线性对流扩散问题的数值方法在渗流力学的应用

对三维非线性对流扩散问题提出一类适合并行计算的二阶迎风分数步差分格式,采用分数步技术,将三维问题化为连续解3个一维问题计算．利用变分形式、能量方法、差分算子乘积交换性、高阶差分算子的分解、微分方程先验估计的理论和技巧,得到收敛性的最佳阶的误差估计．该方法已成功的应用油资源运移聚集渗流力学数值模拟计算、海水入侵预测和防治的工程实践中．     

非线性渗流耦合系统的数值方法及其应用

对非线性二相渗流耦合系统提出修正迎风分数步差分格式,利用变分形式、归纳法假定、高阶差分算子的分解、微分方程的先验估计理论和技巧,得到收敛性的最佳阶误差估计．该方法已成功应用于防治海水入侵主要工程后效预测与调控模式的生产实践中．     

多层非线性渗流耦合系统的特征分数步差分方法

对多层非线性渗流耦合系统提出适合并行计算的特征分数步差分格式, 利用变分形式、能量方法、粗细网格配套、分片双二次插值、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶的l₂误差估计. 该方法已成功的应用到多层油资源评估的生产实际中.     

多层渗流方程组合系统的迎风分数步长差分方法和应用

对多层渗流方程组合系统提出适合并行计算的二阶和一阶两类迎风分数步长差分格式,利用变分形式、能量方法、二维和三维格式的配套、隐显格式的相互结合、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧,对二阶格式得到收敛性的最佳阶的L²误差估计. 对一阶格式亦得到收敛性的L²误差估计. 该方法已成功地应用到多层油资源运移聚集的评估生产实际中,得到了很好的数值模拟结果.     

三维渗流耦合系统动边值问题迎风差分方法的理论和应用

对多层渗流方程耦合系统动边值问题,提出适合并行计算的两类迎风差分格式,利用区域变换、变分形式、能量方法、隐显格式的相互结合、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧,得到收敛性的l~2误差估计.该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中,得到了很好的数值模拟结果.     

非线性渗流耦合系统动边值问题二阶迎风分数步差分方法

对多层非线性渗流方程耦合系统三维动边值问题, 提出适合并行计算的一类二阶迎风分数步差分格式, 利用区域变换、变分形式、能量方法、隐显格式的相互结合、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶l² 误差估计. 该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中, 得到了很好的数值模拟结果.     

三维强化采油渗流耦合系统隐式迎风分数步差分方法的收敛性分析

对强化(化学)采油数值模拟中的渗流耦合系统问题,本文提出隐式迎风分数步差分格式,利用变分形式、能量方法、差分算子乘积交换性理论、高阶差分算子的分解、微分方程先验估计的理论和技巧,得到最佳阶l2模误差估计.该方法已成功应用到强化采油数值模拟的生产实践中,并且得到了很好的数值模拟效果.     

渗流耦合系统动边值问题特征分数步差分方法

对多层渗流方程耦合系统动边值问题,提出适合并行计算的一类特征分数步差分格式,利用区域变换、变分形式、粗细网格配套、乘积型叁二次插值、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧,得到收敛性的最佳阶l2误差估计．该方法已成功地应用到多层油资源运移聚集的资源评估生产实际中,得到了很好的数值模拟结果．     

可压缩两相驱动问题的迎风差分格式及其理论分析

对石油渗流中的可压缩两相渗流驱动问题，提出新的二阶迎风差分格式，利用变分形式，能量方法、微分方程先验估计的理论和技巧，得到最佳阶ι^2误差估计，该方法已应用到海水入侵、防治和油资源评估的数值模拟中。     

The finite difference method for the three-dimensional nonlinear coupled system of dynamics of fluids in porous media

For the three-dimensional coupled system of multilayer dynamics of fluids in porous media, the second-order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in l2 norm are derived to determine the error in the second-order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources.     

A Linearized Spectral-Galerkin Method for Three-Dimensional Riesz-Like Space Fractional Nonlinear Coupled Reaction-Diffusion Equations

In this paper, we establish a novel fractional model arising in the chemical reaction and develop an efficient spectral method for the three-dimensional Riesz-like space fractional nonlinear coupled reaction-diffusion equations. Based on the backward difference method for time stepping and the Legendre-Galerkin spectral method for space discretization, we construct a fully discrete numerical scheme which leads to a linear algebraic system. Then a direct method based on the matrix diagonalization approach is proposed to solve the linear algebraic system, where the cost of the algorithm is of a small multiple of $N^4$ ($N$ is the polynomial degree in each spatial coordinate) flops for each time level. In addition, the stability and convergence analysis are rigorously established. We obtain the optimal error estimate in space, and the results also show that the fully discrete scheme is unconditionally stable and convergent of order one in time. Furthermore, numerical experiments are presented to confirm the theoretical claims. As the applications of the proposed method, the fractional Gray-Scott model is solved to capture the pattern formation with an analysis of the properties of the fractional powers.     

Theory and Application of Characteristic Finite Element Domain Decomposition Procedures for Coupled System of Dynamics of Fluids in Porous Media

For a coupled system of multiplayer dynamics of fluids in porous media,the characteristic finiteelement domain decomposition procedures applicable to parallel arithmetic are put forward.Techniques suchas calculus of variations,domain decomposition,characteristic method,negative norm estimate,energy methodand the theory of prior estimates are adopted.Optimal order estimates in L~2 norm are derived for the error inthe approximate solution.     

The upwind finite difference fractional steps method for nonlinear coupled systems

For nonlinear coupled system of multilayer dynamics of fluids in porous media, a second‐order upwind finite‐difference fractional‐steps scheme applicable to parallel arithmetic are put forward, and two‐ and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, multiplicative commutation rule of difference operators, decomposition of high‐order difference operators, and prior estimates are adopted. Optimal order estimates in l² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of migration‐accumulation of oil resources. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

种间作用的二维非线性动力系统模型及其数值模拟试验研究

系统探讨了种间制约作用的二维非线性动力系统的机理与过程,从生物力学角度扩展了Lotka-Volterra模型,建立了包含种间作用与反作用的二维非线性自治、非自治动力系统新模型,分析了该动力系统平衡点的稳定性、周期解的存在性与稳定性,并对其动力学过程进行了数值模拟试验研究.结果表明两生物种群间作用效率的大小、作用系数与反作用系数变化的同向与异向性和作用力时变等对动力系统稳定性都有一定影响,种间作用效率变大或变小,以及作用系数与反作用系数的异向变化都导致动力系统的不稳定,致使两种生物难以生存,而相互作用力时变则对系统的稳定起促进作用.     

Mathematical Analysis and Numerical Simulation of a Nonlinear Thermoelastic System

Abstract

In this paper, we give a theoretical and numerical analysis of a model for small vertical vibrations of an elastic membrane coupled with a heat equation and subject to linear mixed boundary conditions. We establish the existence, uniqueness, and a uniform decay rate for global solutions to this nonlinear non-local thermoelastic coupled system with linear boundary conditions. Furthermore, we introduced a numerical method based on finite element discretization in a spatial variable and finite difference scheme in time which results in a nonlinear system to be solved by Newton’s method. Numerical experiments for one-dimensional and two-dimensional cases are presented in order to estimate the rate of convergence of the numerical solution that confirm our analysis and show the efficiency of the method.