The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensio...The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.展开更多
In this work,we are concerned with the stability and convergence analysis of the second-order backward difference formula(BDF2)with variable steps for the molecular beam epitaxial model without slope selection.We firs...In this work,we are concerned with the stability and convergence analysis of the second-order backward difference formula(BDF2)with variable steps for the molecular beam epitaxial model without slope selection.We first show that the variable-step BDF2 scheme is convex and uniquely solvable under a weak time-step constraint.Then we show that it preserves an energy dissipation law if the adjacent time-step ratios satisfy r_(k):=τ_(k)/τ_(k-1)<3.561.Moreover,with a novel discrete orthogonal convolution kernels argument and some new estimates on the corresponding positive definite quadratic forms,the L^(2)norm stability and rigorous error estimates are established,under the same step-ratio constraint that ensures the energy stability,i.e.,0<r_(k)<3.561.This is known to be the best result in the literature.We finally adopt an adaptive time-stepping strategy to accelerate the computations of the steady state solution and confirm our theoretical findings by numerical examples.展开更多
The positive definiteness of real quadratic forms with convolution structures plays an important rolein stability analysis for time-stepping schemes for nonlocal operators. In this work, we present a novel analysistoo...The positive definiteness of real quadratic forms with convolution structures plays an important rolein stability analysis for time-stepping schemes for nonlocal operators. In this work, we present a novel analysistool to handle discrete convolution kernels resulting from variable-step approximations for convolution operators.More precisely, for a class of discrete convolution kernels relevant to variable-step L1-type time discretizations, weshow that the associated quadratic form is positive definite under some easy-to-check algebraic conditions. Ourproof is based on an elementary constructing strategy by using the properties of discrete orthogonal convolutionkernels and discrete complementary convolution kernels. To our knowledge, this is the first general result onsimple algebraic conditions for the positive definiteness of variable-step discrete convolution kernels. Using theunified theory, we obtain the stability for some simple nonuniform time-stepping schemes straightforwardly.展开更多
文摘The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.
基金supported by National Natural Science Foundation of China(Grant No.12071216)supported by National Natural Science Foundation of China(Grant No.11731006)+2 种基金the NNW2018-ZT4A06 projectsupported by National Natural Science Foundation of China(Grant Nos.11822111,11688101 and 11731006)the Science Challenge Project(Grant No.TZ2018001)。
文摘In this work,we are concerned with the stability and convergence analysis of the second-order backward difference formula(BDF2)with variable steps for the molecular beam epitaxial model without slope selection.We first show that the variable-step BDF2 scheme is convex and uniquely solvable under a weak time-step constraint.Then we show that it preserves an energy dissipation law if the adjacent time-step ratios satisfy r_(k):=τ_(k)/τ_(k-1)<3.561.Moreover,with a novel discrete orthogonal convolution kernels argument and some new estimates on the corresponding positive definite quadratic forms,the L^(2)norm stability and rigorous error estimates are established,under the same step-ratio constraint that ensures the energy stability,i.e.,0<r_(k)<3.561.This is known to be the best result in the literature.We finally adopt an adaptive time-stepping strategy to accelerate the computations of the steady state solution and confirm our theoretical findings by numerical examples.
基金Hong-Lin Liao was supported by National Natural Science Foundation of China(Grant No.12071216)Tao Tang was supported by Science Challenge Project(Grant No.TZ2018001)+3 种基金National Natural Science Foundation of China(Grants Nos.11731006 and K20911001)Tao Zhou was supported by National Natural Science Foundation of China(Grant No.12288201)Youth Innovation Promotion Association(CAS)Henan Academy of Sciences.
文摘The positive definiteness of real quadratic forms with convolution structures plays an important rolein stability analysis for time-stepping schemes for nonlocal operators. In this work, we present a novel analysistool to handle discrete convolution kernels resulting from variable-step approximations for convolution operators.More precisely, for a class of discrete convolution kernels relevant to variable-step L1-type time discretizations, weshow that the associated quadratic form is positive definite under some easy-to-check algebraic conditions. Ourproof is based on an elementary constructing strategy by using the properties of discrete orthogonal convolutionkernels and discrete complementary convolution kernels. To our knowledge, this is the first general result onsimple algebraic conditions for the positive definiteness of variable-step discrete convolution kernels. Using theunified theory, we obtain the stability for some simple nonuniform time-stepping schemes straightforwardly.
