Need verification for interSettlingFoam for openFoam 2.3.0

vasava · November 26, 2014, 03:00

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.

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

Here is link to the solver.

Navip · July 24, 2015, 04:33

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

vasava · July 29, 2015, 03:20

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.

Navip · July 29, 2015, 03:51

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 ?

November 26, 2014, 03:00	Need verification for interSettlingFoam for openFoam 2.3.0	#1
vasava Senior Member Paritosh Vasava Join Date: Oct 2012 Location: Lappeenranta, Finland Posts: 732 Rep Power: 22	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. In pEqn.H file (line-12) replace *fvc::ddtPhiCorr(rAU, rho, U, phi)* with **fvc::interpolate(rhorAU)fvc::ddtCorr(U, phi). In interSettlingFoam.C replace #include "twoPhaseMixture.H" with #include "incompress ibleTwoPhaseMixture.H". Add #include "fixedFluxPressureFvPatchScalarField.H". In createFields.H (line-33) replace twoPhaseMixture twoPhaseProperties(U, phi); with incompressibleTwoPhaseMixture twoPhaseProperties(U, phi);**. Here is link to the solver.

July 29, 2015, 03:20		#3
vasava Senior Member Paritosh Vasava Join Date: Oct 2012 Location: Lappeenranta, Finland Posts: 732 Rep Power: 22	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. foamF likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Multiple Floating Objects with 6DoF in OF 2.3.0	tyoung	OpenFOAM Running, Solving & CFD	22	Yesterday 00:36
Liquid Evaporation Model Error in OF 2.3.0	brbbhatti	OpenFOAM	11	June 16, 2014 09:40
[snappyHexMesh] Holes in cylinder AMI snappy 2.3.0	pbachant	OpenFOAM Meshing & Mesh Conversion	3	March 14, 2014 21:37
Grid Verification	julie corish	Main CFD Forum	3	January 30, 2003 09:47
Verification data	Achilleas Tsompanos	Main CFD Forum	1	June 19, 2000 12:49

July 24, 2015, 04:33		#2
Navip New Member Xavier Pivan Join Date: May 2015 Posts: 10 Rep Power: 10	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

July 29, 2015, 03:51		#4
Navip New Member Xavier Pivan Join Date: May 2015 Posts: 10 Rep Power: 10	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 ?