A simple characteristic equation solution strategy for deriving the fun- damental analytical solutions of 3D isotropic elasticity is proposed. By calculating the determinant of the differential operator matrix obtaine...A simple characteristic equation solution strategy for deriving the fun- damental analytical solutions of 3D isotropic elasticity is proposed. By calculating the determinant of the differential operator matrix obtained from the governing equations of 3D elasticity, the characteristic equation which the characteristic general solution vectors must satisfy is established. Then, by substitution of the characteristic general solution vectors, which satisfy various reduced characteristic equations, into various reduced ad- joint matrices of the differential operator matrix, the corresponding fundamental analyt- ical solutions for isotropic 3D elasticity, including Boussinesq-Galerkin (B-G) solutions, modified Papkovich-Neuber solutions proposed by Min-zhong WANG (P-N-W), and quasi HU Hai-chang solutions, can be obtained. Furthermore, the independence characters of various fundamental solutions in polynomial form are also discussed in detail. These works provide a basis for constructing complete and independent analytical trial func- tions used in numerical methods.展开更多
The setting work of backup protection using steady-state current is tedious,and mismatches occasionally occur due to the increased proportion of distributed generations(DGs)connected to the power grid.Thus,there is a ...The setting work of backup protection using steady-state current is tedious,and mismatches occasionally occur due to the increased proportion of distributed generations(DGs)connected to the power grid.Thus,there is a practical need to study a backup protection technology that does not require step-by-step setting and can be adaptively coordinated.This paper proposes an action sequence adaptive to fault positions that uses only positive sequence fault component(PSFC)voltage.Considering the influence of DGs,the unified time dial setting can be obtained by selecting specific points.The protection performance is improved by using the adjacent upstream and downstream protections to meet the coordination time interval in the case of metallic faults at the near-and far-ends of the line.Finally,the expression and implementation scheme for inverse-time backup protection(ITBP)based on the unified characteristic equation is given.Simulation results show that this scheme can adapt to DG penetration scenarios and can realize the adaptive coordination of multi-level relays.展开更多
Delay differential systems are widely used in many different fields, It is important to determine the local stability of their equilibria, For systems with' delay dependent parameters, the stability analysis of equil...Delay differential systems are widely used in many different fields, It is important to determine the local stability of their equilibria, For systems with' delay dependent parameters, the stability analysis of equilibria is complicated and difficult In this paper, we shall investigate the ultimate stability of a type of characteristic equation with delay dependent parameters. Our results show that the characteristic equation with delay dependent parameters may be one of ultimately stable, ultimately unstable, and alternate between stable and unstable. Applying our results, the ultimate stability can be often decided directly and need not appeal to mathematic software. Two examples are given in this paper,展开更多
Root of characteristic equation for cylindrical Bessel equation eigenvalue problems on general interval is of great real physical importance at engineering and physical. First, the characteristic equation of cylindric...Root of characteristic equation for cylindrical Bessel equation eigenvalue problems on general interval is of great real physical importance at engineering and physical. First, the characteristic equation of cylindrical Bessel equation eigenvalue problem on general interval is given, second, by mean of compared method, we obtaining roots of characteristic equation with Matlab program is discussed.展开更多
This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using prec...This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.展开更多
Dynamic characteristics of the resonant gyroscope are studied based on the Mathieu equation approximate solution in this paper.The Mathieu equation is used to analyze the parametric resonant characteristics and the ap...Dynamic characteristics of the resonant gyroscope are studied based on the Mathieu equation approximate solution in this paper.The Mathieu equation is used to analyze the parametric resonant characteristics and the approximate output of the resonant gyroscope.The method of small parameter perturbation is used to analyze the approximate solution of the Mathieu equation.The theoretical analysis and the numerical simulations show that the approximate solution of the Mathieu equation is close to the dynamic output characteristics of the resonant gyroscope.The experimental analysis shows that the theoretical curve and the experimental data processing results coincide perfectly,which means that the approximate solution of the Mathieu equation can present the dynamic output characteristic of the resonant gyroscope.The theoretical approach and the experimental results of the Mathieu equation approximate solution are obtained,which provides a reference for the robust design of the resonant gyroscope.展开更多
The approximate equations for self- and mutual-inductance can be further investigated. In addition, the general characteristic equations for the four fundamental elements in three situations are basically completed. T...The approximate equations for self- and mutual-inductance can be further investigated. In addition, the general characteristic equations for the four fundamental elements in three situations are basically completed. This new development not only reveals theoretically the characteristics of the four fundamental electronic elements, but also practically deals with all the media from linearity to nonlinearity. This will meet all the purposes of application. Hence, it is very important and significant to the development of the science of electronic elements.展开更多
A lattice Boltzmann(LB) theory, the analytical characteristic integral(ACI) LB theory, is proposed in this paper.ACI LB theory takes the Bhatnagar–Gross–Krook(BGK)-Boltzmann equation as the exact kinetic equation be...A lattice Boltzmann(LB) theory, the analytical characteristic integral(ACI) LB theory, is proposed in this paper.ACI LB theory takes the Bhatnagar–Gross–Krook(BGK)-Boltzmann equation as the exact kinetic equation behind Navier–Stokes continuum and momentum equations and constructs an LB equation by rigorously integrating the BGK-Boltzmann equation along characteristics. It is a general theory, supporting most existing LB equations including the standard lattice BGK(LBGK) equation inherited from lattice-gas automata, whose theoretical foundation had been questioned. ACI LB theory also indicates that the characteristic parameter of an LB equation is collision number, depicting the particle-interaction intensity in the time span of the LB equation, instead of the traditionally assumed relaxation time, and the over-relaxation time problem is merely a manifestation of the temporal evolution of equilibrium distribution along characteristics under high collision number, irrelevant to particle kinetics. In ACI LB theory, the temporal evolution of equilibrium distribution along characteristics is the determinant of LB method accuracy and we numerically prove this.展开更多
In this paper, the characteristic function method is applied to seek traveling wave solutions of nonlinear partial differential equations in a unified way. We consider the Wu-Zhang equation (which describes (1 + 1)-di...In this paper, the characteristic function method is applied to seek traveling wave solutions of nonlinear partial differential equations in a unified way. We consider the Wu-Zhang equation (which describes (1 + 1)-dimensional disper-sive long wave). The equations governing the wave propagation consist of a pair of non linear partial differential equations. The characteristic function method reduces the system of nonlinear partial differential equations to a system of nonlinear ordinary differential equations which is solved via the shooting method, coupled with Rungekutta scheme. The results include kink-profile solitary wave solutions, periodic wave solutions and rational solutions. As an illustrative example, the properties of some soliton solutions for Wu-Zhang equation are shown by some figures.展开更多
Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow pheno...Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.展开更多
We propose a new way of rewriting the two dimensional Euler equations and derive an original canonical characteristic relation based on the characteristic theory of hyperbolic systems. This relation contains the deriv...We propose a new way of rewriting the two dimensional Euler equations and derive an original canonical characteristic relation based on the characteristic theory of hyperbolic systems. This relation contains the derivatives strictly along the bicharacteristic directions, and can be viewed as the 2D extension of the characteristic relation in 1D case.展开更多
In this paper, we present a differential polynomial characteristic set algorithm for the complete symmetry classification of partial differential equations (PDEs) with some parameters. It can make the solution to th...In this paper, we present a differential polynomial characteristic set algorithm for the complete symmetry classification of partial differential equations (PDEs) with some parameters. It can make the solution to the complete symmetry classification problem for PDEs become direct and systematic. As an illustrative example, the complete potential symmetry classifications of nonlinear and linear wave equations with an arbitrary function parameter are presented. This is a new application of the differential form characteristic set algorithm, i.e., Wu's method, in differential equations.展开更多
In this paper, a high-order finite-volume scheme is presented for the one- dimensional scalar and inviscid Euler conservation laws. The Simpson's quadrature rule is used to achieve high-order accuracy in time. To get...In this paper, a high-order finite-volume scheme is presented for the one- dimensional scalar and inviscid Euler conservation laws. The Simpson's quadrature rule is used to achieve high-order accuracy in time. To get the point value of the Simpson's quadrature, the characteristic theory is used to obtain the positions of the grid points at each sub-time stage along the characteristic curves, and the third-order and fifth-order central weighted essentially non-oscillatory (CWENO) reconstruction is adopted to estimate the cell point values. Several standard one-dimensional examples are used to verify the high-order accuracy, convergence and capability of capturing shock.展开更多
The 1/3 sub-harmonic solution for the Duffing's with damping equation was investigated by using the methods of harmonic balance and numerical integration. The assumed solution is introduced, and the domain of sub-har...The 1/3 sub-harmonic solution for the Duffing's with damping equation was investigated by using the methods of harmonic balance and numerical integration. The assumed solution is introduced, and the domain of sub-harmonic frequencies was found. The asymptotical stability of the subharmonic resonances and the sensitivity of the amplitude responses to the variation of damping coefficient were examined. Then, the subharmonic resonances were analyzed by using the techniques from the general fractal theory. The analysis indicates that the sensitive dimensions of the system time-field responses show sensitivity to the conditions of changed initial perturbation, changed damping coefficient or the amplitude of excitation, thus the sensitive dimension can clearly describe the characteristic of the transient process of the subharmonic resonances.展开更多
The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE s...The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE scheme for the discrete_time along characteristics is presented and error estimates are established.The existence and convergence of MFE solution of the discrete current velocity,elevation of the bottom topography,thickness of fluid column,and mass rate of sediment is demonstrated.展开更多
This paper describes a characteristics-mix finite element method for the computation of incompressible Navi-er-Stokes equations with variable density. We have introduced a mixed scheme which combines a characteristics...This paper describes a characteristics-mix finite element method for the computation of incompressible Navi-er-Stokes equations with variable density. We have introduced a mixed scheme which combines a characteristics finite element scheme for treating the mass conservation equation and a finite element method to deal with the momentum equation and the divergence free constraint. The proposed method has a lot of attractive computational properties: parameter-free, very flexible, and averting the difficulties caused by the original equations. The stability of the method is proved. Finally, several numerical experiments are given to show that this method is efficient for variable density incompressible flows problem.展开更多
A combined approximate scheme is defined for convection-diffusion-reaction equations. This scheme is constructed by two methods. Standard mixed finite element method is used for diffusion term. A second order characte...A combined approximate scheme is defined for convection-diffusion-reaction equations. This scheme is constructed by two methods. Standard mixed finite element method is used for diffusion term. A second order characteristic finite element method is presented to handle the material derivative term, that is, the time derivative term plus the convection term. The stability is proved and the L2-norm error estimates are derived for both the scalar unknown variable and its flux. The scheme is of second order accuracy in time increment, symmetric, and unconditionally stable.展开更多
A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlin...A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlinear systems by two steps. The error analysis shows that the two-grid scheme combined with the characteristic mixed finite-element method can decrease numerical oscillation caused by dominated convections and solve nonlinear advection-dominated diffusion problems efficiently.展开更多
This paper discusses the problem of the extraction of characteristic roots on Bessel-Neumann's mixed equations. It gives the expressions and the evaluation of the minimum root. The advantage of the method has no u...This paper discusses the problem of the extraction of characteristic roots on Bessel-Neumann's mixed equations. It gives the expressions and the evaluation of the minimum root. The advantage of the method has no use for the table of the multi-figure number Bessel function and it does not need computer but can calculate all the characteristic roots The precision of these roots is still high.展开更多
A simple model of chromatographic mechanical mechanism is present, and then a scrics of theoretical chromatographic equations and fundamental Formulae are derived. These theoretical equations and formulae not only res...A simple model of chromatographic mechanical mechanism is present, and then a scrics of theoretical chromatographic equations and fundamental Formulae are derived. These theoretical equations and formulae not only reserve thermodynamic characteristics in the current fundamental chromatographic formulae, but also introduce one or more kinetic parameter, so it is possible to make the macroscopic-control on the effect of kinetic characteristics on chromatographic system.展开更多
基金supported by the National Natural Science Foundation of China (Nos. 10872108 and10876100)the Program for New Century Excellent Talents in University (No. NCET-07-0477)the National Basic Research Programs of China (Nos. 2010CB731503 and 2010CB832701)
文摘A simple characteristic equation solution strategy for deriving the fun- damental analytical solutions of 3D isotropic elasticity is proposed. By calculating the determinant of the differential operator matrix obtained from the governing equations of 3D elasticity, the characteristic equation which the characteristic general solution vectors must satisfy is established. Then, by substitution of the characteristic general solution vectors, which satisfy various reduced characteristic equations, into various reduced ad- joint matrices of the differential operator matrix, the corresponding fundamental analyt- ical solutions for isotropic 3D elasticity, including Boussinesq-Galerkin (B-G) solutions, modified Papkovich-Neuber solutions proposed by Min-zhong WANG (P-N-W), and quasi HU Hai-chang solutions, can be obtained. Furthermore, the independence characters of various fundamental solutions in polynomial form are also discussed in detail. These works provide a basis for constructing complete and independent analytical trial func- tions used in numerical methods.
基金supported in part by the National Natural Science Foundation of China(NSFC-UKRI_EPSRC)(No.52061635105)in part by the Science and Technology Program of SGCC(No.5100-202040327A-0-0-00).
文摘The setting work of backup protection using steady-state current is tedious,and mismatches occasionally occur due to the increased proportion of distributed generations(DGs)connected to the power grid.Thus,there is a practical need to study a backup protection technology that does not require step-by-step setting and can be adaptively coordinated.This paper proposes an action sequence adaptive to fault positions that uses only positive sequence fault component(PSFC)voltage.Considering the influence of DGs,the unified time dial setting can be obtained by selecting specific points.The protection performance is improved by using the adjacent upstream and downstream protections to meet the coordination time interval in the case of metallic faults at the near-and far-ends of the line.Finally,the expression and implementation scheme for inverse-time backup protection(ITBP)based on the unified characteristic equation is given.Simulation results show that this scheme can adapt to DG penetration scenarios and can realize the adaptive coordination of multi-level relays.
基金This work was supported by China Ministry of Science and Technology(2004BA719A01)Postdoctoral Function of China(2005037785).
文摘Delay differential systems are widely used in many different fields, It is important to determine the local stability of their equilibria, For systems with' delay dependent parameters, the stability analysis of equilibria is complicated and difficult In this paper, we shall investigate the ultimate stability of a type of characteristic equation with delay dependent parameters. Our results show that the characteristic equation with delay dependent parameters may be one of ultimately stable, ultimately unstable, and alternate between stable and unstable. Applying our results, the ultimate stability can be often decided directly and need not appeal to mathematic software. Two examples are given in this paper,
文摘Root of characteristic equation for cylindrical Bessel equation eigenvalue problems on general interval is of great real physical importance at engineering and physical. First, the characteristic equation of cylindrical Bessel equation eigenvalue problem on general interval is given, second, by mean of compared method, we obtaining roots of characteristic equation with Matlab program is discussed.
文摘This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.
基金Project supported by the National Natural Science Foundation of China (Grant No. 60927005)the Innovation Foundation of BUAA for Ph. D. Graduates,Chinathe Fundamental Research Funds for the Central Universities,China (Grant No. YWF-10-01-A17)
文摘Dynamic characteristics of the resonant gyroscope are studied based on the Mathieu equation approximate solution in this paper.The Mathieu equation is used to analyze the parametric resonant characteristics and the approximate output of the resonant gyroscope.The method of small parameter perturbation is used to analyze the approximate solution of the Mathieu equation.The theoretical analysis and the numerical simulations show that the approximate solution of the Mathieu equation is close to the dynamic output characteristics of the resonant gyroscope.The experimental analysis shows that the theoretical curve and the experimental data processing results coincide perfectly,which means that the approximate solution of the Mathieu equation can present the dynamic output characteristic of the resonant gyroscope.The theoretical approach and the experimental results of the Mathieu equation approximate solution are obtained,which provides a reference for the robust design of the resonant gyroscope.
文摘The approximate equations for self- and mutual-inductance can be further investigated. In addition, the general characteristic equations for the four fundamental elements in three situations are basically completed. This new development not only reveals theoretically the characteristics of the four fundamental electronic elements, but also practically deals with all the media from linearity to nonlinearity. This will meet all the purposes of application. Hence, it is very important and significant to the development of the science of electronic elements.
基金Project supported by the National Science and Technology Major Project,China(Grant No.2017ZX06002002)
文摘A lattice Boltzmann(LB) theory, the analytical characteristic integral(ACI) LB theory, is proposed in this paper.ACI LB theory takes the Bhatnagar–Gross–Krook(BGK)-Boltzmann equation as the exact kinetic equation behind Navier–Stokes continuum and momentum equations and constructs an LB equation by rigorously integrating the BGK-Boltzmann equation along characteristics. It is a general theory, supporting most existing LB equations including the standard lattice BGK(LBGK) equation inherited from lattice-gas automata, whose theoretical foundation had been questioned. ACI LB theory also indicates that the characteristic parameter of an LB equation is collision number, depicting the particle-interaction intensity in the time span of the LB equation, instead of the traditionally assumed relaxation time, and the over-relaxation time problem is merely a manifestation of the temporal evolution of equilibrium distribution along characteristics under high collision number, irrelevant to particle kinetics. In ACI LB theory, the temporal evolution of equilibrium distribution along characteristics is the determinant of LB method accuracy and we numerically prove this.
文摘In this paper, the characteristic function method is applied to seek traveling wave solutions of nonlinear partial differential equations in a unified way. We consider the Wu-Zhang equation (which describes (1 + 1)-dimensional disper-sive long wave). The equations governing the wave propagation consist of a pair of non linear partial differential equations. The characteristic function method reduces the system of nonlinear partial differential equations to a system of nonlinear ordinary differential equations which is solved via the shooting method, coupled with Rungekutta scheme. The results include kink-profile solitary wave solutions, periodic wave solutions and rational solutions. As an illustrative example, the properties of some soliton solutions for Wu-Zhang equation are shown by some figures.
基金King Mongkut’s University of Technology North Bangkok (KMUTNB)the Office of the Higher Education Commission (OHEC)the National Metal and Materials Technology Center (MTEC) for supporting this research work
文摘Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.
基金Supported by the NNSF of China(10871029)the foundation of LCP(9140C6902020904)
文摘We propose a new way of rewriting the two dimensional Euler equations and derive an original canonical characteristic relation based on the characteristic theory of hyperbolic systems. This relation contains the derivatives strictly along the bicharacteristic directions, and can be viewed as the 2D extension of the characteristic relation in 1D case.
基金supported by the Ph. D. Programs Foundation of Ministry of Education of China(No.20070128001)the Expenditure Budget program of Shanghai Municipal Education Commission (No.2008069)+1 种基金the Innovation Program of Shanghai Municipal Education Commission(No.09YZ239)the Natural Science Foundation of Inner Mongolia (No.200607010103)
文摘In this paper, we present a differential polynomial characteristic set algorithm for the complete symmetry classification of partial differential equations (PDEs) with some parameters. It can make the solution to the complete symmetry classification problem for PDEs become direct and systematic. As an illustrative example, the complete potential symmetry classifications of nonlinear and linear wave equations with an arbitrary function parameter are presented. This is a new application of the differential form characteristic set algorithm, i.e., Wu's method, in differential equations.
基金supported by the National Natural Science Foundation of China (No. 10771134)
文摘In this paper, a high-order finite-volume scheme is presented for the one- dimensional scalar and inviscid Euler conservation laws. The Simpson's quadrature rule is used to achieve high-order accuracy in time. To get the point value of the Simpson's quadrature, the characteristic theory is used to obtain the positions of the grid points at each sub-time stage along the characteristic curves, and the third-order and fifth-order central weighted essentially non-oscillatory (CWENO) reconstruction is adopted to estimate the cell point values. Several standard one-dimensional examples are used to verify the high-order accuracy, convergence and capability of capturing shock.
基金Project supported by the National Natural Science Foundation of China (No.50275024)
文摘The 1/3 sub-harmonic solution for the Duffing's with damping equation was investigated by using the methods of harmonic balance and numerical integration. The assumed solution is introduced, and the domain of sub-harmonic frequencies was found. The asymptotical stability of the subharmonic resonances and the sensitivity of the amplitude responses to the variation of damping coefficient were examined. Then, the subharmonic resonances were analyzed by using the techniques from the general fractal theory. The analysis indicates that the sensitive dimensions of the system time-field responses show sensitivity to the conditions of changed initial perturbation, changed damping coefficient or the amplitude of excitation, thus the sensitive dimension can clearly describe the characteristic of the transient process of the subharmonic resonances.
文摘The mixed finite element(MFE) methods for a shallow water equation system consisting of water dynamics equations,silt transport equation,and the equation of bottom topography change were derived.A fully discrete MFE scheme for the discrete_time along characteristics is presented and error estimates are established.The existence and convergence of MFE solution of the discrete current velocity,elevation of the bottom topography,thickness of fluid column,and mass rate of sediment is demonstrated.
文摘This paper describes a characteristics-mix finite element method for the computation of incompressible Navi-er-Stokes equations with variable density. We have introduced a mixed scheme which combines a characteristics finite element scheme for treating the mass conservation equation and a finite element method to deal with the momentum equation and the divergence free constraint. The proposed method has a lot of attractive computational properties: parameter-free, very flexible, and averting the difficulties caused by the original equations. The stability of the method is proved. Finally, several numerical experiments are given to show that this method is efficient for variable density incompressible flows problem.
文摘A combined approximate scheme is defined for convection-diffusion-reaction equations. This scheme is constructed by two methods. Standard mixed finite element method is used for diffusion term. A second order characteristic finite element method is presented to handle the material derivative term, that is, the time derivative term plus the convection term. The stability is proved and the L2-norm error estimates are derived for both the scalar unknown variable and its flux. The scheme is of second order accuracy in time increment, symmetric, and unconditionally stable.
文摘A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlinear systems by two steps. The error analysis shows that the two-grid scheme combined with the characteristic mixed finite-element method can decrease numerical oscillation caused by dominated convections and solve nonlinear advection-dominated diffusion problems efficiently.
文摘This paper discusses the problem of the extraction of characteristic roots on Bessel-Neumann's mixed equations. It gives the expressions and the evaluation of the minimum root. The advantage of the method has no use for the table of the multi-figure number Bessel function and it does not need computer but can calculate all the characteristic roots The precision of these roots is still high.
文摘A simple model of chromatographic mechanical mechanism is present, and then a scrics of theoretical chromatographic equations and fundamental Formulae are derived. These theoretical equations and formulae not only reserve thermodynamic characteristics in the current fundamental chromatographic formulae, but also introduce one or more kinetic parameter, so it is possible to make the macroscopic-control on the effect of kinetic characteristics on chromatographic system.