richards

    Richards equation test

    Comparison of Richards equation with 2-phase solution Node by node comparison.

    Test Directory: FEHM/fehmpytests/richards

    Example File rich_rlp.dat

    
    NER_DATA March 2008 
    # 11 unsaturated units
    #01_Qal
    #02_Qbt_1g Tshirege
    #03_Qct
    #04_Qbo
    #05_Qbog
    #06_Tpf (upper) 
    #07_Tb4 tholeitic
    #08_Tb4 alkalic
    #09_Tb4 tephra
    #10_Ta
    #11_Tpf (lower)
    #
    #Air water statement
    airwater
    3 
    20.0 0.10
    #end air water statement
    # 100 meter by 322m vadose zone (2-D)
    zone
    1  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    0. 0. 0. 0. -10 -10 -10 -10
    2  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
     -10 -10 -10 -10  -45. -45. -45. -45.
    3  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
     -45. -45. -45. -45. -55. -55. -55. -55.
    4  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
     -55. -55. -55. -55. -115 -115. -115. -115.
    5  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -105 -105. -105. -105. -125. -125. -125. -125.
    6  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -115. -115. -115. -115. -145. -145. -145. -145. 
    7  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -135. -135. -135. -135. -195. -195. -195. -195. 
    8  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -185. -185. -185. -185. -225. -225. -225. -225. 
    9  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -215. -215. -215. -215. -227. -227. -227. -227.
    10  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -217. -217. -217. -217. -237. -237. -237. -237. 
    11  
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -227. -227. -227. -227. -322. -322. -322. -322. 
    
    #begin finite difference grid statement
    fdm
    block
    25 1 322
    0. 0. 0.
    1 4.
    -25 4.
       
    1 1.
    
    1 1.
    -322 1.
    
    #end finite difference grid statement
    #
    #
    #
    #
    #begin intrinsic permeability statement
    perm
    1 0 0 -12 -12 -12
    -1  0 0 4.62E-13      4.62E-13    4.62E-13   #kay ref Apendix A
    -2	0	0	3.68e-13	3.68e-13	3.68e-13
    -3  	0	0	8.82e-13	8.82e-13	8.82e-13
    -4  	0	0	7.25e-13	7.25e-13	7.25e-13
    -5	0	0	1.53e-13	1.53e-13	1.53e-13
    -6  	0	0	4.73e-12	4.73e-12	4.73e-12
    -7  	0	0	2.96e-13	2.96e-13	2.96e-13  #7,8,9 one unit in LANL model
    -8 	0	0	2.96e-13	2.96e-13	2.96e-13
    -9  	0	0	2.96e-13	2.96e-13	2.96e-13
    -10     0       0       1.1489E-13	1.1489E-13      4.4865E-15   # leave as NER values
    -11 	0	0	4.73e-12	4.73e-12	4.73e-12
    
    #end intrinsic permeability statement
    #
    #specify rock or sand properties
    rock  # porosoty models, use  as NER values
    1 0 0 1600 1000 0.2
    -1 0 0 1600 1000 0.43
    -2 0 0 1600 1000 0.48
    -3 0 0 1600 1000 0.53   #placeholder (where did porosity come from?)
    -4 0 0 1600 1000 0.455
    -5 0 0 1600 1000 0.5
    -6 0 0 1600 1000 0.207
    -7 0 0 1600 1000 0.1
    -8 0 0 1600 1000 0.086
    -9 0 0 1600 1000 0.274
    -10 0 0 1600 1000 0.227
    -11 0 0 1600 1000 0.1
    
    #end rock properties
    #
    #
    #start relative permeability statment and capillary pressure
    rlp
    3       0.087   1.      3.85    1.558   5       0.089  Estimate from Kay
    3	0.018	1.	2.22	1.592	2.	0.019	Material 14
    3	0.01	1.	1.52	1.506	2.	0.011	Material 12
    3	0.026	1.	0.66	1.711	2.	0.027	Material 11
    3	0.01	1.	0.081	4.026	2.	0.011	Material 10
    3	0.01	1.	5.	2.68	2.	0.011	Material 6
    3	0.03	1.	5.	1.5	2.	0.040	Material 7 
    3	0.03	1.	5.	1.5	2.	0.040	Material 7 
    3	0.03	1.	5.	1.5	2.	0.040	Material 7
    3 0.001   1. 0.248418  1.689 5 0.002    #left 10 with NER data
    3	0.01	1.	5.	2.68	2.	0.011	Material 6
    
     -1  0 0  1 
     -2  0 0  2 
     -3  0 0  3 
     -4  0 0  4 
     -5  0 0  5 
     -6  0 0  6 
     -7  0 0  7 
     -8  0 0  8 
     -9  0 0  9 
     -10 0 0  10
     -11 0 0  11
    
    #end of intrinsic permeability statment
    #
    #input inital conditions
    #end input initial conditions
    #
    zone   #define top boundary (full boundary)
    21
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -1. -1. -1. -1. 0. 0. 0. 0.
    22    #define top boundary (subset for recharge)
    46. 54. 54. 46. 46. 54. 54. 46.
    0. 0. 1. 1. 0. 0. 1. 1.
    -1. -1. -1. -1. 0. 0. 0. 0.
    31    #define bottom boundary (for disharge)
    0. 100. 100. 0. 0. 100. 100. 0. 
    0. 0. 1. 1. 0. 0. 1. 1.
    -321. -321. -321. -321. -322. -322. -322. -322. 
    
    #
    #input boundary conditions 
    #modify for positive head
    boun
    model 1
    ti
    2 0. 1.e20
    dsw
    -1.60e-3 -1.60e-3
    
    -22 0 0 1
    
    pres
    1 0 0 0.1 .2 2
    -31 0 0 0.1 0.9 2
    
    flow #air pressure = 0.1 on top boundary
    -21 0 0 0.1 -0.1 1.e-12
    
    flow
    -31 0 0 0.1  1 1.e0 
    
    #end boundary conditions
    #
    #solver control options#
    ctrl
    10 1E-10 20 150 gmres
    1 0 0 2
    
    1 3 1 
    10 1.5 1e-8 3650.E6
    0 0
    #
    #end solver control options
    #
    #
    #begin iter statement
    iter
    1.e-10 1.e-10 0.1 -1.e-3 1.2
    0 0 0 20 1000
    #end iter statement
    #
    #
    #time statement
    time
    1e-1 9131.25 100 100 08 01 0
    
    #end time statement
    #
    #node statement
    #see page 8100 for explanation
    node
    bloc
    13 8050 25 
    
    #end node statement
    hist
    tecplot
    mpa
    saturation
    end
    #begin contour statement
    cont
    avs 10000  1.e20
    l 
    p
    s
    endavs
    #end contour statement
    stop
    zone
    50
    -1.e5 1.e5 1.e5 -1.e5 -1.e5 1.e5 1.e5 -1.e5
    -1.e5 -1.e5 1.e5 1.e5 -1.e5 -1.e5 1.e5 1.e5
    -1.e5 -1.e5 -1.e5 -1.e5
    1.e5 1.e5 1.e5 1.e5 
    
    zone
    22    #define top boundary (subset for recharge)
    46. 54. 54. 46. 46. 54. 54. 46.
    0. 0. 1. 1. 0. 0. 1. 1.
    -1. -1. -1. -1. 0. 0. 0. 0.
    
    trac
    0 1 1e-5 1
    1. 1.e8 1.e8 1.e8
    50 1.2 0.1 500.
    1
    1
    1 0 0 1 1e-15 1.0 1.0 1.0
    
    1 0 0 1
    
    1 0 0 0
    
    -22 0 0  -1.0 1. 1.e20
    
    #end tracer statement
    stop