|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Jigs
Join Date: Aug 2017
Location: Pune
Posts: 39
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, |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Pavithran R
Join Date: Mar 2021
Posts: 2
Rep Power: 0 ![]() |
You can keep velocity BC as zeroGradient
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Jigs
Join Date: Aug 2017
Location: Pune
Posts: 39
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 |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Member
Jigs
Join Date: Aug 2017
Location: Pune
Posts: 39
Rep Power: 9 ![]() |
even pressureInletOutletVelocity in U at inlet didn't work. same floating point exception error.
|
|
![]() |
![]() |
![]() |
![]() |
#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 |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wind tunnel flow simulation boundary condition issue | charan3007 | SU2 | 0 | October 21, 2021 08:27 |
CFD analaysis of Pelton turbine | amodpanthee | CFX | 31 | April 19, 2018 18:02 |
Pressure Inlet Boundary Condition Issue | zoeburton1987 | FLUENT | 0 | May 15, 2012 09:20 |
Convective Heat Transfer - Heat Exchanger | Mark | CFX | 6 | November 15, 2004 15:55 |
Pressure Boundary Condition | Matt Umbel | Main CFD Forum | 0 | January 11, 2002 10:06 |