Input Parameters|MOOSE

This class computes the 'Chorin' Predictor equation in fully-discrete (both time and space) form.

Input Parameters

variableThe name of the variable that this Kernel operates on
C++ Type:NonlinearVariableName
Options:
Description:The name of the variable that this Kernel operates on
ux-velocity
C++ Type:std::vector
Options:
Description:x-velocity
component0,1,2 depending on if we are solving the x,y,z component of the Predictor equation
C++ Type:unsigned int
Options:
Description:0,1,2 depending on if we are solving the x,y,z component of the Predictor equation
predictor_typeOne of: OLD, NEW, STAR. Indicates which velocity to use in the predictor.
C++ Type:MooseEnum
Options:OLD NEW STAR
Description:One of: OLD, NEW, STAR. Indicates which velocity to use in the predictor.
u_starstar x-velocity
C++ Type:std::vector
Options:
Description:star x-velocity

Required Parameters

mu_namemuThe name of the dynamic viscosity
Default:mu
C++ Type:MaterialPropertyName
Options:
Description:The name of the dynamic viscosity
rho_namerhoThe name of the density
Default:rho
C++ Type:MaterialPropertyName
Options:
Description:The name of the density
w_starstar z-velocity
C++ Type:std::vector
Options:
Description:star z-velocity
v_starstar y-velocity
C++ Type:std::vector
Options:
Description:star y-velocity
wz-velocity
C++ Type:std::vector
Options:
Description:z-velocity
vy-velocity
C++ Type:std::vector
Options:
Description:y-velocity
displacementsThe displacements
C++ Type:std::vector
Options:
Description:The displacements
blockThe list of block ids (SubdomainID) that this object will be applied
C++ Type:std::vector
Options:
Description:The list of block ids (SubdomainID) that this object will be applied

Optional Parameters

enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
save_inThe name of auxiliary variables to save this Kernel's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector
Options:
Description:The name of auxiliary variables to save this Kernel's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
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
Options:
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.
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Options:
Description:The seed for the master random number generator
diag_save_inThe name of auxiliary variables to save this Kernel's diagonal Jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector
Options:
Description:The name of auxiliary variables to save this Kernel's diagonal Jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Options:
Description:Determines whether this object is calculated using an implicit or explicit form

Advanced Parameters

vector_tagsnontimeThe tag for the vectors this Kernel should fill
Default:nontime
C++ Type:MultiMooseEnum
Options:nontime time
Description:The tag for the vectors this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector
Options:
Description:The extra tags for the vectors this Kernel should fill
matrix_tagssystemThe tag for the matrices this Kernel should fill
Default:system
C++ Type:MultiMooseEnum
Options:nontime system
Description:The tag for the matrices this Kernel should fill
extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector
Options:
Description:The extra tags for the matrices this Kernel should fill

Tagging Parameters

Input Files

modules/navier_stokes/test/tests/ins/lid_driven/lid_driven_chorin.i

modules/navier_stokes/test/tests/ins/lid_driven/lid_driven_chorin.i

[Mesh]
  type = GeneratedMesh
  dim = 2
  xmin = 0
  xmax = 1.0
  ymin = 0
  ymax = 1.0
  nx = 40
  ny = 40
  elem_type = QUAD4
[]

[MeshModifiers]
  [./corner_node]
    type = AddExtraNodeset
    new_boundary = 99
    nodes = '0'
  [../]
[]

[Variables]
  # x-velocity
  [./u]
    order = FIRST
    family = LAGRANGE

    [./InitialCondition]
      type = ConstantIC
      value = 0.0
    [../]
  [../]

  # y-velocity
  [./v]
    order = FIRST
    family = LAGRANGE

    [./InitialCondition]
      type = ConstantIC
      value = 0.0
    [../]
  [../]

  # x-star velocity
  [./u_star]
    order = FIRST
    family = LAGRANGE

    [./InitialCondition]
      type = ConstantIC
      value = 0.0
    [../]
  [../]

  # y-star velocity
  [./v_star]
    order = FIRST
    family = LAGRANGE

    [./InitialCondition]
      type = ConstantIC
      value = 0.0
    [../]
  [../]

  # Pressure
  [./p]
    order = FIRST
    family = LAGRANGE

    [./InitialCondition]
      type = ConstantIC
      value = 0
    [../]
  [../]
