scholarly journals An Eulerian finite element method for PDEs in time-dependent domains

2019 ◽  
Vol 53 (2) ◽  
pp. 585-614 ◽  
Author(s):  
Christoph Lehrenfeld ◽  
Maxim Olshanskii

The paper introduces a new finite element numerical method for the solution of partial differential equations on evolving domains. The approach uses a completely Eulerian description of the domain motion. The physical domain is embedded in a triangulated computational domain and can overlap the time-independent background mesh in an arbitrary way. The numerical method is based on finite difference discretizations of time derivatives and a standard geometrically unfitted finite element method with an additional stabilization term in the spatial domain. The performance and analysis of the method rely on the fundamental extension result in Sobolev spaces for functions defined on bounded domains. This paper includes a complete stability and error analysis, which accounts for discretization errors resulting from finite difference and finite element approximations as well as for geometric errors coming from a possible approximate recovery of the physical domain. Several numerical examples illustrate the theory and demonstrate the practical efficiency of the method.

Mathematics ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 1382
Author(s):  
Denis Spiridonov ◽  
Maria Vasilyeva ◽  
Aleksei Tyrylgin ◽  
Eric T. Chung

In this paper, we present a multiscale model reduction technique for unsaturated filtration problem in fractured porous media using an Online Generalized Multiscale finite element method. The flow problem in unsaturated soils is described by the Richards equation. To approximate fractures we use the Discrete Fracture Model (DFM). Complex geometric features of the computational domain requires the construction of a fine grid that explicitly resolves the heterogeneities such as fractures. This approach leads to systems with a large number of unknowns, which require large computational costs. In order to develop a more efficient numerical scheme, we propose a model reduction procedure based on the Generalized Multiscale Finite element method (GMsFEM). The GMsFEM allows solving such problems on a very coarse grid using basis functions that can capture heterogeneities. In the GMsFEM, there are offline and online stages. In the offline stage, we construct snapshot spaces and solve local spectral problems to obtain multiscale basis functions. These spectral problems are defined in the snapshot space in each local domain. To improve the accuracy of the method, we add online basis functions in the online stage. The construction of the online basis functions is based on the local residuals. The use of online bases will allow us to get a significant improvement in the accuracy of the method. We present results with different number of offline and online multisacle basis functions. We compare all results with reference solution. Our results show that the proposed method is able to achieve high accuracy with a small computational cost.


1983 ◽  
Vol 14 (2) ◽  
pp. 85-92 ◽  
Author(s):  
Tilahun Aberra

The numerical solution of the behaviour of discrete time steps in digital computer analysis of square aquifers containing pumped wells is examined by using the finite element method with a 4 node linear quadrilateral isoparametric surface element. A wide range of time steps are used in the computation. The calculations show that discrete time steps can cause errors and oscillations in the calculations particularly when wells start and stop pumping. Comparison with known results obtained by theoretical and finite difference procedures has been considered. The main objective of this paper is to demonstrate comparison of the finite element and finite difference simulation results over a regular linear 4 node quadrilateral mesh suitable to represent the two numerical schemes with a marked similarity. The dimensionless time drawdown results of the finite element method agreed well with the finite difference and analytical results for small time increment. However, for large time increments, there are from slight to significant oscillations in the results and notable discrepancies are observed in the solutions of the two numerical methods.


Author(s):  
Xuejuan Li ◽  
Ji-Huan He

Purpose The purpose of this paper is to develop an effective numerical algorithm for a gas-melt two-phase flow and use it to simulate a polymer melt filling process. Moreover, the suggested algorithm can deal with the moving interface and discontinuities of unknowns across the interface. Design/methodology/approach The algebraic sub-grid scales-variational multi-scale (ASGS-VMS) finite element method is used to solve the polymer melt filling process. Meanwhile, the time is discretized using the Crank–Nicolson-based split fractional step algorithm to reduce the computational time. The improved level set method is used to capture the melt front interface, and the related equations are discretized by the second-order Taylor–Galerkin scheme in space and the third-order total variation diminishing Runge–Kutta scheme in time. Findings The numerical method is validated by the benchmark problem. Moreover, the viscoelastic polymer melt filling process is investigated in a rectangular cavity. The front interface, pressure field and flow-induced stresses of polymer melt during the filling process are predicted. Overall, this paper presents a VMS method for polymer injection molding. The present numerical method is extremely suitable for two free surface problems. Originality/value For the first time ever, the ASGS-VMS finite element method is performed for the two-phase flow of polymer melt filling process, and an effective numerical method is designed to catch the moving surface.


Sign in / Sign up

Export Citation Format

Share Document