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

Pressure is fluctuating

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 29, 2012, 07:30
Default Pressure is fluctuating
  #1
New Member
 
Sargentini
Join Date: May 2012
Location: Brazil
Posts: 5
Rep Power: 5
Sargentini is on a distinguished road
Dear

I'm trying to simulate a laminar flow using Herschel Bulkley and "non-Newtonian solver."
I adopted nu0 = 1E6, tau0 = 1, k = 0.01 and n = 1 to verify the result gives the analytical value of a Bingham material.
The problem is that pressure is fluctuating too much, although the velocity gain a suitable profile.
Someone please have an idea what might be happening and how to fix this problem.

grateful
Sargentini
Sargentini is offline   Reply With Quote

Old   May 29, 2012, 07:34
Default
  #2
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Maybe you can give some more relevant information, such as discretization and convergence criteria used? Also, which solver are you using?
Bernhard is offline   Reply With Quote

Old   May 29, 2012, 08:08
Default Pressure is fluctuating
  #3
New Member
 
Sargentini
Join Date: May 2012
Location: Brazil
Posts: 5
Rep Power: 5
Sargentini is on a distinguished road
I'm using nonNewtonianIcoFoam.
Follow the settings that have adopted:

fvSolution:
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}
}
PISO
{
nCorrectors 4;
nNonOrthogonalCorrectors 2;
--------------------------------------------------------
fvSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) leastSquares;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
}
laplacianSchemes
{
default none;
laplacian(nu,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p ;
}
-------------------------------------------------------------
p
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 1;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 1;
}
fixedWalls
{
type zeroGradient;
}
frontAndBack
{
type empty;
}

U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0.1 0 0);
}
outlet
{
type zeroGradient;
}
fixedWalls
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;
}
---------------------------------------
Grateful for the support
Sargentini is offline   Reply With Quote

Old   May 29, 2012, 11:52
Default Pressure is fluctuating ,Complementing
  #4
New Member
 
Sargentini
Join Date: May 2012
Location: Brazil
Posts: 5
Rep Power: 5
Sargentini is on a distinguished road
Complementing:

HerschelBulkley.C

);
return

(
min

(
nu0_,
(tau0_ + k_*rtone*pow(tone*sr(), n_))
/(max(sr(), dimensionedScalar ("VSMALL", dimless/dimTime,VSMALL
)))
)
);
}
Sargentini is offline   Reply With Quote

Old   June 2, 2012, 08:44
Default
  #5
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
Hi!

maybe you could upload a residual plot so we can take a look at your "fluctuation"

your running a transient case (1st order euler scheme), plus you have 2 non ortoghonal correctors.. it means the pressure will iterate twice more than velocities..hence you'll c an oscillation which isnt an oscillation..but corrections from the solver. plus it's a transient case.. and you're suppose to have an oscillation like response... now how much it will oscillates it depends on the case..

might have said the obvious but setting things in OF sts makes you forget the basics

hope it helps
__________________
Best Regards
/calim

"Elune will grant us the strength"
calim_cfd is offline   Reply With Quote

Old   June 2, 2012, 15:45
Default
  #6
New Member
 
Sargentini
Join Date: May 2012
Location: Brazil
Posts: 5
Rep Power: 5
Sargentini is on a distinguished road
Hi!

Thank you for your assistance. In fact, as I am a student in both subjects these issues are not so obvious to me.
In this case, I was wondering what the formula written in the file HerschelBulkley.H, which considers nu0 would be more appropriate to the situation with the initial flow tau0 Yield stress, but it is only indicated. How can I use it? The file can be changed? This file is only available for reading.
Sargentini is offline   Reply With Quote

Old   June 2, 2012, 20:25
Default
  #7
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 136
Rep Power: 8
calim_cfd is on a distinguished road
Hi

im not on my desktop right now so i cant search through the files and be more specific.. but i guess you're using them nonNewtonianIcoFoam solver... try to examine the tutorial for this solver and change the viscosity model accordingly i might be of more assistnce once i get to my pc..



ps: c eh br?
__________________
Best Regards
/calim

"Elune will grant us the strength"
calim_cfd 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
Fluctuating Pressure Fields (LES-pisoFoam) fluentfreak OpenFOAM Running, Solving & CFD 23 July 13, 2014 11:07
Pulsatile pressure inlet with pressure outlet a.lynchy FLUENT 3 March 23, 2012 14:45
Fluent natural ventilation pressure boundary condition pierresandre FLUENT 24 November 8, 2011 15:32
UDF to define or adjust pressure??? engahmed FLUENT 0 July 6, 2010 17:19
Neumann pressure BC and velocity field Antech Main CFD Forum 0 April 25, 2006 02:15


All times are GMT -4. The time now is 23:20.