scholarly journals Energy-corrected FEM and explicit time-stepping for parabolic problems

2019 ◽  
Vol 53 (6) ◽  
pp. 1893-1914
Author(s):  
Piotr Swierczynski ◽  
Barbara Wohlmuth

The presence of corners in the computational domain, in general, reduces the regularity of solutions of parabolic problems and diminishes the convergence properties of the finite element approximation introducing a so-called “pollution effect”. Standard remedies based on mesh refinement around the singular corner result in very restrictive stability requirements on the time-step size when explicit time integration is applied. In this article, we introduce and analyse the energy-corrected finite element method for parabolic problems, which works on quasi-uniform meshes, and, based on it, create fast explicit time discretisation. We illustrate these results with extensive numerical investigations not only confirming the theoretical results but also showing the flexibility of the method, which can be applied in the presence of multiple singular corners and a three-dimensional setting. We also propose a fast explicit time-stepping scheme based on a piecewise cubic energy-corrected discretisation in space completed with mass-lumping techniques and numerically verify its efficiency.

Author(s):  
W. Dauksher ◽  
A. F. Emery

The dispersive errors in the finite element solution to the one- and two-dimensional heat equations are examined as a function of element type and size, capacitance matrix formulation, time stepping scheme and time step size.


2016 ◽  
Vol 2016 ◽  
pp. 1-13 ◽  
Author(s):  
S. S. Ravindran

Micropolar fluid model consists of Navier-Stokes equations and microrotational velocity equations describing the dynamics of flows in which microstructure of fluid is important. In this paper, we propose and analyze a decoupled time-stepping algorithm for the evolutionary micropolar flow. The proposed method requires solving only one uncoupled Navier-Stokes and one microrotation subphysics problem per time step. We derive optimal order error estimates in suitable norms without assuming any stability condition or time step size restriction.


2020 ◽  
Vol 20 (2) ◽  
pp. 361-378
Author(s):  
Tamal Pramanick ◽  
Rajen Kumar Sinha

AbstractThe purpose of this paper is to generalize known a priori error estimates of the composite finite element (CFE) approximations of elliptic problems in nonconvex polygonal domains to the time dependent parabolic problems. This is a new class of finite elements which was introduced by [W. Hackbusch and S. A. Sauter, Composite finite elements for the approximation of PDEs on domains with complicated micro-structures, Numer. Math. 75 1997, 4, 447–472] and subsequently modified by [M. Rech, S. A. Sauter and A. Smolianski, Two-scale composite finite element method for Dirichlet problems on complicated domains, Numer. Math. 102 2006, 4, 681–708] for the approximations of stationery problems on complicated domains. The basic idea of the CFE procedure is to work with fewer degrees of freedom by allowing finite element mesh to resolve the domain boundaries and to preserve the asymptotic order convergence on coarse-scale mesh. We analyze both semidiscrete and fully discrete CFE methods for parabolic problems in two-dimensional nonconvex polygonal domains and derive error estimates of order {\mathcal{O}(H^{s}\widehat{\mathrm{Log}}{}^{\frac{s}{2}}(\frac{H}{h}))} and {\mathcal{O}(H^{2s}\widehat{\mathrm{Log}}{}^{s}(\frac{H}{h}))} in the {L^{\infty}(H^{1})}-norm and {L^{\infty}(L^{2})}-norm, respectively. Moreover, for homogeneous equations, error estimates are derived for nonsmooth initial data. Numerical results are presented to support the theoretical rates of convergence.


2012 ◽  
Vol 12 (3) ◽  
pp. 789-806 ◽  
Author(s):  
P. Sváček ◽  
J. Horáček

AbstractIn this paper the numerical method for solution of an aeroelastic model describing the interactions of air flow with vocal folds is described. The flow is modelled by the incompressible Navier-Stokes equations spatially discretized with the aid of the stabilized finite element method. The motion of the computational domain is treated with the aid of the Arbitrary Lagrangian Eulerian method. The structure dynamics is replaced by a mechanically equivalent system with the two degrees of freedom governed by a system of ordinary differential equations and discretized in time with the aid of an implicit multistep method and strongly coupled with the flow model. The influence of inlet/outlet boundary conditions is studied and the numerical analysis is performed and compared to the related results from literature.


