A POSTERIORI ERROR ESTIMATION OF THE NEW MIXED ELEMENT SCHEMES FOR SECOND ORDER ELLIPTIC PROBLEM ON ANISOTROPIC MESHES

This paper presents a posteriori residual error estimator for the new mixed el-ement scheme for second order elliptic problem on anisotropic meshes. The reliability and efficiency of our estimator are established without any regularity assumption on the mesh.

A MIXED FINITE ELEMENT AND CHARACTERISTIC MIXED FINITE ELEMENT FOR INCOMPRESSIBLE MISCIBLE DARCY-FORCHHEIMER DISPLACEMENT AND NUMERICAL ANALYSIS

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.

UNSTEADY/STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.

定常Navier-Stokes方程的一个二阶非协调三角型混合元格式（英文）

In this paper, a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure. The convergence analysis is presented and optimal error estimates of both broken H^1-norm and L^2-norm for velocity as well as the L^2-norm for the pressure are derived.

A NEW CHARACTERISTIC EXPANDED MIXED METHOD FOR SOBOLEV EQUATION WITH CONVECTION TERM

In this paper,a new numerical method based on a new expanded mixed scheme and the characteristic method is developed and discussed for Sobolev equation with convection term.The hyperbolic part d(x)∂u/∂t+c(x,t)·∇u is handled by the characteristic method and the diffusion term∇·(a(x,t)∇u+b(x,t)∇ut)is approximated by the new expanded mixed method,whose gradient belongs to the simple square integrable(L^(2)(Ω))^(2)space instead of the classical H(div;Ω)space.For a priori error estimates,some important lemmas based on the new expanded mixed projection are introduced.An optimal priori error estimates in L^(2)-norm for the scalar unknown u and a priori error estimates in(L^(2))^(2)-norm for its gradientλ,and its fluxσ(the coefficients times the negative gradient)are derived.In particular,an optimal priori error estimate in H1-norm for the scalar unknown u is obtained.

ON RESIDUAL-BASED A POSTERIORI ERROR ESTIMATORS FOR LOWEST-ORDER RAVIART-THOMAS ELEMENT APPROXIMATION TO CONVECTION-DIFFUSION-REACTION EQUATIONS

A new technique of residual-type a posteriori error analysis is developed for the lowest- order Raviart-Thomas mixed finite element discretizations of convection-diffusion-reaction equations in two- or three-dimension. Both centered mixed scheme and upwind-weighted mixed scheme are considered. The a posteriori error estimators, derived for the stress variable error plus scalar displacement error in L_2-norm, can be directly computed with the solutions of the mixed schemes without any additional cost, and are proven to be reliable. Local efficiency dependent on local variations in coefficients is obtained without any saturation assumption, and holds from the cases where convection or reaction is not present to convection- or reaction-dominated problems. The main tools of the analysis are the postprocessed approximation of scalar displacement, abstract error estimates, and the property of modified Oswald interpolation. Numerical experiments are carried out to support our theoretical results and to show the competitive behavior of the proposed posteriori error estimates.

Radiative decays J/ψ → η (') γ in perturbative QCD

With the recent investigations of the g＊g＊-η（＇） transition form factor and mixing scheme, we present an updated study of the radiative decays J/ψ→η（＇）γ in perturbative QCD. The decays are taken as a test ground for the g＊g＊-η（＇） transition form factors and the η-η＇ mixing scheme. The form factors are found to be working for glunic η＇ production and the mixing angle is constrained to be φ = 35.1°±1 0.8°.

IMPACT OF DIABATIC PROCESSES IN AGCM ON 4-DIMENSIONAL VARIATIONAL DATA ASSIMILATION

The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to 'observation' (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by 'on-off' switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.