lge: add input files for three sizes (0119a1de) · Commits · William D. Fullmer / gtsh

liu.ge.wang/build.19.3.gnu/usr1.f

+1 −1

Original line number	Diff line number	Diff line
		@@ -65,7 +65,7 @@
		!
		! data analysis & printing
		8882 FORMAT (4(E16.8, 2x))
		8885 FORMAT (4(E16.8, 2x))
		8885 FORMAT (5(E16.8, 2x))

		! Normalize vertical velocity variables to Vk' = Vk - jy
		! ie ' variables are relative to the zero net flux plane

liu.ge.wang/large/mfix.dat

+13 −14

Original line number	Diff line number	Diff line
		@@ -5,10 +5,10 @@
		RUN_TYPE = 'new' ! new, restart_1
		UNITS = 'SI'
		TIME = 0.0 ! start time
		TSTOP = 0.127 ! stop time
		DT = 1.0d-05 ! time step
		DT_MAX = 1.0d-03 ! max time-step
		DT_MIN = 1.0d-08 ! min time-step
		TSTOP = 5.08d-01 ! stop time
		DT = 5.08d-06 ! time step
		DT_MAX = 5.08d-06 ! max time-step
		DT_MIN = 2.54d-09 ! min time-step

		MOMENTUM_X_EQ(0) = .TRUE. ! solve gas x-mom eq, granular: use .FALSE.
		MOMENTUM_Y_EQ(0) = .TRUE. ! solve gas y-mom eq, granular: use .FALSE.
		@@ -22,7 +22,6 @@
		SPECIES_EQ(0) = .FALSE. ! do not solve gas species eq
		SPECIES_EQ(1) = .FALSE. ! do not solve solid species equations

		!19.3 GRANULAR_ENERGY = .TRUE. ! solve the granular energy transport equation
		KT_TYPE = 'GTSH' ! kinetic theory, granular: GD_99
		FRICTION_MODEL = 'SCHAEFFER' ! Schaeffer frictional viscosity
		BLENDING_FUNCTION = 'TANH_BLEND' ! hyperbolic tangent blending
		@@ -51,16 +50,16 @@
		! --------------------------- Geometry Section ---------------------------------------

		COORDINATES = 'cartesian' !use cartesian coordinates
		XLENGTH = 2.88d-3 !length of domain in x direction
		YLENGTH = 11.52d-3 !
		XLENGTH = 11.52d-3 !length of domain in x direction
		YLENGTH = 46.08d-3 !
		ZLENGTH = 0.48d-3 !
		IMAX = 48 !number of cells in x direction (i)
		JMAX = 192 !
		IMAX = 192 !number of cells in x direction (i)
		JMAX = 768 !
		KMAX = 8 !

		CYCLIC_X = .TRUE. !use periodic boundaries in all directions
		CYCLIC_Y_PD = .TRUE.
		DELP_Y = 170.6648831999999913d0 !pressure drop in Y-dir
		DELP_Y = 682.6595327999999654d0 !pressure drop in Y-dir
		CYCLIC_Z = .TRUE.

		GRAVITY = 98.0d+0 !acts in -ve y direction
		@@ -89,9 +88,9 @@
		!----------------------- Initial Conditions Section ----------------------------------

		IC_X_w(1) = 0.0d0 !x boundary location west
		IC_X_e(1) = 2.88d-3 !x boundary location east
		IC_X_e(1) = 11.52d-3 !x boundary location east
		IC_Y_s(1) = 0.0d0
		IC_Y_n(1) = 11.52d-3
		IC_Y_n(1) = 46.08d-3
		IC_Z_b(1) = 0.0d0
		IC_Z_t(1) = 0.48d-3
		! initial values in the region
		@@ -112,7 +111,7 @@
		NLOG = 100000 ! write logfile R.LOG every 100 time steps
		FULL_LOG = .T. ! display residuals on screen
		CALL_USR = .T. ! call the user functions
		USR_DT = 1. !2.54d-4
		USR_DT = 2.54d-4

		! SPX_DT values determine how often SPx files are written (in seconds).
		! SP1 SP2 SP3 SP4 SP5 SP6 SP7 SP8 SP9 SPA SPB
		@@ -125,7 +124,7 @@

		! The decomposition in I J and K directions for a Distributed Memory Parallel machine

		NODESI = 4 NODESJ = 10 NODESK = 1
		NODESI = 16 NODESJ = 64 NODESK = 1

		! eor EOR EOF eof

liu.ge.wang/large/run_joule2_general.q

+2 −2

