Add New Example and Doc

COM.f03

@@ -63,3 +63,65 @@ integer Function sgn(x)
 end Function sgn
 
 end module COM
+
+! Following is a list of identifiers in module COM with their meanings
+! Grid                Data structure for surface multi-block structured mesh
+! Force               Data structure for driven force and collection property
+! Energy              Data structure for heat transfer related coefficients
+! Flux                Data structure for flow flux of water film
+! Imp                 Data structure for matrix coefficients and source terms of an implicit water film flow equation
+! Icecoordinate       Data structure for thickness of water and ice (unknown variables) and coordinates of ice shape
+! Boundata            Data structure for boundary field
+! topo                Data structure for topology information
+! block               Patches of the surface structured mesh
+! maxl                Maximum iteration steps for solving the linear equations in an implicit time marching
+! timestep            Current time steps in a time marching
+! timeout             Total time steps in a time marching
+! pi                  Circular constant
+! g                   Gravitational acceleration
+! R                   Universal gas constant
+! Ma                  Molecular weight of air
+! rhoa                Density of air
+! rhow                Density of water
+! rhoi                Density of ice
+! mua                Viscosity of air
+! muw                Viscosity of water
+! ka                  Thermal conductivity of air
+! kw                  Thermal conductivity of water
+! ki                  Thermal conductivity of ice
+! ca                  Specific heat at constant pressure of air
+! cw                  Specific heat of water
+! Pr                  Prandtl number
+! Prt                 Turbulent Prandtl number
+! P0                  Pressure of free stream
+! Wf                  Velocity of free stream
+! epsi                Heat emissivity
+! sigmar              Stefan constant
+! lwc                 Liquid Water Content
+! mvd                 Mean Volume Diameter
+! Ta                  Temperature of air
+! Td                  Temperature of droplets
+! Ts                  Temperature of substract
+! Tf                  Fusion point of ice
+! Lf                  Latent heat of fusion
+! Lv                  Latent heat of evaporation
+! c                   Scaling factor
+! hp,bp               Precursor thickness of water and ice
+! dt                  Time step
+! t                   Total time
+! time                Current time
+! err                 Minimum residuals for solving the linear equations in an implicit time marching
+! alpha               Empirical factor of thermal conduction in ice layer
+! surtempcontrol      Whether using an inner iteration to calulate the surface temperature
+! solutioncontrol     Whether using explicit or implicit time marching method
+! initialtimecontrol  Whether initializing from a clean surface
+! discretecontrol     The discrete method for flow flux at cell interface
+! icecoupled          Whether coupling with the icing equation
+! advancemethod       The pseudo-steady method
+! fluidproperty       The air property
+! Scpt                Whether using a script file
+! filename            File name
+! In,Jn               Dimensions of the grid nodes in each direction of one patch
+! Ig,Jg               Dimensions of the grid cells in each direction of one patch
+! formn,formc         Formatted strings
+! sgn()               Sign function

Makefile

@@ -2,7 +2,9 @@ vpath %.o obj
 vpath %.exe bin
 FC=gfortran
 FCFLAGS=-W -fdefault-real-8 -fdefault-double-8 -ffree-line-length-none -O3
+#FCFLAGS+=-ID:\Libsrc\hypre-2.11.1\src\hypre\dll\include
 LDFLAGS=-O3 -s
+#LDFLAGS+=-LD:\Libsrc\hypre-2.11.1\src\hypre\dll\lib -lHYPRE
 src:=$(sort $(wildcard *.f03))
 objects:=$(src:.f03=.o)
 exec:=ICE3D.exe
@@ -16,11 +18,11 @@ rd:=rm -r -f
 del:=rm -f
 cp:=cp -f
 %.o:%.f03 COM.mod
-	$(FC) -c $(FCFLAGS) -I$(odir) $< -o $(addprefix $(oprefix),$@)
+	$(FC) $(FCFLAGS) -I$(odir) -c $< -o $(addprefix $(oprefix),$@)
 %.mod:%.f03|$(odir)
-	$(FC) -c $(FCFLAGS) -J$(odir) $< -o $(addprefix $(oprefix),$(subst .mod,.o,$@))
+	$(FC) $(FCFLAGS) -J$(odir) -c $< -o $(addprefix $(oprefix),$(subst .mod,.o,$@))
 $(exec):$(objects)|$(bdir)
-	$(FC) $(LDFLAGS) $(addprefix $(oprefix),$(objects)) -o $(addprefix $(bprefix),$(exec))
+	$(FC) $(addprefix $(oprefix),$(objects)) $(LDFLAGS) -o $(addprefix $(bprefix),$(exec))
 $(odir):
 	-$(md) $(odir)
 $(bdir):

cases/GLC305/Script.txt

@@ -3,9 +3,9 @@ GLC305surf.xyz
 Give name of surface grids file:
 GLC305surf.xyz
 Give name of air flow property file:
-wallflux(GLC305).dat
+wallflux_GLC305.dat
 Give name of drop property file:
-drop(GLC305).dat
+drop_GLC305.dat
 Set air fluid property(com/incom):
 com
 Select march method(singl/multi):

cases/Sphere/Script.txt

@@ -3,9 +3,9 @@ spheresurf.xyz
 Give name of surface grids file:
 spheresurf.xyz
 Give name of air flow property file:
-wallflux(sphere).dat
+wallflux_sphere.dat
 Give name of drop property file:
-drop(sphere).dat
+drop_sphere.dat
 Set air fluid property(com/incom):
 com
 Select march method(singl/multi):

cases/WindBody/Script.txt

@@ -0,0 +1,46 @@
+Give name of initial surface grids file:
+wingbodysurf.xyz
+Give name of surface grids file:
+wingbodysurf.xyz
+Give name of air flow property file:
+wallflux_wingbody.dat
+Give name of drop property file:
+drop_wingbody.dat
+Set air fluid property(com/incom):
+com
+Select march method(singl/multi):
+singl
+Include Ice Accretion(Y/N)?
+Y
+Input grid scale factor:
+1.0
+Input liquid water content(kg/m3):
+6.8e-4
+Input mean volume diameter(m):
+2e-5
+Input air far field tempreture(K):
+262.04
+Input drop tempreture(K):
+262.04
+Input initial solid surface tempreture(K):
+262.04
+Input air far field velocity(m/s):
+75
+Give time step(s):
+5e-3
+Give total march time(s):
+180
+Start march from clean surface(Y/N)?
+Y
+Give name of initial data file:(pre film height(m):)
+5e-9
+Start inner surface tempreture iteration(Y/N)?
+N
+Select solution method(explicit/implicit)?
+implicit
+Select discretization scheme(upwind/TVD)?
+upwind
+Give max step of inner implicit iteration:
+100
+Give residual of inner implicit iteration:
+1e-8