INSFVNoSlipWallBC between subdomains

[j-bowhay]

Hi

I am trying to model fluid flow through two concentric pipes (one pipe inside of another). To do this I have split my mesh into two subdomain and I am trying to apply a no slip BC between the two subdomains as per the diagram bellow.

However when I run the model the bc between the two subdomains doesn't appear to be applied as I am only getting one boundary layer developing.

I have also tried without using separate subdomains but the same issue still persists.

Thanks in advance for any thoughts on how to achieve this.
Jake

Input File

# 2d siumulation of a water through a pipe.

mu=1e-3
rho=997.0
Re_inner=1000.0
pipe_length=1
inner_radius=0.1
outer_radius=0.1
u_inner_inlet=${fparse (mu * Re_inner)/(2 * inner_radius * rho)}
advected_interp_method='average'
velocity_interp_method='average'

[GlobalParams]
  rhie_chow_user_object = 'rc'
[]

[UserObjects]
  [rc]
    type = INSFVRhieChowInterpolator
    u = u
    v = v
    pressure = pressure
    block = '0 1'
  []
[]

[Mesh]
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    xmin = 0
    xmax = ${pipe_length}
    ymin = 0
    ymax = ${fparse inner_radius + outer_radius}
    nx = 100
    ny = 150
  []
  [id_outer]
    type = ParsedSubdomainMeshGenerator
    input = gen
    combinatorial_geometry = 'y > ${inner_radius}'
    block_id = 1
  []
  [delete_sideset]
    type = BoundaryDeletionGenerator
    input = id_outer
    boundary_names = 'left right top bottom'
  []
  [name_outer_top]
    type = ParsedGenerateSideset
    input = delete_sideset
    combinatorial_geometry = 'y = ${fparse inner_radius + outer_radius}'
    replace = true
    new_sideset_name = 'outer_top'
  []
  [name_outer_left]
    type = ParsedGenerateSideset
    input = name_outer_top
    combinatorial_geometry = 'y > ${inner_radius} & x = 0'
    replace = true
    new_sideset_name = 'outer_left'
  []
  [name_outer_right]
    type = ParsedGenerateSideset
    input = name_outer_left
    combinatorial_geometry = 'y > ${inner_radius} & x = ${pipe_length}'
    replace = true
    new_sideset_name = 'outer_right'
  []
  [middle_interface]
    type = SideSetsBetweenSubdomainsGenerator
    input = name_outer_right
    primary_block = '0'
    paired_block = '1'
    new_boundary = 'middle_interface'
  []
  [name_inner_bottom]
    type = ParsedGenerateSideset
    input = middle_interface
    combinatorial_geometry = 'y = 0'
    replace = true
    new_sideset_name = 'inner_bottom'
  []
  [name_inner_left]
    type = ParsedGenerateSideset
    input = name_inner_bottom
    combinatorial_geometry = 'y < ${inner_radius} & x = 0'
    replace = true
    new_sideset_name = 'inner_left'
  []
  [name_inner_right]
    type = ParsedGenerateSideset
    input = name_inner_left
    combinatorial_geometry = 'y < ${inner_radius} & x = ${pipe_length}'
    replace = true
    new_sideset_name = 'inner_right'
  []
[]

[Problem]
  fv_bcs_integrity_check = true
  coord_type = 'RZ'
  rz_coord_axis = x
[]

[Variables]
  [u]
    type = INSFVVelocityVariable
    initial_condition = ${u_inner_inlet}
    block = '0 1'
  []
  [v]
    type = INSFVVelocityVariable
    initial_condition = 1e-15
    block = '0 1'
  []
  [pressure]
    type = INSFVPressureVariable
    block = '0 1'
  []
[]

