CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

Boundary Condition: Antisymmetric periodicity

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 4, 2011, 09:16
Default Boundary Condition: Antisymmetric periodicity
  #1
New Member
 
Jacob Jensen
Join Date: Mar 2011
Posts: 1
Rep Power: 0
Jacob Jensen is on a distinguished road
Has anyone come across boundary conditions in OpenFoam that impose anti-symmetric periodicity …

An anti-symmetric periodic boundary condition is just like periodic boundary conditions only swapped or mirrored … i.e.

u(x, y)=u(x+L, ymax-y)
v(x, y)=-v(x+L, ymax-y)
p(x, y)=p(x+L, ymax-y)

rather then (pure periodicity):

u(x, y)=u(x+L, y)
v(x, y)=v(x+L, y)
p(x, y)=p(x+L, y)

here u =velocity in x, v = velocity in y, p=pressure, x and y coordinates where y[0, ymax].

The boundary condition is described in “ON THE IMPLEMENTATION OF SYMMETRIC AND ANTISYMMETRIC PERIODIC BOUNDARY CONDITIONS
FOR INCOMPRESSIBLE FLOW” by GUUS SEGAL, KEES VUIK AND KEES KASSELS.

Last edited by Jacob Jensen; March 7, 2011 at 08:40.
Jacob Jensen is offline   Reply With Quote

Old   August 17, 2014, 05:35
Default
  #2
New Member
 
Oleg Rogozin
Join Date: Aug 2014
Location: Moscow, Russia
Posts: 1
Rep Power: 0
olegrog is on a distinguished road
Cyclic boundary conditions can be used in this case.
olegrog is offline   Reply With Quote

Old   December 1, 2014, 12:28
Default
  #3
New Member
 
Fredo
Join Date: Apr 2014
Posts: 3
Rep Power: 11
Fredo is on a distinguished road
Quote:
Originally Posted by olegrog View Post
Cyclic boundary conditions can be used in this case.
Hello,

yes they can, but then you need a cyclic geometry.

Let me explain my case :
I have a cubic domain which describe a porous media (liquid and solid). Only the liquid phase is meshed. Since the media is pretty random, the faces of my cube do not have the same mesh.
To apply a cyclic boundary conditions (because that's what I want), I mirror the mesh in the x, y and finally z directions.
I have then a cube 8 times bigger than my initial cube with all the same faces.
I can then use cyclic boundary.

But, I compute 7 times something I am not really interessted in.
If I could apply antisymmetric periodicity, that would save a lot a time.

Do you have any suggestions ?

Fred.
Fredo is offline   Reply With Quote

Old   December 1, 2014, 16:32
Default
  #4
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
Hi Fredo,

What I do in this case: I had a thin layer of cells surrounding my domain. Like that I can impose cyclic boundary conditions.

The two main problems when mirroring the domain: 1) it significantly increases the number of cells, 2) it biases the anisotropy (primordial for permeability estimation!).
Cyp is offline   Reply With Quote

Old   December 2, 2014, 03:29
Default
  #5
New Member
 
Fredo
Join Date: Apr 2014
Posts: 3
Rep Power: 11
Fredo is on a distinguished road
Hi Cyp,

thank you for this quick reply. Your solution seems very relevant ! Did you quantify the errors induced by this thin layer. My media is not very percolating so I fear that all the fluid goes through this thin layer.

I will try and keep you updated.

Fred.
Fredo is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Domain Imbalance HMR CFX 5 October 10, 2016 05:57
How exactly the "pressure outlet" bdry condition compute properties on the boundary? yating9901 FLUENT 3 June 28, 2010 12:26
Axis Boundary Condition..what is it? CFDtoy FLUENT 6 February 13, 2007 05:51
How to set boundary condition in Fluent for the fo Peiyong FLUENT 1 November 10, 2006 11:44
How to resolve boundary condition problem? sam FLUENT 2 July 20, 2003 02:19


All times are GMT -4. The time now is 00:43.