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

pressure spikes due to addition/deletion of layers

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 25, 2011, 08:45
Default pressure spikes due to addition/deletion of layers
  #1
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
Hello All
I was running icoDyMFoam on the tutorial case movingConeTopo. I added one line to icoDyMFoam to write out the maximum pressure in the domain at every time step. I found that almost every time there is a cell addition/deletion, there is a pressure spike in the domain. Has anyone else seen these spikes? and has anyone found a solution?
Thanks
Kshitij
kdneroorkar is offline   Reply With Quote

Old   October 26, 2011, 03:48
Default
  #2
Senior Member
 
Karl-Johan Nogenmyr
Join Date: Mar 2009
Location: Linköping
Posts: 249
Rep Power: 12
kalle is on a distinguished road
Hi,

Which version are you using? The implementation of the movingCone(or Body)TopoMesh is altered in later revisions of the 1.6-ext version. In the newer (and parallel) version, a lot of lines are commented out in Foam::movingBodyTopoFvMesh::update() . These lines affects how meshPhi (in fvc::makeAbsolute(phi, U) and fvc::makeRelative(phi, U)) is determined (if I understood correctly). If meshPhi is wrong, I would guess it could result in pressure peaks. Do you have any possibility to repeat your check on a previous, non-parallel version of the movingCone/Body case?

Kalle
kalle is offline   Reply With Quote

Old   October 27, 2011, 09:48
Default
  #3
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
Hi Kalle
Thanks for your post. We updated to the latest version of 1.6-ext to see what difference it made, and we found that the spikes were considerably smaller, but still existed.
if you/ anyone else would like to test this, you can just add the following functionObject to the controlDict and run the movingConeTopo test case with icoDyMFoam.
Thanks a lot for your help
Kshitij


functions{
fieldMinMax
{
type fieldMinMax;
functionObjectLibs ("libfieldFunctionObjects.so");
log yes;
outputControl timeStep;
outputInterval 1;
mode magnitude;
fields
(
U
p
);
}
}
kdneroorkar is offline   Reply With Quote

Old   October 28, 2011, 03:35
Default
  #4
Senior Member
 
Karl-Johan Nogenmyr
Join Date: Mar 2009
Location: Linköping
Posts: 249
Rep Power: 12
kalle is on a distinguished road
Hi,

Thank you for bringing this up. I tried the case, and there are spikes as you say. Pressure is initially fluctuating around ~4e-4 to ~1e-4 as layerAR is working. I tried a few modifications to the movingBodyTopoFvMesh-class without success. Right now I do not have any clue how to get on. To solve this issue, I think one should find out how this adding layer functionality is handled within the finite volume method. This paper may provide some help:

http://web.univ-ubs.fr/limatb/EG2M/D...pers/01178.pdf see especially section 2.3

I also need to fix this for my own work, so I will likely get back to it when I have time.

Regards,
Kalle
kalle is offline   Reply With Quote

Old   October 28, 2011, 08:58
Default
  #5
Member
 
kshitij neroorkar
Join Date: Mar 2009
Location: Michigan, USA
Posts: 32
Rep Power: 8
kdneroorkar is on a distinguished road
Hi Kalle
Thanks for looking into this and thanks for the paper. Please let me know if you find anything.
Kshitij
kdneroorkar is offline   Reply With Quote

Old   November 29, 2014, 13:33
Default
  #6
New Member
 
James
Join Date: Feb 2014
Location: London
Posts: 13
Rep Power: 3
jimteb is on a distinguished road
Hi,

I can see this thread is quite old, but I would like to know if this problem was ever resolved?

I am having the same issue but with pimpleDyMFoam with layer addition/removal and translation (foam-extend 3.1). I have tried to resolve this by mostly messing around with the following settings, but with no real success:

nCorrectors - increase
nNonOrthogonalCorrectors - increase
nOuterCorrectors - increase
relaxation factor for pressure - reduce
use of smoother & preconditioner with solver
minimum layer thickness in dynamicMeshDict - increase
Gauss upwind for div schemes and Guess linear limited for laplacian schemes zeroGradient for all pressures at boundary patches
relTol for all solvers set to 0

Any advice would be greatly appreciated

Thanks
jimteb 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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Pressure Outlet Guage pressure Mohsin FLUENT 34 January 12, 2012 06:45
Pressure Condition sidd CD-adapco 0 April 2, 2007 09:31
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11


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