Approximate Formulation and Numerical Solution for Hypersingular Boundary Integral Equations in Plane Elasticity

Based on the fact that the singular boundary integrals in the sense of Cauchy principal value can be represented approximately by the mean values of two companion nearly singular boundary integrals, a vary general approach was developed in the paper. In the approach, the approximate formulation before discretization was constructed to cope with the difficulties encountered in the corner treatment in the formulations of hypersingular boundary integral equations. This makes it possible to solve the hypersingular boundary integral equation numerically in a non regularized form and in a local manner by using conforming C 0 quadratic boundary elements and standard Gaussian quadratures similar to those employed in the conventional displacement BIE formulations. The approximate formulation is very convenient to use because the corner information is comprised naturally in the representations of those approximate integrals. Numerical examples in plane elasticity show that with the present approach, the compatible or better results can be achieved in comparison with those of the conventional BIE formulations.

THE STABILITY AND CONVERGENCE OF THE FINITE ANALYTIC METHOD FOR THE NUMERICAL SOLUTION OF CONVECTIVE DIFFUSION EQUATION

In this paper we make a close study of the finite analytic method by means of the maximum principles in differential equations and give the proof of the stability and convergence of the finite analytic method.

Lie group analysis, numerical and non-traveling wave solutions for the (2+1)-dimensional diffusion–advection equation with variable coefficients

In this paper, the variable-coefficient diffusion-advection （DA） equation, which arises in modeling various physical phenomena, is studied by the Lie symmetry approach. The similarity reductions are derived by determining the complete sets of point symmetries of this equation, and then exact and numerical solutions are reported for the reduced second-order nonlinear ordinary differential equations. Further, an extended （Gl/G）-expansion method is applied to the DA equation to construct some new non-traveling wave solutions.

Qualitative Properties and Numerical Solution of the Kolmogorov-Fisher Type Biological Population Task with Double Nonlinear Diffusion

In the present work we study the global solvability of the Kolmogorov-Fisher type biological population task with double nonlinear diffusion and qualitative properties of the solution of the task based on the self-similar analysis. In additional, in this paper we consider the model of two competing population with dual nonlinear cross-diffusion.

Numerical Simulation of Diffusion Type Traffic Flow Model Using Second-Order Lax-Wendroff Scheme Based on Exponential Velocity Density Function

In order to control traffic congestion, many mathematical models have been used for several decades. In this paper, we study diffusion-type traffic flow model based on exponential velocity density relation, which provides a non-linear second-order parabolic partial differential equation. The analytical solution of the diffusion-type traffic flow model is very complicated to approximate the initial density of the Cauchy problem as a function of x from given data and it may cause a huge error. For the complexity of the analytical solution, the numerical solution is performed by implementing an explicit upwind, explicitly centered, and second-order Lax-Wendroff scheme for the numerical solution. From the comparison of relative error among these three schemes, it is observed that Lax-Wendroff scheme gives less error than the explicit upwind and explicit centered difference scheme. The numerical, analytical analysis and comparative result discussion bring out the fact that the Lax-Wendroff scheme with exponential velocity-density relation of diffusion type traffic flow model is suitable for the congested area and shows a better fit in traffic-congested regions.

On the Construction of Analytic-Numerical Approximations for a Class of Coupled Differential Models in Engineering

In this paper, a method to construct an analytic-numerical solution for homogeneous parabolic coupled systems with homogeneous boundary conditions of the type ut = Auxx, A1u(o,t) + B1ux(o,t) = 0, A2u(1,t) + B2ux(1,t) = 0, ot>0, u (x,0) = f(x), where A is a positive stable matrix and A1, B1, B1, B2,? ?are arbitrary matrices for which the block matrix is non-singular, is proposed.

STRONG COMPACTNESS OF APPROXIMATE SOLUTIONS TO DEGENERATE ELLIPTIC-HYPERBOLIC EQUATIONS WITH DISCONTINUOUS FLUX FUNCTION

