mofem/html/sorption__cube_8jou_source.html

reset

set duplicate block elements on


brick x 0.0025 y 0.00041 z 0.00041


# Set block with volume with gel material

block 1 volume all

block 1 name 'GEL'


# Set Dirichelt boundary condition for solvent concentration

block 2 surface 2 3 4

block 2 name 'CHEMICAL_LOAD_1'

block 2 attribute count 1

block 2 attribute index 1 6.624729054813950e-18


# Set solvent flux, i.e.

block 4 surface 1 5 6

block 4 name 'FLUX_CHEMICAL_LOAD_2'

block 4 attribute count 1

block 4 attribute index 1 0


# Kinematic boundary condition for mechanical field

create displacement on surface 6 dof 1 fix 0


# Kinematic boundary condition for mechanical field

create displacement on surface 5 dof 2 fix 0


# Kinematic boundary condition for mechanical field

create displacement on surface 1 dof 3 fix 0


# Make a mesh

volume all scheme Tetmesh

volume all size auto factor 7

mesh volume all


# Set block 3 and set 10 node tetrahedrons in that block

block 3 tet all

block 3 element type TETRA10