Warning: Declaration of syntax_plugin_mathpublish::handle($match, $state, $pos, &$handler) should be compatible with DokuWiki_Syntax_Plugin::handle($match, $state, $pos, Doku_Handler $handler) in /home/np29546/public_html/elmerice/wiki/lib/plugins/mathpublish/syntax.php on line 29

Warning: Declaration of syntax_plugin_mathpublish::render($mode, &$R, $data) should be compatible with DokuWiki_Syntax_Plugin::render($format, Doku_Renderer $renderer, $data) in /home/np29546/public_html/elmerice/wiki/lib/plugins/mathpublish/syntax.php on line 29

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/lib/plugins/mathpublish/syntax.php:29) in /home/np29546/public_html/elmerice/wiki/inc/common.php on line 2045

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/lib/plugins/mathpublish/syntax.php:29) in /home/np29546/public_html/elmerice/wiki/inc/actions.php on line 38
userfunctions:coulomb [Elmer/Ice Wiki]

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
userfunctions:coulomb [2012/11/15 04:45]
gag [SIF contents]
userfunctions:coulomb [2015/12/02 12:24]
gag [General Description]
Line 3: Line 3:
  
 ==== General Informations ==== ==== General Informations ====
-  * **USF Fortran File:​** ​//USF_Sliding.f90// +  * **USF Fortran File:​** ​''​USF_Sliding.f90''​ 
-  * **USF Name:​** ​//Friction_Coulomb// +  * **USF Name:​** ​''​Friction_Coulomb''​ 
-  * **Required Input Variable(s):​** A //Flow Solution// in //Flow Solution Name////Normal Vector////Stress//+  * **Required Input Variable(s):​** A ''​Flow Solution'' ​in ''​Flow Solution Name''​''​Normal Vector''​''​Stress''​ or the ''​Effective Pressure''​ variable.
  
  
 ==== General Description ==== ==== General Description ====
-The file //USF_Sliding.f90// contains ​two user functions to apply non-linear friction at the base of glacier.  ​+The file ''​USF_Sliding.f90'' ​contains ​three user functions to apply non-linear friction at the base of glacier.  ​
  
-The first user function (//Sliding_Weertman//) is a non-linear Weertman-type friction law and is described [[:​userfunctions:​weertman|here]]. The second user function (//Friction_Coulomb//) is a non-linear water pressure dependant friction law, as proposed by Schoof (2005) and Gagliardini et al. (2007), and is presented in this page.+The first user function (''​Sliding_Weertman''​) is a non-linear Weertman-type friction law and is described [[:​userfunctions:​weertman|here]]. The second user function (''​Friction_Coulomb''​) is a non-linear water pressure dependant friction law, as proposed by Schoof (2005) and Gagliardini et al. (2007), and is presented in this page. The third user function (''​Sliding_Budd''​) is described [[:​userfunctions:​budd|here]] and is from Budd et al 1984 (Annals of Glaciology 5, page 29-36)
   ​   ​
  
Line 23: Line 23:
 The //Slip Coefficient//​ in Elmer is then given as\\ The //Slip Coefficient//​ in Elmer is then given as\\
 <​m> ​ C.N {[{ {chi . {u_b}^{-n} }/ {(1 + a . chi^q)} }]}^{1/​n} ​ </​m>​\\ <​m> ​ C.N {[{ {chi . {u_b}^{-n} }/ {(1 + a . chi^q)} }]}^{1/​n} ​ </​m>​\\
-When  <​m>​u_b < u_{t0}</​m>,​ <​m>​u_b</​m>​ in the previous is replaced by <​m>​u_{t0}</​m>​.+When  <​m>​u_b < u_{t0}</​m>,​ <​m>​u_b</​m>​ in the previous ​equation ​is replaced by <​m>​u_{t0}</​m>​.
  
  
 The parameters to be given are:\\ The parameters to be given are:\\
-  * Friction Law Sliding Coefficient ​     -> <​m>​A_s</​m>​ +  * ''​Friction Law Sliding Coefficient'' ​     -> <​m>​A_s</​m>​ 
-  * Friction Law Post-Peak Exponent ​      -> <m>q >= 1</​m>​ +  * ''​Friction Law Post-Peak Exponent'' ​      -> <m>q >= 1</​m>​ 
-  * Friction Law Maximum Value            -> <​m>​C</​m>​ ~ max bed slope +  * ''​Friction Law Maximum Value'' ​           -> <​m>​C</​m>​ ~ max bed slope 
-  * Friction Law Exponent ​       -> m = (n Glen's law) +  * ''​Friction Law Exponent'' ​      -> m = (n Glen's law) 
-  * Friction Law Linear Velocity -> <​m>​u_{t0}</​m>  ​+  * ''​Friction Law Linear Velocity'' ​-> <​m>​u_{t0}</​m>  ​
  
-The effective pressure <m>N = -sigma_{nn} -p_w</​m>​ where <​m>​sigma_{nn}</​m>​ is the normal Cauchy stress and <​m>​p_w</​m>​ the water pressure. ​Hereat time t+dt, the normal Cauchy stress is estimated from the stress computed at time t. The water pressure is prescribed as an //External Pressure// (Positive ​- Compressive convention). ​+The effective pressure ​is defined as <m>N = -sigma_{nn} -p_w</​m>​,  ​where <​m>​sigma_{nn}</​m>​ is the normal Cauchy stress and <​m>​p_w</​m>​ the water pressure. ​If a variable ''​Effective Pressure''​ existsit is used to evaluate directly <​m>​N</​m>​. Else, the normal Cauchy stress is estimated from the stress computed at previous timestep. The water pressure is prescribed as an ''​External Pressure'' ​(Negative ​- Compressive convention, and therefore '​External Pressure'​ should be equal to the opposite of the water pressure in the sif). 
       ​       ​
    
Line 47: Line 47:
   Flow Force BC = Logical True   Flow Force BC = Logical True
   ​   ​
-  !! Water pressure given through ​this parameter (Positive ​= Compressive)+  !! Water pressure given through ​the Stokes '​External Pressure' ​parameter ​ 
 +  !! (Negative ​= Compressive)
   External Pressure = Equals Water Pressure   External Pressure = Equals Water Pressure
        
Line 53: Line 54:
   ​   ​
   Slip Coefficient 2 =  Variable Coordinate 1   Slip Coefficient 2 =  Variable Coordinate 1
-    Real Procedure "./​USF_Sliding" "​Friction_Coulomb"​+    Real Procedure "ElmerIceUSF" "​Friction_Coulomb"​
   Slip Coefficient 3 =  Variable Coordinate 1   Slip Coefficient 3 =  Variable Coordinate 1
-    Real Procedure "./​USF_Sliding" "​Friction_Coulomb"​+    Real Procedure "ElmerIceUSF" "​Friction_Coulomb"​
     ​     ​
   !! Parameters needed for the Coulomb Friction Law   !! Parameters needed for the Coulomb Friction Law
Line 68: Line 69:
  
 ==== Examples ==== ==== Examples ====
-An example ​of the usage of the user function ​//Friction_Weertman// TODO +The Coulomb friction law is tested in ''​[ELMER_TRUNK]/​elmerice/​Tests/​Friction_Coulomb''​ with a direct input of the effective pressure and ''​[ELMER_TRUNK]/elmerice/Tests/Friction_Coulomb_Pw''​ with the effective pressure computed from the stress and a prescribed water pressure.  ​
  
  
userfunctions/coulomb.txt · Last modified: 2015/12/02 12:24 by gag
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0