perturbU in twoPhaseEulerFoam
 

Hi Foamers,

I want to check that the twoPhaseEulerFoam is reliable, and one of the test case is that I set the domain as one fluid (i.e. the beta=1 everywhere).

I modified the code to include pressure gradient which can give a designed Re_tau, and If I'm not wrong, when I turn on the k-epsilon model, the wall-function is automatically on. I set the y_plus=50 at the first grid.

OK, Maybe too much information. Now the problem is that when I want to use perturbU to introduce some perturbation, I get stuck. The following is the Perturb.C file. Attachment 7436

I think the problem is Ub.write(), since if I let it print out Ub[celli], it goes well, however, I just can't have the new Ub file in 0/, so what should I do to let the Ub.write() work?

Any hint will be highly appreciated!