Although it intuitively should work to set-up a pressure drop problem by just prescribing differing pressures, the finite element method also implicitly prescribes homogeneous Neumann (flux) boundary conditions on non-Dirichlet (fixed value) boundaries. In the case of prescribed pressure boundaries this means that the momentum equations will also have an implicit condition like "(-p*I + ux_i)*n = 0" affecting the pressure, which must be negated to get the expected solution. Alternatively one can prescribe periodic boundaries with a given pressure drop.

Unfortunately, both these approaches require some non-standard modifications of the fea data struct which currently can not be done in the FEATool GUI. However, there is an m-file script model example showing how to do this on the MATLAB command line (the full code listing can be found in the

*examples/ex_navierstokes14.m* subdirectory of the FEATool installation directory):

MATLAB m-file script example for flow in a channel driven by a pressure difference
Leon Robers wrote

I'm struggling with setting up a 2D CFD (Navier Stokes) Simulation where the flow is induced by a pressure difference.

Let's assume we're having a perfect 2D pipe and we want to apply a certain pressure at one boundary surface and the other surface to be at ambient pressure.

Independent of material properties (I tried it with air and water) and of the applied pressures, the solution will not converge. I even run a simulation with an inflow velocity as BC, then computed the pressure, and in a second simulation, I took exactly this pressure as BC. It still didn't converge...