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

Create Matrix for pressure correction

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

LinkBack Thread Tools Display Modes
Old   September 20, 2011, 15:48
Default Create Matrix for pressure correction
New Member
Join Date: Dec 2009
Location: Austria
Posts: 13
Rep Power: 10
nep101 is on a distinguished road
Hey Everybody,

i try to create a Matrix UEqn for the correction Loop.
The Equatin is:
U and k are vectors and T and Ra are scalars.

Is there a way to create a Matrix like
fvVectorMatrix UEqn
- Ra*T*k
but this doesnt work

Thanks for Help!
nep101 is offline   Reply With Quote

Old   November 7, 2011, 00:59
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 358
Rep Power: 17
mprinkey will become famous soon enough
You need to use fvm::Sp(U) (or fvm::S for the linear source term, so something like:

fvVectorMatrix UEqn
fvm::Sp(U) - Ra*T*k

should work. Note that this is just a diagonal matrix.
mprinkey 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
how to set periodic boundary conditions Ganesh FLUENT 13 January 22, 2014 05:11
Include list of points Hikachu OpenFOAM Native Meshers: blockMesh 0 June 20, 2011 09:03
Actuator disk model audrich FLUENT 0 September 21, 2009 07:06
Where's the singularity/mesh flaw? audrich FLUENT 3 August 4, 2009 01:07
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55

All times are GMT -4. The time now is 18:52.