channel

Author: aroccon

channel/input.inp

@@ -2,17 +2,17 @@
 0           ! iteration to start from (fresh start 0 otherwise different)
 200000        ! Final time step
 1000         ! Saving frequency
-4.0         ! Streamwise domain lenght Lx
-4.0         ! Spanwise domain width Ly
+6.28         ! Streamwise domain lenght Lx
+3.14         ! Spanwise domain width Ly
 2.0         ! Wall-normal domain height Lz=2h
 0.01         ! Stretching factor quasi~0=uniform 1.5 mid-stretching and 3.0 extreme (do not use zero)
 1           ! inflow condition (0=zero, 1=perturbed 3D field, 2=read from fields)
 0           ! phase-field initial condition (0=circular drop, 1=read from fields)
 0           ! theta initial condition (0=uniform, 1=read from fields)
-0.0001      ! dt for the simulation
-0.0066666   ! mu (viscosity)
+0.0005      ! dt for the simulation
+0.00666   ! mu (viscosity)
 1.000       ! density
-0.0        ! Pressure gradient along x (negative if flow along x) 
+-1.0        ! Pressure gradient along x (negative if flow along x) 
 0.0         ! Pressure gradient along y (negative if flow along y) 
 0.5         ! Initial drop radius for phase-field
 0.1d0       ! Surface tension

channel/module.f90

@@ -1,7 +1,7 @@
 module param
     integer, parameter :: nx=200
-    integer, parameter :: ny=200
-    integer, parameter :: nz=100
+    integer, parameter :: ny=100
+    integer, parameter :: nz=200
     double precision :: pi, rhoi, twopi
     double precision :: lx, dx, dxi, ddxi 
     double precision :: ly, dy, dyi, ddyi