adv-diff-sn-e3v-case.sif

Header
  CHECK KEYWORDS Warn
  Mesh DB "." "Mesh_1"
  !Mesh DB "." "Mesh_2"
  Include Path ""
  Results Directory ""
End

Simulation
  Max Output Level = 4
  Coordinate System = Cartesian
  Coordinate Mapping(3) = 1 2 3
  Coordinate Scaling = 1.0e-4
  Simulation Type = Transient
  Steady State Max Iterations = 10
  Output Intervals = 1
  Timestepping Method = BDF
  BDF Order = 1
  Timestep intervals = 3000
  Timestep Sizes = .005
  Solver Input File = adv-diff-sn-e3v-case.sif
  Post File = case.vtu
  Use Mesh Names = True
End

Constants
  Gravity(4) = 0 -1 0 9.82
  Stefan Boltzmann = 5.67e-08
  Permittivity of Vacuum = 8.8542e-12
  Boltzmann Constant = 1.3807e-23
  Unit Charge = 1.602e-19
End

Body 1
  !Target Bodies(1) = 1
  Name = "electrolyte"
  Equation = 1
  Material = 1
  Body Force = 1
  Initial condition = 1
End

Solver 5
  Equation = SaveScalars
  Filename = case.dat
  Coefficient 1 = Concentration Diffusivity
  Procedure = "SaveData" "SaveScalars"
  Operator 1 = diffusive flux
  Variable 1 = Concentration
  Variable 2 = Concentration
  Operator 2 = min
  Operator 3 = max
  Operator 4 = mean
  Exec Solver = String After Timestep
End

Solver 1
  Equation = Static Current Conduction
  Procedure = "StatCurrentSolve" "StatCurrentSolver"
  Variable = -dofs 1 Potential
  Exec Solver = Always
  Stabilize = True
  Bubbles = False
  Lumped Mass Matrix = False
  Optimize Bandwidth = True
  Steady State Convergence Tolerance = 1.0e-5
  Nonlinear System Convergence Tolerance = 1.0e-8
  Nonlinear System Max Iterations = 20
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Newton After Tolerance = 1.0e-3
  Nonlinear System Relaxation Factor = 1
  Linear System Solver = Direct
  Linear System Direct Method = Umfpack
End

Solver 4
  Equation = Navier-Stokes
  Procedure = "FlowSolve" "FlowSolver"
  Variable = Flow Solution[Velocity:2 Pressure:1]
  Exec Solver = Always
  Stabilize = True
  Bubbles = False
  Lumped Mass Matrix = False
  Optimize Bandwidth = True
  Steady State Convergence Tolerance = 1.0e-5
  Nonlinear System Convergence Tolerance = 1.0e-8
  Nonlinear System Max Iterations = 20
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Newton After Tolerance = 1.0e-3
  Nonlinear System Relaxation Factor = 1
  Linear System Solver = Iterative
  Linear System Iterative Method = BiCGStab
  Linear System Max Iterations = 500
  Linear System Convergence Tolerance = 1.0e-8
  Linear System Preconditioning = ILU0
  Linear System ILUT Tolerance = 1.0e-3
  Linear System Abort Not Converged = False
  Linear System Residual Output = 1
  Linear System Precondition Recompute = 1
End

Solver 2
  Equation = Advection Diffusion Equation
  Variable = -dofs 1 Concentration
  Procedure = "AdvectionDiffusion" "AdvectionDiffusionSolver"
  Exec Solver = Always
  Stabilize = True
  Bubbles = False
  Lumped Mass Matrix = False
  Optimize Bandwidth = True
  Steady State Convergence Tolerance = 1.0e-5
  Nonlinear System Convergence Tolerance = 1.0e-8
  Nonlinear System Max Iterations = 20
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Newton After Tolerance = 1.0e-3
  Nonlinear System Relaxation Factor = 1
  Linear System Solver = Direct
  Linear System Direct Method = Umfpack
End

Solver 3
  Equation = Heat Equation
  Variable = -dofs 1 Temperature
  Procedure = "HeatSolve" "HeatSolver"
  Exec Solver = Always
  Stabilize = True
  Bubbles = False
  Lumped Mass Matrix = False
  Optimize Bandwidth = True
  Steady State Convergence Tolerance = 1.0e-5
  Nonlinear System Convergence Tolerance = 1.0e-8
  Nonlinear System Max Iterations = 20
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Newton After Tolerance = 1.0e-3
  Nonlinear System Relaxation Factor = 1
  Linear System Solver = Iterative
  Linear System Iterative Method = BiCGStab
  Linear System Max Iterations = 500
  Linear System Convergence Tolerance = 1.0e-8
  Linear System Preconditioning = ILU0
  Linear System ILUT Tolerance = 1.0e-3
  Linear System Abort Not Converged = False
  Linear System Residual Output = 1
  Linear System Precondition Recompute = 1
End

Equation 1
  Name = "Electrolyte_Calculation"
  Convection = Computed
  Concentration Units = Absolute Mass
  Concentration Convection = Computed
  Active Solvers(5) = 5 1 4 2 3
End

Material 1
  Name = "Sn solution"
  Concentration Diffusivity = 9.8e-11 !5.085e-10
  Reference Temperature = 298.15
  Viscosity = .001
  Heat expansion Coefficient = 1e-4
  Electric Conductivity = 0.024
  Heat Conductivity = .6
  Heat Capacity = 4184
  Density = 1000
  Concentration Ion Charge = 2
End

Body Force 1
  Name = "Natural convection"
  Boussinesq = True
  Joule Heat = True
End

Initial Condition 1
  Name = "Initial concentration in mol/L" !SI unit is kg/m^3
  Concentration = 0.0
  Temperature = 298
End

Boundary Condition 1
  !Target Boundaries(1) = 4 !Anode
  Name = "left"
  Potential = 3.0
  !Current Density BC = Logical True
  !Current Density = Real 31.714 !150 A/m^2 Determination from Butler Volmer kinetics  
  ! Ekilik2000, Amorim2007, Steigerwald1962, He2014, Potentiostat measurements
  Temperature = 298 !303
  Concentration =  7.5e-6 !mol/L He2014 2.113e-3 kg/m^3 measured from ion chromatography (not used)
  Noslip wall BC = True
  Save Scalars = True
End

Boundary Condition 2
  !Target Boundaries(1) = 1
  Name = "bottom"
  Noslip wall BC = True
End

Boundary Condition 3
  !Target Boundaries(1) = 3
  Name = "top"
  Velocity 2 = 0
End

Boundary Condition 4
  !Target Boundaries(1) = 2 !Cathode
  Name = "right"
  Potential = 0
  !Current Density BC = Logical True
  !Current Density = Real -0.7816 !A/m^2 Determination from Butler Volmer kinetics 
  Temperature = 298 !293
  Concentration = 0.0
  Noslip wall BC = True
  Save Scalars = True
End