CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Add a custom body force to the momentum equation, f=f(x,y,z)

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 5, 2012, 05:17
Default Add a custom body force to the momentum equation, f=f(x,y,z)
  #1
Member
 
Jeong Kim
Join Date: Feb 2010
Posts: 39
Rep Power: 6
enoch is on a distinguished road
I wanted to add a force to the alpha-phase moment equation.
This is a funciton of x, y and z location, which force is not dependent on time. For an exmaple, customF = customF(x,y,z) where i,j and k are indexing counters. Below is the folmula expression that I'd like to use.

if(x(i)>x0 or x(i)<xl )
customF(i) = a*x(i-1)*y(j)*z(k)
+ a*x(i)*y(j)*z(k)
+ a*x(i+1)*y(j)*z(k);
else
customF(i) = 0;

The vector customF is composed of three components (x,y,z). How can I write this force in the OF?
enoch is offline   Reply With Quote

Old   March 5, 2012, 15:44
Default
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,889
Rep Power: 38
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by enoch View Post
I wanted to add a force to the alpha-phase moment equation.
This is a funciton of x, y and z location, which force is not dependent on time. For an exmaple, customF = customF(x,y,z) where i,j and k are indexing counters. Below is the folmula expression that I'd like to use.

if(x(i)>x0 or x(i)<xl )
customF(i) = a*x(i-1)*y(j)*z(k)
+ a*x(i)*y(j)*z(k)
+ a*x(i+1)*y(j)*z(k);
else
customF(i) = 0;

The vector customF is composed of three components (x,y,z). How can I write this force in the OF?
You might want to have a look at this (the thread discusses various approaches to a similar problem, but my favourite is - for obvious reasons - expressionSource from swak4Foam):
transient source term
gschaider is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
how to add source-term in momentum equation for interFoam? anon_g OpenFOAM 9 October 18, 2011 13:47
Force can not converge colopolo CFX 13 October 4, 2011 23:03
Question on the discretization of momentum equation in icoFoam MPJ OpenFOAM 3 October 4, 2011 10:44
momentum equation (Normal to wall_ kk Main CFD Forum 0 July 3, 2006 00:10
Force on cylinders through momentum balance Mauro Oliveira Main CFD Forum 0 April 14, 2000 09:27


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