Under a non-degeneracy condition on the nonlinearities we show that sequences of approximate entropy solutions of mixed elliptic-hyperbolic equations are strongly precompact in the general case of a Caratheodory flux vector. The proofs are based on deriving localization principles for H-measures associated to sequences of measurevalued functions. This main result implies existence of solutions to degenerate parabolic convection-diffusion equations with discontinuous flux. Moreover, it provides a framework in which one can prove convergence of various types of approximate solutions, such as those generated by the vanishing viscosity method and numerical schemes.

Studies on fluid model for numerical simulation of gas discharges in color plasma displays

The fluid models of gas discharge in alternating current plasma display panel (AC PDP) cell are discussed. From the Boltzmann equation, the hydrodynamic equations are derived, but this model consumes much computa- tional time for simulation. The drift-diffusion approximation model and the local field approximation model are ob- tained to simplify the numerical computation, and the approximation conditions of these two models are discussed in detail. The drift-diffusion approximation model gives more satisfactory result for PDP simulation, and the expression of energy balance equation is given completely in this model.

Numerical Approximation of Stochastic Theta Method for Random Periodic Solution of Stochastic Differential Equations

In this paper,we make use of stochastic theta method to study the existence of the numerical approximation of random periodic solution.We prove that the error between the exact random periodic solution and the approximated one is at the 1/4 order time step in mean sense when the initial time tends to∞.

Approximate derivative-dependent functional variable separation for quasi-linear diffusion equations with a weak source

By using the approximate derivative-dependent functional variable separation approach, we study the quasi-linear diffusion equations with a weak source ut = （A（u）Ux）x ＋ eB（u, Ux）. A complete classification of these perturbed equations which admit approximate derivative-dependent functional separable solutions is listed. As a consequence, some approxi- mate solutions to the resulting perturbed equations are constructed via examples.

An Approximate Riemann Solver for Advection-Diffusion Based on the Generalized Riemann Problem

We construct an approximate Riemann solver for scalar advection-diffusion equations with piecewise polynomial initial data.The objective is to handle advection and diffusion simultaneously to reduce the inherent numerical diffusion produced by the usual advection flux calculations.The approximate solution is based on the weak formulation of the Riemann problem and is solved within a space-time discontinuous Galerkin approach with two subregions.The novel generalized Riemann solver produces piecewise polynomial solutions of the Riemann problem.In conjunction with a recovery polynomial,the Riemann solver is then applied to define the numerical flux within a finite volume method.Numerical results for a piecewise linear and a piecewise parabolic approximation are shown.These results indicate a reduction in numerical dissipation compared with the conventional separated flux calculation of advection and diffusion.Also,it is shown that using the proposed solver only in the vicinity of discontinuities gives way to an accurate and efficient finite volume scheme.

Numerical investigation of Dufour and Soret effects on unsteadyMHD natural convection flow past vertical plate embedded innon-Darcy porous medium

The Dufour and Soret effects on the unsteady twodimensional magnetonyaro dynamics （MHD） doublediffusive free convective flow of an electrically conducting fluid past a vertical plate embedded in a nonDarcy porous medium are investigated numeri cally. The governing nonlinear dimensionless equations are solved by an implicit finite difference scheme of the CrankNicolson type with a tridiagonal matrix manipulation. The effects of various parameters entering into the problem on the unsteady dimension less velocity, temperature, and concentration profiles are studied in detail. Furthermore, the time variation of the skin friction coefficient, the Nusselt number, and the Sherwood number is presented and analyzed. The results show that the unsteady velocity, tem perature, and concentration profiles are substantially influenced by the Dufour and Soret effects. When the Dufour number increases or the Soret number decreases, both the skin friction and the Sherwood number decrease, while the Nusselt number increases. It is found that, when the magnetic parameter increases, the velocity and the temperature decrease in the boundary layer.

Time-frequency analysis of Li solid-phase diffusion in spherical active particles under typical discharge modes

