Recently we have developed the optimal local truncation error method (OLTEM) for scalar PDEs on irregular domains and unfitted Cartesian meshes. Here, OLTEM is extended to a much more general case of a system of PDEs for the 2‐D time‐independent elasticity equations on irregular domains. Compact 9‐point uniform and nonuniform stencils (with the computational costs of linear finite elements) are used with OLTEM. The stencil coefficients are assumed to be unknown and are calculated by the minimization of the local truncation error. It is shown that the second order of accuracy is the maximum possible accuracy for 9‐point stencils independent of the numerical technique used for their derivations. The special treatment of the Neumann boundary conditions has been developed that does not increase the size of the stencils. The numerical examples are in agreement with the theoretical findings. They also show that due to the minimization of the local truncation error, OLTEM is much more accurate than linear finite elements and than quadratic finite elements (up to engineering accuracy of 0.1%–1%) at the same numbers of degrees of freedom. Due to the computational efficiency and trivial unfitted Cartesian meshes that are independent of irregular domains, the proposed technique with no remeshing for the shape change of irregular domains will be effective for many engineering applications.
- Award ID(s):
- 1935452
- NSF-PAR ID:
- 10346121
- Date Published:
- Journal Name:
- International Journal of Numerical Methods for Heat & Fluid Flow
- Volume:
- 32
- Issue:
- 8
- ISSN:
- 0961-5539
- Page Range / eLocation ID:
- 2719 to 2749
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Abstract The optimal local truncation error method (OLTEM) has been developed for 2‐D time‐dependent elasticity problems on irregular domains and trivial unfitted Cartesian meshes. Compact nine‐point uniform and nonuniform stencils (similar to those for linear finite elements on uniform meshes) are used with OLTEM. The stencil coefficients are assumed to be unknown and are calculated by the minimization of the local truncation error. It is shown that the second order of accuracy is the maximum possible accuracy for nine‐point stencils independent of the numerical technique used for their derivations. The special treatment of the Neumann boundary conditions has been developed that does not increase the size of the stencils. The cases of the nondiagonal and diagonal mass matrices are considered for OLTEM. The results of numerical examples are in agreement with the theoretical findings. They also show that due to the minimization of the local truncation error, OLTEM with the nondiagonal mass matrix is much more accurate than linear finite elements and than quadratic and cubic finite elements (up to the engineering accuracy of 0.1%–1%) at the same numbers of degrees of freedom. The proposed numerical technique can be efficiently used for many engineering applications including geomechanics.
-
The primal variational formulation of the fourth-order Cahn-Hilliard equation requires C1-continuous finite element discretizations, e.g., in the context of isogeometric analysis. In this paper, we explore the variational imposition of essential boundary conditions that arise from the thermodynamic derivation of the Cahn-Hilliard equation in primal variables. Our formulation is based on the symmetric variant of Nitsche's method, does not introduce additional degrees of freedom and is shown to be variationally consistent. In contrast to strong enforcement, the new boundary condition formulation can be naturally applied to any mapped isogeometric parametrization of any polynomial degree. In addition, it preserves full accuracy, including higher-order rates of convergence, which we illustrate for boundary-fitted discretizations of several benchmark tests in one, two and three dimensions. Unfitted Cartesian B-spline meshes constitute an effective alternative to boundary-fitted isogeometric parametrizations for constructing C1-continuous discretizations, in particular for complex geometries. We combine our variational boundary condition formulation with unfitted Cartesian B-spline meshes and the finite cell method to simulate chemical phase segregation in a composite electrode. This example, involving coupling of chemical fields with mechanical stresses on complex domains and coupling of different materials across complex interfaces, demonstrates the flexibility of variational boundary conditions in the context of higher-order unfitted isogeometric discretizations.more » « less
-
We present three new semi-Lagrangian methods based on radial basis function (RBF) interpolation for numerically simulating transport on a sphere. The methods are mesh-free and are formulated entirely in Cartesian coordinates, thus avoiding any irregular clustering of nodes at artificial boundaries on the sphere and naturally bypassing any apparent artificial singularities associated with surface-based coordinate systems. For problems involving tracer transport in a given velocity field, the semi-Lagrangian framework allows these new methods to avoid the use of any stabilization terms (such as hyperviscosity) during time-integration, thus reducing the number of parameters that have to be tuned. The three new methods are based on interpolation using 1) global RBFs, 2) local RBF stencils, and 3) RBF partition of unity. For the latter two of these methods, we find that it is crucial to include some low degree spherical harmonics in the interpolants. Standard test cases consisting of solid body rotation and deformational flow are used to compare and contrast the methods in terms of their accuracy, efficiency, conservation properties, and dissipation/dispersion errors. For global RBFs, spectral spatial convergence is observed for smooth solutions on quasi-uniform nodes, while high-order accuracy is observed for the local RBF stencil and partition of unity approaches.more » « less
-
Immersed finite element methods are designed to solve interface problems on interface- unfitted meshes. However, most of the study, especially analysis, is mainly limited to the two-dimension case. In this paper, we provide an a priori analysis for the trilinear immersed finite element method to solve three-dimensional elliptic interface problems on Cartesian grids consisting of cuboids. We establish the trace and inverse inequalities for trilinear IFE functions for interface elements with arbitrary interface-cutting configuration. Optimal a priori error estimates are rigorously proved in both energy and L2 norms, with the constant in the error bound independent of the interface location and its dependence on coefficient contrast explicitly specified. Numerical examples are provided not only to verify our theoretical results but also to demonstrate the applicability of this IFE method in tackling some real-world 3D interface models.more » « less