inputs/cellsolver.parset

+
+
+#path for logfile
+#outputPath = "../../outputs/output.txt"
+#outputPath = "/home/klaus/Desktop/DUNE/dune-microstructure/outputs/output.txt"
+#outputPath = "/home/klaus/Desktop/DUNE/dune-microstructure/outputs"
+
+#############################################
+#  Debug Output
+#############################################
+
+#print_debug = true    #default = false
+
+
+
+#############################################
+#  Cell Domain
+#############################################
+# Domain 1: (-1/2, 1/2)^3  , Domain 2 : [0,1)^2 x (-1/2, 1/2)
+
+cellDomain = 1
+
+#############################################
+#  Grid parameters
+#############################################
+
+
+#######################################################################
+## numLevels : Number of Levels on which solution is computed. starting with a 2x2x2 cube mesh.
+## {start,finish} computes on all grid from 2^(start) to 2^finish refinement
+########################################################################
+
+#numLevels =  1 3   # computes all levels from first to second entry
+numLevels = 3 3     # computes all levels from first to second entry
+#numLevels = 1 6
+
+
+
+
+
+#Elements_Cell = 20 20 20	               # number elements in each direction (y1 y2 x3) 
+#nElements_Cell = 30 30 30	
+#nElements_Cell = 30 30 30	
+#nElements_Cell = 50 50 50	
+#nElements_Cell = 100 100 2
+#nElements_Cell = 100 100 100  // does not work
+#nElements_Cell = 10 10 10
+#nElements_Cell = 2 2 2
+#nElements_Cell = 4 4 4
+#nElements_Cell = 8 8 8
+#nElements_Cell = 16 16 16
+#nElements_Cell = 32 32 32
+#nElements_Cell = 64 64 64 
+
+
+#gamma=50.0
+#gamma=1.0
+gamma=0.01
+
+#############################################
+#  Material parameters
+#############################################
+
+write_materialFunctions = true   # VTK mu-functions , lambda-functions
+
+beta = 2.0    # ratio between material parameters mu1 & mu2 .... beta = 1.0 corresponds to homogeneous case 
+
+mu1 = 10.0
+lambda1 = 0.0 
+#lambda1 = 20.0 
+#lambda1 = 20.0 
+#lambda1 = 5.0 
+
+####  material_implementation("analytical_Example") ? 
+
+
+material_implementation = "analytical_Example"
+#material_implementation ="isotropic_bilayer"
+
+#material_implementation ="matrix_material_circles"
+#material_implementation ="matrix_material_squares"
+
+#material_implementation = "circle_fiber"    #TEST
+#material_implementation = "square_fiber"    #TEST
+
+#############################################
+#  Prestrain parameters
+#############################################
+
+write_prestrainFunctions = true  # VTK norm of B ,
+
+rho1 = 1.0
+alpha = 5.0    # ratio between prestrain parameters rho1 & rho2 
+
+#theta = 0.3    # volume fraction   #default = 1.0/4.0
+#theta = 0.25   # volume fraction
+#theta = 0.75   # volume fraction
+
+prestrainType = "analytical_Example"
+#prestrainType = "isotropic_bilayer"
+
+------------Matrix Material --------------
+
+#prestrainType = "matrix_material_circles"
+#prestrainType = "matrix_material_squares"
+
+nF = 8   #number of Fibers (in each Layer)
+#rF = 0.05  #Fiber radius    max-fiber-radius =  (width/(2.0*nF)
+------------------------------------------
+
+width = 1.0
+height = 1.0
+
+
+
+
+# Prestrain Types: 
+
+#1 Isotropic Pressure
+#    Func2Tensor B1_ = [this] (const Domain& x) {              // ISOTROPIC PRESSURE
+#                if (abs(x[0]) > (theta/2)  && x[2] > 0)
+#                    return MatrixRT{{p1, 0.0 , 0.0}, {0.0, p1, 0.0}, {0.0, 0.0, p1}};
+#                if (abs(x[0]) < (theta/2)  && x[2] < 0)
+#                    return MatrixRT{{p2, 0.0 , 0.0}, {0.0, p2, 0.0}, {0.0, 0.0, p2}};
+#                else
+#                    return MatrixRT{{0.0, 0.0 , 0.0}, {0.0, 0.0, 0.0}, {0.0, 0.0, 0.0}};
+#            };
+
+
+#############################################
+#  Assembly options
+#############################################
+
+set_IntegralZero = true
+#set_IntegralZero = false
+
+#arbitraryLocalIndex = 7
+#arbitraryElementNumber = 3
+
+#arbitraryLocalIndex = 0
+#arbitraryElementNumber = 0
+
+
+#############################################
+#  Solver Type
+#############################################
+# 1: CG-Solver  # 2: GMRES   # 3: QR
+
+Solvertype = 1
+
+#write_corrector_phi1 = false
+#write_corrector_phi2 = false
+#write_corrector_phi3 = false
+#write_corrector_phi1 = true
+#write_corrector_phi2 = true
+#write_corrector_phi3 = true
+
+write_L2Error = true
+#write_IntegralMean = true
+
+#############################################
+#  Define Analytic Solutions
+#############################################
+
+#b1 = (-(theta/4.0)*mu1*mu2)/(theta*mu1+(1.0-theta)*mu2)
+
+
+
+
+
+
+
+

src/CMakeLists.txt

+add_executable("dune-microstructure" dune-microstructure.cc)
+target_link_dune_default_libraries("dune-microstructure")
+set(CMAKE_BUILD_TYPE Debug)