The contraction geometry: boundary conditions
struct base {
string sys_coord = omega.coordinate_system_name();
}
};
struct u_upstream: base {
};
space Xh (omega, approx,
"vector");
Xh.block ("wall");
Xh.block ("upstream");
Xh[1].block ("axis");
Xh[1].block ("downstream");
return Xh;
}
geo omega = Xh.get_geo();
string approx =
"P" +
itos(Xh.degree());
space Wh (omega[
"upstream"], approx);
uh[0][
"upstream"] =
interpolate (Wh, u_upstream(omega));
return uh;
}
space Ph (omega, approx);
Ph.block("upstream");
Ph.block("wall");
Ph.block("axis");
return Ph;
}
struct psi_upstream: base {
} else {
}
}
};
geo omega = Ph.get_geo();
space Wh (omega[
"upstream"], Ph.get_approx());
psi_upstream psi_up (omega);
psi_h["wall"] = 0;
psi_h["axis"] = -1;
return psi_h;
}
};