2013 ◽  
Vol 13 (2) ◽  
pp. 386-410 ◽  
Author(s):  
Björn Sjögreen ◽  
Jeffrey W. Banks

AbstractWe consider multi-physics computations where the Navier-Stokes equations of compressible fluid flow on some parts of the computational domain are coupled to the equations of elasticity on other parts of the computational domain. The different subdomains are separated by well-defined interfaces. We consider time accurate computations resolving all time scales. For such computations, explicit time stepping is very efficient. We address the issue of discrete interface conditions between the two domains of different physics that do not lead to instability, or to a significant reduction of the stable time step size. Finding such interface conditions is non-trivial.We discretize the problem with high order centered difference approximations with summation by parts boundary closure. We derive L2 stable interface conditions for the linearized one dimensional discretized problem. Furthermore, we generalize the interface conditions to the full non-linear equations and numerically demonstrate their stable and accurate performance on a simple model problem. The energy stable interface conditions derived here through symmetrization of the equations contain the interface conditions derived through normal mode analysis by Banks and Sjögreen in [8] as a special case.


2014 ◽  
Vol 7 (6) ◽  
pp. 3017-3035 ◽  
Author(s):  
P. A. Ullrich

Abstract. This paper presents a novel nodal finite-element method for either continuous and discontinuous elements, as applied to the 2-D shallow-water equations on the cubed sphere. The cornerstone of this method is the construction of a robust derivative operator that can be applied to compute discrete derivatives even over a discontinuous function space. A key advantage of the robust derivative is that it can be applied to partial differential equations in either a conservative or a non-conservative form. However, it is also shown that discontinuous penalization is required to recover the correct order of accuracy for discontinuous elements. Two versions with discontinuous elements are examined, using either the g1 and g2 flux correction function for distribution of boundary fluxes and penalty across nodal points. Scalar and vector hyperviscosity (HV) operators valid for both continuous and discontinuous elements are also derived for stabilization and removal of grid-scale noise. This method is validated using four standard shallow-water test cases, including geostrophically balanced flow, a mountain-induced Rossby wave train, the Rossby–Haurwitz wave and a barotropic instability. The results show that although the discontinuous basis requires a smaller time step size than that required for continuous elements, the method exhibits better stability and accuracy properties in the absence of hyperviscosity.


1973 ◽  
Vol 95 (1) ◽  
pp. 126-129 ◽  
Author(s):  
R. E. Beckett ◽  
S.-C. Chu

By use of an implicit iteration technique, the finite-element method applied to the heat-conduction problems of solids is no longer restricted to the linear heat-flux boundary conditions, but is extended to include nonlinear radiation–convection boundary conditions. The variation of surface temperatures within each time increment is taken into account; hence a rather large time-step size can be assigned to obtain transient heat-conduction solutions without introducing instability in the surface temperature of a body.


2018 ◽  
Vol 10 (1) ◽  
pp. 168781401774807 ◽  
Author(s):  
Yonghu Wang ◽  
Dongwei Shu ◽  
Yusaku Fujii ◽  
Akihiro Takita ◽  
Tsuneaki Ishima ◽  
...  

In order to precisely measure water impact loads of a spherical structure vertically dropping onto a calm water surface, a new validity check of the analysis using the levitation mass method experiment is proposed. The main feature of levitation mass method experiment is to obtain a better estimation of early water impact loads through the application of Doppler effect. Experimental results of different heights are verified based on the Assessment Index and are in comparison with the classical experimental data for validation purpose. It shows that the levitation mass method measurement is useful and effective to obtain the water impact loads for the crashworthiness analysis. Besides, early water impact hydrodynamic behaviors are simulated based on the nonlinear explicit finite element method, together with application of a multi-material arbitrary Lagrangian–Eulerian solver. A penalty coupling algorithm is utilized to realize fluid–structure interaction between the spherical body and fluids. Convergence studies are performed to construct the proper finite element model by the comparison with experimental results, where mesh sensitivity, contact stiffness, and time-step size parametric studies are thoroughly investigated. The comparisons between experimental and numerical results show good consistency by the prediction of the water impact coefficients on the structure.


Sign in / Sign up

Export Citation Format

Share Document