[Data]
; Periodicity of the small scale oscillating functions
Eta= 0.001
; length of one edge of the domain for square/cubic micro problem
Delta= 0.001
; Choice of predefined magnetic permeability
Mu= 2
; Choice of predefined electric permittivity
Eps= 2
; Choice of predefined grid
Grid= 1
; Choice of predefined right-hand side
RHS= 2
; Choice of analytical solution (has to match Mu, Eps, Grid and RHS)
ExactSolution= 2
; Coefficients for the cavity solution (Ez is given by Ex and Ey)
ExEy = 1,0
; ExactSolution = 4 is the sum of the cavity solution and the polynomial
; solution weighted with the following two factors. The right-hand side is only
; PolynomialFactor times polynomial right-hand side.
CavityFactor = 1
PolynomialFactor = 1
; Coefficients for the Polynomial solution (for RHS=3 or ExactSolution=3)
; (e.g. 1, 2.5, 3 <-> x^2 + 2.5*x + 3)
Coefficients = 1,0,0,0,1,0,0,0,1,0,0,0,1
[Solving]
HMM_Solution= true
HOM_Solution= false
ODE_Solution= false
EX_Solution= true
[Space_Discretisation]
; Determine the macro refinement levels that should be run through
; (seperate with ',')
Macro_Refinement= 4,8,16,32,64
; Determine the micro refinement levels that should be run through
; (seperate with ',')
Micro_Refinement= 2,4,8,16,32
; Determine the polynomial space of Nedelec's elements for the macro
; discretization (seperate with ',')
Macro_Pol_Deg = 1
; Determine the polynomial space or Lagrangian elements for the micro
; discretization (seperate with ',')
Micro_Pol_Deg = 1
[Assembling]
; Reuse local mass matrix in HMM scheme (does only work properly if the
; material parameters Mu and Eps are periodic and if the edge length of the
; macro mesh cells are multiples of Delta
Reuse_Local= true
; Evaluate slow variable only at macro quadrature point
Collocation= true
; Take periodic (Eta==Delta) or homogeneous Dirichlet boundary conditions for
; micro problems.
Periodic_BC= true
; The mean value constraint assures the
; mean of the solution over the whole micro domain to be
; zero. This yields the existence of a unique
; solution. However, this means that all dofs are
; coupled and therefore the sparsity pattern
; created by DoFTools::make_sparsity_pattern is
; huge. Numerical experiments have shown, that
; for periodic boundary conditions with eta=delta,
; the CG solver also converges without this
; constraint. As we are only interested in the
; gradient of the solution, knowing the solution
; up to a constant is enough.
Mean_Value_Constraint = false
[Time_Integration]
; Time integrator that will be used (Leapfrog, ImplicitMidpoint, CrankNicolson)
Time_Integrator= Leapfrog
; Final time of the integration
End_Time= 1
; Number of timesteps, that will be used for the time integration
; (seperate with ',')
Timesteps= 100
; Time integrator used for the reference solution (Leapfrog, ImplicitMidpoint, CrankNicolson)
ODE_Time_Integrator = Leapfrog
; Number of timesteps for reference solution
ODE_Timesteps = 1000
[Output]
; Output of the solutions (*.vtu) every Steps timesteps
; (Steps= 1 for output of all timesteps)
Solution= false
Steps= 50
; Compute error between HMM and effective solution
Error_HMM_HOM= true
; Compute error between HMM and analytical solution
Error_HMM_ANA= true
; Compute error between effective and analytical solution
Error_HOM_ANA= true
; Compute error between HMM and reference solution
Error_HMM_ODE= false
; Compute error between effective and reference solution
Error_HOM_ODE= false
; Compute error between reference and effective solution
Error_ODE_ANA= false
; Write errors to console for every timestep
Error_Console= false
[Testing]
; For Sweep_Delta= false, all computations are done with Delta given above as
; edge length of the sampling domain.
; Sweep_Delta=true allows for variation of the size of the sampling domain
; corresponding to Sweep_Div and Sweep_Range. The values used for Delta
; are Delta * (1 + iter / Sweep_Div), for iter from 1 to Sweep_Range. As the
; Micro_Refinement is varied in the same way, such that the mesh width of
; the local problems is invariant, Sweep_Div should be a divisor of
; Micro_Refinement.
Sweep_Delta = false
Sweep_Div = 10
Sweep_Range = 7