CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Anisotropic, implicit momentum source term

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

LinkBack Thread Tools Display Modes
Old   June 11, 2013, 06:58
Default Anisotropic, implicit momentum source term
New Member
Chris Meyer
Join Date: Mar 2009
Location: Cape Town, Western Province, South Africa
Posts: 5
Rep Power: 10
chris_j_meyer is on a distinguished road
Hi All

Will not bore you with the details, but I am adding a momentum source term that is anisotropic. To do this explicitly is quite easy as a volVectorField with the sources can just be added to the U-equation. For an implicit formulation the following can be used and added to the U-equation:

fvm::Sp(a, U)

However, "a" can at best be a volScalarField and hence the source will be locally isotropic. Is there any way in which I can add different implicit source terms to the momentum equations in the three coordinate directions?

chris_j_meyer is offline   Reply With Quote

Old   June 11, 2013, 10:56
Senior Member
Join Date: Dec 2011
Location: Leuven, Belgium
Posts: 297
Rep Power: 15
Lieven will become famous soon enough
Hi Chris,

Have a look at the implementation of the porous media source term. I think this is what you are looking for.


Lieven is offline   Reply With Quote

Old   July 18, 2013, 18:42
New Member
Join Date: Jul 2009
Posts: 16
Rep Power: 10
ski is on a distinguished road
I'd like to do this too..

Issue is that it appears to be using Jacobi method to solve, which tends to give poor convergence rates. Anyone have any experience on this?
ski is offline   Reply With Quote

Old   July 25, 2016, 16:11
Question Did anyone find a way to program an implicit anisotropic source???
New Member
Mateus Palharini Schwalbert
Join Date: May 2011
Posts: 21
Rep Power: 8
mateusps is on a distinguished road
Hello y'all,

I have the same problem: how to program an anisotropic implicit source in the momentum equation?

fvm::Sp simply does not accept a tensor as coefficient (fvm::Sp(tensor,vector) does not seem to be implemented).

Did anyone find a way?

mateusps is offline   Reply With Quote


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
swak4foam building problem GGerber OpenFOAM Installation 54 April 24, 2015 16:02
build problem swak4Foam OF 2.2.0 mcathela OpenFOAM Installation 14 April 23, 2013 13:59
How does fluent handles the momentum source term? QBeast FLUENT 0 April 22, 2013 14:12
pisoFoam compiling error with OF 1.7.1 on MAC OSX Greg Givogue OpenFOAM Programming & Development 3 March 4, 2011 18:18
Momentum source UDF Matthew Brannock FLUENT 5 May 3, 2001 21:18

All times are GMT -4. The time now is 03:16.