Merge pull request #288 from loganoz/newmobility

update mobility definition

Author: AbbBallout

Solver/src/libs/physics/cahnhilliard/PhysicsStorage_CH.f90

@@ -7,7 +7,7 @@ Module Physics_CHKeywordsModule
 !        Required arguments
 !        ******************
 !
-         character(len=KEYWORD_LENGTH), parameter    :: TCH_KEY               = "chemical characteristic time (s)"
+         character(len=KEYWORD_LENGTH), parameter    :: TCH_KEY               = "chem characteristic time (s)"
          character(len=KEYWORD_LENGTH), parameter    :: INTERFACE_WIDTH_KEY   = "interface width (m)"
          character(len=KEYWORD_LENGTH), parameter    :: INTERFACE_TENSION_KEY = "interface tension (n/m)"
          CHARACTER(LEN=KEYWORD_LENGTH), DIMENSION(3) :: physics_CHKeywords = [INTERFACE_WIDTH_KEY, &
@@ -84,7 +84,7 @@ SUBROUTINE ConstructPhysicsStorage_CH( controlVariables, Lref, timeRef, pRef, su
       multiphase_ % sigma_wDim = controlVariables % DoublePrecisionValueForKey(INTERFACE_TENSION_KEY)
 
       if ( .not. almostEqual(multiphase_ % tCH_wDim, 0.0_RP) ) then
-         multiphase_ % M0_wDim = POW2(Lref)*multiphase_ % eps_wDim / (multiphase_ % tCH_wDim * multiphase_ % sigma_wDim)
+         multiphase_ % M0_wDim = POW2(Lref)*(multiphase_ % eps_wDim**2) / (multiphase_ % tCH_wDim * multiphase_ % sigma_wDim)
       else
          multiphase_ % M0_wDim = 0.0_RP
       end if
@@ -95,7 +95,7 @@ SUBROUTINE ConstructPhysicsStorage_CH( controlVariables, Lref, timeRef, pRef, su
       multiphase_ % sigma = multiphase_ % sigma_wDim / pRef
 
       if ( .not. almostEqual(multiphase_ % tCH, 0.0_RP) ) then
-         multiphase_ % M0 = multiphase_ % eps / (multiphase_ % tCH * multiphase_ % sigma)
+         multiphase_ % M0 = (multiphase_ % eps**2) / (multiphase_ % tCH * multiphase_ % sigma)
       else
          multiphase_ % M0 = 0.0_RP
       end if

Solver/test/CahnHilliard/TJokisaari/TJokisaari.control

@@ -22,7 +22,7 @@ time integration = imex
 reference length (m) = 1.0d0
 interface width (m) = 3.1623d0
 interface tension (N/m) = 1.0d0
-chemical characteristic time (s) = 0.31623
+chem characteristic time (s) = 1.000014129
 dt                    = 1.0d-3
 simulation type       = time-accurate
 final time            = 20.0d3

Solver/test/Multiphase/ActuatorLineInterpolation/Actuator.control

@@ -21,7 +21,7 @@ fluid 2 sound speed square (m/s) = 116136.222222222
 fluid 2 density (kg/m^3) = 1.22145354210477
 fluid 2 viscosity (Pa.s) = 1.831943887d-5
 
-chemical characteristic time (s) = 1000000
+chem characteristic time (s) = 40000
 interface width (m) = 0.04
 interface tension (N/m) = 72.8e-3
 

Solver/test/Multiphase/Convergence/Convergence.control

@@ -22,7 +22,7 @@
 
 interface width (m) = 0.70711
 interface tension (N/m) = 6.236E-3
-chemical characteristic time (s) = 1000.0
+chem characteristic time (s) = 707.11
 artificial sound speed square (m/s) = 1000.0
 !------------------------- Discretization:-
                         Polynomial order = 5

Solver/test/Multiphase/Convergence_cs/Convergence_cs.control

@@ -21,7 +21,7 @@
 
 interface width (m) = 0.70711
 interface tension (N/m) = 6.236E-3
-chemical characteristic time (s) = 1000.0
+chem characteristic time (s) = 707.11      !1000
 !artificial sound speed square (m/s) = 1000.0
 fluid 1 sound speed square (m/s) = 1000.0   
 fluid 2 sound speed square (m/s) = 2000.0

Solver/test/Multiphase/Convergence_cs/SETUP/ProblemFile.f90

@@ -559,7 +559,7 @@ SUBROUTINE UserDefinedFinalize(mesh, time, iter, maxResidual &
                                           7.4217194511428453E-17_RP, & 
                                           6.1637135702849238E-04_RP]
             integer        :: i, j,k, eID 
-            CHARACTER(LEN=29)                  :: testName           = "Multiphase convergence non constant sound speed"
+            CHARACTER(LEN=49)                  :: testName           = "Multiphase convergence non constant sound speed"
             real(kind=RP)  :: error(5)
             TYPE(FTAssertionsManager), POINTER :: sharedManager
             LOGICAL                            :: success

Solver/test/Multiphase/Convergence_cs/convergence_cs.control

@@ -19,7 +19,7 @@
                 fluid 1 viscosity (Pa.s) = 0.005
                 fluid 2 viscosity (Pa.s) = 0.01
 ! MAINTAINED TO COMPARE!!        mobility = 0.18856
-        chemical characteristic time (s) = 900.0
+        chem characteristic time (s) = 38.16
                      interface width (m) = 0.0424
                  Interface tension (N/m) = 0.00025
        artificial compressibility factor = 1000.0

Solver/test/Multiphase/EntropyConservingTest/EntropyConservingTest.control

@@ -50,7 +50,7 @@ Averaging = Skew-symmetric 2
 Viscous discretization = BR1   ! IP/BR2
 Cahn-Hilliard discretization = BR1
 
-chemical characteristic time (s) = 1e+16  !10000
+chem characteristic time (s) = 1e+14  !10000
 interface width (m) = 0.01 
 interface tension (N/m) = 1e-16   !72.8e-3
 

Solver/test/Multiphase/MixedRK/mixedRK.control

@@ -44,7 +44,7 @@ Lambda stabilization = 1.0
 Viscous discretization = BR1 
 Cahn-Hilliard discretization = BR1
 
-chemical characteristic time (s) = 400 !400 !1.0d12 !1000.0  !50000
+chem characteristic time (s) = 100 !1.0d12 !1000.0  !50000
 interface width (m) = 0.25 
 interface tension (N/m) = 0.0728 !1.0d-12 !0.01