Skip to content
Snippets Groups Projects
Commit a9efa2ab authored by Paul Dechamps's avatar Paul Dechamps :speech_balloon:
Browse files

(format) Fix format in boundaryLayer

parent 8de8978e
No related branches found
No related tags found
No related merge requests found
......@@ -280,9 +280,9 @@ std::vector<std::vector<bldouble>> BoundaryLayer::evalGradwrtX()
for (size_t i = 1; i < getnNodes(); ++i)
{
bldouble val = 1 / rhoe[i] *
(rhoe[i] * vt[i] * deltaStar[i] -
rhoe[i - 1] * vt[i - 1] * deltaStar[i - 1]) /
((nodes[i]->xi - nodes[i - 1]->xi) * (nodes[i]->xi - nodes[i - 1]->xi));
(rhoe[i] * vt[i] * deltaStar[i] -
rhoe[i - 1] * vt[i - 1] * deltaStar[i - 1]) /
((nodes[i]->xi - nodes[i - 1]->xi) * (nodes[i]->xi - nodes[i - 1]->xi));
gradVe_loc[i - 1][i] = -val;
gradVe_loc[i - 1][i - 1] = val;
}
......@@ -335,9 +335,9 @@ std::vector<std::vector<bldouble>> BoundaryLayer::evalGradwrtY()
for (size_t i = 1; i < getnNodes(); ++i)
{
bldouble val = 1 / rhoe[i] *
(rhoe[i] * vt[i] * deltaStar[i] -
rhoe[i - 1] * vt[i - 1] * deltaStar[i - 1]) /
((nodes[i]->xi - nodes[i - 1]->xi) * (nodes[i]->xi - nodes[i - 1]->xi));
(rhoe[i] * vt[i] * deltaStar[i] -
rhoe[i - 1] * vt[i - 1] * deltaStar[i - 1]) /
((nodes[i]->xi - nodes[i - 1]->xi) * (nodes[i]->xi - nodes[i - 1]->xi));
gradVe_loc[i - 1][i] = -val;
gradVe_loc[i - 1][i - 1] = val;
}
......
......@@ -21,9 +21,9 @@ class BLAST_API BoundaryLayer : public fwk::wSharedObject
private:
unsigned int const nVar =
5; ///< Number of variables of the partial differential problem.
5; ///< Number of variables of the partial differential problem.
bldouble nCrit = 9.0; ///< Critical amplification factor.
std::string name; ///< Name of the region.
std::string name; ///< Name of the region.
public:
Closures *closSolver; ///< Closure relations class.
......@@ -33,22 +33,17 @@ public:
std::vector<int> no;
// Boundary layer variables.
std::vector<bldouble> u; ///< Solution vector (theta, H, N, ue, Ct).
std::vector<bldouble>
Ret; ///< Reynolds number based on the momentum thickness (theta).
std::vector<bldouble> cd; ///< Local dissipation coefficient.
std::vector<bldouble> cf; ///< Local friction coefficient.
std::vector<bldouble>
Hstar; ///< Kinetic energy shape parameter (thetaStar/theta).
std::vector<bldouble>
Hstar2; ///< Density shape parameter (deltaStarStar/theta).
std::vector<bldouble> Hk; ///< Kinematic shape parameter (int(1-u/ue d_eta)).
std::vector<bldouble>
ctEq; ///< Equilibrium shear stress coefficient (turbulent BL).
std::vector<bldouble> us; ///< Equivalent normalized wall slip velocity.
std::vector<bldouble> delta; ///< Boundary layer thickness.
std::vector<bldouble>
deltaStar; ///< Displacement thickness (int(1-rho*u/rhoe*ue d_eta)).
std::vector<bldouble> u; ///< Solution vector (theta, H, N, ue, Ct).
std::vector<bldouble> Ret; //< Reynolds number based on the momentum thickness (theta).
std::vector<bldouble> cd; ///< Local dissipation coefficient.
std::vector<bldouble> cf; ///< Local friction coefficient.
std::vector<bldouble> Hstar; ///< Kinetic energy shape parameter (thetaStar/theta).
std::vector<bldouble> Hstar2; ///< Density shape parameter (deltaStarStar/theta).
std::vector<bldouble> Hk; ///< Kinematic shape parameter (int(1-u/ue d_eta)).
std::vector<bldouble> ctEq; ///< Equilibrium shear stress coefficient (turbulent BL).
std::vector<bldouble> us; ///< Equivalent normalized wall slip velocity.
std::vector<bldouble> delta; ///< Boundary layer thickness.
std::vector<bldouble> deltaStar; ///< Displacement thickness (int(1-rho*u/rhoe*ue d_eta)).
std::vector<bldouble> Me; ///< Mach number.
std::vector<bldouble> vt; ///< Velocity.
......@@ -63,7 +58,7 @@ public:
private:
bldouble xtrF; ///< Forced transition location in the local frame of reference
///< (x/c).
///< (x/c).
bldouble xtr; ///< Transition location in the local frame of reference (x/c).
void reset();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment