Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions

doi:10.1016/j.jcp.2017.12.037

Title: Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions

Abstract

A discrete De Rham complex enables compatible, structure-preserving discretizations for a broad range of partial differential equations problems. Such discretizations can correctly reproduce the physics of interface problems, provided the grid conforms to the interface. However, large deformations, complex geometries, and evolving interfaces makes generation of such grids difficult. We develop and demonstrate two formally equivalent approaches that, for a given background mesh, dynamically construct an interface-conforming discrete De Rham complex. Both approaches start by dividing cut elements into interface-conforming subelements but differ in how they build the finite element basis on these subelements. The first approach discards the existing non-conforming basis of the parent element and replaces it by a dynamic set of degrees of freedom of the same kind. The second approach defines the interface-conforming degrees of freedom on the subelements as superpositions of the basis functions of the parent element. These approaches generalize the Conformal Decomposition Finite Element Method (CDFEM) and the extended finite element method with algebraic constraints (XFEM-AC), respectively, across the De Rham complex.

Authors:

Kramer, Richard M. J. ^[1]; Siefert, Christopher M. ^[1]; Voth, Thomas E. ^[1]; Bochev, Pavel B. ^[1]

Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Publication Date:: 2018-01-02

Research Org.:: Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Sponsoring Org.:: USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1469640

Report Number(s):: SAND2018-9818J
Journal ID: ISSN 0021-9991; 667673

Grant/Contract Number:: AC04-94AL85000

Resource Type:: Journal Article: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 359; Journal Issue: C; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS

Citation Formats


                    Kramer, Richard M. J., Siefert, Christopher M., Voth, Thomas E., and Bochev, Pavel B. Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions.  United States: N. p., 2018. 
        Web.  doi:10.1016/j.jcp.2017.12.037.

Copy to clipboard


                    Kramer, Richard M. J., Siefert, Christopher M., Voth, Thomas E., & Bochev, Pavel B. Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions.  United States.  doi:10.1016/j.jcp.2017.12.037.

Copy to clipboard


                    Kramer, Richard M. J., Siefert, Christopher M., Voth, Thomas E., and Bochev, Pavel B. Tue .  
        "Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions".  United States.  doi:10.1016/j.jcp.2017.12.037.  https://www.osti.gov/servlets/purl/1469640.

