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
V0.6.0 container subdomain multiple variable General	1	172	September 26, 2023
Non-matching coordinate ordering for CG1 functions defined on submeshes dolfinx mesh , dolfinx	2	43	November 22, 2024
Subdomain in dolfinx dolfinx	1	1174	July 4, 2022
Mixed-dimensional branch - Error: Unable to set given (local) rows to identity matrix Legacy DOLFIN	2	129	May 18, 2024
Submesh from parallel mesh dolfinx	9	83	September 5, 2024

Coupling of two subdomains with new mixed-dimensional assembly

Related topics