scholarly journals Implementing model reduction to the JuliaFEM platform

2018 ◽  
Vol 51 (1) ◽  
pp. 36-54 ◽  
Author(s):  
Marja Liisa Rapo ◽  
Jukka Aho ◽  
Hannu Koivurova ◽  
Tero Frondelius
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Model Reduction ◽  
Open Source ◽  
Degrees Of Freedom ◽  
Mass Matrices ◽  
Dynamic Condensation ◽  
Dynamic Analyses ◽  
Reduction Methods ◽  
Large Model

JuliaFEM is an open source finite element method solver written in the Julia language. This paper presents an implementation of two common model reduction methods: the Guyan reduction and the Craig-Bampton method. The goal was to implement these algorithms to the JuliaFEM platform and demonstrate that the code works correctly. This paper first describes the JuliaFEM concept briefly after which it presents the theory of model reduction, and finally, it demonstrates the implemented functions in an example model. This paper includes instructions for using the implemented algorithms, and reference the code itself in GitHub. The reduced stiness and mass matrices give the same results in both static and dynamic analyses as the original matrices, which proves that the code works correctly. The code runs smoothly on relatively large model of 12.6 million degrees of freedom. In future, damping could be included in the dynamic condensation now that it has been shown to work.

scholarly journals Higher Order Multiscale Finite Element Method for Heat Transfer Modeling

Materials ◽  
2021 ◽  
Vol 14 (14) ◽  
pp. 3827
Author(s):  
Marek Klimczak ◽  
Witold Cecot
Keyword(s):  
Heat Transfer ◽  
Finite Element Method ◽  
Finite Element ◽  
Degrees Of Freedom ◽  
Higher Order ◽  
Shape Functions ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Steady State Heat ◽  
Steady State Heat Transfer

In this paper, we present a new approach to model the steady-state heat transfer in heterogeneous materials. The multiscale finite element method (MsFEM) is improved and used to solve this problem. MsFEM is a fast and flexible method for upscaling. Its numerical efficiency is based on the natural parallelization of the main computations and their further simplifications due to the numerical nature of the problem. The approach does not require the distinct separation of scales, which makes its applicability to the numerical modeling of the composites very broad. Our novelty relies on modifications to the standard higher-order shape functions, which are then applied to the steady-state heat transfer problem. To the best of our knowledge, MsFEM (based on the special shape function assessment) has not been previously used for an approximation order higher than p = 2, with the hierarchical shape functions applied and non-periodic domains, in this problem. Some numerical results are presented and compared with the standard direct finite-element solutions. The first test shows the performance of higher-order MsFEM for the asphalt concrete sample which is subject to heating. The second test is the challenging problem of metal foam analysis. The thermal conductivity of air and aluminum differ by several orders of magnitude, which is typically very difficult for the upscaling methods. A very good agreement between our upscaled and reference results was observed, together with a significant reduction in the number of degrees of freedom. The error analysis and the p-convergence of the method are also presented. The latter is studied in terms of both the number of degrees of freedom and the computational time.

scholarly journals An Online Generalized Multiscale Finite Element Method for Unsaturated Filtration Problem in Fractured Media

Mathematics ◽  
2021 ◽  
Vol 9 (12) ◽  
pp. 1382
Author(s):  
Denis Spiridonov ◽  
Maria Vasilyeva ◽  
Aleksei Tyrylgin ◽  
Eric T. Chung
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Model Reduction ◽  
Basis Functions ◽  
Computational Domain ◽  
Filtration Problem ◽  
Spectral Problems ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Element Method

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.

Application of Extended Finite Element Method (XFEM) to Stress Intensity Factor Calculations

2015 ◽  
Author(s):  
Do-Jun Shim ◽  
Mohammed Uddin ◽  
Sureshkumar Kalyanam ◽  
Frederick Brust ◽  
Bruce Young
Keyword(s):  
Finite Element Method ◽  
Stress Intensity Factor ◽  
Finite Element ◽  
Stress Intensity ◽  
Intensity Factor ◽  
Degrees Of Freedom ◽  
Extended Finite Element Method ◽  
Partition Of Unity ◽  
Extended Finite Element ◽  
Element Method