Original line number	Diff line number	Diff line
		@@ -5,7 +5,7 @@
		#SBATCH --qos=normal

		##Submission
		#SBATCH --nodes=1
		#SBATCH --nodes=26
		#SBATCH --job-name="gtsh"
		#SBATCH --output=job.out

		@@ -14,4 +14,4 @@ module load gnu/6.5.0
		module load openmpi/3.1.3_gnu6.5

		##Run the program
		mpirun -np 40 ./mfixsolver > screen.txt
		mpirun -np 1024 ./mfixsolver > screen.txt

liu.ge.wang/med/calc_dp_ufs.f

0 → 100755

+49 −0

Original line number	Diff line number	Diff line
		PROGRAM calc_DP
		!---------------------------------------------------------------------------------
		!
		! Purpose:
		! Calculate the body force and pressure drop for the gas-solid
		! sedimentation problems
		!
		!---------------------------------------------------------------------------------
		IMPLICIT NONE
		!---------------------------------------------------------------------------------
		! Local variables
		!---------------------------------------------------------------------------------
		INTEGER :: i, j
		DOUBLE PRECISION :: rho_g, rho_s, Drho, rrho, mu_g
		DOUBLE PRECISION :: phi, diap, Ar, grav, Ly, DP
		!---------------------------------------------------------------------------------
		! set variables
		!---------------------------------------------------------------------------------
		rho_s = 1.50d+3
		rho_g = 1.30d+0
		diap = 75.00d-6
		mu_g = 1.80d-5
		phi = 0.10d+0
		Ly = 23.04d-3
		grav = 98.00d+0
		!---------------------------------------------------------------------------------
		! read in
		!---------------------------------------------------------------------------------
		!---------------------------------------------------------------------------------
		! Calc other variables
		! Ar = rho_gDrhogravdiap3 / mu_g*2
		!---------------------------------------------------------------------------------
		rrho = rho_s/rho_g
		Drho = rho_s - rho_g
		Ar = rho_gDrhodiap3/mu_g2
		DP = (phirho_s + (1.0d0 - phi)rho_g)gravLy
		!--------------------------------------------------------------------------------
		! Write out
		!---------------------------------------------------------------------------------
		6 FORMAT (A9,F12.6)
		16 FORMAT (A9,F24.16)
		WRITE(,) ' '
		WRITE(*,16) ' DP = ', DP
		WRITE(,) ' '
		!
		!
		!
		END PROGRAM calc_DP

liu.ge.wang/med/mfix.dat

0 → 100755

+130 −0

