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

ajout cas k pour Re = 500

parent f625f29d
No related branches found
No related tags found
No related merge requests found
......@@ -160,7 +160,7 @@ int main(int argc, char *argv[])
double c = 0.;
c = finite_volumes_diffusion.conservatif(unknowns);
/*
// Ecriture des valeurs initiales dans un fichier
const Kokkos::View<const Rd*> xj = mesh_data.xj();
......@@ -257,7 +257,7 @@ int main(int argc, char *argv[])
}
}
diff.close();
*/
while((t<tmax) and (iteration<itermax)) {
......@@ -325,7 +325,7 @@ int main(int argc, char *argv[])
++iteration;
std::cout << "temps t : " << t << std::endl;
/*
// ECRITURE DANS UN FICHIER
if ((std::fmod(t,0.01) < 0.0001) or (t == tmax)) {
......@@ -544,7 +544,7 @@ int main(int argc, char *argv[])
riffout.close();
}
*/
// ENTROPY TEST
//finite_volumes_diffusion.entropie(unknowns);
......
......@@ -279,7 +279,8 @@ public:
m_kj[j]=0.;
} else {
if (xj[j][0]<0.9){
m_kj[j]=0.05;
// m_kj[j]=0.05;
m_kj[j] = 0.0028; // Re = 500
} else {
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