By using the revision of the momentum for a particle moving with high velocity and by investigating the famous Bucherer's experiment of an electron deflecting with high velocity in the electromagnetic fields in 19...By using the revision of the momentum for a particle moving with high velocity and by investigating the famous Bucherer's experiment of an electron deflecting with high velocity in the electromagnetic fields in 1908, the paper determines that mass of the electron with high velocity is still to observe the law of conservation of mass.展开更多
The purpose is to reestablish the coupled conservation laws, the local conservation equations and the jump conditions of mass and inertia for polar continuum theories. In this connection the new material derivatives o...The purpose is to reestablish the coupled conservation laws, the local conservation equations and the jump conditions of mass and inertia for polar continuum theories. In this connection the new material derivatives of the deformation gradient, the line element, the surface element and the volume element were derived and the generalized Reynolds transport theorem was presented. Combining these conservation laws of mass and inertia with the balance laws of momentum, angular momentum and energy derived in our previous papers of this series, a rather complete system of coupled basic laws and principles for polar continuum theories is constituted on the whole. From this system the coupled nonlocal balance equations of mass, inertia, momentum, angular momentum and energy may be obtained by the usual localization.展开更多
The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differenti...The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.展开更多
The physical model is described by a seepage coupled system for simulating numerically three-dimensional chemical oil recovery, whose mathematical description includes three equations to interpret main concepts. The p...The physical model is described by a seepage coupled system for simulating numerically three-dimensional chemical oil recovery, whose mathematical description includes three equations to interpret main concepts. The pressure equation is a nonlinear parabolic equation, the concentration is defined by a convection-diffusion equation and the saturations of different components are stated by nonlinear convection-diffusion equations. The transport pressure appears in the concentration equation and saturation equations in the form of Darcy velocity, and controls their processes. The flow equation is solved by the conservative mixed volume element and the accuracy is improved one order for approximating Darcy velocity. The method of characteristic mixed volume element is applied to solve the concentration, where the diffusion is discretized by a mixed volume element method and the convection is treated by the method of characteristics. The characteristics can confirm strong computational stability at sharp fronts and it can avoid numerical dispersion and nonphysical oscillation. The scheme can adopt a large step while its numerical results have small time-truncation error and high order of accuracy. The mixed volume element method has the law of conservation on every element for the diffusion and it can obtain numerical solutions of the concentration and adjoint vectors. It is most important in numerical simulation to ensure the physical conservative nature. The saturation different components are obtained by the method of characteristic fractional step difference. The computational work is shortened greatly by decomposing a three-dimensional problem into three successive one-dimensional problems and it is completed easily by using the algorithm of speedup. Using the theory and technique of a priori estimates of differential equations, we derive an optimal second order estimates in 12 norm. Numerical examples are given to show the effectiveness and practicability and the method is testified as a powerful tool to solve the important problems.展开更多
In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and t...In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order.The concentration equation is solved by the method of mixed finite element,where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method.The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation.In actual computations the characteristics adopts a large time step without any loss of accuracy.The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux.In order to derive the optimal 3/2-order error estimate in L^(2) norm,a post-processing technique is included in the approximation to the scalar unknowns.Numerical experiments are illustrated finally to validate theoretical analysis and efficiency.This method can be used to solve such an important problem.展开更多
In this paper,the authors discuss a three-dimensional problem of the semiconductor device type involved its mathematical description,numerical simulation and theoretical analysis.Two important factors,heat and magneti...In this paper,the authors discuss a three-dimensional problem of the semiconductor device type involved its mathematical description,numerical simulation and theoretical analysis.Two important factors,heat and magnetic influences are involved.The mathematical model is formulated by four nonlinear partial differential equations(PDEs),determining four major physical variables.The influences of magnetic fields are supposed to be weak,and the strength is parallel to the z-axis.The elliptic equation is treated by a block-centered method,and the law of conservation is preserved.The computational accuracy is improved one order.Other equations are convection-dominated,thus are approximated by upwind block-centered differences.Upwind difference can eliminate numerical dispersion and nonphysical oscillation.The diffusion is approximated by the block-centered difference,while the convection term is treated by upwind approximation.Furthermore,the unknowns and adjoint functions are computed at the same time.These characters play important roles in numerical computations of conductor device problems.Using the theories of priori analysis such as energy estimates,the principle of duality and mathematical inductions,an optimal estimates result is obtained.Then a composite numerical method is shown for solving this problem.展开更多
In this paper, we propose a class of stable finite difference schemes for the initial-boundary value problem of the Cahn-Hilliard equation. These schemes are proved to inherit the total mass conservation and energy di...In this paper, we propose a class of stable finite difference schemes for the initial-boundary value problem of the Cahn-Hilliard equation. These schemes are proved to inherit the total mass conservation and energy dissipation in the discrete level. The dissipation of the total energy implies boundness of the numerical solutions in the discrete H1 norm. This in turn implies boundedness of the numerical solutions in the maximum norm and hence the stability of the difference schemes. Unique existence of the numerical solutions is proved by the fixed-point theorem. Convergence rate of the class of finite difference schemes is proved to be O(h2 + r2) with time step T and mesh size h. An efficient iterative algorithm for solving these nonlinear schemes is proposed and discussed in detail.展开更多
In this paper the authors discuss a numerical simulation problem of three-dimensional compressible contamination treatment from nuclear waste. The mathematical model, a nonlinear convection-diffusion system of four PD...In this paper the authors discuss a numerical simulation problem of three-dimensional compressible contamination treatment from nuclear waste. The mathematical model, a nonlinear convection-diffusion system of four PDEs, determines four major physical unknowns: the pressure, the concentrations of brine and radionuclide, and the temperature. The pressure is solved by a conservative mixed finite volume element method, and the computational accuracy is improved for Darcy velocity. Other unknowns are computed by a composite scheme of upwind approximation and mixed finite volume element. Numerical dispersion and nonphysical oscillation are eliminated, and the convection-dominated diffusion problems are solved well with high order computational accuracy. The mixed finite volume element is conservative locally, and get the objective functions and their adjoint vector functions simultaneously. The conservation nature is an important character in numerical simulation of underground fluid. Fractional step difference is introduced to solve the concentrations of radionuclide factors, and the computational work is shortened significantly by decomposing a three-dimensional problem into three successive one-dimensional problems. By the theory and technique of a priori estimates of differential equations, we derive an optimal order estimates in L2norm. Finally, numerical examples show the effectiveness and practicability for some actual problems.展开更多
A physically based numerical approach is presented for modeling multiphase flow and transport processes in fractured rock.In particular,a general framework model is discussed for dealing with fracture-matrix interacti...A physically based numerical approach is presented for modeling multiphase flow and transport processes in fractured rock.In particular,a general framework model is discussed for dealing with fracture-matrix interactions,which is applicable to both continuum and discrete fracture conceptualization.The numerical modeling approach is based on a general multiple-continuum concept,suitable for modeling any types of fractured reservoirs,including double-,triple-,and other multiplecontinuum conceptual models.In addition,a new,physically correct numerical scheme is discussed to calculate multiphase flow between fractures and the matrix,using continuity of capillary pressure at the fracture-matrix interface.The proposed general modeling methodology is verified in special cases using analytical solutions and laboratory experimental data,and demonstrated for its application in modeling flow through fractured vuggy reservoirs.展开更多
In this paper the authors discuss the numerical simulation problem of threedimensional compressible contamination treatment from nuclear waste.The mathematical model is defined by an initial-boundary nonlinear convect...In this paper the authors discuss the numerical simulation problem of threedimensional compressible contamination treatment from nuclear waste.The mathematical model is defined by an initial-boundary nonlinear convection-diffusion system of four partial differential equations:a parabolic equation for the pressure,two convection-diffusion equations for the concentrations of brine and radionuclide and a heat conduction equation for the temperature.The pressure appears within the concentration equations and heat conduction equation,and the Darcy velocity controls the concentrations and the temperature.The pressure is solved by the conservative mixed volume element method,and the order of the accuracy is improved by the Darcy velocity.The concentration of brine and temperature are computed by the upwind mixed volume element method on a changing mesh,where the diffusion is discretized by a mixed volume element and the convection is treated by an upwind scheme.The composite method can solve the convection-dominated diffusion problems well because it eliminates numerical dispersion and nonphysical oscillation and has high order computational accuracy.The mixed volume element has the local conservation of mass and energy,and it can obtain the brine and temperature and their adjoint vector functions simultaneously.The conservation nature plays an important role in numerical simulation of underground fluid.The concentrations of radionuclide factors are solved by the method of upwind fractional step difference and the computational work is decreased by decomposing a three-dimensional problem into three successive one-dimensional problems and using the method of speedup.By the theory and technique of a priori estimates of differential equations,we derive an optimal order result in L^(2) norm.Numerical examples are given to show the effectiveness and practicability and the composite method is testified as a powerful tool to solve the well-known actual problem.展开更多
This article describes a number of velocity-based moving mesh numerical methods for multidimensional nonlinear time-dependent partial differential equations(PDEs).It consists of a short historical review followed by a...This article describes a number of velocity-based moving mesh numerical methods for multidimensional nonlinear time-dependent partial differential equations(PDEs).It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation.Finite element algorithms are derived for both mass-conserving and non mass-conserving problems,and results shown for a number of multidimensional nonlinear test problems,including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem.Further applications and extensions are referenced.展开更多
文摘By using the revision of the momentum for a particle moving with high velocity and by investigating the famous Bucherer's experiment of an electron deflecting with high velocity in the electromagnetic fields in 1908, the paper determines that mass of the electron with high velocity is still to observe the law of conservation of mass.
文摘The purpose is to reestablish the coupled conservation laws, the local conservation equations and the jump conditions of mass and inertia for polar continuum theories. In this connection the new material derivatives of the deformation gradient, the line element, the surface element and the volume element were derived and the generalized Reynolds transport theorem was presented. Combining these conservation laws of mass and inertia with the balance laws of momentum, angular momentum and energy derived in our previous papers of this series, a rather complete system of coupled basic laws and principles for polar continuum theories is constituted on the whole. From this system the coupled nonlocal balance equations of mass, inertia, momentum, angular momentum and energy may be obtained by the usual localization.
基金supported the Natural Science Foundation of Shandong Province(ZR2016AM08)Natural Science Foundation of Hunan Province(2018JJ2028)National Natural Science Foundation of China(11871312).
文摘The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.
基金Supported by the National Natural Science Foundation of China(11101124 and 11271231)Natural Science Foundation of Shandong Province(ZR2016AM08)National Tackling Key Problems Program(2011ZX05052,2011ZX05011-004)
文摘The physical model is described by a seepage coupled system for simulating numerically three-dimensional chemical oil recovery, whose mathematical description includes three equations to interpret main concepts. The pressure equation is a nonlinear parabolic equation, the concentration is defined by a convection-diffusion equation and the saturations of different components are stated by nonlinear convection-diffusion equations. The transport pressure appears in the concentration equation and saturation equations in the form of Darcy velocity, and controls their processes. The flow equation is solved by the conservative mixed volume element and the accuracy is improved one order for approximating Darcy velocity. The method of characteristic mixed volume element is applied to solve the concentration, where the diffusion is discretized by a mixed volume element method and the convection is treated by the method of characteristics. The characteristics can confirm strong computational stability at sharp fronts and it can avoid numerical dispersion and nonphysical oscillation. The scheme can adopt a large step while its numerical results have small time-truncation error and high order of accuracy. The mixed volume element method has the law of conservation on every element for the diffusion and it can obtain numerical solutions of the concentration and adjoint vectors. It is most important in numerical simulation to ensure the physical conservative nature. The saturation different components are obtained by the method of characteristic fractional step difference. The computational work is shortened greatly by decomposing a three-dimensional problem into three successive one-dimensional problems and it is completed easily by using the algorithm of speedup. Using the theory and technique of a priori estimates of differential equations, we derive an optimal second order estimates in 12 norm. Numerical examples are given to show the effectiveness and practicability and the method is testified as a powerful tool to solve the important problems.
基金supported by the Natural ScienceFoundation of Shandong Province(ZR2021MA019)。
文摘In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order.The concentration equation is solved by the method of mixed finite element,where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method.The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation.In actual computations the characteristics adopts a large time step without any loss of accuracy.The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux.In order to derive the optimal 3/2-order error estimate in L^(2) norm,a post-processing technique is included in the approximation to the scalar unknowns.Numerical experiments are illustrated finally to validate theoretical analysis and efficiency.This method can be used to solve such an important problem.
基金supported by the National Natural Science Foundation of China(Grant No.11871312)Natural Science Foundation of Shandong Province(Grant No.ZR2021MA019).
文摘In this paper,the authors discuss a three-dimensional problem of the semiconductor device type involved its mathematical description,numerical simulation and theoretical analysis.Two important factors,heat and magnetic influences are involved.The mathematical model is formulated by four nonlinear partial differential equations(PDEs),determining four major physical variables.The influences of magnetic fields are supposed to be weak,and the strength is parallel to the z-axis.The elliptic equation is treated by a block-centered method,and the law of conservation is preserved.The computational accuracy is improved one order.Other equations are convection-dominated,thus are approximated by upwind block-centered differences.Upwind difference can eliminate numerical dispersion and nonphysical oscillation.The diffusion is approximated by the block-centered difference,while the convection term is treated by upwind approximation.Furthermore,the unknowns and adjoint functions are computed at the same time.These characters play important roles in numerical computations of conductor device problems.Using the theories of priori analysis such as energy estimates,the principle of duality and mathematical inductions,an optimal estimates result is obtained.Then a composite numerical method is shown for solving this problem.
基金Supported by National Natural Science Foundation of China(Nos.11201239,11571181)
文摘In this paper, we propose a class of stable finite difference schemes for the initial-boundary value problem of the Cahn-Hilliard equation. These schemes are proved to inherit the total mass conservation and energy dissipation in the discrete level. The dissipation of the total energy implies boundness of the numerical solutions in the discrete H1 norm. This in turn implies boundedness of the numerical solutions in the maximum norm and hence the stability of the difference schemes. Unique existence of the numerical solutions is proved by the fixed-point theorem. Convergence rate of the class of finite difference schemes is proved to be O(h2 + r2) with time step T and mesh size h. An efficient iterative algorithm for solving these nonlinear schemes is proposed and discussed in detail.
基金supported by the Natural Science Foundation of Shangdong Province (Grant No.ZR2021MA019)Natural Science Foundation of Hunan Province (Grant No.2018JJ2028)。
文摘In this paper the authors discuss a numerical simulation problem of three-dimensional compressible contamination treatment from nuclear waste. The mathematical model, a nonlinear convection-diffusion system of four PDEs, determines four major physical unknowns: the pressure, the concentrations of brine and radionuclide, and the temperature. The pressure is solved by a conservative mixed finite volume element method, and the computational accuracy is improved for Darcy velocity. Other unknowns are computed by a composite scheme of upwind approximation and mixed finite volume element. Numerical dispersion and nonphysical oscillation are eliminated, and the convection-dominated diffusion problems are solved well with high order computational accuracy. The mixed finite volume element is conservative locally, and get the objective functions and their adjoint vector functions simultaneously. The conservation nature is an important character in numerical simulation of underground fluid. Fractional step difference is introduced to solve the concentrations of radionuclide factors, and the computational work is shortened significantly by decomposing a three-dimensional problem into three successive one-dimensional problems. By the theory and technique of a priori estimates of differential equations, we derive an optimal order estimates in L2norm. Finally, numerical examples show the effectiveness and practicability for some actual problems.
文摘A physically based numerical approach is presented for modeling multiphase flow and transport processes in fractured rock.In particular,a general framework model is discussed for dealing with fracture-matrix interactions,which is applicable to both continuum and discrete fracture conceptualization.The numerical modeling approach is based on a general multiple-continuum concept,suitable for modeling any types of fractured reservoirs,including double-,triple-,and other multiplecontinuum conceptual models.In addition,a new,physically correct numerical scheme is discussed to calculate multiphase flow between fractures and the matrix,using continuity of capillary pressure at the fracture-matrix interface.The proposed general modeling methodology is verified in special cases using analytical solutions and laboratory experimental data,and demonstrated for its application in modeling flow through fractured vuggy reservoirs.
基金The authors express their deep appreciation to Prof.J.Douglas Jr,Prof.R.E.Ewing,and Prof.L.S.Jiang for their many helpful suggestions in the series research on numerical simulation of energy sciences.Also,the project is supported by NSAF(Grant No.U1430101)Natural Science Foundation of Shandong Province(Grant No.ZR2016AM08)National Tackling Key Problems Program(Grant Nos.2011ZX05052,2011ZX05011-004,20050200069).
文摘In this paper the authors discuss the numerical simulation problem of threedimensional compressible contamination treatment from nuclear waste.The mathematical model is defined by an initial-boundary nonlinear convection-diffusion system of four partial differential equations:a parabolic equation for the pressure,two convection-diffusion equations for the concentrations of brine and radionuclide and a heat conduction equation for the temperature.The pressure appears within the concentration equations and heat conduction equation,and the Darcy velocity controls the concentrations and the temperature.The pressure is solved by the conservative mixed volume element method,and the order of the accuracy is improved by the Darcy velocity.The concentration of brine and temperature are computed by the upwind mixed volume element method on a changing mesh,where the diffusion is discretized by a mixed volume element and the convection is treated by an upwind scheme.The composite method can solve the convection-dominated diffusion problems well because it eliminates numerical dispersion and nonphysical oscillation and has high order computational accuracy.The mixed volume element has the local conservation of mass and energy,and it can obtain the brine and temperature and their adjoint vector functions simultaneously.The conservation nature plays an important role in numerical simulation of underground fluid.The concentrations of radionuclide factors are solved by the method of upwind fractional step difference and the computational work is decreased by decomposing a three-dimensional problem into three successive one-dimensional problems and using the method of speedup.By the theory and technique of a priori estimates of differential equations,we derive an optimal order result in L^(2) norm.Numerical examples are given to show the effectiveness and practicability and the composite method is testified as a powerful tool to solve the well-known actual problem.
文摘This article describes a number of velocity-based moving mesh numerical methods for multidimensional nonlinear time-dependent partial differential equations(PDEs).It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation.Finite element algorithms are derived for both mass-conserving and non mass-conserving problems,and results shown for a number of multidimensional nonlinear test problems,including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem.Further applications and extensions are referenced.