For a generalized quasi-Newtonian flow, a new stabilized method focused on the low-order velocity-pressure pairs, (bi)linear/(bi)linear and (bi)linear/constant element, is presented. The pressure projection stab...For a generalized quasi-Newtonian flow, a new stabilized method focused on the low-order velocity-pressure pairs, (bi)linear/(bi)linear and (bi)linear/constant element, is presented. The pressure projection stabilized method is extended from Stokes problems to quasi-Newtonian flow problems. The theoretical framework developed here yields an estimate bound, which measures error in the approximate velocity in the W 1,r(Ω) norm and that of the pressure in the L r' (Ω) (1/r + 1/r' = 1). The power law model and the Carreau model are special ones of the quasi-Newtonian flow problem discussed in this paper. Moreover, a residual-based posterior bound is given. Numerical experiments are presented to confirm the theoretical results.展开更多
In this paper,we introduce and analyze an augmented mixed discontinuous Galerkin(MDG)method for a class of quasi-Newtonian Stokes flows.In the mixed formulation,the unknowns are strain rate,stress and velocity,which a...In this paper,we introduce and analyze an augmented mixed discontinuous Galerkin(MDG)method for a class of quasi-Newtonian Stokes flows.In the mixed formulation,the unknowns are strain rate,stress and velocity,which are approximated by a discontinuous piecewise polynomial triplet ■for k≥0.Here,the discontinuous piecewise polynomial function spaces for the field of strain rate and the stress field are designed to be symmetric.In addition,the pressure is easily recovered through simple postprocessing.For the benefit of the analysis,we enrich the MDG scheme with the constitutive equation relating the stress and the strain rate,so that the well-posedness of the augmented formulation is obtained by a nonlinear functional analysis.For k≥0,we get the optimal convergence order for the stress in broken ■(div)-norm and velocity in L^(2)-norm.Furthermore,the error estimates of the strain rate and the stress in-norm,and the pressure in L^(2)-norm are optimal under certain conditions.Finally,several numerical examples are given to show the performance of the augmented MDG method and verify the theoretical results.Numerical evidence is provided to show that the orders of convergence are sharp.展开更多
This paper proposes a weak Galerkin finite element method to solve incompressible quasi-Newtonian Stokes equations. We use piecewise polynomials of degrees k + 1(k 0) and k for the velocity and pressure in the interio...This paper proposes a weak Galerkin finite element method to solve incompressible quasi-Newtonian Stokes equations. We use piecewise polynomials of degrees k + 1(k 0) and k for the velocity and pressure in the interior of elements, respectively, and piecewise polynomials of degrees k and k + 1 for the boundary parts of the velocity and pressure, respectively. Wellposedness of the discrete scheme is established. The method yields a globally divergence-free velocity approximation. Optimal priori error estimates are derived for the velocity gradient and pressure approximations. Numerical results are provided to confirm the theoretical results.展开更多
Bayesian model averaging(BMA) is a recently proposed statistical method for calibrating forecast ensembles from numerical weather models.However,successful implementation of BMA requires accurate estimates of the weig...Bayesian model averaging(BMA) is a recently proposed statistical method for calibrating forecast ensembles from numerical weather models.However,successful implementation of BMA requires accurate estimates of the weights and variances of the individual competing models in the ensemble.Two methods,namely the Expectation-Maximization(EM) and the Markov Chain Monte Carlo(MCMC) algorithms,are widely used for BMA model training.Both methods have their own respective strengths and weaknesses.In this paper,we first modify the BMA log-likelihood function with the aim of removing the addi-tional limitation that requires that the BMA weights add to one,and then use a limited memory quasi-Newtonian algorithm for solving the nonlinear optimization problem,thereby formulating a new approach for BMA(referred to as BMA-BFGS).Several groups of multi-model soil moisture simulation experiments from three land surface models show that the performance of BMA-BFGS is similar to the MCMC method in terms of simulation accuracy,and that both are superior to the EM algo-rithm.On the other hand,the computational cost of the BMA-BFGS algorithm is substantially less than for MCMC and is al-most equivalent to that for EM.展开更多
基金Project supported by the Key Technology Research and Development Program of Sichuan Province of China(No.05GG006-006-2)
文摘For a generalized quasi-Newtonian flow, a new stabilized method focused on the low-order velocity-pressure pairs, (bi)linear/(bi)linear and (bi)linear/constant element, is presented. The pressure projection stabilized method is extended from Stokes problems to quasi-Newtonian flow problems. The theoretical framework developed here yields an estimate bound, which measures error in the approximate velocity in the W 1,r(Ω) norm and that of the pressure in the L r' (Ω) (1/r + 1/r' = 1). The power law model and the Carreau model are special ones of the quasi-Newtonian flow problem discussed in this paper. Moreover, a residual-based posterior bound is given. Numerical experiments are presented to confirm the theoretical results.
基金supported by the National Natural Science Foundation of China(Grant No.12171383)the National Natural Science Foundation of China(Grant No.11971377).
文摘In this paper,we introduce and analyze an augmented mixed discontinuous Galerkin(MDG)method for a class of quasi-Newtonian Stokes flows.In the mixed formulation,the unknowns are strain rate,stress and velocity,which are approximated by a discontinuous piecewise polynomial triplet ■for k≥0.Here,the discontinuous piecewise polynomial function spaces for the field of strain rate and the stress field are designed to be symmetric.In addition,the pressure is easily recovered through simple postprocessing.For the benefit of the analysis,we enrich the MDG scheme with the constitutive equation relating the stress and the strain rate,so that the well-posedness of the augmented formulation is obtained by a nonlinear functional analysis.For k≥0,we get the optimal convergence order for the stress in broken ■(div)-norm and velocity in L^(2)-norm.Furthermore,the error estimates of the strain rate and the stress in-norm,and the pressure in L^(2)-norm are optimal under certain conditions.Finally,several numerical examples are given to show the performance of the augmented MDG method and verify the theoretical results.Numerical evidence is provided to show that the orders of convergence are sharp.
基金supported by Major Research Plan of National Natural Science Foundation of China (Grant No. 91430105)
文摘This paper proposes a weak Galerkin finite element method to solve incompressible quasi-Newtonian Stokes equations. We use piecewise polynomials of degrees k + 1(k 0) and k for the velocity and pressure in the interior of elements, respectively, and piecewise polynomials of degrees k and k + 1 for the boundary parts of the velocity and pressure, respectively. Wellposedness of the discrete scheme is established. The method yields a globally divergence-free velocity approximation. Optimal priori error estimates are derived for the velocity gradient and pressure approximations. Numerical results are provided to confirm the theoretical results.
基金supported by National Basic Research Program of China (Grant No. 2010CB428403)National Natural Science Foundation of China (Grant No.41075076)Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No.KZCX2-EW-QN207)
文摘Bayesian model averaging(BMA) is a recently proposed statistical method for calibrating forecast ensembles from numerical weather models.However,successful implementation of BMA requires accurate estimates of the weights and variances of the individual competing models in the ensemble.Two methods,namely the Expectation-Maximization(EM) and the Markov Chain Monte Carlo(MCMC) algorithms,are widely used for BMA model training.Both methods have their own respective strengths and weaknesses.In this paper,we first modify the BMA log-likelihood function with the aim of removing the addi-tional limitation that requires that the BMA weights add to one,and then use a limited memory quasi-Newtonian algorithm for solving the nonlinear optimization problem,thereby formulating a new approach for BMA(referred to as BMA-BFGS).Several groups of multi-model soil moisture simulation experiments from three land surface models show that the performance of BMA-BFGS is similar to the MCMC method in terms of simulation accuracy,and that both are superior to the EM algo-rithm.On the other hand,the computational cost of the BMA-BFGS algorithm is substantially less than for MCMC and is al-most equivalent to that for EM.