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

UEqn.A() and UEqn.H()

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

LinkBack Thread Tools Display Modes
Old   May 25, 2010, 04:31
Default UEqn.A() and UEqn.H()
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 16
l_r_mcglashan will become famous soon enough
Dear all,

I've been messing about with the momentum equations in bubbleFoam to try and get an implementation of the product elimination algorithm so that the drag forces are solved for fully implicitly.

If I have an fvVectorMatrix, UEqn;
UEqn =
+ fvm::div(phib, U)
should I not be able to reconstruct it as follows:

UEqnReconstructed = fvm::Sp(UEqn.A(),U) - UEqn.H();
If I do this UEqnReconstructed does not appear to equal UEqn, as the simulation results on the bubbleColumn tutorial are different. Is the above correct?
Laurence R. McGlashan :: Website
l_r_mcglashan 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
rasInterFoam Channel Flowlooking for working case nico765 OpenFOAM Running, Solving & CFD 4 March 24, 2009 01:49
A stupid question luckyluke OpenFOAM Running, Solving & CFD 14 August 13, 2007 04:25
SIMPLE for viscous NavierStokes almanstoetter OpenFOAM Running, Solving & CFD 2 March 10, 2006 00:23

All times are GMT -4. The time now is 17:38.