|
[Sponsors] |
July 18, 2022, 03:39 |
Pressure gradient flows boundary condition
|
#1 |
Member
Jigs
Join Date: Sep 2017
Location: Pune
Posts: 42
Rep Power: 9 |
Hi Team,
I am trying to simulate a problem with pressure condition at inlet and outlet. The same problem works great with fixed value at inlet velocity but doesn't converge for pressure inputs at inlet and outlet. Let's say I am trying to simulate with pressure inlet at 1 psi and atmospheric pressure at outlet. Used fixedValue for inlet and outlet for pressure condition. (I did divide the density with the pressure to apply pressure gradient here as input value) - Hope this part is correct. Do correct me if I am wrong. Now what i am confused is with the velocity (U) at inlet and outlet. Should I keep zeroGradient or any other type i should be using. Regards, |
|
July 20, 2022, 15:28 |
|
#2 |
New Member
Pavithran R
Join Date: Mar 2021
Posts: 2
Rep Power: 0 |
You can keep velocity BC as zeroGradient
|
|
July 21, 2022, 02:11 |
|
#3 |
Member
Jigs
Join Date: Sep 2017
Location: Pune
Posts: 42
Rep Power: 9 |
I did try that and it isn't converging. The simulation stops after few steps throwing below error.
Code:
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0835649, No Iterations 9 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0832669, No Iterations 9 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0962039, No Iterations 6 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:? #6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? #7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? #8 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:? #9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam" #10 Foam::fvMatrix<double>::solve() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam" #11 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam" #12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #13 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/simpleFoam" Floating point exception |
|
July 21, 2022, 02:54 |
|
#4 |
Member
Jigs
Join Date: Sep 2017
Location: Pune
Posts: 42
Rep Power: 9 |
even pressureInletOutletVelocity in U at inlet didn't work. same floating point exception error.
|
|
July 21, 2022, 14:30 |
|
#5 |
New Member
Pavithran R
Join Date: Mar 2021
Posts: 2
Rep Power: 0 |
Check Schems in your case and try using relaxation factors. If you using tet mesh means use nNonOrthogonalCorrectors in fvSolution
|
|
Tags |
boundary condition, pressure gradient |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wind tunnel flow simulation boundary condition issue | charan3007 | SU2 | 0 | October 21, 2021 09:27 |
CFD analaysis of Pelton turbine | amodpanthee | CFX | 31 | April 19, 2018 19:02 |
Pressure Inlet Boundary Condition Issue | zoeburton1987 | FLUENT | 0 | May 15, 2012 10:20 |
Convective Heat Transfer - Heat Exchanger | Mark | CFX | 6 | November 15, 2004 16:55 |
Pressure Boundary Condition | Matt Umbel | Main CFD Forum | 0 | January 11, 2002 11:06 |