Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFinite difference approach to fourth-order linear boundary-value problems
93
0
0.0
(
0
)
Added by
Matania Ben Artzi
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Discrete approximations to the equation begin{equation*} L_{cont}u = u^{(4)} + D(x) u^{(3)} + A(x) u^{(2)} + (A(x)+H(x)) u^{(1)} + B(x) u = f, ; xin[0,1] end{equation*} are considered. This is an extension of the Sturm-Liouville case $D(x)equiv H(x)equiv 0$ [ M. Ben-Artzi, J.-P. Croisille, D. Fishelov and R. Katzir, Discrete fourth-order Sturm-Liouville problems, IMA J. Numer. Anal. {bf 38} (2018), 1485-1522. doi: 10.1093/imanum/drx038] to the non-self-adjoint setting. The natural boundary conditions in the Sturm-Liouville case are the values of the function and its derivative. The inclusion of a third-order discrete derivative entails a revision of the underlying discrete functional calculus. This revision forces evaluations of accurate discrete approximations to the boundary values of the second, third and fourth order derivatives. The resulting functional calculus provides the discrete analogs of the fundamental Sobolev properties--compactness and coercivity. It allows to obtain a general convergence theorem of the discrete approximations to the exact solution. Some representative numerical examples are presented.
rate research
Read More
Many practical problems occur due to the boundary value problem. This paper evaluates the finite element solution of the boundary value problem of Poissons equation and proposes a novel a posteriori local error estimation based on the Hypercircle method. Compared to the existing literature on qualitative error estimation, the proposed error estimation provides an explicit and sharp bound for the approximation error in the subdomain of interest and is applicable to problems without the $H^2$ regularity. The efficiency of the proposed method is demonstrated by numerical experiments for both convex and non-convex 2D domains.
We develop a numerical method for solving the boundary value problem of The Linear Seventh Ordinary Boundary Value Problem by using seventh degree B-Spline function. Formulation is based on particular terms of order of seventh order boundary value problem. We obtain Septic B-Spline formulation and the Collocation B-spline Method is formulated as an approximation solution. We apply the presented method to solve an example of seventh-order boundary value problem which the results show that there is an agreement between approximate solutions and exact solutions. Resulting low absolute errors show that the presented numerical method is effective for solving high order boundary value problems. Finally, a general conclusion has been included.
In this paper, we develop a high order residual distribution (RD) method for solving steady state conservation laws in a novel Hermite weighted essentially non-oscillatory (HWENO) framework recently developed in [23]. In particular, we design a high order HWENO reconstructions for the integrals of source term and fluxes based on the point values of the solution and its spatial derivatives, and the principles of residual distribution schemes are adapted to obtain steady state solutions. The proposed novel HWENO framework enjoys two advantages. First, compared with the traditional HWENO framework, the proposed methods do not need to introduce additional auxiliary equations to update the derivatives of the unknown function, and compute them from the current value and the old spatial derivatives. This approach saves the computational storage and CPU time, which greatly improves the computational efficiency of the traditional HWENO framework. Second, compared with the traditional WENO method, reconstruction stencil of the HWENO methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller at the same grid. Thus, it is also a compact scheme when we design the higher order accuracy, compared with that in [11] Chou and Shu proposed. Extensive numerical experiments for one and two-dimensional scalar and systems problems confirm the high order accuracy and good quality of our scheme.
We develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.
In this paper, we present a unified analysis of the superconvergence property for a large class of mixed discontinuous Galerkin methods. This analysis applies to both the Poisson equation and linear elasticity problems with symmetric stress formulations. Based on this result, some locally postprocess schemes are employed to improve the accuracy of displacement by order min(k+1, 2) if polynomials of degree k are employed for displacement. Some numerical experiments are carried out to validate the theoretical results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
