Hello,
Has anyone solved 3D Taylor-Green Vortex problem with a periodic boundary condition using FeniCSx?
Thanks!
You can consider the code at:
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()show original
anti_periodic_rt.py
from script import create_periodic_mesh, transfer_meshtags_to_periodic_mesh
from mpi4py import MPI
import dolfinx.fem.petsc
import dolfinx.nls.petsc
import dolfinx
import numpy as np
import ufl
from petsc4py import PETSc
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(show original
There are more than three files. show original