I tried to create define Dirichlet BC as a spatially varying vector. To do so, I defined a SpatialCoordinate and wrapped the function is ufl.as_vector. The MWE from dolfinx import mesh domain = mesh.create_unit_square(MPI.COMM_WORLD, 8, 8, mesh.CellType.quadrilateral) V = functionspace(domain, ("L…

Interpolate the u_D into the function space of the bc. This is for instance covered in: Diffusion of a Gaussian function — FEniCSx tutorial

I created a function and interpolated it into the function space, but I got an error. This is a mwe domain = mesh.create_unit_square(MPI.COMM_WORLD, 8, 8, mesh.CellType.quadrilateral) V = functionspace(domain, ("Lagrange", 1, (domain.geometry.dim,))) #vector function space x = ufl.SpatialCoordinat…

Your function rect is wrong. See dolfinx/python/demo/demo_navier-stokes.py at main · FEniCS/dolfinx · GitHub for an example on how to define it correctly. (Hint: do not get confused by the newlines, that is a vector expression with two components, look for commas to determine where the first compon…

Thank you. I also found another problem: instead of [image] oszachter: bc_up = fem.dirichletbc(u_D, dofs_up, V) I should have written bc_up = fem.dirichletbc(u_D, dofs_up) (without the vector space)

Spatially varying BC

Errors

oszachter February 14, 2024, 3:27pm 3

Thank you!

I also found this helpful

Topic		Replies	Views
How to set a dirichlet boundary condition, which depends on a Function? variational formulation	2	484	July 1, 2022
Boundary condition value must have a dtype attribute General boundary-conditions , dolfinx	3	215	January 11, 2024
Time-dependent BC dolfinx boundary-conditions	5	213	December 19, 2023
How to impose spatially varying dirichletbc boundary conditions for one direction in a linear elasticity problem? dolfinx boundary-conditions , dolfinx	6	311	April 5, 2024
Problems about setting boundary conditions with two different ways dolfinx	7	247	September 28, 2023

Spatially varying BC

Related topics