Li transient concentration distribution in spherical active material particles can affect the maximum power density and the safe operating regime of the electric vehicles(EVs). On one hand, the quasiexact/exact solution obtained in the time/frequency domain is time-consuming and just as a reference value for approximate solutions;on the other hand, calculation errors and application range of approximate solutions not only rely on approximate algorithms but also on discharge modes. For the purpose to track the transient dynamics for Li solid-phase diffusion in spherical active particles with a tolerable error range and for a wide applicable range, it is necessary to choose optimal approximate algorithms in terms of discharge modes and the nature of active material particles. In this study, approximation methods,such as diffusion length method, polynomial profile approximation method, Padé approximation method,pseudo steady state method, eigenfunction-based Galerkin collocation method, and separation of variables method for solving Li solid-phase diffusion in spherical active particles are compared from calculation fundamentals to algorithm implementation. Furthermore, these approximate solutions are quantitatively compared to the quasi-exact/exact solution in the time/frequency domain under typical discharge modes, i.e., start-up, slow-down, and speed-up. The results obtained from the viewpoint of time-frequency analysis offer a theoretical foundation on how to track Li transient concentration profile in spherical active particles with a high precision and for a wide application range. In turn, optimal solutions of Li solid diffusion equations for spherical active particles can improve the reliability in predicting safe operating regime and estimating maximum power for automotive batteries.

Maximum modulus principle estimates for one dimensional fractional diffusion equation

We present scheme I for solving one-dimensional fractional diffusion equation with variable coefficients based on the maximum modulus principle and two Grunwald approxima- tions. Scheme II is obtained by using classic Crank-Nicolson approximations in order to improve the time convergence. Schemes are proved to be unconditionally stable and second-order accuracy in spatial grid size for the problem with order of fractional derivative belonging to [（√17- 1）/2, 2] using the maximum modulus principle. A numerical example is given to confirm the theoretical analysis result.

New Analytical Study of the Effects Thermo-Diffusion, Diffusion-Thermo and Chemical Reaction of Viscous Fluid on Magneto Hydrodynamics Flow in Divergent and Convergent Channels

In this paper, the magneto hydrodynamic (MHD) flow of viscous fluid in a channel with non-parallel plates is studied. The governing partial differential equation was transformed into a system of dimensionless non-similar coupled ordinary differential equation. The transformed conservations equations were solved by using new algorithm. Basically, this new algorithm depends mainly on the Taylor expansion application with the coefficients of power series resulting from integrating the order differential equation. Results obtained from new algorithm are compared with the results of numerical Range-Kutta fourth-order algorithm with help of the shooting algorithm. The comparison revealed that the resulting solutions were excellent agreement. Thermo-diffusion and diffusion-thermo effects were investigated to analyze the behavior of temperature and concentration profile. Also the influences of the first order chemical reaction and the rate of mass and heat transfer were studied. The computed analytical solution result for the velocity, temperature and concentration distribution with the effect of various important dimensionless parameters was analyzed and discussed graphically.

MOVING BOUNDARY PROBLEM FOR DIFFUSION RELEASE OF DRUG FROM A CYLINDER POLYMERIC MATRIX

An approximate analytical solution of moving boundary problem for diffusion release of drug from a cylinder polymeric matrix was obtained by use of refined integral method. The release kinetics has been analyzed for non-erodible matrices with perfect sink condition. The formulas of the moving boundary and the fractional drug release were given. The moving boundary and the fractional drug release have been calculated at various drug loading levels, mid the calculated results were in good agreement with those of experiments. The comparison of the moving boundary in spherical, cylinder, planar matrices has been completed. An approximate formula for estimating the available release time was presented. These results are useful for the clinic experiments. This investigation provides a new theoretical tool for studying the diffusion release of drug from a cylinder polymeric matrix and designing the controlled released drug.

Numerical Solutions for Optimal Control of Stochastic Kolmogorov Systems

This work is concerned with controlled stochastic Kolmogorov systems.Such systems have received much attention recently owing to the wide range of applications in biology and ecology.Starting with the basic premise that the underlying system has an optimal control,this paper is devoted to designing numerical methods for approximation.Different from the existing literature on numerical methods for stochastic controls,the Kolmogorov systems take values in the first quadrant.That is,each component of the state is nonnegative.The work is designing an appropriate discrete-time controlled Markov chain to be in line with(locally consistent)the controlled diffusion.The authors demonstrate that the Kushner and Dupuis Markov chain approximation method still works.Convergence of the numerical scheme is proved under suitable conditions.

