Need verification for interSettlingFoam for openFoam 2.3.0
 

Recently I came across the interSettlingFoam solver by Pedram Ramin (Chalmers University). It was made for version 2.2.0. After few hick-ups I managed to compile it for openFoam 2.3.0 (except unused variable ‘maxDeltaT’ warning).

Could some one please check and verify if the modifications are correct.
Following are the modifications.

1. In pEqn.H file (line-12) replace fvc::ddtPhiCorr(rAU, rho, U, phi) with fvc::interpolate(rho*rAU)*fvc::ddtCorr(U, phi).
2. In interSettlingFoam.C replace #include "twoPhaseMixture.H" with #include "incompress ibleTwoPhaseMixture.H". Add #include "fixedFluxPressureFvPatchScalarField.H".
3. In createFields.H (line-33) replace twoPhaseMixture twoPhaseProperties(U, phi); with incompressibleTwoPhaseMixture twoPhaseProperties(U, phi);.

Here is link to the solver.

Hi,

I downloaded your solver and I've been successful to compile it with the maxdeltaT error that you mention.

I tried to make it run with one of my initial condition I used with settling but it looks like the AlphaS dramatically increase instead of settling out.

Have you been sucessful to get good results with this solver ?

Many thanks,

Xavier

Hello Xavier, thanks for your interest in my post. It has been a while since I last had a look at this solver. I tried the dahl tutorial from settlingFoam solver. It seemed to run fine (Again it has been a while). You can run and have a look. Here is the link.

Hi Vasava. I also use settlingFoam to produce a gravity current from Bonnecaze et al 1993 - I use this solver instead of interFoam because I need to include the settling.

But my current doesn't extend further enough so I thought it could be a good idea to mix the interFOAM and settling FOAM solver.

I think the problem might be due to the missing advection term in the Drift equation (i.e AlphaEqn.H). I reckon I have to implement or include the advectiveFvPatchFields.H somewhere - most likely in AlphaEqn.H. By luck, did you check something like this ?