From 34bef42527ef6ee575da3c06482ccd4d5f63c203 Mon Sep 17 00:00:00 2001
From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr>
Date: Tue, 12 Jun 2018 17:04:59 +0200
Subject: [PATCH] ajout cas k pour Re = 500
---
src/main.cpp | 8 ++++----
src/scheme/FiniteVolumesEulerUnknowns.hpp | 3 ++-
2 files changed, 6 insertions(+), 5 deletions(-)
diff --git a/src/main.cpp b/src/main.cpp
index 0a83b28d9..6325e25b9 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -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);
diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp
index b1abbde58..98e5859ec 100644
--- a/src/scheme/FiniteVolumesEulerUnknowns.hpp
+++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp
@@ -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. ;
}
--
GitLab