# Pressure is fluctuating

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

 May 29, 2012, 07:30 Pressure is fluctuating #1 New Member   Sargentini Join Date: May 2012 Location: Brazil Posts: 5 Rep Power: 6 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

 May 29, 2012, 07:34 #2 Senior Member   Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 14 Maybe you can give some more relevant information, such as discretization and convergence criteria used? Also, which solver are you using?

 May 29, 2012, 08:08 Pressure is fluctuating #3 New Member   Sargentini Join Date: May 2012 Location: Brazil Posts: 5 Rep Power: 6 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

 May 29, 2012, 11:52 Pressure is fluctuating ,Complementing #4 New Member   Sargentini Join Date: May 2012 Location: Brazil Posts: 5 Rep Power: 6 Complementing: HerschelBulkley.C );return ( min ( nu0_, (tau0_ + k_*rtone*pow(tone*sr(), n_)) /(max(sr(), dimensionedScalar ("VSMALL", dimless/dimTime,VSMALL ))) ) ); }

 June 2, 2012, 08:44 #5 Senior Member     mauricio Join Date: Jun 2011 Posts: 141 Rep Power: 9 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"

 June 2, 2012, 15:45 #6 New Member   Sargentini Join Date: May 2012 Location: Brazil Posts: 5 Rep Power: 6 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.

 June 2, 2012, 20:25 #7 Senior Member     mauricio Join Date: Jun 2011 Posts: 141 Rep Power: 9 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"

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post fluentfreak OpenFOAM Running, Solving & CFD 24 August 30, 2016 15:08 a.lynchy FLUENT 3 March 23, 2012 14:45 pierresandre FLUENT 24 November 8, 2011 15:32 engahmed FLUENT 0 July 6, 2010 17:19 Antech Main CFD Forum 0 April 25, 2006 02:15

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