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

verification convergence de k_n et mise en commentaire des fichiers-films

parent a9da2bd1
No related branches found
No related tags found
No related merge requests found
...@@ -160,7 +160,7 @@ int main(int argc, char *argv[]) ...@@ -160,7 +160,7 @@ int main(int argc, char *argv[])
double c = 0.; double c = 0.;
c = finite_volumes_diffusion.conservatif(unknowns); c = finite_volumes_diffusion.conservatif(unknowns);
/*
// Ecriture des valeurs initiales dans un fichier // Ecriture des valeurs initiales dans un fichier
const Kokkos::View<const Rd*> xj = mesh_data.xj(); const Kokkos::View<const Rd*> xj = mesh_data.xj();
...@@ -249,16 +249,15 @@ int main(int argc, char *argv[]) ...@@ -249,16 +249,15 @@ 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();
*/
while((t<tmax) and (iteration<itermax)) { while((t<tmax) and (iteration<itermax)) {
......
...@@ -287,6 +287,9 @@ public: ...@@ -287,6 +287,9 @@ public:
// Re = 28 // Re = 28
// m_kj[j]=0.05; // m_kj[j]=0.05;
// Re = 100
//m_kj[j] = 0.014;
// Re = 500 // Re = 500
//m_kj[j] = 0.0028; //m_kj[j] = 0.0028;
...@@ -296,6 +299,12 @@ public: ...@@ -296,6 +299,12 @@ public:
// Re = 2000 // Re = 2000
//m_kj[j] = 0.0007; //m_kj[j] = 0.0007;
// Re = 10 000
m_kj[j] = 0.00014;
// Re = 100 000
//m_kj[j] = 0.000014;
} else { } else {
m_kj[j]=0. ; m_kj[j]=0. ;
} }
...@@ -305,6 +314,7 @@ public: ...@@ -305,6 +314,7 @@ public:
// k regulier // k regulier
int n = 1.; int n = 1.;
m_kj[j] = std::exp(1.)*std::exp(-1./(1.-( (xj[j][0]-(0.7+0.1/n)) / (0.1/n) )*( (xj[j][0]-(0.7+0.1/n)) / (0.1/n) ))) * (xj[j][0]>0.7)*(xj[j][0]<0.7+0.1/n) + std::exp(1.)*std::exp(-1./(1.-( (xj[j][0]-(0.9-0.1/n)) / (0.1/n) )*( (xj[j][0]-(0.9-0.1/n)) / (0.1/n) ))) * (xj[j][0]>0.9-0.1/n)*(xj[j][0]<0.9) + (xj[j][0]>0.7+0.1/n)*(xj[j][0]<0.9-0.1/n); m_kj[j] = std::exp(1.)*std::exp(-1./(1.-( (xj[j][0]-(0.7+0.1/n)) / (0.1/n) )*( (xj[j][0]-(0.7+0.1/n)) / (0.1/n) ))) * (xj[j][0]>0.7)*(xj[j][0]<0.7+0.1/n) + std::exp(1.)*std::exp(-1./(1.-( (xj[j][0]-(0.9-0.1/n)) / (0.1/n) )*( (xj[j][0]-(0.9-0.1/n)) / (0.1/n) ))) * (xj[j][0]>0.9-0.1/n)*(xj[j][0]<0.9) + (xj[j][0]>0.7+0.1/n)*(xj[j][0]<0.9-0.1/n);
m_kj[j] = 0.014*m_kj[j];
}); });
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment