The contraction geometry: boundary conditions.
The contraction geometry: boundary conditions
struct base {
string sys_coord = omega.coordinate_system_name();
}
};
struct u_upstream: base {
u_upstream (
geo omega) : base(omega) {}
};
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" + to_string(Xh.degree());
space Wh (omega[
"upstream"], approx);
uh[0]["upstream"] = lazy_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 {
psi_upstream (
geo omega) : base(omega) {}
} else {
}
}
};
geo omega = Ph.get_geo();
space Wh (omega[
"upstream"], Ph.get_approx());
psi_upstream psi_up (omega);
psi_h["upstream"] = lazy_interpolate (Wh, psi_up);
psi_h["wall"] = 0;
psi_h["axis"] = -1;
return psi_h;
}
};
see the Float page for the full documentation
see the field page for the full documentation
see the geo page for the full documentation
see the point page for the full documentation
see the space page for the full documentation
Float operator()(const point &x) const
Float operator()(const point &x) const
static field streamf_field(space Ph)
static space velocity_space(geo omega, string approx)
static space streamf_space(geo omega, string approx)
static field velocity_field(space Xh)