Юлианский язык FEniCS периодическое граничное условие - PullRequest
Новая
       31

Юлианский язык FEniCS периодическое граничное условие

0 голосов
/ 29 октября 2018

Я хотел бы применить периодическое граничное условие в FEniCS для языка Julia, но все примеры, которые я нашел, находятся на C ++ или Python. Как создать периодическое граничное условие, используя Юлию? Это кажется трудным, потому что у Юлии нет занятий. вот минимальный пример: 

using FEniCS

using PyCall

length=2.2

height=0.41

channel = Rectangle(Point([0.0, 0.0]), Point([length, height]))

domain = channel

mesh = generate_mesh(domain, 64)

# insert function here for PeriodicBoundarycondition

Q = FunctionSpace(mesh, "P", 1,constrained_domain=#the function that i am looking for)

1 Ответ

0 голосов
/ 29 октября 2018

Я посмотрел код julia, который составляет FEniCS, пример периодического граничного условия на странице fenics, и некоторые из моих старых кодов python для fenics, и это вдохновило меня на написание этого: 

    using FEniCS
using PyCall
@pyimport fenics
py"""
from dolfin import *
from mshr import *
length=2.2

height=0.41
channel = Rectangle(Point([0.0, 0.0]), Point([length, height]))

domain = channel

mesh = generate_mesh(domain, 64)
subdomains = MeshFunction("size_t", mesh, 1)
subdomains.set_all(0)
class Wall(SubDomain):
    def inside(self,x,on_boundary):
        return (near(x[1],height) or near(x[1],height)) and on_boundary
wall=Wall()
class PeriodicBoundary(SubDomain):

    # Left boundary is "target domain" G
    def inside(self, x, on_boundary):
        return bool(x[0] < DOLFIN_EPS and x[0] > -DOLFIN_EPS and on_boundary)

    # Map right boundary (H) to left boundary (G)
    def map(self, x, y):
        y[0] = x[0] - length
        y[1] = x[1]
pbc=PeriodicBoundary()
"""
Q=FunctionSpace(fenics.VectorFunctionSpace(py"mesh", "P", 1,constrained_domain=py"pbc"))

решение не является оптимальным, поскольку оно просто делает все это в python, но я думаю, что мне придется с этим смириться.

...