共查询到20条相似文献,搜索用时 15 毫秒
1.
The two-grid method is studied for solving a two-dimensional second-order nonlinear hyperbolic equation using finite volume element method. The method is based on two different finite element spaces defined on one coarse grid with grid size H and one fine grid with grid size h, respectively. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. A prior error estimate in the H1-norm is proved to be O(h+H3|lnH|) for the two-grid semidiscrete finite volume element method. With these proposed techniques, solving such a large class of second-order nonlinear hyperbolic equations will not be much more difficult than solving one single linearized equation. Finally, a numerical example is presented to validate the usefulness and efficiency of the method. 相似文献
2.
Robert Eymard Thierry Gallouït Raphaèle Herbin Anthony Michel 《Numerische Mathematik》2002,92(1):41-82
Summary. One approximates the entropy weak solution u of a nonlinear parabolic degenerate equation by a piecewise constant function using a discretization in space and time and a finite volume scheme. The convergence of to u is shown as the size of the space and time steps tend to zero. In a first step, estimates on are used to prove the convergence, up to a subsequence, of to a measure valued entropy solution (called here an entropy process solution). A result of uniqueness of the entropy process
solution is proved, yielding the strong convergence of to{\it u}. Some on a model equation are shown.
Received September 27, 2000 / Published online October 17, 2001 相似文献
3.
To solve spatially semidiscrete approximative solution of a class of semilinear hyperbolic equations, the finite element method (FEM) with interpolated coefficients is discussed. By use of semidiscrete finite element for linear problem as comparative function, the error estimate in L∞-norm is derived by the nonlinear argument in Chen [Structure theory of superconvergence of finite elements, Hunan Press of Science and Technology, Changsha, 2001 (in Chinese)]. This indicates that convergence of FEMs with interpolated coefficients for a semilinear equation is similar to that of classical FEMs. 相似文献
4.
In this paper, high-resolution finite volume schemes are combined with an adaptive mesh technique inspired by multiresolution analysis to improve the computational efficiency for two-dimensional hyperbolic conservation laws. The method is conservative. Moreover, it is stable which is proven numerically in this paper. The computational grid is dynamically adapted so that higher spatial resolution is automatically allocated to regions where strong gradients are observed. Using this proposed scheme, we compute several two-dimensional model problems and a compressive rate ranging from about 5–10 is observed in all simulations. 相似文献
5.
Two-grid methods are studied for solving a two dimensional nonlinear parabolic equation using finite volume element method. The methods are based on one coarse-grid space and one fine-grid space. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine-grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. The two-grid methods achieve asymptotically optimal approximation as long as the mesh sizes satisfy h=O(H3|lnH|). As a result, solving such a large class of nonlinear parabolic equations will not be much more difficult than solving one single linearized equation. 相似文献
6.
In this paper we present an analysis of a numerical method for a degenerate partial differential equation, called the Black–Scholes
equation, governing American and European option pricing. The method is based on a fitted finite volume spatial discretization
and an implicit time stepping technique. The analysis is performed within the framework of the vertical method of lines, where
the spatial discretization is formulated as a Petrov–Galerkin finite element method with each basis function of the trial
space being determined by a set of two-point boundary value problems. We establish the stability and an error bound for the
solutions of the fully discretized system. Numerical results are presented to validate the theoretical results. 相似文献
7.
In this article, a class of nonlinear evolution equations – reaction–diffusion equations with time delay – is studied. By
combining the domain decomposition technique and the finite difference method, the results for the existence, convergence
and the stability of the numerical solution are obtained in the case of subdomain overlap and when the time-space is completely
discretized.
This revised version was published online in June 2006 with corrections to the Cover Date. 相似文献
8.
Polina Vinogradova 《Journal of Computational and Applied Mathematics》2009,231(1):1-10
This article investigates the projection-difference method for a Cauchy problem for a linear operator-differential equation with a leading self-adjoint operator A(t) and a subordinate linear operator K(t) in Hilbert space. This method leads to the solution of a system of linear algebraic equations on each time level; moreover, the projection subspaces are linear spans of eigenvectors of an operator similar to A(t). The convergence estimates are obtained. The application of the developed method for solving the initial boundary value problem is given. 相似文献
9.
This paper addresses the finite element method for the two-dimensional time-dependent Schrödinger equation on an infinite strip by using artificial boundary conditions. We first reduce the original problem into an initial-boundary value problem in a bounded domain by introducing a transparent boundary condition, then fully discretize this reduced problem by applying the Crank-Nicolson scheme in time and a bilinear or quadratic finite element approximation in space. This scheme, by a rigorous analysis, has been proved to be unconditionally stable and convergent, and its convergence order has also been obtained. Finally, two numerical examples are given to verify the accuracy of the scheme. 相似文献
10.
Kenneth Hvistendahl Karlsen Siddhartha Mishra Nils Henrik Risebro 《Numerische Mathematik》2009,111(4):559-589
We consider non-strictly hyperbolic systems of conservation laws in triangular form, which arise in applications like three-phase
flows in porous media. We device simple and efficient finite volume schemes of Godunov type for these systems that exploit
the triangular structure. We prove that the finite volume schemes converge to weak solutions as the discretization parameters
tend to zero. Some numerical examples are presented, one of which is related to flows in porous media.
The research of K. H. Karlsen was supported by an Outstanding Young Investigators Award from the Research Council of Norway. 相似文献
11.
We study the convergence of a finite volume scheme for the linear advection equation with a Lipschitz divergence-free speed
in R
d
. We prove a h
1/2-error estimate in the L
∞(0,t;L
1)-norm for BV data. This result was expected from numerical experiments and is optimal. 相似文献
12.
Julien Vovelle 《Numerische Mathematik》2002,90(3):563-596
Summary. This paper is devoted to the study of the finite volume methods used in the discretization of conservation laws defined on
bounded domains. General assumptions are made on the data: the initial condition and the boundary condition are supposed to
be measurable bounded functions. Using a generalized notion of solution to the continuous problem (namely the notion of entropy
process solution, see [9]) and a uniqueness result on this solution, we prove that the numerical solution converges to the
entropy weak solution of the continuous problem in for every . This also yields a new proof of the existence of an entropy weak solution.
Received May 18, 2000 / Revised version received November 21, 2000 / Published online June 7, 2001 相似文献
13.
Francis Filbet 《Numerische Mathematik》2006,104(4):457-488
A finite volume method is presented to discretize the Patlak–Keller–Segel (PKS) modeling chemosensitive movements. First, we prove existence and uniqueness of a numerical solution to the proposed scheme. Then, we give a priori estimates and establish a threshold on the initial mass, for which we show that the numerical approximation converges to the solution to the PKS system when the initial mass is lower than this threshold. Numerical simulations are performed to verify accuracy and the properties of the scheme. Finally, in the last section we investigate blow-up of the solution for large mass. 相似文献
14.
In this paper, we propose a new method to compute the numerical flux of a finite volume scheme, used for the approximation
of the solution of the nonlinear partial differential equation ut+div(qf(u))−ΔΦ(u)=0 in a 1D, 2D or 3D domain. The function Φ is supposed to be strictly increasing, but some values s such that Φ′(s)=0 can exist. The method is based on the solution, at each interface between two control volumes, of the nonlinear elliptic
two point boundary value problem (qf(υ)+(Φ(υ))′)′=0 with Dirichlet boundary conditions given by the values of the discrete approximation in both control volumes. We prove
the existence of a solution to this two point boundary value problem. We show that the expression for the numerical flux can
be yielded without referring to this solution. Furthermore, we prove that the so designed finite volume scheme has the expected
stability properties and that its solution converges to the weak solution of the continuous problem. Numerical results show
the increase of accuracy due to the use of this scheme, compared to some other schemes. 相似文献
15.
16.
The continuous Galerkin finite element method for linear delay-differential equation with several terms is studied. Adding some lower terms in the remainder of orthogonal expansion in an element so that the remainder satisfies more orthogonal condition in the element, and obtain a desired superclose function to finite element solution, thus the superconvergence of p -degree finite element approximate solution on (p+1)-order Lobatto points is derived. 相似文献
17.
In this work, we introduce the convergence analysis of the recently developed finite volume scheme to solve a pure aggregation population balance equation that is of substantial interest in many areas such as chemical engineering, aerosol physics, astrophysics, polymer science, pharmaceutical sciences, and mathematical biology. The notion of the finite volume scheme is to conserve total mass of the particles in the system by introducing weight in the formulation. The consistency of the finite volume scheme is also analyzed thoroughly as it is an influential factor. The convergence study of the numerical scheme shows second order convergence on uniform, nonuniform smooth (geometric) as well as on locally uniform meshes independent of the aggregation kernel. Moreover, the first‐order convergence is shown when the finite volume scheme is implemented on oscillatory and random meshes. In order to check the accuracy, the numerical experimental order of convergence is also computed for the physically relevant as well as analytically tractable kernels and validated against its analytical results. 相似文献
18.
Summary. A new characteristic finite element scheme is presented for It is of second order accuracy in time increment, symmetric,
and unconditionally stable. Optimal error estimates are proved in the framework of -theory. Numerical results are presented for two examples, which show the advantage of the scheme.
Received November 22, 2000 / Revised version received July 11, 2001 / Published online October 17, 2001 相似文献
19.
Daniel Bouche Jean-Michel Ghidaglia 《Journal of Computational and Applied Mathematics》2011,235(18):5394-5410
In this paper we estimate the error of upwind first order finite volume schemes applied to scalar conservation laws. As a first step, we consider standard upwind and flux finite volume scheme discretization of a linear equation with space variable coefficients in conservation form. We prove that, in spite of their lack of consistency, both schemes lead to a first order error estimate. As a final step, we prove a similar estimate for the nonlinear case. Our proofs rely on the notion of geometric corrector, introduced in our previous paper by Bouche et al. (2005) [24] in the context of constant coefficient linear advection equations. 相似文献
20.
An adaptive semi-Lagrangian scheme for solving the Cauchy problem associated to the periodic 1+1-dimensional Vlasov-Poisson
system in the two- dimensional phase space is proposed and analyzed. A key feature of our method is the accurate evolution
of the adaptive mesh from one time step to the next one, based on a rigorous analysis of the local regularity and how it gets
transported by the numerical flow. The accuracy of the scheme is monitored by a prescribed tolerance parameter ε which represents the local interpolation error at each time step, in the L
∞ metric. The numerical solutions are proved to converge in L
∞ towards the exact ones as ε and Δt tend to zero provided the initial data is Lipschitz and has a finite total curvature, or in other words, that it belongs
to . The rate of convergence is , which should be compared to the results of Besse who recently established in (SIAM J Numer Anal 42(1):350–382, 2004) similar
rates for a uniform semi-Lagrangian scheme, but requiring that the initial data are in . Several numerical tests illustrate the effectiveness of our approach for generating the optimal adaptive discretizations. 相似文献