Coupling of two subdomains with new mixed-dimensional assembly

dokken · November 20, 2024, 7:28pm

I would consider

and the demo I made in

gist.github.com

https://gist.github.com/jorgensd/8a5c32f491195e838f5863ca88b27bce

dg.py

from mpi4py import MPI
import dolfinx
import dolfinx.fem.petsc
import dolfinx.nls.petsc
import numpy
import argparse
from ufl import (
    Circumradius,
    FacetNormal,
    SpatialCoordinate,

This file has been truncated. show original

problem.py

# SPDX-License-Identifier: MIT

from mpi4py import MPI
import dolfinx
import dolfinx.fem.petsc
import ufl
import numpy as np
from petsc4py import PETSc

This file has been truncated. show original

which couples two poisson problems discretised with different finite element spaces over a common boundary.

I have a similar implementation for RTC at

gist.github.com

https://gist.github.com/jorgensd/befb28b02de55d254ffab6adcbc3261c

rtc.py

import dolfinx.fem.petsc
from mpi4py import MPI

import dolfinx
import numpy as np
import ufl
from petsc4py import PETSc

import numpy.typing as npt

This file has been truncated. show original

rtc_dg.py

from mpi4py import MPI

import dolfinx
import ufl
import numpy as np

mesh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD, 10, 10, cell_type=dolfinx.mesh.CellType.triangle)


V = dolfinx.fem.functionspace(mesh, ("DG", 1))

This file has been truncated. show original

Topic		Replies	Views
Discontinuity at interface using mixed domains variational formulation dolfinx	12	447	September 13, 2024
V0.6.0 container subdomain multiple variable General	1	172	September 26, 2023
Integrating over interface between two functionspaces dolfinx	11	207	December 25, 2024
Solve poisson equation with three submeshes dolfinx	12	309	May 8, 2025
Mixed finite element for two fields one defined on whole mesh and the other defined on submesh General	11	234	July 2, 2024

Coupling of two subdomains with new mixed-dimensional assembly

Related topics