Remove the two darcy test cases that were there in 1d and replace them with versions that do 1,2 and 3d.
Change the sign of the ct_rhs add to in Divergence_CV.F90 for incompressible flow.
Adjust the darcy test that was there to test this so that it runs and the 1,2 and 3d version use weak velocity BC.
Add a similar test that uses two phase but one with a large absorption such as to be the solid porous media. This tests the multiphase ct_rhs add to sign change.
the absorption term is included in the pressure correction (being a necessity as it is the only term in the matrix equation)
24
24
25
the geometry is 3d (although this is a 1d problem) and one time step is done (as nothing depends on time)
26
27
to get this to work properly a weak inlet velocity boundary condition and a strong zero pressure outlet boundary condition was found necessary</comment>
25
the geometry is 3d (although this is a 1d problem) and one time step is done (as nothing depends on time)</comment>