区域注浆扰动下渗流场−化学场演化及耦合作用

在淮北煤田普遍采用地面定向钻高压注浆技术治理煤系底板岩溶水害。注浆后渗流场补径排路径发生变化,且注浆实施中注浆高压以及浆液析水会在一定时期内改变受注目标含水层地下水渗流场和化学场的分布特征。以淮北煤田桃园煤矿为研究对象,利用Feflow软件,构建了区域注浆扰动下目标层渗流场与化学场(简称“双场”)演化模型,探讨了“双场”耦合机制,选择受环境影响较小的Cl−为模拟因子,开展了“双场”耦合作用下溶质运移影响因素分析,在此基础上构建了区域注浆扰动下溶质扩散预测模型。研究结果表明:识别验证后渗流场模型中,模拟水位与实测水位误差小于3 m的点占71.9%,模拟效果较好;识别验证后化学场模型中,Ⅱ4采区溶质运移模拟效果较好,Ⅱ1采区模拟值偏低约14.4 mg/L,误差约为实际值的6.6%,低于10%,总体模拟效果较好。渗透系数相对较大的Ⅱ4采区和Ⅱ2采区Cl−较容易运移扩散;而渗透系数较小的Ⅱ1采区岩石致密,渗透性较差,模型运行50 a期间,其基本以“滞水”状态存在。基于注浆扰动下Ⅱ4采区局部渗流场演化特征,认为溶质运移主要受渗透系数、弥散度、水力梯度、渗透流速、注浆时间、浆液相对密度等参数控制,并发现在注浆结束后18~22 a Cl−质量浓度达到峰值,随后Cl−质量浓度开始降低,约在40 a后达到区域注浆之前的平衡状态。基于多次设参运行获得数据,建立了“双场”耦合作用下溶质扩散的预测模型,误差率小于10%的数据占比达81.4%,说明所建立的溶质扩散预测模型基本可靠。研究可为区域注浆扰动下受注目标含水层水环境演化研究以及煤矿水害预测预警提供科学依据,具有重要的理论和实践意义。

Java Simulation of Au Diffusion in Si Affected by Vacancies and Self-Interstitials: Partial Differential Equations

A Java program in a GUI environment has been developed for the numerical solution of basic partial differential equations and applied to Au diffusion in Si affected by vacancies and self-interstitials. Text fields of selected parameters for the calculation are set on the display, and the calculation starts by checking the start button after putting values on the text fields. The calculated results are plotted immediately after the finish of the calculation as the concentration profiles of substitutional Au, interstitial Au, vacancies and self-interstitials, and their diffusion can be presented immediately, resulting in the identification of the diffusion mechanism. By changing the values of the text fields, new results can be represented immediately. The diffusion of Au in Si can be simulated correctly and easily by this program. Results from the program for one set of conditions are shown, including images produced on the display.

一类新型自适应反扩散近似Riemann求解器及其应用

对于包含激波、剪切层等复杂结构的流动问题,为了精确模拟剪切层等精细结构,且保证激波计算的稳定性,必须采用低耗散且强鲁棒的数值通量方法。传统的HLL近似Riemann求解器的耗散性较大,Roe、HLLEM和HLLC等近似Riemann求解器在计算某些含有强激波的物理问题时会出现非物理解,容易导致不稳定。针对这一问题,本文在Riemann求解器中通过合理设计反扩散矩阵,发展了一类具有自适应反扩散的新型Riemann求解器,并将其应用到高阶加权紧致格式,实现了高阶精度求解。通过典型数值算例验证了新型方法的计算精度和稳定性,结果表明本文提出的新型自适应反扩散Riemann求解器克服了传统Riemann求解器的缺陷,既能准确识别剪切层等精细结构,又能保证激波解的稳定性。