From aab22b2b3331760f22e3c11d9af52c2267a9d899 Mon Sep 17 00:00:00 2001
From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr>
Date: Tue, 26 Jun 2018 14:09:07 +0200
Subject: [PATCH] =?UTF-8?q?correction=20CL=20pour=20flux=20Bl,=20mais=20pa?=
=?UTF-8?q?s=20de=20changement=20pour=20les=20r=C3=A9sultats?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
src/scheme/FiniteVolumesDiffusion.hpp | 13 ++++---------
1 file changed, 4 insertions(+), 9 deletions(-)
diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp
index 49d246485..8fa7ce14e 100644
--- a/src/scheme/FiniteVolumesDiffusion.hpp
+++ b/src/scheme/FiniteVolumesDiffusion.hpp
@@ -240,20 +240,11 @@ private:
// Conditions aux bords
- // Diffusion pure
- /*
int cell_here = face_cells(0,0);
m_Bl(0) = (nuL(0) + nuj(cell_here))*(1./(2*Vl(0)))*(Tj(cell_here) - TL(0));
cell_here = face_cells(m_mesh.numberOfFaces()-1,0);
m_Bl(m_mesh.numberOfFaces()-1) = -(nuR(0) + nuj(cell_here))*(1/(2.*Vl(m_mesh.numberOfFaces()-1)))*(Tj(cell_here) - TR(0));
- */
-
- // Kidder
-
- double h = std::sqrt(1. - (t*t)/(50./9.));
- m_Bl(0) = ((1. + x0[0][0])*3.*x0[0][0])/(100.*h*h*h*h);
- m_Bl(m_mesh.numberOfFaces()-1) = ((1. + xmax[0][0])*3.*xmax[0][0])/(100.*h*h*h*h);
return m_Bl ;
}
@@ -420,6 +411,10 @@ public:
// Diffusion pure
//TR(0) = 2-0.5*pi*pi*(std::exp(-2.*t)-1.);
+
+ // Kidder
+ TL(0) = (1./(100*h*h))*((3.*x0[0][0]*x0[0][0])/(h*h) + 100.);
+ TR(0) = (1./(100*h*h))*((3.*xmax[0][0]*xmax[0][0])/(h*h) + 100.);
// Calcule les flux
computeExplicitFluxes(uj, Cjr, kj, uL, uR, kL, kR, Tj, nuj, TL, TR, nuL, nuR, t);
--
GitLab