all_good

Author: aroccon

channel/input.inp

@@ -1,21 +1,21 @@
 0           ! fresh start=0 restart =1
 0           ! iteration to start from (fresh start 0 otherwise different)
-200000        ! Final time step
-1000         ! Saving frequency
-6.28         ! Streamwise domain lenght Lx
-3.14         ! Spanwise domain width Ly
+200000      ! Final time step
+1000        ! Saving frequency
+6.28        ! Streamwise domain lenght Lx
+3.14        ! Spanwise domain width Ly
 2.0         ! Wall-normal domain height Lz=2h
-1.6         ! Stretching factor quasi~0=uniform 1.5 mid-stretching and 3.0 extreme (do not use zero)
+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.00169   ! mu (viscosity)
+0.00002     ! dt for the simulation
+0.00555     ! mu (viscosity)
 1.000       ! density
 -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
-1.0         ! Ratio eps/dx
+0.51         ! Ratio eps/dx
 0.00666666  ! kappa (thermal diffusivity)
 0.0         ! alphag for buoyancy

channel/main.f90

@@ -54,7 +54,7 @@ program main
 !real(kind=8), parameter :: beta(3)   = (/ 0.0d0, -0.122243120495896d0, -0.377756879504104d0 /)
 
 ! Enable or disable phase field 
-#define phiflag 1
+#define phiflag 0
 ! Enable or disable temperature field
 #define thetaflag 0
 
@@ -1013,7 +1013,7 @@ program main
           enddo
       enddo
    enddo
-   write(*,*) "Max divergence after correction ", maxdiv
+   !write(*,*) "Max divergence after correction ", maxdiv
 
    call MPI_Allreduce(umax,gumax,1,MPI_DOUBLE_PRECISION,MPI_MAX,MPI_COMM_WORLD, ierr)
    call MPI_Allreduce(vmax,gvmax,1,MPI_DOUBLE_PRECISION,MPI_MAX,MPI_COMM_WORLD, ierr)