scholarly journals A Numerical Method to Improve the Representativeness of Real Microstructure Cut-Outs Applied in Finite Element Simulations

Crystals ◽  
2021 ◽  
Vol 11 (4) ◽  
pp. 382
Author(s):  
Yanling Schneider ◽  
Werner Wasserbäch ◽  
Siegfried Schmauder ◽  
Zhangjian Zhou ◽  
Reiner Zielke ◽  
...  

To improve the representativeness of a real microstructural cut-out for modeling purposes, a numerical method named as “boundary pixel color alteration (BPCA)” is presented to modify measured 2D microstructure cut-outs. Its physical background is related to the phase growth. For the application, the precondition is that the representativeness of the microstructure is already satisfied to a certain extent. This method resolves the problem that the phase composition of a small cut-out can have a large discrepancy to the real one. The main idea is to change the pixel color among neighboring pixels belonging to different phases. Our process simultaneously maintains most of the characteristics of the original morphology and is applicable for nearly all kinds of multi-phase or polycrystalline metallic alloys, as well. From our axisymmetric finite element (FE) simulations (ABAQUS ) applied with 2D real microstructures, it shows that the volume ratios of microstructural phases, as a function of the structure position to the symmetric axis, converge to phase area ratios in the 2D cut-out, even though the axisymmetric element volume is position dependent. A mathematical proof provides the reason for the aforementioned convergence. As examples to achieve real compositions and to numerically prove the aforementioned convergence, four different materials including multiphase polycrystals are implemented. An improvement of the predicted FE result is presented for the application of a modified microstructure (with a higher representativeness) compared to the original one.

2017 ◽  
Vol 14 (01) ◽  
pp. 1750003 ◽  
Author(s):  
Yuching Wu ◽  
Jianzhuang Xiao

In this study, a multi-scale finite element method was proposed to solve two linear scale-coupling stochastic elliptic PDE problems, a tightly stretched wire and flow through porous media. At microscopic level, the main idea was to form coarse-scale equations with a prescribed analytic form that may differ from the underlying fine-scale equations. The relevant stochastic homogenization theory was proposed to model the effective global material coefficient matrix. At the macroscopic level, the Karhunen–Loeve decomposition was coupled with a Polynomial Chaos expansion in conjunction with a Galerkin projection to achieve an efficient implementation of the randomness into the solution procedure. Various stochastic methods were used to plug the microscopic cell to the global system. Strategy and relevant algorithms were developed to boost computational efficiency and to break the curse of dimension. The results of numerical examples were shown consistent with ones from literature. It indicates that the proposed numerical method can act as a paradigm for general stochastic partial differential equations involving multi-scale stochastic data. After some modification, the proposed numerical method could be extended to diverse scientific disciplines such as geophysics, material science, biological systems, chemical physics, oceanography, and astrophysics, etc.


2020 ◽  
Vol 85 (2) ◽  
Author(s):  
R. Abgrall ◽  
J. Nordström ◽  
P. Öffner ◽  
S. Tokareva

AbstractIn the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.


Sign in / Sign up

Export Citation Format

Share Document