Pre-correction (called per substep after the integrator step): enforce the inflow Dirichlet and convective outflow BCs, then apply an exact uniform mass correction to the outlet, following Incompact3d's pre_correc.
Steps: 1. cx = uxmaxgdt/dx from the post-integration field. 2. Convective outflow at x=nx and Dirichlet inflow at x=1 for u,v,w (no correction yet): u(nx) = u(nx) - cx(u(nx) - u(nx-1)). 3. Measure the inlet mean (from the prescribed Dirichlet plane) and the outlet mean AFTER the convective update, then shift the outlet plane of u uniformly so mean(u_outlet) = mean(u_inlet) exactly (Incompact3d: bxxn = bxxn - ut + ut1). Applied to u only.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(case_wind_turbine_t) | :: | self | ||||
| class(field_t), | intent(inout) | :: | u | |||
| class(field_t), | intent(inout) | :: | v | |||
| class(field_t), | intent(inout) | :: | w |