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

correction CL pour flux Bl, mais pas de changement pour les résultats

parent 1ee6b885
No related branches found
No related tags found
No related merge requests found
......@@ -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 ;
}
......@@ -421,6 +412,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);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment