WIP: Mixed Boundary Conditions

[1]:

import dune
import dune.grid
from dune.fem.space import lagrange

[2]:

from ufl import SpatialCoordinate, as_vector

gridView = dune.grid.structuredGrid([-1, -1], [1, 1], [100, 100])
space = lagrange(gridView, dimRange=1)
x = SpatialCoordinate(dune.ufl.cell(2))

[3]:

import ddfem
from ddfem.geometry import Ball

b1 = Ball(0.5, [0.2, 0])
b2 = Ball(0.5, [-0.2, 0])
omega = b1 | b2
omega.epsilon = 1e-1

[4]:



class Model:
    boundary = {
        b1: ddfem.boundary.BndValue(as_vector([0])),
        b2: ddfem.boundary.BndValue(as_vector([1])),
    }

[5]:

uh = space.function(name="uh")

ddfem.plot.plotBndExt(omega, Model, uh, [10, 10], gridLines=None)
_images/boundary_nb_5_0.jpg

This page was generated from the notebook boundary_nb.ipynb and is part of the tutorial for the ddfem package using the dune-fem python bindings DOI