Copy to clipboard


                    
@article{osti_1469640,

  title        = {Formulation and computation of dynamic, interface-compatible Whitney complexes in three dimensions},

  author       = {Kramer, Richard M. J. and Siefert, Christopher M. and Voth, Thomas E. and Bochev, Pavel B.},

  abstractNote = {A discrete De Rham complex enables compatible, structure-preserving discretizations for a broad range of partial differential equations problems. Such discretizations can correctly reproduce the physics of interface problems, provided the grid conforms to the interface. However, large deformations, complex geometries, and evolving interfaces makes generation of such grids difficult. We develop and demonstrate two formally equivalent approaches that, for a given background mesh, dynamically construct an interface-conforming discrete De Rham complex. Both approaches start by dividing cut elements into interface-conforming subelements but differ in how they build the finite element basis on these subelements. The first approach discards the existing non-conforming basis of the parent element and replaces it by a dynamic set of degrees of freedom of the same kind. The second approach defines the interface-conforming degrees of freedom on the subelements as superpositions of the basis functions of the parent element. These approaches generalize the Conformal Decomposition Finite Element Method (CDFEM) and the extended finite element method with algebraic constraints (XFEM-AC), respectively, across the De Rham complex.},

  doi          = {10.1016/j.jcp.2017.12.037},

  journal      = {Journal of Computational Physics},

  issn         = {0021-9991},

  number       = C,

  volume       = 359,

  place        = {United States},

  year         = {2018},

  month        = {1}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

DOI: 10.1016/j.jcp.2017.12.037

Other availability

Search WorldCat to find libraries that may hold this journal

Figures / Tables:

Table 1.: Summary of the differential and boundary operators for the model PDEs, the spaces for their weak variational forms (22), the Whitney spaces necessary for the compatible spatial discretization (23) of the weak forms with boundary conditions (20) and the mesh entities containing the degrees of freedom of thesemore »

All figures and tables (12 total)

Save / Share:

Export Metadata

Save to My Library

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

Similar records in OSTI.GOV collections:

Geophysical subsurface imaging and interface identification.

Technical Report Pendley, Kevin ; Bochev, Pavel Blagoveston ; Day, David Minot ; ...

Electromagnetic induction is a classic geophysical exploration method designed for subsurface characterization--in particular, sensing the presence of geologic heterogeneities and fluids such as groundwater and hydrocarbons. Several approaches to the computational problems associated with predicting and interpreting electromagnetic phenomena in and around the earth are addressed herein. Publications resulting from the project include [31]. To obtain accurate and physically meaningful numerical simulations of natural phenomena, computational algorithms should operate in discrete settings that reflect the structure of governing mathematical models. In section 2, the extension of algebraic multigrid methods for the time domain eddy current equations to the frequency domainmore »« less
DOI: 10.2172/875615

Full Text Available
Space-Time Approximation with Sparse Grids

Journal Article Griebel, M ; Oeltz, D ; Vassilevski, P S - SIAM Journal on Scientific Computing

In this article we introduce approximation spaces for parabolic problems which are based on the tensor product construction of a multiscale basis in space and a multiscale basis in time. Proper truncation then leads to so-called space-time sparse grid spaces. For a uniform discretization of the spatial space of dimension d with O(N{sup d}) degrees of freedom, these spaces involve for d > 1 also only O(N{sup d}) degrees of freedom for the discretization of the whole space-time problem. But they provide the same approximation rate as classical space-time Finite Element spaces which need O(N{sup d+1}) degrees of freedoms. Thismore »« less
Full Text Available
A computational model for molten corium spreading and solidification

Journal Article Kucala, Alec ; Rao, Rekha ; Erickson, Lindsay - Computers and Fluids

We report that when the core is breached during a severe nuclear accident, a molten mixture of nuclear fuel, cladding, and structural supports is discharged from the reactor vessel. This molten mixture of ceramic and metal is often referred to as “corium”. Predicting the flow and solidification of corium poses challenges for numerical models due to the presence of large Peclet numbers when convective transport dominates the physics. Here, we utilize a control volume finite-element method (CVEM) discretization to stabilize the advection dominated flow and heat transport. This CVFEM approach is coupled with the conformal decomposition finite-element method (CDFEM), whichmore »« less
DOI: 10.1016/j.compfluid.2018.11.009
An interface tracking model for droplet electrocoalescence.

Technical Report Erickson, Lindsay Crowl

This report describes an Early Career Laboratory Directed Research and Development (LDRD) project to develop an interface tracking model for droplet electrocoalescence. Many fluid-based technologies rely on electrical fields to control the motion of droplets, e.g. microfluidic devices for high-speed droplet sorting, solution separation for chemical detectors, and purification of biodiesel fuel. Precise control over droplets is crucial to these applications. However, electric fields can induce complex and unpredictable fluid dynamics. Recent experiments (Ristenpart et al. 2009) have demonstrated that oppositely charged droplets bounce rather than coalesce in the presence of strong electric fields. A transient aqueous bridge forms betweenmore »« less
DOI: 10.2172/1096257

Full Text Available
Cutting Tetrahedra by Node Identifiers

Technical Report Kramer, Richard Michael Jack

This report briefly outlines an algorithm for dividing a tetrahedron intersected by a planar interface into conforming sub-tetrahedra. The problem of conformal decomposition of tetrahedral meshes arises in enriched finite element methods; in particular, we are concerned with the Conformal Decomposition Finite Element Method (CDFEM) and variants of the eXtended Finite Element Method (XFEM). The algorithm presented is based on the paper How to Subdivide Pyramids, Prisms and Hexahedra into Tetrahedra by Dompierre, Labbe, Vallet, and Camarero (1999), and here is applied and extended to the problem of fully defining and tracking all geometric features of the sub-tetrahedra generated whenmore »« less
DOI: 10.2172/1459088

Full Text Available

Similar Records