- per1first perpendicular direction to propagation
C++ Type:std::vector<VariableName>
Controllable:No
Description:first perpendicular direction to propagation
- per2second perpendicular direction to propagation
C++ Type:std::vector<VariableName>
Controllable:No
Description:second perpendicular direction to propagation
- variableThe name of the variable that this object applies to
C++ Type:AuxVariableName
Controllable:No
Description:The name of the variable that this object applies to
Birefringence
Computes the difference between refractive indices (birefringence).
Overview
Computes the difference between refractive indices and where and denote user-defined directions. Since this is an AuxKernel, these quantites are computed locally and can modulate in the simulation box if there are inhomogeneous fields (i.e. electric or elastic) that cause variations.
Example Input File Syntax
Input Parameters
- blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
- boundaryThe list of boundaries (ids or names) from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies
- check_boundary_restrictedTrueWhether to check for multiple element sides on the boundary in the case of a boundary restricted, element aux variable. Setting this to false will allow contribution to a single element's elemental value(s) from multiple boundary sides on the same element (example: when the restricted boundary exists on two or more sides of an element, such as at a corner of a mesh
Default:True
C++ Type:bool
Controllable:No
Description:Whether to check for multiple element sides on the boundary in the case of a boundary restricted, element aux variable. Setting this to false will allow contribution to a single element's elemental value(s) from multiple boundary sides on the same element (example: when the restricted boundary exists on two or more sides of an element, such as at a corner of a mesh
- execute_onLINEAR TIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, ALWAYS, PRE_DISPLACE.
Default:LINEAR TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, ALWAYS, PRE_DISPLACE
Controllable:No
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM, ALWAYS, PRE_DISPLACE.
- prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
Optional Parameters
- control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
- enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.
- seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Controllable:No
Description:The seed for the master random number generator
- use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Advanced Parameters
Input Files
(test/tests/electrooptics/electrooptic_test.i)
[Mesh]
[file]
type = FileMeshGenerator
file = 6grains.e
[]
[./add_side_sets]
input = file
type = SideSetsFromNormalsGenerator
normals = '1 0 0
0 0 1
0 0 -1'
new_boundary = 'left side1 side2'
[../]
[]
[GlobalParams]
displacements = 'disp_x disp_y disp_z'
n_a = 2.91
n_b = 2.91
n_g = 3.27
[]
[Variables]
[./disp_x]
order = FIRST
family = LAGRANGE
[../]
[./disp_y]
order = FIRST
family = LAGRANGE
[../]
[./disp_z]
order = FIRST
family = LAGRANGE
[../]
[./potential_E_int]
order = FIRST
family = LAGRANGE
block = '1 2 3 4 5 6'
[../]
[]
[AuxVariables]
[./stress_xx_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_yy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_xy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_xz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_zz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_yz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xx_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_yy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_zz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_yz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_1]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_2]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_3]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_bire_12]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_bire_23]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_bire_13]
order = CONSTANT
family = MONOMIAL
[../]
[./n_1]
order = CONSTANT
family = MONOMIAL
[../]
[./n_2]
order = CONSTANT
family = MONOMIAL
[../]
[./n_3]
order = CONSTANT
family = MONOMIAL
[../]
[]
[AuxKernels]
[./matl_e11]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 0
variable = strain_xx_elastic
execute_on = 'timestep_end'
[../]
[./matl_e12]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 1
variable = strain_xy_elastic
execute_on = 'timestep_end'
[../]
[./matl_e13]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 2
variable = strain_xz_elastic
execute_on = 'timestep_end'
[../]
[./matl_e22]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 1
index_j = 1
variable = strain_yy_elastic
execute_on = 'timestep_end'
[../]
[./matl_e23]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 1
index_j = 2
variable = strain_yz_elastic
execute_on = 'timestep_end'
[../]
[./matl_e33]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 2
index_j = 2
variable = strain_zz_elastic
execute_on = 'timestep_end'
[../]
[./matl_s11]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 0
variable = stress_xx_elastic
execute_on = 'timestep_end'
[../]
[./matl_s12]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 1
variable = stress_xy_elastic
execute_on = 'timestep_end'
[../]
[./matl_s13]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 2
variable = stress_xz_elastic
execute_on = 'timestep_end'
[../]
[./matl_s22]
type = RankTwoAux
rank_two_tensor = stress
index_i = 1
index_j = 1
variable = stress_yy_elastic
execute_on = 'timestep_end'
[../]
[./matl_s23]
type = RankTwoAux
rank_two_tensor = stress
index_i = 1
index_j = 2
variable = stress_yz_elastic
execute_on = 'timestep_end'
[../]
[./matl_s33]
type = RankTwoAux
rank_two_tensor = stress
index_i = 2
index_j = 2
variable = stress_zz_elastic
execute_on = 'timestep_end'
[../]
[./dn_s1]
type = ChangeInRefractiveIndexElectro
index_i = 0
index_j = 0
index_k = 0
index_l = 0
variable = dn_1
execute_on = 'timestep_end'
[../]
[./dn_s2]
type = ChangeInRefractiveIndexElectro
index_i = 1
index_j = 1
index_k = 1
index_l = 1
variable = dn_2
execute_on = 'timestep_end'
[../]
[./dn_s3]
type = ChangeInRefractiveIndexElectro
index_i = 2
index_j = 2
index_k = 2
index_l = 2
variable = dn_3
execute_on = 'timestep_end'
[../]
[./dn_bire_s12]
type = Birefringence
variable = dn_bire_12
per1 = n_1
per2 = n_2
execute_on = 'timestep_end'
[../]
[./dn_bire_s23]
type = Birefringence
variable = dn_bire_23
per1 = n_2
per2 = n_3
execute_on = 'timestep_end'
[../]
[./dn_bire_s13]
type = Birefringence
variable = dn_bire_13
per1 = n_1
per2 = n_3
execute_on = 'timestep_end'
[../]
[./n_1_c]
type = RefractiveIndex
variable = n_1
electro = true
index_j = 0
index_k = 0
var1 = dn_1
execute_on = 'timestep_end'
[../]
[./n_2_c]
type = RefractiveIndex
variable = n_2
electro = true
index_j = 1
index_k = 1
var1 = dn_2
execute_on = 'timestep_end'
[../]
[./n_3_c]
type = RefractiveIndex
variable = n_3
electro = true
index_j = 2
index_k = 2
var1 = dn_3
execute_on = 'timestep_end'
[../]
[]
################################################
# Block list: #
# #
# No 99? #
# #
################################################
[Materials]
[./eigen_strain_zz] #Use for stress-free strain (ie epitaxial)
type = ComputeEigenstrain
block = '1 2 3 4 5 6'
# eigen_base = 'exx exy exz eyx eyy eyz ezx ezy ezz'
eigen_base = '1 0 0 0 1 0 0 0 0'
eigenstrain_name = eigenstrain
prefactor = 0.0
[../]
[./elasticity_tensor_1]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 0.0
euler_angle_2 = 15.0
euler_angle_3 = 65.0
block = '1'
[../]
[./strain_1]
type = ComputeSmallStrain
block = '1'
eigenstrain_names = eigenstrain
[../]
[./stress_1]
type = ComputeLinearElasticStress
block = '1'
[../]
[./photoelastic_tensor_1]
type = ComputeElectroopticTensor
fill_method = general #For HgS (32)
# Fill Method (see Nye page 113)
# r111, r112, r113, r121, r122, r123, r131, r132, r133,
# r211, r212, r213, r221, r222, r223, r231, r232, r233,
# r311, r312, r313, r321, r322, r323, r331, r332, r333,
r_ijk ='0.0031 0 0 0 -0.0031 0 0 -0.0015 0 0 -0.0031 0 -0.0031 0 0 0.0015 0 0 0 -0.0015 0 0.0015 0 0 0 0 0'
euler_angle_1 = 0.0
euler_angle_2 = 15.0
euler_angle_3 = 65.0
block = '1'
[../]
[./beta_tensor_1]
type = ComputeIndicatrix
block = '1'
euler_angle_1 = 0.0
euler_angle_2 = 15.0
euler_angle_3 = 65.0
[../]
[./delta_beta_tensor_1]
type = ComputeDeltaIndicatrixElectro
block = '1'
potential_E_int = potential_E_int
[../]
[./elasticity_tensor_2]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 12.0
euler_angle_2 = -15.0
euler_angle_3 = 65.0
block = '2'
[../]
[./strain_2]
type = ComputeSmallStrain
block = '2'
eigenstrain_names = eigenstrain
[../]
[./stress_2]
type = ComputeLinearElasticStress
block = '2'
[../]
[./photoelastic_tensor_2]
type = ComputeElectroopticTensor
fill_method = general #For HgS (32)
# Fill Method (see Nye page 113)
# r111, r112, r113, r121, r122, r123, r131, r132, r133,
# r211, r212, r213, r221, r222, r223, r231, r232, r233,
# r311, r312, r313, r321, r322, r323, r331, r332, r333,
r_ijk ='0.0031 0 0 0 -0.0031 0 0 -0.0015 0 0 -0.0031 0 -0.0031 0 0 0.0015 0 0 0 -0.0015 0 0.0015 0 0 0 0 0'
euler_angle_1 = 12.0
euler_angle_2 = -15.0
euler_angle_3 = 65.0
block = '2'
[../]
[./beta_tensor_2]
type = ComputeIndicatrix
block = '2'
euler_angle_1 = 12.0
euler_angle_2 = -15.0
euler_angle_3 = 65.0
[../]
[./delta_beta_tensor_2]
type = ComputeDeltaIndicatrixElectro
block = '2'
potential_E_int = potential_E_int
[../]
[./elasticity_tensor_3]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 1.0
euler_angle_2 = 50.0
euler_angle_3 = -132.0
block = '3'
[../]
[./strain_3]
type = ComputeSmallStrain
block = '3'
eigenstrain_names = eigenstrain
[../]
[./stress_3]
type = ComputeLinearElasticStress
block = '3'
[../]
[./photoelastic_tensor_3]
type = ComputeElectroopticTensor
fill_method = general #For HgS (32)
# Fill Method (see Nye page 113)
# r111, r112, r113, r121, r122, r123, r131, r132, r133,
# r211, r212, r213, r221, r222, r223, r231, r232, r233,
# r311, r312, r313, r321, r322, r323, r331, r332, r333,
r_ijk ='0.0031 0 0 0 -0.0031 0 0 -0.0015 0 0 -0.0031 0 -0.0031 0 0 0.0015 0 0 0 -0.0015 0 0.0015 0 0 0 0 0'
euler_angle_1 = 1.0
euler_angle_2 = 50.0
euler_angle_3 = -132.0
block = '3'
[../]
[./beta_tensor_3]
type = ComputeIndicatrix
block = '3'
euler_angle_1 = 1.0
euler_angle_2 = 50.0
euler_angle_3 = -132.0
[../]
[./delta_beta_tensor_3]
type = ComputeDeltaIndicatrixElectro
block = '3'
potential_E_int = potential_E_int
[../]
[./elasticity_tensor_4]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 100.0
euler_angle_2 = 0.0
euler_angle_3 = 0.0
block = '4'
[../]
[./strain_4]
type = ComputeSmallStrain
block = '4'
eigenstrain_names = eigenstrain
[../]
[./stress_4]
type = ComputeLinearElasticStress
block = '4'
[../]
[./photoelastic_tensor_4]
type = ComputeElectroopticTensor
fill_method = general #For HgS (32)
# Fill Method (see Nye page 113)
# r111, r112, r113, r121, r122, r123, r131, r132, r133,
# r211, r212, r213, r221, r222, r223, r231, r232, r233,
# r311, r312, r313, r321, r322, r323, r331, r332, r333,
r_ijk ='0.0031 0 0 0 -0.0031 0 0 -0.0015 0 0 -0.0031 0 -0.0031 0 0 0.0015 0 0 0 -0.0015 0 0.0015 0 0 0 0 0'
euler_angle_1 = 100.0
euler_angle_2 = 0.0
euler_angle_3 = 0.0
block = '4'
[../]
[./beta_tensor_4]
type = ComputeIndicatrix
block = '4'
euler_angle_1 = 100.0
euler_angle_2 = 0.0
euler_angle_3 = 0.0
[../]
[./delta_beta_tensor_4]
type = ComputeDeltaIndicatrixElectro
block = '4'
potential_E_int = potential_E_int
[../]
[./elasticity_tensor_5]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 200.0
euler_angle_2 = 33.0
euler_angle_3 = 5.0
block = '5'
[../]
[./strain_5]
type = ComputeSmallStrain
block = '5'
eigenstrain_names = eigenstrain
[../]
[./stress_5]
type = ComputeLinearElasticStress
block = '5'
[../]
[./photoelastic_tensor_5]
type = ComputeElectroopticTensor
fill_method = general #For HgS (32)
# Fill Method (see Nye page 113)
# r111, r112, r113, r121, r122, r123, r131, r132, r133,
# r211, r212, r213, r221, r222, r223, r231, r232, r233,
# r311, r312, r313, r321, r322, r323, r331, r332, r333,
r_ijk ='0.0031 0 0 0 -0.0031 0 0 -0.0015 0 0 -0.0031 0 -0.0031 0 0 0.0015 0 0 0 -0.0015 0 0.0015 0 0 0 0 0'
euler_angle_1 = 200.0
euler_angle_2 = 33.0
euler_angle_3 = 5.0
block = '5'
[../]
[./beta_tensor_5]
type = ComputeIndicatrix
block = '5'
euler_angle_1 = 200.0
euler_angle_2 = 33.0
euler_angle_3 = 5.0
[../]
[./delta_beta_tensor_5]
type = ComputeDeltaIndicatrixElectro
block = '5'
potential_E_int = potential_E_int
[../]
[./elasticity_tensor_6]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 110.0
euler_angle_2 = 0.0
euler_angle_3 = 100.0
block = '6'
[../]
[./strain_6]
type = ComputeSmallStrain
block = '6'
eigenstrain_names = eigenstrain
[../]
[./stress_6]
type = ComputeLinearElasticStress
block = '6'
[../]
[./photoelastic_tensor_6]
type = ComputeElectroopticTensor
fill_method = general #For HgS (32)
# Fill Method (see Nye page 113)
# r111, r112, r113, r121, r122, r123, r131, r132, r133,
# r211, r212, r213, r221, r222, r223, r231, r232, r233,
# r311, r312, r313, r321, r322, r323, r331, r332, r333,
r_ijk ='0.0031 0 0 0 -0.0031 0 0 -0.0015 0 0 -0.0031 0 -0.0031 0 0 0.0015 0 0 0 -0.0015 0 0.0015 0 0 0 0 0'
euler_angle_1 = 110.0
euler_angle_2 = 0.0
euler_angle_3 = 100.0
block = '6'
[../]
[./beta_tensor_6]
type = ComputeIndicatrix
block = '6'
euler_angle_1 = 110.0
euler_angle_2 = 0.0
euler_angle_3 = 100.0
[../]
[./delta_beta_tensor_6]
type = ComputeDeltaIndicatrixElectro
block = '6'
potential_E_int = potential_E_int
[../]
[./permitivitty_1]
type = GenericConstantMaterial
prop_names = 'permittivity'
prop_values = '1.0'
[../]
[]
[Kernels]
#Elastic problem
[./TensorMechanics]
#This is an action block
[../]
[./ElectroStats]
type = Electrostatics
variable = potential_E_int
block = '1 2 3 4 5 6'
[../]
[]
[BCs]
[./pot_top]
type = DirichletBC
variable = 'potential_E_int'
value = -10
boundary = 'side1'
[../]
[./pot_bottom]
type = DirichletBC
variable = 'potential_E_int'
value = 0.0
boundary = 'side2'
[../]
[./z_bottom]
type = DirichletBC
variable = 'disp_z'
value = 0.0
boundary = 'side2'
[../]
[./y_bottom]
type = DirichletBC
variable = 'disp_y'
value = 0.0
boundary = 'side2'
[../]
[./x_bottom]
type = DirichletBC
variable = 'disp_z'
value = 0.0
boundary = 'side2'
[../]
[]
[Preconditioning]
[./smp]
type = SMP
full = true
petsc_options_iname = '-ksp_gmres_restart -snes_atol -snes_rtol -ksp_rtol -pc_type -pc_hypre_type'
petsc_options_value = ' 250 1e-10 1e-8 1e-6 hypre boomeramg '
[../]
[]
[Executioner]
type = Steady
solve_type = 'NEWTON' #"PJFNK, JFNK, NEWTON"
[]
[Outputs]
print_linear_residuals = false
[./out]
type = Exodus
execute_on = 'timestep_end'
file_base = out_6grain_electro
elemental_as_nodal = true
[../]
[]
(test/tests/photoelasticity/photoelastic_test.i)
[Mesh]
[file]
type = FileMeshGenerator
file = 6grains.e
[]
[./add_side_sets]
input = file
type = SideSetsFromNormalsGenerator
normals = '1 0 0
0 0 1
0 0 -1'
new_boundary = 'left side1 side2'
[../]
[]
[GlobalParams]
displacements = 'disp_x disp_y disp_z'
n_a = 2.437
n_b = 2.437
n_g = 2.365
[]
[Variables]
[./disp_x]
order = FIRST
family = LAGRANGE
[../]
[./disp_y]
order = FIRST
family = LAGRANGE
[../]
[./disp_z]
order = FIRST
family = LAGRANGE
[../]
[]
[AuxVariables]
[./stress_xx_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_yy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_xy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_xz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_zz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_yz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xx_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_yy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_zz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_yz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_1]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_2]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_3]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_bire_12]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_bire_23]
order = CONSTANT
family = MONOMIAL
[../]
[./dn_bire_13]
order = CONSTANT
family = MONOMIAL
[../]
[./n_1]
order = CONSTANT
family = MONOMIAL
[../]
[./n_2]
order = CONSTANT
family = MONOMIAL
[../]
[./n_3]
order = CONSTANT
family = MONOMIAL
[../]
[]
[AuxKernels]
[./matl_e11]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 0
variable = strain_xx_elastic
execute_on = 'timestep_end'
[../]
[./matl_e12]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 1
variable = strain_xy_elastic
execute_on = 'timestep_end'
[../]
[./matl_e13]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 2
variable = strain_xz_elastic
execute_on = 'timestep_end'
[../]
[./matl_e22]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 1
index_j = 1
variable = strain_yy_elastic
execute_on = 'timestep_end'
[../]
[./matl_e23]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 1
index_j = 2
variable = strain_yz_elastic
execute_on = 'timestep_end'
[../]
[./matl_e33]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 2
index_j = 2
variable = strain_zz_elastic
execute_on = 'timestep_end'
[../]
[./matl_s11]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 0
variable = stress_xx_elastic
execute_on = 'timestep_end'
[../]
[./matl_s12]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 1
variable = stress_xy_elastic
execute_on = 'timestep_end'
[../]
[./matl_s13]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 2
variable = stress_xz_elastic
execute_on = 'timestep_end'
[../]
[./matl_s22]
type = RankTwoAux
rank_two_tensor = stress
index_i = 1
index_j = 1
variable = stress_yy_elastic
execute_on = 'timestep_end'
[../]
[./matl_s23]
type = RankTwoAux
rank_two_tensor = stress
index_i = 1
index_j = 2
variable = stress_yz_elastic
execute_on = 'timestep_end'
[../]
[./matl_s33]
type = RankTwoAux
rank_two_tensor = stress
index_i = 2
index_j = 2
variable = stress_zz_elastic
execute_on = 'timestep_end'
[../]
[./dn_s1]
type = ChangeInRefractiveIndex
index_i = 0
index_j = 0
index_k = 0
index_l = 0
variable = dn_1
execute_on = 'timestep_end'
[../]
[./dn_s2]
type = ChangeInRefractiveIndex
index_i = 1
index_j = 1
index_k = 1
index_l = 1
variable = dn_2
execute_on = 'timestep_end'
[../]
[./dn_s3]
type = ChangeInRefractiveIndex
index_i = 2
index_j = 2
index_k = 2
index_l = 2
variable = dn_3
execute_on = 'timestep_end'
[../]
[./dn_bire_s12]
type = Birefringence
variable = dn_bire_12
per1 = n_1
per2 = n_2
execute_on = 'timestep_end'
[../]
[./dn_bire_s23]
type = Birefringence
variable = dn_bire_23
per1 = n_2
per2 = n_3
execute_on = 'timestep_end'
[../]
[./dn_bire_s13]
type = Birefringence
variable = dn_bire_13
per1 = n_1
per2 = n_3
execute_on = 'timestep_end'
[../]
[./n_1_c]
type = RefractiveIndex
variable = n_1
elasto = true
index_j = 0
index_k = 0
var1 = dn_1
execute_on = 'timestep_end'
[../]
[./n_2_c]
type = RefractiveIndex
variable = n_2
elasto = true
index_j = 1
index_k = 1
var1 = dn_2
execute_on = 'timestep_end'
[../]
[./n_3_c]
type = RefractiveIndex
variable = n_3
elasto = true
index_j = 2
index_k = 2
var1 = dn_3
execute_on = 'timestep_end'
[../]
[]
################################################
# Block list: #
# #
# No 99? #
# #
################################################
[Materials]
[./eigen_strain_zz] #Use for stress-free strain (ie epitaxial)
type = ComputeEigenstrain
block = '1 2 3 4 5 6'
# eigen_base = 'exx exy exz eyx eyy eyz ezx ezy ezz'
eigen_base = '1 0 0 0 1 0 0 0 0'
eigenstrain_name = eigenstrain
prefactor = 0.0
[../]
[./elasticity_tensor_1]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 0.0
euler_angle_2 = 15.0
euler_angle_3 = 65.0
block = '1'
[../]
[./strain_1]
type = ComputeSmallStrain
block = '1'
eigenstrain_names = eigenstrain
[../]
[./stress_1]
type = ComputeLinearElasticStress
block = '1'
[../]
[./photoelastic_tensor_1]
type = ComputeElastoopticTensor
fill_method = symmetric21 #BTO is not symmetric21 FIX!!
# Use BaTiO3, crystal symmetry P4mm.
P_mnkl = '0.5 0.106 0.2 0.0 0.0 0.0 0.5 0.2 0.0 0.0 0.0 0.77 0.0 0.0 0.0 1.0 0.0 0.0 1.0 0.0 0.1'
euler_angle_1 = 0.0
euler_angle_2 = 15.0
euler_angle_3 = 65.0
block = '1'
[../]
[./beta_tensor_1]
type = ComputeIndicatrix
block = '1'
euler_angle_1 = 0.0
euler_angle_2 = 15.0
euler_angle_3 = 65.0
[../]
[./delta_beta_tensor_1]
type = ComputeDeltaIndicatrix
block = '1'
[../]
[./elasticity_tensor_2]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 12.0
euler_angle_2 = -15.0
euler_angle_3 = 65.0
block = '2'
[../]
[./strain_2]
type = ComputeSmallStrain
block = '2'
eigenstrain_names = eigenstrain
[../]
[./stress_2]
type = ComputeLinearElasticStress
block = '2'
[../]
[./photoelastic_tensor_2]
type = ComputeElastoopticTensor
fill_method = symmetric21 #BTO is not symmetric21 FIX!!
# Use BaTiO3, crystal symmetry P4mm.
P_mnkl = '0.5 0.106 0.2 0.0 0.0 0.0 0.5 0.2 0.0 0.0 0.0 0.77 0.0 0.0 0.0 1.0 0.0 0.0 1.0 0.0 0.1'
euler_angle_1 = 12.0
euler_angle_2 = -15.0
euler_angle_3 = 65.0
block = '2'
[../]
[./beta_tensor_2]
type = ComputeIndicatrix
block = '2'
euler_angle_1 = 12.0
euler_angle_2 = -15.0
euler_angle_3 = 65.0
[../]
[./delta_beta_tensor_2]
type = ComputeDeltaIndicatrix
block = '2'
[../]
[./elasticity_tensor_3]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 1.0
euler_angle_2 = 50.0
euler_angle_3 = -132.0
block = '3'
[../]
[./strain_3]
type = ComputeSmallStrain
block = '3'
eigenstrain_names = eigenstrain
[../]
[./stress_3]
type = ComputeLinearElasticStress
block = '3'
[../]
[./photoelastic_tensor_3]
type = ComputeElastoopticTensor
fill_method = symmetric21 #BTO is not symmetric21 FIX!!
# Use BaTiO3, crystal symmetry P4mm.
P_mnkl = '0.5 0.106 0.2 0.0 0.0 0.0 0.5 0.2 0.0 0.0 0.0 0.77 0.0 0.0 0.0 1.0 0.0 0.0 1.0 0.0 0.1'
euler_angle_1 = 1.0
euler_angle_2 = 50.0
euler_angle_3 = -132.0
block = '3'
[../]
[./beta_tensor_3]
type = ComputeIndicatrix
block = '3'
euler_angle_1 = 1.0
euler_angle_2 = 50.0
euler_angle_3 = -132.0
[../]
[./delta_beta_tensor_3]
type = ComputeDeltaIndicatrix
block = '3'
[../]
[./elasticity_tensor_4]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 100.0
euler_angle_2 = 0.0
euler_angle_3 = 0.0
block = '4'
[../]
[./strain_4]
type = ComputeSmallStrain
block = '4'
eigenstrain_names = eigenstrain
[../]
[./stress_4]
type = ComputeLinearElasticStress
block = '4'
[../]
[./photoelastic_tensor_4]
type = ComputeElastoopticTensor
fill_method = symmetric21 #BTO is not symmetric21 FIX!!
# Use BaTiO3, crystal symmetry P4mm.
P_mnkl = '0.5 0.106 0.2 0.0 0.0 0.0 0.5 0.2 0.0 0.0 0.0 0.77 0.0 0.0 0.0 1.0 0.0 0.0 1.0 0.0 0.1'
euler_angle_1 = 100.0
euler_angle_2 = 0.0
euler_angle_3 = 0.0
block = '4'
[../]
[./beta_tensor_4]
type = ComputeIndicatrix
block = '4'
euler_angle_1 = 100.0
euler_angle_2 = 0.0
euler_angle_3 = 0.0
[../]
[./delta_beta_tensor_4]
type = ComputeDeltaIndicatrix
block = '4'
[../]
[./elasticity_tensor_5]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 200.0
euler_angle_2 = 33.0
euler_angle_3 = 5.0
block = '5'
[../]
[./strain_5]
type = ComputeSmallStrain
block = '5'
eigenstrain_names = eigenstrain
[../]
[./stress_5]
type = ComputeLinearElasticStress
block = '5'
[../]
[./photoelastic_tensor_5]
type = ComputeElastoopticTensor
fill_method = symmetric21 #BTO is not symmetric21 FIX!!
# Use BaTiO3, crystal symmetry P4mm.
P_mnkl = '0.5 0.106 0.2 0.0 0.0 0.0 0.5 0.2 0.0 0.0 0.0 0.77 0.0 0.0 0.0 1.0 0.0 0.0 1.0 0.0 0.1'
euler_angle_1 = 200.0
euler_angle_2 = 33.0
euler_angle_3 = 5.0
block = '5'
[../]
[./beta_tensor_5]
type = ComputeIndicatrix
block = '5'
euler_angle_1 = 200.0
euler_angle_2 = 33.0
euler_angle_3 = 5.0
[../]
[./delta_beta_tensor_5]
type = ComputeDeltaIndicatrix
block = '5'
[../]
[./elasticity_tensor_6]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 110.0
euler_angle_2 = 0.0
euler_angle_3 = 100.0
block = '6'
[../]
[./strain_6]
type = ComputeSmallStrain
block = '6'
eigenstrain_names = eigenstrain
[../]
[./stress_6]
type = ComputeLinearElasticStress
block = '6'
[../]
[./photoelastic_tensor_6]
type = ComputeElastoopticTensor
fill_method = symmetric21 #BTO is not symmetric21 FIX!!
# Use BaTiO3, crystal symmetry P4mm.
P_mnkl = '0.5 0.106 0.2 0.0 0.0 0.0 0.5 0.2 0.0 0.0 0.0 0.77 0.0 0.0 0.0 1.0 0.0 0.0 1.0 0.0 0.1'
euler_angle_1 = 110.0
euler_angle_2 = 0.0
euler_angle_3 = 100.0
block = '6'
[../]
[./beta_tensor_6]
type = ComputeIndicatrix
block = '6'
euler_angle_1 = 110.0
euler_angle_2 = 0.0
euler_angle_3 = 100.0
[../]
[./delta_beta_tensor_6]
type = ComputeDeltaIndicatrix
block = '6'
[../]
[]
[Kernels]
#Elastic problem
[./TensorMechanics]
#This is an action block
[../]
[]
[BCs]
[./center_disp_z_top]
type = DirichletBC
variable = 'disp_z'
value = -0.1
boundary = 'side1'
[../]
[./center_disp_z_bottom]
type = DirichletBC
variable = 'disp_z'
value = 0.0
boundary = 'side2'
[../]
[./center_disp_x_bottom]
type = DirichletBC
variable = 'disp_x'
value = 0.0
boundary = 'side2'
[../]
[./center_disp_y_bottom]
type = DirichletBC
variable = 'disp_y'
value = 0.0
boundary = 'side2'
[../]
[]
[Preconditioning]
[./smp]
type = SMP
full = true
petsc_options_iname = '-ksp_gmres_restart -snes_atol -snes_rtol -ksp_rtol -pc_type -pc_hypre_type'
petsc_options_value = ' 250 1e-10 1e-8 1e-6 hypre boomeramg '
[../]
[]
[Executioner]
type = Steady
solve_type = 'NEWTON' #"PJFNK, JFNK, NEWTON"
[]
[Outputs]
print_linear_residuals = false
[./out]
type = Exodus
execute_on = 'timestep_end'
file_base = out_6grain_photo
elemental_as_nodal = true
[../]
[]
(examples/other/Al2O3_photo_elastic.i)
[Mesh]
file = exodus_disk_r8_h1_photo.e
uniform_refine = 1
[]
[MeshModifiers]
[./centernodeset_1]
type = AddExtraNodeset
new_boundary = 'center_node_1'
coord = '0.0 0.0 -0.5'
[../]
[./centernodeset_2]
type = AddExtraNodeset
new_boundary = 'center_node_2'
coord = '0.0 0.0 0.5'
[../]
[./centernodeset_3]
type = AddExtraNodeset
new_boundary = 'center_node_3'
coord = '0.0 0.0 0.166667'
[../]
[./centernodeset_4]
type = AddExtraNodeset
new_boundary = 'center_node_4'
coord = '0.0 0.0 -0.166667'
[../]
[]
[GlobalParams]
displacements = 'disp_x disp_y disp_z'
[]
[Variables]
[./disp_x]
order = FIRST
family = LAGRANGE
[../]
[./disp_y]
order = FIRST
family = LAGRANGE
[../]
[./disp_z]
order = FIRST
family = LAGRANGE
[../]
[]
[AuxVariables]
[./stress_xx_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_yy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_xy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_xz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_zz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./stress_yz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xx_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_yy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xy_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_xz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_zz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./strain_yz_elastic]
order = CONSTANT
family = MONOMIAL
[../]
[./n_o_ij]
order = CONSTANT
family = MONOMIAL
[../]
[./n_e_ij]
order = CONSTANT
family = MONOMIAL
[../]
[./dn]
order = CONSTANT
family = MONOMIAL
[../]
[]
[AuxKernels]
[./matl_e11]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 0
variable = strain_xx_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_e12]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 1
variable = strain_xy_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_e13]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 0
index_j = 2
variable = strain_xz_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_e22]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 1
index_j = 1
variable = strain_yy_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_e23]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 1
index_j = 2
variable = strain_yz_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_e33]
type = RankTwoAux
rank_two_tensor = elastic_strain
index_i = 2
index_j = 2
variable = strain_zz_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_s11]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 0
variable = stress_xx_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_s12]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 1
variable = stress_xy_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_s13]
type = RankTwoAux
rank_two_tensor = stress
index_i = 0
index_j = 2
variable = stress_xz_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_s22]
type = RankTwoAux
rank_two_tensor = stress
index_i = 1
index_j = 1
variable = stress_yy_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_s23]
type = RankTwoAux
rank_two_tensor = stress
index_i = 1
index_j = 2
variable = stress_yz_elastic
execute_on = 'timestep_end'
block = '1 2 3 4'
[../]
[./matl_s33]
type = RankTwoAux
rank_two_tensor = stress
index_i = 2
index_j = 2
variable = stress_zz_elastic
block = '1 2 3 4'
execute_on = 'timestep_end'
[../]
[./ref_deps_00_ij]
type = RefractiveIndex
index_one = 0
index_two = 0
variable = n_o_ij
execute_on = 'timestep_end'
[../]
[./ref_deps_22_ij]
type = RefractiveIndex
index_one = 2
index_two = 2
variable = n_e_ij
execute_on = 'timestep_end'
[../]
[./dn_r]
type = Birefringence
variable = dn
n_o = n_o_ij
n_e = n_e_ij
execute_on = 'timestep_end'
[../]
[]
[Materials]
[./eigen_strain_zz] #Use for stress-free strain (ie epitaxial)
type = ComputeEigenstrain
block = '1 2 3 4'
# eigen_base = 'exx exy exz eyx eyy eyz ezx ezy ezz'
eigen_base = '1 0 0 0 1 0 0 0 0'
eigenstrain_name = eigenstrain
prefactor = 0.0
[../]
[./elasticity_tensor_1]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
block = '2 3'
[../]
[./strain_1]
type = ComputeSmallStrain
block = '2 3'
eigenstrain_names = eigenstrain
[../]
[./stress_1]
type = ComputeLinearElasticStress
block = '2 3'
[../]
[./photoelastic_tensor_1]
type = ComputePhotostrictiveTensor
fill_method = symmetric9
#S11 S21 S31 S12 S22 S32 S13 S23 S33
P_mnkl = '0.5 0.106 0.20 0.50 0.106 0.5 0.77 0.77 0.77'
block = '2 3'
[../]
[./beta_tensor_1]
type = ComputeBetaTensor
block = '2 3'
[../]
[./indicatrix_1]
type = ComputeIndicatrix
block = '2 3'
#from Handbook of Optics Chp. 17
refractive_index_bulk_ordinary = 2.437
refractive_index_bulk_extraordinary = 2.365
[../]
[./elasticity_tensor_2]
type = ComputeElasticityTensor
fill_method = symmetric9
#BaTiO3 from MaterialsProject
C_ijkl = '260.06 105.79 76.90 260.06 105.79 260.06 81.57 81.57 116.28'
euler_angle_1 = 0.0 # 45.0
euler_angle_2 = 54.74 #37.0
euler_angle_3 = 45.0 #-76.0 #ne = -2.162e-3, no = -4.279e-4, bf = -1.734e-3
block = '1 4'
[../]
[./strain_2]
type = ComputeSmallStrain
block = '1 4'
eigenstrain_names = eigenstrain
[../]
[./stress_2]
type = ComputeLinearElasticStress
block = '1 4'
[../]
[./photoelastic_tensor_2]
type = ComputePhotostrictiveTensor
fill_method = symmetric9
P_mnkl = '0.5 0.106 0.20 0.50 0.106 0.5 0.77 0.77 0.77'
block = '1 4'
[../]
[./beta_tensor_2]
type = ComputeBetaTensor
block = '1 4'
[../]
[./indicatrix_2]
type = ComputeIndicatrix
block = '1 4'
#from BTO paper
refractive_index_bulk_ordinary = 2.412
refractive_index_bulk_extraordinary = 2.360
euler_angle_1 = 0.0 # 45.0
euler_angle_2 = 54.74 #37.0
euler_angle_3 = 45.0 #-76.0
[../]
[]
[Kernels]
#Elastic problem
[./TensorMechanics]
#This is an action block
[../]
[]
[BCs]
[./center_disp_z_top]
type = DirichletBC
variable = 'disp_z'
value = -0.0025
boundary = '5'
[../]
[./center_disp_z_bottom]
type = DirichletBC
variable = 'disp_z'
value = 0.0025
boundary = '6'
[../]
[]
[Preconditioning]
[./smp]
type = SMP
full = true
petsc_options_iname = '-ksp_gmres_restart -snes_atol -snes_rtol -ksp_rtol -pc_type -pc_hypre_type'
petsc_options_value = ' 250 1e-10 1e-8 1e-6 hypre boomeramg '
[../]
[]
[Executioner]
type = Steady
solve_type = 'NEWTON' #"PJFNK, JFNK, NEWTON"
[]
[Outputs]
print_linear_residuals = false
print_perf_log = true
[./out]
type = Exodus
execute_on = 'timestep_end'
file_base = out_disk
elemental_as_nodal = true
[../]
[]