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

UEqn.A() and UEqn.H()

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 25, 2010, 05:31
Default UEqn.A() and UEqn.H()
  #1
Senior Member
 
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 23
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;
Code:
UEqn =
fvm::ddt(U)
+ fvm::div(phib, U)
should I not be able to reconstruct it as follows:

Code:
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

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
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 05: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 01:53.