[MattS]

I posted this in the main OpenFOAM forum and got no replies, I am hoping it was just the wrong spot for the question:

I am currently trying to use the new wave toolbox in OpenFoam-dev with compressibleMultiphaseInterFoam. (I have previously used waves2Foam with OF2.1)

I want to use compressibleMultiphaseInterFoam to be able to solve for the temperature field in addition to the velocity field. The simulation is a simple channel flow with a water phase and an overlying air phase, think numerical wave tank.

I am setting the waves with setWaves and the boundary condition for the liquid phase (phase.water) appears to work fine, straightforward:

inlet
    {
        type            waveAlpha;
        U                U;
        inletOutlet    true;

    }

I am also using setFields with fieldToCell based on alpha.water to initialize the rho for the air and water phases and that appears ok.

However, I am having trouble with the inlet boundary condition for the air phase. After some research (see other thread OpenFoam-dev waveAlpha for multiple phases) I settled on using the inlet BC for the air phase (phase.air) as follows:

inlet
    {
        type            waveAlpha;
        U                U;
        inletOutlet    true;
        liquid           no;
 }

I am initializing the waves by running
setWaves -alpha alpha.water
then
setWaves -alpha alpha.air

But with this, while everything appears to initialize correctly, as the simulation goes on, the inlet switches back to the water phase (instead of air), essentially my alpha.air and my alpha.water at the inlet have the same value at t > 0 with the above inlet boundary conditions for alpha.water and alpha.air, respectively.

I would so much appreciate any leads on how to troubleshoot this. Could this maybe be a bug in the waveAlpha condition?

Thank you!

(I am not allowed to post a picture, sorry)