CFD Online Discussion Forums - Pressure gradient flows boundary condition

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)

- - Pressure gradient flows boundary condition (https://www.cfd-online.com/Forums/openfoam/244019-pressure-gradient-flows-boundary-condition.html)

Pressure gradient flows boundary condition

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,

You can keep velocity BC as zeroGradient

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

even pressureInletOutletVelocity in U at inlet didn't work. same floating point exception error.

Check Schems in your case and try using relaxation factors. If you using tet mesh means use nNonOrthogonalCorrectors in fvSolution