From d8034f7d837652317fdd3f36e15ca7243bd635dd Mon Sep 17 00:00:00 2001
From: Adrien Magni <Adrien.Magni@gmail.com>
Date: Mon, 20 Jun 2011 14:39:03 +0000
Subject: [PATCH] shpere dans un canal
---
NavierStokes3D-Penalization/advection_mod.f90 | 4 ++--
NavierStokes3D-Penalization/main.f90 | 2 +-
NavierStokes3D-Penalization/parameter | 6 +++---
3 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/NavierStokes3D-Penalization/advection_mod.f90 b/NavierStokes3D-Penalization/advection_mod.f90
index 1838de551..88c87c21f 100644
--- a/NavierStokes3D-Penalization/advection_mod.f90
+++ b/NavierStokes3D-Penalization/advection_mod.f90
@@ -260,8 +260,8 @@ contains
r2=kx**2+ky**2+kz**2
coef=-fac*cmplx(0.,1.)/(2.*pi*r2)
- !uinfx=1/((yh-yb)*(zf-zd))!debit/surface du domaine numérique
- uinfx=(4.*(1./0.4844)**2)/((yh-yb)*(zf-zd))!debit/surface du domaine numérique
+
+ uinfx=(4.*(1./0.4844)**2)/((yh-yb)*(zf-zd))!debit*surface domaine physique/surface du domaine numérique
!if ((time>=3.).and.(time<=4.)) then
! uinfy=sin(pi*(time-3.))
diff --git a/NavierStokes3D-Penalization/main.f90 b/NavierStokes3D-Penalization/main.f90
index cff4bcf59..c7d086765 100644
--- a/NavierStokes3D-Penalization/main.f90
+++ b/NavierStokes3D-Penalization/main.f90
@@ -66,7 +66,7 @@ program NS
!dt=0.0125
!dt=0.001
!dt=min(cfl*dx/maxv,cfl pour diffusion)
- dt=0.05
+ dt=0.005
!dt=0.01
!dt=3./8.*dx**2/nu
diff --git a/NavierStokes3D-Penalization/parameter b/NavierStokes3D-Penalization/parameter
index dd53e8214..33f2ab81a 100644
--- a/NavierStokes3D-Penalization/parameter
+++ b/NavierStokes3D-Penalization/parameter
@@ -1,12 +1,12 @@
-2.06441 2.06441 -3. 3. -2.06441 2.06441 ! xg xd yb yh zd zf
-33 33 33 33 65 129 257!nx ny nz
+129 129 129 33 65 129 257!nx ny nz
30. ! TFin
1.D-8 !cutoff
4 ! type bloc
2 !longeur bloc
"test.vtk" ! sortie fin vtk
1.D6 !lambda
-7.496251874D-03 (Re=133.4) 4.74833808D-03 (Re=210.6) !nu=1/Re
+4.74833808D-03 (Re=210.6) 7.496251874D-03 (Re=133.4) !nu=1/Re
"RES/test/omg_" "RES/test/vg_" 2000000 !vtk
-"RES/test/coupe1.res" "RES/test/test.drag" !coupe vitesse drag/lift
+"RES/test/test.res" "RES/test/test.drag" !coupe vitesse drag/lift
0 200000 "RES/test/sauv.vtk" !num_suite (0non 1 oui) num_ite name_relance
\ No newline at end of file
--
GitLab