Original line number	Diff line number	Diff line
		! ------------------------- RUN CONTROL SECTION -------------------------------

		RUN_NAME = 'R'
		DESCRIPTION = 'LWG Problem: phi = 0.10'
		RUN_TYPE = 'new' ! new, restart_1
		UNITS = 'SI'
		TIME = 0.0 ! start time
		TSTOP = 5.08d-01 ! stop time
		DT = 5.08d-06 ! time step
		DT_MAX = 5.08d-06 ! max time-step
		DT_MIN = 2.54d-09 ! min time-step

		MOMENTUM_X_EQ(0) = .TRUE. ! solve gas x-mom eq, granular: use .FALSE.
		MOMENTUM_Y_EQ(0) = .TRUE. ! solve gas y-mom eq, granular: use .FALSE.
		MOMENTUM_Z_EQ(0) = .TRUE. ! solve gas z-mom eq, granular: use .FALSE.

		MOMENTUM_X_EQ(1) = .TRUE.
		MOMENTUM_Y_EQ(1) = .TRUE.
		MOMENTUM_Z_EQ(1) = .TRUE.

		ENERGY_EQ = .FALSE. ! do not solve energy eq
		SPECIES_EQ(0) = .FALSE. ! do not solve gas species eq
		SPECIES_EQ(1) = .FALSE. ! do not solve solid species equations

		KT_TYPE = 'GTSH' ! kinetic theory, granular: GD_99
		FRICTION_MODEL = 'SCHAEFFER' ! Schaeffer frictional viscosity
		BLENDING_FUNCTION = 'TANH_BLEND' ! hyperbolic tangent blending


		! --------------------------- NUMERICAL PARAMETERS SECTION --------------------------

		NORM_S = 0 ! normalize continuity equation residual
		NORM_G = 0
		MAX_NIT = 100 !xx 100 ! number of total iterations allowed
		TOL_RESID = 2e-04 !xx 1e-08 ! total tolerance
		TOL_RESID_Th = 2e-04 !xx 1e-08 ! Tolerance for granular temperature
		!
		!0 =fou, 2 =superbee, 6 =MUSCL, 8 =minmod
		!
		DISCRETIZE = 9*3 ! smart
		# DEF_COR = .TRUE. ! use deferred correction method for discretization
		LEQ_IT = 9*500 ! set number of iterations in linear solver to 500 for all nine base equation types
		LEQ_TOL = 9*1D-10 ! set tolerance for all nine equation types
		! P_g phi U_k V_k W_k enrgy spcs Thta k-e
		UR_FAC = 1.0 1.0 1.0 1.0 1.0 1.0 1.0 0.8 1.0
		LEQ_PC = 9*'NONE'



		! --------------------------- Geometry Section ---------------------------------------

		COORDINATES = 'cartesian' !use cartesian coordinates
		XLENGTH = 5.76d-3 !length of domain in x direction
		YLENGTH = 23.04d-3 !
		ZLENGTH = 0.48d-3 !
		IMAX = 96 !number of cells in x direction (i)
		JMAX = 384 !
		KMAX = 8 !

		CYCLIC_X = .TRUE. !use periodic boundaries in all directions
		CYCLIC_Y_PD = .TRUE.
		DELP_Y = 341.3297663999999827d0 !pressure drop in Y-dir
		CYCLIC_Z = .TRUE.

		GRAVITY = 98.0d+0 !acts in -ve y direction


		! --------------------------- Gas-phase Section--------------------------------------

		RO_g0 = 1.3d+0 ! constant density, granular: use 0.0d0 for vaccuum
		MU_g0 = 1.8d-5 !constant gas viscosity
		P_ref = 0
		NMAX_g = 0
		NMAX_s = 0


		! ------------------------- Solids-phase Section ------------------------------------

		DRAG_TYPE = 'BVK' !use wen-yu drag correlation WEN_YU, KOCH_HILL
		RO_s0 = 1.5d+3 !solids density
		D_p0(1) = 75.0d-6 !particle diameter
		e = 0.8d0 !restitution coefficient
		e_w = 1 !particle-wall restitution (not used)
		Phi = 0 !angle of internal friction
		EP_star = 0.35644 !void fraction at minimum fluidization


		!----------------------- Initial Conditions Section ----------------------------------

		IC_X_w(1) = 0.0d0 !x boundary location west
		IC_X_e(1) = 5.76d-3 !x boundary location east
		IC_Y_s(1) = 0.0d0
		IC_Y_n(1) = 23.04d-3
		IC_Z_b(1) = 0.0d0
		IC_Z_t(1) = 0.48d-3
		! initial values in the region
		IC_EP_g(1) = 0.90d0 ! gas void fraction
		IC_U_g(1) = 0.0d0 ! x-dir gas velocity
		IC_V_g(1) = 0.0d0 ! y-dir gas velocity
		IC_W_g(1) = 0.0d0 ! z-dir gas velocity

		IC_U_s(1,1) = 0.0d0 !initial solid velocity x
		IC_V_s(1,1) = 0.0d0 !y
		IC_W_s(1,1) = 0.0d0 !z

		IC_Theta_m(1,1) = 1.0d-04 !initial granular temperature

		! --------------------------- Output Control -----------------------------------
		OUT_DT = 100 ! write text file .OUT every 10s
		RES_DT = 1e-1 ! write binary restart file .RES every 10s
		NLOG = 100000 ! write logfile R.LOG every 100 time steps
		FULL_LOG = .T. ! display residuals on screen
		CALL_USR = .T. ! call the user functions
		USR_DT = 2.54d-4

		! SPX_DT values determine how often SPx files are written (in seconds).
		! SP1 SP2 SP3 SP4 SP5 SP6 SP7 SP8 SP9 SPA SPB
		! EP_g P_g U_g U_s ROP_s T_g X_g Theta usr scalars Rxn k
		! P_star V_g V_s T_s X_s epsilon
		! W_g W_s
		! SPX_DT = 2.5e-4 1e9 2.5e-4 2.5e-4 1e9 1e9 1e9 2.5e-4 1e9 1e9 1e9
		SPX_DT = 1e9 1e9 1e9 1d9 1e9 1e9 1e9 1e9 1e9 1e9 1e9
		RESID_STRING = 'P0' 'U0' 'V0' 'W0' 'U1' 'V1' 'W1' 'G1' !'R1'

		! The decomposition in I J and K directions for a Distributed Memory Parallel machine

		NODESI = 8 NODESJ = 32 NODESK = 1

		! eor EOR EOF eof