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

Memory leak due to source term???

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 5, 2011, 10:38
Default Memory leak due to source term???
  #1
RoE
New Member
 
Roland Engberg
Join Date: Jan 2011
Posts: 14
Rep Power: 6
RoE is on a distinguished road
Hi!

Simulating a rising bubble I've implemented a moving reference frame similar to the approach outlined in Henrik Rusche's PhD thesis.
To take the moving reference frame into account, I've added a "source" term to the UEqn:

rhoaF = rho*aF; // density*acceleration (vector)

fvVectorMatrix UEqn

(
fvm::ddt(rho, U)
+ rhoaF //---------------------------acceleration term!
+ fvm::div(rhoPhi, U)
- fvm::laplacian(eta, U)
- (fvc::grad(U) & fvc::grad(eta))
);


rhoaF is a volVectorField. I've also tried adding rho*aF directly to the UEqn, which seems to make no difference.

At the begining of a simulation a bool is read from the solution dictionary to decide whether the reference frame has to be moved. If the reference frame is static, aF is set to (0 0 0) so UEqn remains unchanged.

Now if I simulate with a static frame, everything is fine. When I change to the moving reference frame mode, the results are reasonable, but working memory rises constantly. I've run my code with valgrind but couldn't find any errors. Moreover, there are no additional variables declared when I run the code with moving reference frame.

At the moment I'm pretty perplexed... Any advice would be welcome!

Thanks!
Roland
RoE is offline   Reply With Quote

Old   January 10, 2012, 01:37
Default OpenFOAM Memory Leak?
  #2
New Member
 
Jie Zhang
Join Date: Sep 2010
Location: FL, U.S.
Posts: 20
Rep Power: 6
jiez is on a distinguished road
I'm simulating a LES case using pisoFoam with a source term added. The memory requirement increased constantly from 2G to 8G after I ran 4E4 steps. I don't know what is the reason. May I have some advice from you?
Appreciate your help in advance.

Jie
jiez 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
transient source term strakakl OpenFOAM 38 November 19, 2013 02:18
Problem of SOURCE term gradient in UDS wind Fluent UDF and Scheme Programming 5 June 21, 2013 05:39
How to write source term into scalar Fiel JimKnopf OpenFOAM Programming & Development 0 March 23, 2011 06:59
Compiling gmshFoam with OpenFOAM-1.5 BlGene Open Source Meshers: Gmsh, Netgen, CGNS, ... 10 August 6, 2009 04:26
Adding a momentum source term segersson OpenFOAM Running, Solving & CFD 5 March 3, 2006 00:06


All times are GMT -4. The time now is 04:31.