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
Last revision Both sides next revision
solvers:fabric [2012/11/27 05:03]
gag [General Description]
solvers:fabric [2014/10/07 13:22]
gag [Discontinuous Galerkin Fabric Solver]
Line 1: Line 1:
-===== Discontinuous Galerkin Fabric Solver =====+===== Fabric Evolution - DG Fabric Solver =====
  
 ==== General Informations ==== ==== General Informations ====
Line 26: Line 26:
  
  
-  Exported Variable 1 = Temperature        !Define Temperature Mandatory!!+  Exported Variable 1 = Temperature        ! Define Temperature Mandatory!
   Exported Variable 1 DOFS = Integer 1   Exported Variable 1 DOFS = Integer 1
  
-  Exported Variable 2 = Fabric             !Define Fabric Variable !!Mandatory if Isotropic=False +  Exported Variable 2 = Fabric             ! Define Fabric Variable 
-  Exported Variable 2 DOFS = Integer 5+  Exported Variable 2 DOFS = Integer 5     ! Mandatory if Isotropic=False
  
-  Exported Variable 3 =  StrainRate        ! Compute SR 
-  Exported Variable 3 DOFS = Integer 6     !4 in 2D  6 in 3D (11,22,33,12,23,31) 
  
-  Exported Variable 4 =  DeviatoricStress  Compute Stresses +If non-linearity introduced using deviatoric stress second invariant  
-  Exported Variable 4 DOFS Integer 6     !4 in 2D  6 in 3D  (11,22,33,12,23,31)+  Procedure "ElmerIceSolvers" "AIFlowSolver_nlS2" 
 +End
  
-  Exported Variable 4 =  Spin              Compute Spin +! Fabric solver itself 
-  Exported Variable 4 DOFS Integer 3     !in 2D  3 in 3D (12,23,31)+Solver 2 
 +  Equation Fabric 
 +  Variable = -nooutput Compfab    dumy variable 
 +  Variable DOFs = 1               !FabricSolver compute each variable independently 
 +                                  !Picard Type iterations
  
-! If non-linearity introduced using deviatoric stress second invariant  +  Procedure = "ElmerIceSolvers" "FabricSolver
-  Procedure = "AIFlowSolve_nlS2" "AIFlowSolver_nlS2+  Discontinuous Galerkin = Logical True 
-   +End
-! If non-linearity introduced using strain-rate second invariant    +
-!  Procedure = "AIFlowSolve_nlD2" "AIFlowSolver_nlD2"+
 </code> </code>
  
  
 <code> <code>
-! Body Force 
-Body Force 1 
-  AIFlow Force 1 = Real 0.0 
-  AIFlow Force 1 = Real 0.0 
-  AIFlow Force 3 = Real -0.00899  ! body force, i.e. gravity * density 
-End 
- 
 ! Material ! Material
 Material 1 Material 1
-!!!!! For AIFlows... +!!!! For Fabric Solver 
-  Powerlaw Exponent = Real 3.0         ! sqrt(tr(S^2/2))^n if AIFlow_nlS2 sqrt(tr(2D^2))^(1/n-1) if  AIFlow_nlD2 +  Interaction Parameter = Real 00 => Fabric Evolution function of Strain-rates  
-  Min Second Invariant Real 1.0e-10  ! Min value for the second invariant of strain-rates +                                  ! 1 => Fabric Evolution function of dev stresses 
-  Reference Temperature = Real -10.0   T0 (Celsius)! +                                  !If not defined set to the default value given in the Viscosity File 
-  Fluidity Parameter = Real 20.        ! Bn(T0) +                                   
-  Limit Temperature = Real -5.0        ! TL  (Celsius)! +  Diffusion Parameter Real 0.   ! Diffusion term. By default set to 0 if not defined
-  Activation Energy 1 = Real 7.8e04    ! Joule/mol for T&lt;TL +
-  Activation Energy 2 = Real 7.8e04    Joule/mol for T&gt;TL +
- +
-  Viscosity File = FILE "040010010.Va" +
- +
-  Isotropic Logical False   !If set to true Glen flow law (no need to define Fabric)+
 End End
  
Line 81: Line 69:
   Fabric 4 = Real 0.               !a2_23   Fabric 4 = Real 0.               !a2_23
   Fabric 5 = Real 0.               !a2_13   Fabric 5 = Real 0.               !a2_13
- +  ...
-  AIFlow 1 = Real 0.0              ! u +
-  AIFlow 2 = Real 0.0              ! v +
-  AIFlow 3 = Real 0.0              ! w +
-  AIFlow 4 = Real 0.0              ! p+
 End End
  
Line 91: Line 75:
 Boundary Condition 1 Boundary Condition 1
   Target Boundaries = 1   Target Boundaries = 1
-!dirichlet condition for velocity +Dirichlet conditions for Fabric  
-   AIFlow 1 = Real 0.0 +only if inflow boundary condition, no condition for outflow 
-   AIFlow 2 = Real 0.0 +  Fabric 1 = Real 0.33333333333333 !a2_11 
-End +  Fabric 2 = Real 0.33333333333333 !a2_22 
- +  Fabric 3 = Real 0.               !a2_12 
- +  Fabric 4 = Real 0.               !a2_23 
-Boundary Condition 2 +  Fabric 5 = Real 0.               !a2_13
-  Target Boundaries = 2 +
-Neuman condition for AIFlow +
-  Normal force = Real 0.0    force along normal +
-  Force 1 = Real 0.0        force along x +
-  Force 2 = Real 0.0        force along y +
-  Force 3 = Real 0.0        force along z +
- +
-  AIFlow Slip Coeff 1 = Real 0.0   Slip coeff. +
 End End
 </code> </code>
Line 112: Line 87:
  
 ==== Examples ==== ==== Examples ====
-Download an example using the FlowDepth Solver. TODO+Download an example using the ''Fabric'' Solver. TODO
  
 ==== References ==== ==== References ====
solvers/fabric.txt · Last modified: 2014/10/07 13:23 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