[junewookee]

Hello,

I tried to simulate the water droplet impinging onto flat surface using interFoam solver.

I countered some problems.

I know that interFoam deals with fluid flow only. So defining solid zone is hard for me.

I modified UEqn.H shown below.

SL(solid layer) is defined as volScalarField in createFields.H.
If SL is 1, then cell is considered as solid cell. So there's no velocity values in solid cell.
If SL is 0, then the cell is considered as fluid cell. So momentum equation will be calculated.

@ UEqn.H

fvVectorMatrix UEqn
(
fvm::ddt(rho, U) * pos (0.01 - mag(fvc::interpolate(SL)))
+ fvm::div(rhoPhi, U) * pos (0.01 - mag(fvc::interpolate(SL)))
- fvm::laplacian(muEff, U) * pos (0.01 - mag(fvc::interpolate(SL)))
- (fvc::grad(U) & fvc::grad(muEff)) * pos (0.01 - mag(fvc::interpolate(SL)))
);

UEqn.relax();

if (pimple.momentumPredictor())
{
solve
(
UEqn
==
fvc::reconstruct
(
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1) * pos (0.01 - mag(fvc::interpolate(SL)))
- ghf*fvc::snGrad(rho) * pos (0.01 - mag(fvc::interpolate(SL)))
- fvc::snGrad(p_rgh) * pos (0.01 - mag(fvc::interpolate(SL)))
) * mesh.magSf()
)
);
}

When i compile the modified solver, the error message came out like this.
'no match for 'operator*' ..... '
at the first transient term of momentum equation (bold and italic fonts)

I want to know why this is wrong.

Thank you.

Sincerely.