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

Modified pisoFoam for pitzDaily - solver blowing up !!

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 12, 2016, 13:23
Default Modified pisoFoam for pitzDaily - solver blowing up !!
  #1
Member
 
Suman Chakraborty
Join Date: Sep 2014
Location: Mumbai, India
Posts: 42
Rep Power: 5
suman91 is on a distinguished road
Hi there,

Am having some problem in solving flow over backward facing step in OpenFOAM. I am using LES and Smagorinsky Model. At start it all goes well but after some time utau values at the wall starts increasing and the deltaT of simulation goes on decreasing (till if falls to 10^-12 or so order). I cannot detect what is the problem. I have modified the pisoFoam solver so impose stress boundary condition using the 'fixedGradient' feature of OpenFOAM.

Am attaching my required files here. Please have a look and suggest where am I going wrong. Will be waiting for a reply. Thanks in advance.

solver.tar.gz

TestCase.tar.gz
suman91 is offline   Reply With Quote

Old   October 13, 2016, 02:43
Default
  #2
Member
 
Suman Chakraborty
Join Date: Sep 2014
Location: Mumbai, India
Posts: 42
Rep Power: 5
suman91 is on a distinguished road
Hi there,

I think the problems are coming with the pressure solver settings. At t=0.01 pressure values are about -120 at the inlet and then increases linearly to close to zero at the outlet. But at T=0.05, it falls down to -4000 approx. Obviously in an incompressible flow, the absolute value of pressure does not matter, it is the gradient that matters. But what I am not understanding are such low values. Considering inlet velocity has a magnitude of 10 in the x direction (which is max amongst all 3 directions), a value of -50 to -70 may still be understandable but -4000, am not getting at all. Not understanding why the pressure field is behaving such. Maybe that is triggering the problem.
suman91 is offline   Reply With Quote

Old   October 26, 2016, 09:32
Default
  #3
Member
 
Suman Chakraborty
Join Date: Sep 2014
Location: Mumbai, India
Posts: 42
Rep Power: 5
suman91 is on a distinguished road
Hi,

Back again with an update on this same issue. So to start with, setting the gradient at the wall was not so much of a good idea it seems. Somehow the wall normal velocity was being calculated to be a non zero value at the wall. This might be one of the reasons for the solver to blow up in that manner. The modified code I am using now, uses fixedValue (0 0 0) as a boundary condition at the wall and in the modified solver, I calculate the velocity at the wall using the gradient and set the velocity. The solver now runs properly and does not blow up.

The only issue am facing with this is, the reynolds stress when calculated comes out to be very less; flow is not becoming turbulent only. Am wondering if it's an issue with the gradient schemes being used or not. If anyone has any views on this, it will be really helpful.

So the question is, what is the best gradient scheme to be used for solving incompressible turbulent flow (using LES) in openFOAM?
suman91 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
New Solver for a modified pUCoupledFoam noiseMaker OpenFOAM 0 August 2, 2016 06:14
How to add a new scheme to a solver and run a case with the modified solver? SINA.KHAYYAM OpenFOAM Programming & Development 0 January 12, 2016 12:41
Star cd es-ice solver error ernarasimman STAR-CD 2 September 12, 2014 00:01
Error with modified interFoam solver robertmaier9 OpenFOAM 3 September 30, 2011 06:26
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25


All times are GMT -4. The time now is 15:21.