Does DOLFINX_MPC have the support for Nedelec element

WjjSteve · January 28, 2025, 8:33am

Hello, I’m wondering that can I use DOLFINX_MPC to impose periodic BC for Nedelec elements? I find this Issue: Github repo, but don’t find any update in it. So I guess that this is not possible so far, is that true?

Thanks in advance.

sbhasan · January 28, 2025, 5:34pm

It is possible with some kinks. Please refer to the following:

I was able to use the method for 2D geometries and did not encounter the anomaly indicated in the thread. But that somehow meant I did not stumble upon the corner case?

Also look down for an alternate approach by @dokken that makes the mesh itself periodic.

dokken · January 29, 2025, 12:28pm

The implementation with the periodic mesh referenced in the above post has been extended to parallel in:

gist.github.com

https://gist.github.com/jorgensd/e0e610fd17f73337178bc12788fa17ba#file-script-py

README.md

# TODO

1.  Create interface that reads in the original global index from MSH directly
```python
from mpi4py import MPI
import dolfinx.io
import numpy.typing as npt
import gmsh

gmsh.initialize()

This file has been truncated. show original

create_mesh.py

import gmsh
import numpy as np
import argparse
from pathlib import Path

parser = argparse.ArgumentParser(description="Create a mesh for a box in a channel")
parser.add_argument(
    "--res", type=float, default=0.05, help="Resolution of the mesh (at inlet)"
)
parser.add_argument(

This file has been truncated. show original

mwe_failing.py

from dolfinx import fem
from dolfinx.fem.petsc import LinearProblem
from dolfinx.io.gmshio import model_to_mesh
from dolfinx.io import VTXWriter
import ufl
import gmsh
import dolfinx
from petsc4py import PETSc
from basix.ufl import element
from mpi4py import MPI

This file has been truncated. show original

There are more than three files. show original

It will likely be incorporated in scifem in the foreseeable future.

Topic		Replies	Views
Periodic conditions for vector valued spaces dolfinx boundary-conditions	12	265	December 22, 2024
Periodic boundary constraints for vector-valued functionspaces in dolfinx_mpc, when V.sub(i) is insufficient dolfinx	0	230	August 26, 2023
Periodic boundary condition with non-periodic mesh using dolfinx_mpc General	1	463	July 21, 2023
Periodic boundary condition with DG and MPC General	6	382	December 6, 2023
Periodic boundary conditions Homogenization dolfinx boundary-conditions , dolfinx_mpc	4	127	November 29, 2024

Does DOLFINX_MPC have the support for Nedelec element

Related topics