[foamF]

Hi, everyone. I read a good example of effluent disposal mixing zone into ocean in the following link:

https://www.hatarilabs.com/ih-en/mod...nfoam-tutorial

In this example, interMixingFoam is worked with 3 phases - air, effluent and sea water; where effluent and sea water are the miscible liquids with different physical properties, e.g. density.

In the example, a post processing python script is used with paraView to compute the spread of pollutant in the mixture of effluent and sea water. It appears that the effect of advection (due to bulk fluid movement) and diffusion (due to the concentration gradient) of pollutants in the mixture of effluent and sea water are not considered in the model.* It looks a good approximation for near-field that advection and diffusion of pollutants may be negligible in initial dilution.* But, the effect might be significant in far-field.*

I am thinking the feasibility of using passive scalar function object in controlDict.* However, to my understanding, the utility works only for a single phase (which needs to be specified in the function object), whilst the movement of pollutants across two fluid phases (sea water and effluent) does not work.*

Any experts can give me some hints of using the standard solvers / utilities for the passive scalar transport across two miscible phases?