Axisymmetric2D3DSolutionFunction

Function for reading a 2D axisymmetric solution from file and mapping it to a 3D Cartesian model

The 2D solution is likely to be the output of a 2D-RZ calculation, which we want to compare to a full 3D model. This is useful for assessing the validity of the 2D-RZ geometric approximation.

The axis of symmetry for the original 2D axisymmetric calculation and for mapping this 2D axisymmetric calculation into the 3D space can both be specified.

Example input syntax

In this example, three Axisymmetric2D3DSolutionFunction are used to load results for a 2D axisymmetric simulation in a 3D mechanics simulation. These results are then used in the BCs block to impose a displacement on a boundary, and in a temperature variable using a FunctionAux.

[Functions]
  [./soln_func_temp]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'temp'
  [../]
  [./soln_func_disp_x]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'disp_x disp_y'
    component = 0
  [../]
  [./soln_func_disp_y]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'disp_x disp_y'
    component = 1
  [../]
  [./soln_func_disp_z]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'disp_x disp_y'
    component = 2
  [../]
[]

Input Parameters

solutionThe SolutionUserObject to extract data from.
C++ Type:UserObjectName
Unit:(no unit assumed)
Controllable:No
Description:The SolutionUserObject to extract data from.

Required Parameters

2d_axis_point10 0 0Start point for axis of symmetry for the 2d model
Default:0 0 0
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Start point for axis of symmetry for the 2d model
2d_axis_point20 1 0End point for axis of symmetry for the 2d model
Default:0 1 0
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:End point for axis of symmetry for the 2d model
3d_axis_point10 0 0Start point for axis of symmetry for the 3d model
Default:0 0 0
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Start point for axis of symmetry for the 3d model
3d_axis_point20 1 0End point for axis of symmetry for the 3d model
Default:0 1 0
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:End point for axis of symmetry for the 3d model
add_factor0Add this value (b) to the solution (x): ax+b, where a is the 'scale_factor'
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Add this value (b) to the solution (x): ax+b, where a is the 'scale_factor'
axial_dimension_ratio1Ratio of the axial dimension in the 3d model to that in the 2d model. Optionally permits the 3d model to be larger than the 2d model in that dimension, and scales vector solutions in that direction by this factor.
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Ratio of the axial dimension in the 3d model to that in the 2d model. Optionally permits the 3d model to be larger than the 2d model in that dimension, and scales vector solutions in that direction by this factor.
componentComponent of the variable to be computed if it is a vector
C++ Type:unsigned int
Unit:(no unit assumed)
Controllable:No
Description:Component of the variable to be computed if it is a vector
from_variablesThe names of the variables in the file that are to be extracted, in x, y order if they are vector components
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
Controllable:No
Description:The names of the variables in the file that are to be extracted, in x, y order if they are vector components
scale_factor1Scale factor (a) to be applied to the solution (x): ax+b, where b is the 'add_factor'
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scale factor (a) to be applied to the solution (x): ax+b, where b is the 'add_factor'

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
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
Unit:(no unit assumed)
Controllable:No
Description:Set the enabled status of the MooseObject.

Advanced Parameters

(moose/modules/combined/test/tests/axisymmetric_2d3d_solution_function/3dy.i)

[GlobalParams]
  order = FIRST
  family = LAGRANGE
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  file = 3dy.e
[]

[Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]
  [./disp_z]
  [../]
[]

[AuxVariables]
  [./temp]
  [../]
  [./hoop_stress]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[UserObjects]
  [./soln]
    type = SolutionUserObject
    mesh = 2d_out.e
    system_variables = 'disp_x disp_y temp'
  [../]
[]

[Functions]
  [./soln_func_temp]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'temp'
  [../]
  [./soln_func_disp_x]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'disp_x disp_y'
    component = 0
  [../]
  [./soln_func_disp_y]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'disp_x disp_y'
    component = 1
  [../]
  [./soln_func_disp_z]
    type = Axisymmetric2D3DSolutionFunction
    solution = soln
    from_variables = 'disp_x disp_y'
    component = 2
  [../]
[]

[Physics/SolidMechanics/QuasiStatic]
  [./all]
    volumetric_locking_correction = true
    add_variables  = true
    incremental = true
    strain = FINITE
    eigenstrain_names = thermal_expansion
    generate_output = 'stress_xx stress_yy stress_zz vonmises_stress hydrostatic_stress'
  [../]
[]

[AuxKernels]
  [./t_soln_aux]
    type = FunctionAux
    variable = temp
    block = '1 2'
    function = soln_func_temp
  [../]
  [./hoop_stress]
    type = RankTwoScalarAux
    rank_two_tensor = stress
    variable = hoop_stress
    scalar_type = HoopStress
    execute_on = timestep_end
  [../]
[]

[BCs]
  [./x_soln_bc]
    type = FunctionDirichletBC
    variable = disp_x
    preset = false
    boundary = '1 2'
    function = soln_func_disp_x
  [../]
  [./y_soln_bc]
    type = FunctionDirichletBC
    variable = disp_y
    preset = false
    boundary = '1 2'
    function = soln_func_disp_y
  [../]
  [./z_soln_bc]
    type = FunctionDirichletBC
    variable = disp_z
    preset = false
    boundary = '1 2'
    function = soln_func_disp_z
  [../]
[]

[Materials]
  [./elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    block = '1 2'
    youngs_modulus = 193.05e9
    poissons_ratio = 0.3
  [../]

  [./stress]
    type = ComputeFiniteStrainElasticStress
    block = '1 2'
  [../]

  [./thermal_expansion]
    type = ComputeThermalExpansionEigenstrain
    block = '1 2'
    thermal_expansion_coeff = 13e-6
    stress_free_temperature = 295.00
    temperature = temp
    eigenstrain_name = thermal_expansion
  [../]

  [./density]
    type = Density
    block = '1'
    density = 8000.0
  [../]
[]

[Executioner]
  type = Transient
  solve_type = 'PJFNK'

  petsc_options = '-ksp_snes_ew'
  petsc_options_iname = '-ksp_gmres_restart -pc_type -pc_hypre_type -pc_hypre_boomeramg_max_iter'
  petsc_options_value = ' 201                hypre    boomeramg      4'
  line_search = 'none'
  l_max_its = 25
  nl_max_its = 20
  nl_rel_tol = 1e-10
  l_tol = 1e-2

  start_time = 0.0
  dt = 1
  end_time = 1
  dtmin = 1
[]

[Outputs]
  file_base = 3dy_out
  exodus = true
  [./console]
    type = Console
    max_rows = 25
  [../]
[]

Example input syntax
Input Parameters