I added velocity in PrepareLattice like this :
”
std::vector<T> Vg( 2,T(velocity) );
AnalyticalConst2D<T,T> constantVelocity( Vg);
// add second fluid (tank)
SmoothIndicatorCuboid2D<T,T> cuboid( {T(nx2)/T(2),T(nx)/T(2)}, T(nx2),T(nx+2), converter.getPhysLength(alpha) );
AnalyticalIdentity2D<T,T> rho( one );
AnalyticalIdentity2D<T,T> phi( one – cuboid – cuboid );
“