[]



[Kernels]
  [./x_chorin_predictor]
    type = INSChorinPredictor
    variable = u_star
    u = u
    v = v
    u_star = u_star
    v_star = v_star
    component = 0
    predictor_type = 'new'
  [../]

  [./y_chorin_predictor]
    type = INSChorinPredictor
    variable = v_star
    u = u
    v = v
    u_star = u_star
    v_star = v_star
    component = 1
    predictor_type = 'new'
  [../]

  [./x_chorin_corrector]
    type = INSChorinCorrector
    variable = u
    u_star = u_star
    v_star = v_star
    p = p
    component = 0
  [../]

  [./y_chorin_corrector]
    type = INSChorinCorrector
    variable = v
    u_star = u_star
    v_star = v_star
    p = p
    component = 1
  [../]

  [./chorin_pressure_poisson]
    type = INSChorinPressurePoisson
    variable = p
    u_star = u_star
    v_star = v_star
  [../]
[]

[BCs]
  [./u_no_slip]
    type = DirichletBC
    variable = u
    boundary = 'bottom right left'
    value = 0.0
  [../]

  [./u_lid]
    type = DirichletBC
    variable = u
    boundary = 'top'
    value = 100.0
  [../]

  [./v_no_slip]
    type = DirichletBC
    variable = v
    boundary = 'bottom right top left'
    value = 0.0
  [../]

  # Make u_star satsify all the same variables as the real velocity.
  [./u_star_no_slip]
    type = DirichletBC
    variable = u_star
    boundary = 'bottom right left'
    value = 0.0
  [../]

  [./u_star_lid]
    type = DirichletBC
    variable = u_star
    boundary = 'top'
    value = 100.0
  [../]

  [./v_star_no_slip]
    type = DirichletBC
    variable = v_star
    boundary = 'bottom right top left'
    value = 0.0
  [../]

  # With solid walls everywhere, we specify dp/dn=0, i.e the
  # "natural BC" for pressure.  Technically the problem still
  # solves without pinning the pressure somewhere, but the pressure
  # bounces around a lot during the solve, possibly because of
  # the addition of arbitrary constants.
  [./pressure_pin]
    type = DirichletBC
    variable = p
    boundary = '99'
    value = 0
  [../]
[]

[Materials]
  [./const]
    type = GenericConstantMaterial
    block = 0
    # rho = 1000    # kg/m^3
    # mu = 0.798e-3 # Pa-s at 30C
    # cp = 4.179e3  # J/kg-K at 30C
    # k = 0.58      # W/m-K at ?C

    # Dummy parameters
    prop_names = 'rho mu cp k'
    prop_values = '1  1  1  1'
  [../]
[]

[Preconditioning]
#active = 'FDP_Newton'
#active = 'SMP_PJFNK'
active = 'SMP_Newton'

[./FDP_Newton]
type = FDP
full = true

solve_type = 'NEWTON'

#petsc_options_iname = '-mat_fd_coloring_err'
#petsc_options_value = '1.e-10'
[../]

# For some reason, nonlinear convergence with JFNK is poor, but it
# seems to be OK for SMP_Newton.  This may indicate a a scaling issue
# in the JFNK case....
[./SMP_PJFNK]
  type = SMP
  full = true

  #Preconditioned JFNK (default)
  solve_type = 'PJFNK'


[../]

[./SMP_Newton]
  type = SMP
  full = true

  solve_type = 'NEWTON'

[../]
[]


[Executioner]
  type = Transient
  # Note: the explicit case with lid velocity = 100 and a 40x40 was unstable
  # for dt=1.e-4, even though the restriction should be dt < dx/|u| = 1/4000 = 2.5e-4
  #
  dt = 1.e-3
  dtmin = 1.e-6
  petsc_options_iname = '-ksp_gmres_restart '
  petsc_options_value = '300                '

  line_search = 'none'

  nl_rel_tol = 1e-5
  nl_max_its = 6
  l_tol = 1e-6
  l_max_its = 300
  start_time = 0.0
  num_steps = 5
[]




[Outputs]
  file_base = lid_driven_chorin_out
  exodus = true
[]

Input Parameters
Input Files

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools

Input Parameters

Required Parameters

Optional Parameters

Advanced Parameters

Tagging Parameters

Input Files