The extended finite element method (XFEM) is an extension of the conventional finite element method based on the concept of partition of unity. In this method, the presence of a crack is ensured by the special enriched functions in conjunction with additional degrees of freedom. This approach also removes the requirement for explicitly defining the crack front or specifying the virtual crack extension direction when evaluating the contour integral. In this paper, stress intensity factors (SIF) for various crack types in plates and pipes were calculated using the XFEM embedded in ABAQUS. These results were compared against handbook solutions, results from conventional finite element method, and results obtained from finite element alternating method (FEAM). Based on these results, applicability of the ABAQUS XFEM to stress intensity factor calculations was investigated. Discussions are provided on the advantages and limitations of the XFEM.

scholarly journals Adaptive modelling of variably saturated seepage problems

2021 ◽  
Author(s):  
B Ashby ◽  
C Bortolozo ◽  
A Lukyanov ◽  
T Pryer
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Degrees Of Freedom ◽  
A Priori ◽  
Water Flux ◽  
Posteriori Error ◽  
Adaptive Finite Element Method ◽  
Posteriori Error Estimate ◽  
Adaptive Finite Element ◽  
Element Method

Summary In this article, we present a goal-oriented adaptive finite element method for a class of subsurface flow problems in porous media, which exhibit seepage faces. We focus on a representative case of the steady state flows governed by a nonlinear Darcy–Buckingham law with physical constraints on subsurface-atmosphere boundaries. This leads to the formulation of the problem as a variational inequality. The solutions to this problem are investigated using an adaptive finite element method based on a dual-weighted a posteriori error estimate, derived with the aim of reducing error in a specific target quantity. The quantity of interest is chosen as volumetric water flux across the seepage face, and therefore depends on an a priori unknown free boundary. We apply our method to challenging numerical examples as well as specific case studies, from which this research originates, illustrating the major difficulties that arise in practical situations. We summarise extensive numerical results that clearly demonstrate the designed method produces rapid error reduction measured against the number of degrees of freedom.

scholarly journals Quasi-3-D spectral wavelet method for a thermal quench simulation

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Jonas Bundschuh ◽  
Laura A. M. D’Angelo ◽  
Herbert De Gersem
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Spectral Method ◽  
Degrees Of Freedom ◽  
Hot Spot ◽  
Longitudinal Direction ◽  
Particle Accelerators ◽  
Wavelet Basis ◽  
Element Simulation ◽  
Element Method

AbstractThe finite element method is widely used in simulations of various fields. However, when considering domains whose extent differs strongly in different spatial directions a finite element simulation becomes computationally very expensive due to the large number of degrees of freedom. An example of such a domain are the cables inside of the magnets of particle accelerators. For translationally invariant domains, this work proposes a quasi-3-D method. Thereby, a 2-D finite element method with a nodal basis in the cross-section is combined with a spectral method with a wavelet basis in the longitudinal direction. Furthermore, a spectral method with a wavelet basis and an adaptive and time-dependent resolution is presented. All methods are verified. As an example the hot-spot propagation due to a quench in Rutherford cables is simulated successfully.

scholarly journals Nonconforming Finite Element Method Applied to the Driven Cavity Problem

2017 ◽  
Vol 21 (4) ◽  
pp. 1012-1038 ◽  
Author(s):  
Roktaek Lim ◽  
Dongwoo Sheen
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Degrees Of Freedom ◽  
Nonconforming Finite Element ◽  
Element Space ◽  
Nonconforming Finite Element Method ◽  
Driven Cavity ◽  
Dirichlet Data ◽  
Minimal Modification ◽  
Element Method

AbstractA cheapest stable nonconforming finite element method is presented for solving the incompressible flow in a square cavity without smoothing the corner singularities. The stable cheapest nonconforming finite element pair based on P1×P0 on rectangularmeshes [29] is employed with a minimal modification of the discontinuous Dirichlet data on the top boundary, where is the finite element space of piecewise constant pressures with the globally one-dimensional checker-board pattern subspace eliminated. The proposed Stokes elements have the least number of degrees of freedom compared to those of known stable Stokes elements. Three accuracy indications for our elements are analyzed and numerically verified. Also, various numerous computational results obtained by using our proposed element show excellent accuracy.

Sign in / Sign up

Export Citation Format

Share Document