Software Required
========================
- foam-extend-4.1


Features
========================
This submission includes two solvers i) biotFoam and ii) pressureBuildupFoam, as well as examples  to validate each solvers.
In addition, two boundaries conditions are included: i) propagatingPressureWave and ii) tractionDisplacement

biotFoam
~~~~~~~~

biotFoam is a freely available OpenFoam package developed by Roenby [1].  It is transient segregated finite-volume solver of the Biot equations for linear-elastic, small-strain deformation of a solid skeleton, coupled with pore water flow and pressure governed by Darcy's law.

pressureBuildupFoam
~~~~~~~~~~~~~~~~~~~
pressureBuildupFoam is an incremental development on biotFoam. It adds an additional equation to solve for the pore pressure buildup. Also, the liquefaction criteria is implemented in this solver to predict the onset of liquefaction. 


propagaringPressureWave
~~~~~~~~~~~~~~~~~~~~~~~
propagaringPressureWave is a boundary condition from the biotFoam package to simulate the effect of waves on the seabed by means of the pressure due to the progressive waves acting at the mudline.

tractionDisplacement
~~~~~~~~~~~~~~~~~~~~

tractionDisplacemnt is also a boundary condition from biotFoam package, which contains a fixed displacement gradient boundary due to fixed traction, which is a very common boundary condition for stress analysis problem.


File structure
==============

Structure
~~~~~~~~~
	-src
		-boundaries
		        -propagatingPressureWave
		        -tractionDsiplacemnet   
		       
		-solvers
		        -biotFoam
		        -pressureBuildupFoam   
	-run
		-biotFoam
			-Hsu_and_Jeng
               -pressureBuilupFoam
			-SummerExample_b
			-SummerExample_p
			-stoneColumn_b
			-stoneColumn_p



Results
========    

bioFoam - Hsu_and_Jeng  case
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
	- Allrun script is used to run the case
	- Comparison.m file produces the comaparison of the numerical results of the pore pressure and the shear stress against the analytical solutions   
	
pressureBuilupFoam - cases
~~~~~~~~~~~~~~~~~~~~~~~~~~
	- In pressureBuildupFoam directory, the subscript 'b' indicates that biotFoam has to be solved and 'p' indicates pressureBuildupFoam has to be solved.
        - There are scripts runScript_stoneColumn.m and runScript_SumerExample.m which runs the two different example under pressureBuildupFoam
        - Note that for the examples in pressureBuildupFoam directory, biotFoam need to be solved first to calculate shear stress and later pressureBuildupFoam to calculate accumulated pore pressure and residual liquefaction 
	- post_criteria.m file in Sumer_example_p directory produces the numerical results of the accumulated pore pressure comaparing
	 against the analytical solutions 
	- Comparison.m file in stoneColumn_p directory delivers the comaprison of the numerical results of accumulated pore pressure of Sumer_Example case and stoneColumn cases.
	- The script plot_p_tau.m delivers the comparison of the numerical results of  pore pressure and the shear stress for the SumerExample and the stoneColumn case
	


[1] J. Roenby, Openfoam library for soil mechanics, https://github.com/roenby/soilFoam, 2013.

