Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA numerical implementation of the unified Fokas transform for evolution problems on a finite interval
101
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present the numerical solution of two-point boundary value problems for a third order linear PDE, representing a linear evolution in one space dimension. The difficulty of this problem is in the numerical imposition of the boundary conditions, and to our knowledge, no such computations exist. Instead of computing the evolution numerically, we evaluate the solution representation formula obtained by the unified transform of Fokas. This representation involves complex line integrals, but in order to evaluate these integrals numerically, it is necessary to deform the integration contours using appropriate deformation mappings. We formulate a strategy to implement effectively this deformation, which allows us to obtain accurate numerical results.
rate research
Read More
In this work, we describe a simple finite element approach that is able to resolve weak discontinuities in interface problems accurately. The approach is based on a fixed patch mesh consisting of quadrilaterals, that will stay unchanged independent of the position of the interface. Inside the patches we refine once more, either in eight triangles or in four quadrilaterals, in such a way that the interface is locally resolved. The resulting finite element approach can be considered a fitted finite element approach. In our practical implementation, we do not construct this fitted mesh, however. Instead, the local degrees of freedom are included in a parametric way in the finite element space, or to be more precise in the local mappings between a reference patch and the physical patches. We describe the implementation in the open source C++ finite element library deal.II in detail and present two numerical examples to illustrate the performance of the approach. Finally, detailed studies of the behavior of iterative linear solvers complement this work.
In this paper we develop a class of efficient Galerkin boundary element methods for the solution of two-dimensional exterior single-scattering problems. Our approach is based upon construction of Galerkin approximation spaces confined to the asymptotic behaviour of the solution through a certain direct sum of appropriate function spaces weighted by the oscillations in the incident field of radiation. Specifically, the function spaces in the illuminated/shadow regions and the shadow boundaries are simply algebraic polynomials whereas those in the transition regions are generated utilizing novel, yet simple, emph{frequency dependent changes of variables perfectly matched with the boundary layers of the amplitude} in these regions. While, on the one hand, we rigorously verify for smooth convex obstacles that these methods require only an $mathcal{O}left( k^{epsilon} right)$ increase in the number of degrees of freedom to maintain any given accuracy independent of frequency, and on the other hand, remaining in the realm of smooth obstacles they are applicable in more general single-scattering configurations. The most distinctive property of our algorithms is their emph{remarkable success} in approximating the solution in the shadow region when compared with the algorithms available in the literature.
The scaled boundary finite element method (SBFEM) is a semi-analytical computational scheme, which is based on the characteristics of the finite element method (FEM) and combines the advantages of the boundary element method (BEM). This paper integrates the scaled boundary finite element method (SBFEM) and the polygonal mesh technique into a new approach to solving the steady-state and transient seepage problems. The proposed method is implemented in Abaqus using a user-defined element (UEL). The detailed implementations of the procedure, defining the UEL element, updating the RHS and AMATRX, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from seepage are solved to validate the proposed implementation. Results show that the polygonal element of PS-SBFEM has a higher accuracy rate than the standard FEM element in the same element size. For the transient problems, the results between PS-SBFEM and the FEM are in excellent agreement. Furthermore, the PS-SBFEM with quadtree meshes shows a good effect for solving complex geometric in the seepage problem. Hence, the proposed method is robust accurate for solving the steady-state and transient seepage problems. The developed UEL source code and the associated input files can be downloaded from https://github.com/yangyLab/PS-SBFEM.
We study the large time behaviour of the solution of linear dispersive partial differential equations posed on a finite interval, when at least one of the prescribed boundary conditions is time periodic. We use the Q equation approach, pioneered in Fokas & Lenells 2012 and applied to linear problems on the half-line in Fokas & van der Weele 2021, to characterise necessary conditions for the solution of such problem to be periodic, at least in an asymptotic sense. We then fully describe the periodicity properties of the solution in three important illustrative examples, recovering known results for the second-order cases and establishing new results for the third order case.
We derive a posteriori error estimates in the $L_infty((0,T];L_infty(Omega))$ norm for approximations of solutions to linear para bolic equations. Using the elliptic reconstruction technique introduced by Makridakis and Nochetto and heat kernel estimates for linear parabolic pr oblems, we first prove a posteriori bounds in the maximum norm for semidiscrete finite element approximations. We then establish a posteriori bounds for a fully discrete backward Euler finite element approximation. The elliptic reconstruction technique greatly simplifies our development by allow ing the straightforward combination of heat kernel estimates with existing elliptic maximum norm error estimators.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
