Skip to content
Snippets Groups Projects
Commit ce8e30bf authored by Fanny CHOPOT's avatar Fanny CHOPOT
Browse files

ajout cas k regulier

parent e4898ac5
Branches
No related tags found
No related merge requests found
...@@ -249,12 +249,13 @@ int main(int argc, char *argv[]) ...@@ -249,12 +249,13 @@ int main(int argc, char *argv[])
std::ofstream diff("diffinter"); std::ofstream diff("diffinter");
diff.precision(5); diff.precision(5);
for (size_t j=0; j<mesh.numberOfCells(); ++j) { for (size_t j=0; j<mesh.numberOfCells(); ++j) {
//diff << std::fixed << xj[j][0] << ' ' << std::fixed << kj[j] << '\n'; diff << std::fixed << xj[j][0] << ' ' << std::fixed << kj[j] << '\n';
if (kj[j]>0.) { /* if (kj[j]>0.) {
diff << std::fixed << xj[j][0] << ' ' << std::fixed << 4. << '\n'; diff << std::fixed << xj[j][0] << ' ' << std::fixed << 4. << '\n';
} else { } else {
diff << std::fixed << xj[j][0] << ' ' << std::fixed << -0.1 << '\n'; diff << std::fixed << xj[j][0] << ' ' << std::fixed << -0.1 << '\n';
} }
*/
} }
diff.close(); diff.close();
......
...@@ -275,10 +275,15 @@ public: ...@@ -275,10 +275,15 @@ public:
//m_kj[j] = xj[j][0]; //m_kj[j] = xj[j][0];
//m_kj[j] = 0.5; //m_kj[j] = 0.5;
// Sod
double haut = 0.0028; // Re = 500
if (xj[j][0]<0.7) { if (xj[j][0]<0.7) {
m_kj[j]=0.; m_kj[j]=0.;
} else { } else {
if (xj[j][0]<0.9){ if (xj[j][0]<0.9){
// k non regulier
// Re = 28 // Re = 28
// m_kj[j]=0.05; // m_kj[j]=0.05;
// Re = 500 // Re = 500
...@@ -286,7 +291,12 @@ public: ...@@ -286,7 +291,12 @@ public:
// Re = 1000 // Re = 1000
//m_kj[j] = 0.0014; //m_kj[j] = 0.0014;
// Re = 2000 // Re = 2000
m_kj[j] = 0.0007; //m_kj[j] = 0.0007;
// k regulier
m_kj[j] = (haut*std::exp(100))*std::exp(-1./(0.1*0.1-(xj[j][0]-0.8)*(xj[j][0]-0.8)));
} else { } else {
m_kj[j]=0. ; m_kj[j]=0. ;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment