bodyforce

Test greater flexibility for specifying body forces

Compares results generated by three available syntax inputs to an old output files known to be correct. This is stress_bodyforce in VV Test Suite

Syntax tested:

bodyforce
  
---
bodyforce force
1	0	0	0	0	-0.0245
  
---
bodyforce acceleration
1	0	0	0	0	-9.8
  

Test Directory: FEHM/fehmpytests/bodyforce

Example File bodyforce1.dat

# input.dat
# ----------------------------SOLUTION TYPE----------------------------
sol
1	-1	
# -----------------------CONTOUR OUTPUT REQUESTS-----------------------
cont
surf	1000	1e+30	time	
xyz
stress
end
# -------------------INITIAL TEMPERATURE AND PRESSURE-------------------
pres
1	0	0	5.0	20.0	1	

# ------------------------------GENERATORS------------------------------
zone
994		#ZMAX
-0.1	1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	
1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	-0.1	
0.1	0.1	0.1	0.1	-50.0	-50.0	-50.0	-50.0	

flow
-994	0	0	0.1	-30.0	1000000.0	

# -----------------------------PERMEABILITY-----------------------------
perm
1	0	0	1e-14	1e-14	1e-14	

# -------------------------MATERIAL PARAMETERS-------------------------
rock
1	0	0	2500.0	1000.0	0.1	

# --------------------------ROCK CONDUCTIVITY--------------------------
cond
1	0	0	2.5	2.5	2.5	

# -----------------------TIME STEPPING PARAMETERS-----------------------
time
1.0	10.0	1	1	

# --------------------SIMULATION CONTROL PARAMETERS--------------------
ctrl
10	1e-05	8	24	gmre	
1	0	0	2	

1	3	1.0	
7	1.5	1e-05	30.0	
0	0	
# --------------------------SOLVER PARAMETERS--------------------------
iter
1e-05	1e-05	0.001	-1e-05	1.1	
0	0	0	0	3600	
# ----------------------------STRESS MODULE----------------------------
strs
1	-3	
# ----------------------------------------------------------------------
bodyforce
initcalc
fem
# ----------------------------------------------------------------------
zone
994		#ZMAX
-0.1	1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	
1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	-0.1	
0.1	0.1	0.1	0.1	-50.0	-50.0	-50.0	-50.0	

stressboun
-994	0	0	0	-3	

zone
995		#ZMIN
-0.1	1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	
1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	-0.1	
-950.0	-950.0	-950.0	-950.0	-1000.1	-1000.1	-1000.1	-1000.1	

stressboun
-995	0	0	0	3	

zone
997		#YMIN
-0.1	1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	
0.5	0.5	-0.1	-0.1	0.5	0.5	-0.1	-0.1	
0.1	0.1	0.1	0.1	-1000.1	-1000.1	-1000.1	-1000.1	

stressboun
-997	0	0	0	2	

zone
999		#XMIN
-0.1	0.5	0.5	-0.1	-0.1	0.5	0.5	-0.1	
1.1	1.1	-0.1	-0.1	1.1	1.1	-0.1	-0.1	
0.1	0.1	0.1	0.1	-1000.1	-1000.1	-1000.1	-1000.1	

stressboun
-999	0	0	0	1	

# ----------------------------------------------------------------------
elastic
1	0	0	10000.0	0.25	

tolerance
-0.01
stressend
stop