DDtU.H 225 Bytes
Newer Older
shunbo's avatar
shunbo committed
1
2
3
4
5
6
7
8
9
for (phaseModel& phase : fluid.phases())
{
    phase.DDtU() =
        fvc::ddt(phase.U())
      + fvc::div(phase.phi(), phase.U())
      - fvc::div(phase.phi())*phase.U();

    MRF.addAcceleration(phase.U(), phase.DDtU());
}