DOLFIN

Is it possible to strongly apply no-normal flow on boundaries that are not aligned with an axis?

Asked by Patrick Farrell

Is it possible to strongly impose the condition

u . n = 0

in a general way in dolfin? It's easy if the boundary is aligned with the axes, but is it currently possible if the boundary is not aligned?

Question information

Language:
English Edit question
Status:
Answered
For:
DOLFIN Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Johan Jansson (jjan) said : 		#1

Hi Patrick,

In the Unicorn application component you have the class SlipBC which encapsulates a strong slip boundary condition (u . n = 0):

http://launchpad.net/unicorn

Usage example:

http://bazaar.launchpad.net/~unicorn-core/unicorn/hpc/view/head:/ucsolver/icns/cylinder/main.cpp

Best,
  Johan

Revision history for this message
Anders Logg (logg) said : 		#2

On Sat, Apr 28, 2012 at 03:25:51PM -0000, Patrick Farrell wrote:
> New question #195121 on DOLFIN:
> https://answers.launchpad.net/dolfin/+question/195121
>
> Is it possible to strongly impose the condition
>
> u . n = 0
>
> in a general way in dolfin? It's easy if the boundary is aligned with the axes, but is it currently possible if the boundary is not aligned?

No, not unless you either use BDM (or RT) elements, or implement the
bc weakly.

A new class NoslipBC, or a more general class NormalBC (of which
NoslipBC could be a subclass) would be a welcome contribution.

--
Anders

Can you help with this problem?

Provide an answer of your own, or ask Patrick Farrell for more information if necessary.

To post a message you must log in.