[FVKernels]
  [mass]
    type = INSFVMassAdvection
    variable = pressure
    advected_interp_method = 'average'#${advected_interp_method}
    velocity_interp_method = 'average'#${velocity_interp_method}
    rho = ${rho}
    block = '0 1'
  []

  [u_advection]
    type = INSFVMomentumAdvection
    variable = u
    advected_interp_method = 'average'#${advected_interp_method}
    velocity_interp_method = 'average'#${velocity_interp_method}
    rho = ${rho}
    momentum_component = 'x'
    block = '0 1'
  []
  [u_viscosity]
    type = INSFVMomentumDiffusion
    variable = u
    mu = ${mu}
    momentum_component = 'x'
    block = '0 1'
  []
  [u_pressure]
    type = INSFVMomentumPressure
    variable = u
    momentum_component = 'x'
    pressure = pressure
    block = '0 1'
  []

  [v_advection]
    type = INSFVMomentumAdvection
    variable = v
    advected_interp_method = 'average'#${advected_interp_method}
    velocity_interp_method = 'average'#${velocity_interp_method}
    rho = ${rho}
    momentum_component = 'y'
    block = '0 1'
  []
  [v_viscosity]
    type = INSFVMomentumDiffusion
    variable = v
    mu = ${mu}
    momentum_component = 'y'
    block = '0 1'
  []
  [v_pressure]
    type = INSFVMomentumPressure
    variable = v
    momentum_component = 'y'
    pressure = pressure
    block = '0 1'
  []
[]

[FVBCs]
  [outer-inlet-u]
    type = INSFVInletVelocityBC
    boundary = 'outer_left'
    variable = u
    function = ${u_inner_inlet}
  []
  [outer-inlet-v]
    type = INSFVInletVelocityBC
    boundary = 'outer_left'
    variable = v
    function = 0
  []
  [outer-top-no-slip-wall-u]
    type = INSFVNoSlipWallBC
    boundary = 'outer_top'
    variable = u
    function = 0
  []
  [outer-top-no-slip-wall-v]
    type = INSFVNoSlipWallBC
    boundary = 'outer_top'
    variable = v
    function = 0
  []
  [outer-outlet-p]
    type = INSFVOutletPressureBC
    boundary = 'outer_right'
    variable = pressure
    function = 1e-15
  []
  [no-slip-middle-u]
    type = INSFVNoSlipWallBC
    boundary = 'middle_interface'
    variable = u
    function = 0
  []
  [no-slip-middle-v]
    type = INSFVNoSlipWallBC
    boundary = 'middle_interface'
    variable = v
    function = 0
  [] 
  [inner-inlet-u]
    type = INSFVInletVelocityBC
    boundary = 'inner_left'
    variable = u
    function = ${u_inner_inlet}
  []
  [inner-inlet-v]
    type = INSFVInletVelocityBC
    boundary = 'inner_left'
    variable = v
    function = 0
  []
  [inner-outlet-p]
    type = INSFVOutletPressureBC
    boundary = 'inner_right'
    variable = pressure
    function = 1e-15
  []
  [axis-u]
    type = INSFVSymmetryVelocityBC
    boundary = 'inner_bottom'
    variable = u
    u = u
    v = v
    mu = ${mu}
    momentum_component = x
  []
  [axis-v]
    type = INSFVSymmetryVelocityBC
    boundary = 'inner_bottom'
    variable = v
    u = u
    v = v
    mu = ${mu}
    momentum_component = y
  []
  [axis-p]
    type = INSFVSymmetryPressureBC
    boundary = 'inner_bottom'
    variable = pressure
  []
[]

[Executioner]
  type = Steady
  solve_type = 'NEWTON'
  petsc_options_iname = '-pc_type -ksp_gmres_restart -sub_pc_type -sub_pc_factor_shift_type'
  petsc_options_value = 'asm      100                lu           NONZERO'
  line_search = 'none'
  nl_max_its = 50       
  l_max_its = 50
  nl_rel_tol = 1e-12
  automatic_scaling = true
[]

[Outputs]
  exodus = true
[]

[GiudGiud]

Hello

We do not support boundary conditions inside the flow domains. You will have to either:

• cut the domain by deleting an element or generating two separate meshes (you can keep them in the same solve using the CombinerGenerator
  https://mooseframework.inl.gov/source/meshgenerators/CombinerGenerator.html
• use different flow variables (u,v,P) in both blocks.
• use my branch where I'm working on interface kernels and a boundaries_to_avoid parameter to restrict the execution of kernels on some sidesets. But I havent figured out everything out for this yet, and it's being done for restricting flow with some sort of diode behavior, not for adding a wall (though it's very similar)
  https://github.com/giudgiud/moose/tree/PR_wall_flow_diode

btw average velocity and advected quantity interpolation are unstable for fluid flow, so you may want to use Rhie Chow for velocity and first order upwind or try our newly merged second order schemes for the advected quantity.

Guillaume

[j-bowhay]

Thank you for the helpful response