Problems with own CFD
 

Hello everyone!
I am creating my own solver of unsteady Stokes flow in 2D(pressure at the inlet rises linearly and stop after few seconds). At first I solve and then computed pressure is put into Stokes equations. It bases on finite difference method and uses Conjugate Gradient Method to solve the pressure equation. Next the stokes equations are solved using ADI algorithm and CGM. First step of time is okey, because the solution is zero(so it is put during initiazation). Unfortunately in next step when pressure at the inlet rises pressure don't want to converge. It is linear problem and theoritically it should converge(matrix from laplace equation is symmetricall and PDM) in some cases slowly but still it should converge. Unfortunately in my case calculations are diverging. I have checked everything but I wasn't able to find an error. Have anyone got some ideas how to find where is the problem ?

I don't understand the Laplace equation for the pressure...

If take a divergence of Stokes Equation(incompressible), you will get Laplace equation for pressure.

be carefull, numerically the diffusive terms in the momentum equation produce a source term into the pressure equation, that is Div (1/Re Div Grad v). The bc for the normal derivative of the pressure must ensure the compatibility

Problem found, my matrix turn out to be not symmetricall. Thank you for advice.

solver controle
 

hi everyone
I have a question concerning the difference between conservative and aggresive length scale option for simulation

Dpm
 

Hi Dears

could you please help me about solid particle injection in liquid phase? particles have injected to inlet but why they don't escape from outlet?