immersed_boundaries_markers Module

Treatment of surface markers for direct-forcing immersed boundary methods


Uses

  • module~~immersed_boundaries_markers~~UsesGraph module~immersed_boundaries_markers immersed_boundaries_markers module~leapblock leapBlock module~immersed_boundaries_markers->module~leapblock module~leapdiffop leapDiffOp module~immersed_boundaries_markers->module~leapdiffop module~leapeulerian leapEulerian module~immersed_boundaries_markers->module~leapeulerian module~leaphypre leapHypre module~immersed_boundaries_markers->module~leaphypre module~leapio leapIO module~immersed_boundaries_markers->module~leapio module~leapkinds leapKinds module~immersed_boundaries_markers->module~leapkinds module~leaplagrangian leapLagrangian module~immersed_boundaries_markers->module~leaplagrangian module~leapmonitor leapMonitor module~immersed_boundaries_markers->module~leapmonitor module~leapparallel leapParallel module~immersed_boundaries_markers->module~leapparallel module~leapparser leapParser module~immersed_boundaries_markers->module~leapparser module~leaptimer leapTimer module~immersed_boundaries_markers->module~leaptimer module~leapblock->module~leapkinds module~leapblock->module~leapparallel mpi_f08 mpi_f08 module~leapblock->mpi_f08 module~leapdiffop->module~leapblock module~leapdiffop->module~leapeulerian module~leapdiffop->module~leapkinds module~leapdiffop->module~leapparallel module~leapbc leapBC module~leapdiffop->module~leapbc module~leapeulerian->module~leapblock module~leapeulerian->module~leapio module~leapeulerian->module~leapkinds module~leapeulerian->module~leapparallel module~leapeulerian->module~leapparser module~leapio_mpi leapIO_mpi module~leapeulerian->module~leapio_mpi module~leaputils leapUtils module~leapeulerian->module~leaputils module~leapeulerian->mpi_f08 module~leaphypre->module~leapblock module~leaphypre->module~leapeulerian module~leaphypre->module~leapkinds module~leaphypre->module~leapparallel iso_c_binding iso_c_binding module~leaphypre->iso_c_binding module~leaphypre->mpi_f08 module~leapio_h5hut leapIO_h5hut module~leapio->module~leapio_h5hut module~leapio_hdf5 leapIO_hdf5 module~leapio->module~leapio_hdf5 module~leapio->module~leapio_mpi module~leapio_silo leapIO_silo module~leapio->module~leapio_silo iso_fortran_env iso_fortran_env module~leapkinds->iso_fortran_env module~leaplagrangian->module~leapblock module~leaplagrangian->module~leapio module~leaplagrangian->module~leapkinds module~leaplagrangian->module~leapparallel module~leaplagrangian->module~leapparser module~leaplagrangian->module~leaputils module~leaplagrangian->mpi_f08 module~leapmonitor->module~leapkinds module~leapmonitor->module~leapparallel module~leapmonitor->module~leaputils module~leapparallel->module~leapkinds module~leapparallel->module~leapparser module~leapparallel->mpi_f08 module~leapparser->module~leapkinds module~leapcli leapCli module~leapparser->module~leapcli module~leaptimer->module~leapkinds module~leaptimer->module~leapparser module~leapbc->module~leapblock module~leapbc->module~leapeulerian module~leapbc->module~leapkinds module~leapbc->module~leapparallel module~leapbc->module~leaputils module~leapbc->mpi_f08 module~leapcli->module~leapkinds module~leapio_h5hut->module~leapkinds module~leapio_h5hut->module~leapparallel h5hut h5hut module~leapio_h5hut->h5hut module~leapio_hdf5->module~leapkinds module~leapio_hdf5->module~leapparallel module~leapio_hdf5->module~leaputils hdf5 hdf5 module~leapio_hdf5->hdf5 module~leapio_mpi->module~leapblock module~leapio_mpi->module~leapkinds module~leapio_mpi->module~leapparallel module~leapio_mpi->mpi_f08 module~leapio_silo->module~leapblock module~leapio_silo->module~leapkinds module~leapio_silo->module~leapparallel module~leapio_silo->module~leaputils module~leapio_silo->mpi_f08 module~leaputils->module~leapkinds module~leaputils->mpi_f08

Used by

  • module~~immersed_boundaries_markers~~UsedByGraph module~immersed_boundaries_markers immersed_boundaries_markers module~immersed_boundaries immersed_boundaries module~immersed_boundaries->module~immersed_boundaries_markers module~immersed_boundaries_solids immersed_boundaries_solids module~immersed_boundaries->module~immersed_boundaries_solids module~immersed_boundaries_solids->module~immersed_boundaries_markers module~cdifs cdifs module~cdifs->module~immersed_boundaries module~particles_resolved particles_resolved module~cdifs->module~particles_resolved module~grans grans module~grans->module~immersed_boundaries module~grans_cases grans_cases module~grans->module~grans_cases module~grans->module~particles_resolved module~grans_cases->module~immersed_boundaries module~grans_cases->module~particles_resolved module~particles_resolved->module~immersed_boundaries proc~add_markers add_markers proc~add_markers->module~immersed_boundaries proc~add_markers->module~particles_resolved proc~add_markers~2 add_markers proc~add_markers~2->module~immersed_boundaries proc~add_markers~2->module~particles_resolved proc~cdifs_cylinder_ib cdifs_cylinder_IB proc~cdifs_cylinder_ib->module~immersed_boundaries proc~cdifs_rebound_ib cdifs_rebound_IB proc~cdifs_rebound_ib->module~immersed_boundaries proc~cdifs_settling_sphere_ib cdifs_settling_sphere_IB proc~cdifs_settling_sphere_ib->module~immersed_boundaries module~cdifs_advancesolution_smod cdifs_AdvanceSolution_smod module~cdifs_advancesolution_smod->module~cdifs module~cdifs_monitor_smod cdifs_Monitor_smod module~cdifs_monitor_smod->module~cdifs module~cdifs_preparesolver_smod cdifs_PrepareSolver_smod module~cdifs_preparesolver_smod->module~cdifs module~cdifs_writeoutputdata_smod cdifs_WriteOutputData_smod module~cdifs_writeoutputdata_smod->module~cdifs module~cdifs_writerestartdata_smod cdifs_WriteRestartData_smod module~cdifs_writerestartdata_smod->module~cdifs module~grans_advancesolution_smod grans_AdvanceSolution_smod module~grans_advancesolution_smod->module~grans module~grans_module_smod grans_module_smod module~grans_module_smod->module~grans module~grans_preparesolver_smod grans_PrepareSolver_smod module~grans_preparesolver_smod->module~grans module~grans_writeoutputdata_smod grans_WriteOutputData_smod module~grans_writeoutputdata_smod->module~grans module~grans_writerestartdata_smod grans_WriteRestartData_smod module~grans_writerestartdata_smod->module~grans module~headon_collision_smod headon_collision_smod module~headon_collision_smod->module~grans_cases module~homogeneous_cooling_smod homogeneous_cooling_smod module~homogeneous_cooling_smod->module~grans_cases module~rotating_drum_smod rotating_drum_smod module~rotating_drum_smod->module~grans_cases proc~cdifs_rebound_respart cdifs_rebound_ResPart proc~cdifs_rebound_respart->module~particles_resolved proc~cdifs_settling_sphere_respart cdifs_settling_sphere_ResPart proc~cdifs_settling_sphere_respart->module~particles_resolved program~main main program~main->module~cdifs program~main->module~grans

Derived Types

type, public, extends(lagrangian_obj) ::  marker_obj

An extended Lagrangian object that represents a marker placed at the centroid of a triangular element on the surface of an immersed body

Components

Type Visibility Attributes Name Initial
real(kind=WP), public :: SA

Sufrace area
integer, public :: c(3)

nearest cell
real(kind=WP), public :: f(3)

marker forcing
integer(kind=leapI8), public :: id

Identifying number (inactive if <0)
real(kind=WP), public :: n(3)

normal at the marker
real(kind=WP), public :: p(3)

position
real(kind=WP), public :: pold(3)

Old position
integer, public :: s

A tag
integer, public :: t

Another tag
real(kind=WP), public :: v(3)

marker velocity
real(kind=WP), public :: vold(3)

Old velocity

Type-Bound Procedures

procedure, public :: Extrapolate => lagrangian_obj_Extrapolate
procedure, public :: Info => marker_obj_Info
procedure, public :: Interpolate => lagrangian_obj_Interpolate
procedure, public :: Locate => lagrangian_obj_Locate
procedure, public :: assign => marker_obj_assign
generic, public :: assignment(=) => assign

type, public, extends(lagrangian_set) ::  marker_set

A collection of IB markers

Components

Type Visibility Attributes Name Initial
integer, public :: MPI_SIZE = 44

MPI size
type(MPI_Datatype), public :: MPI_TYPE

MPI variable type
type(block_obj), public, pointer :: block => null()

Associated block structure
type(block_obj), public :: cblock

Collision block
integer, public :: count = 0

Total count across all MPI ranks
integer, public :: count_ = 0

Local count for this rank
integer, public, allocatable :: count_proc(:) Read more…
procedure(kernel_1D), public, nopass, pointer :: g1ex => int_g1_triangle

1D kernel used in extrapolations
procedure(kernel_1D), public, nopass, pointer :: g1in => g1_triangle

1D kernel used in interpolations
real(kind=wp), public :: l_filter

Half filter width
type(monitor_set), public, pointer :: monitors => null()

Monitors to print to stdout and files
character(len=:), public, allocatable :: name

Name of the Lagrangian set
type(sllist_obj), public, allocatable :: neighbors(:,:,:)

Singly linked list of neighbors
integer, public, allocatable :: objincell(:,:,:)

Number of objects in this list
type(op_obj), public, pointer :: op => null()

operators object
logical, public :: overwrite = .true.

Switch to overwrite IO files
class(lagrangian_obj), public, allocatable :: p(:)

Array of Lagrangian_obj or any extended type
type(parallel_obj), public, pointer :: parallel => null()

Associated parallel structure
type(parser_obj), public, pointer :: parser => null()

Parser for input file
character(len=str64), public :: read_file

file to read
class(lagrangian_obj), public, allocatable :: sample

Sample used in allocation of polymorphic data
integer, public :: stib = 3

Stencil size for filtering
type(timer_obj), public, pointer :: timer => null()

Timer utility
character(len=str64), public :: write_file

file to write

Type-Bound Procedures

procedure, public :: ApplyPeriodicity => lagrangian_set_ApplyPeriodicity
procedure, public :: CoM => marker_set_CoM
procedure, public :: Communicate => lagrangian_set_Communicate
procedure, public :: ComputeSolidVolFrac => marker_set_ComputeSolidVolFrac
procedure, public :: CreateMPIType => lagrangian_set_CreateMPIType
procedure, public :: Filter => marker_set_Filter
generic, public :: Finalize => lagrangian_set_Final
procedure, public :: FreeNeighborList => lagrangian_set_FreeNeighborList
procedure, public :: GetIBForcing => marker_set_GetIBForcing
procedure, public :: GetOwnerRankByBlock => lagrangian_set_GetOwnerRankByBlock
procedure, public :: GetReadFileName => lagrangian_set_GetReadFileName
procedure, public :: GetWriteFileName => lagrangian_set_GetWriteFileName
procedure, public :: Info => lagrangian_set_Info
generic, public :: Initialize => lagrangian_set_Init
procedure, public :: LoadSTL => marker_set_LoadSTL
procedure, public :: Localize => lagrangian_set_Localize
procedure, public :: Prepare => marker_set_Prepare
procedure, public :: Read => marker_set_Read
procedure, public :: Recycle => lagrangian_set_Recycle
procedure, public :: Resize => lagrangian_set_Resize
procedure, public :: SetFilterKernel => lagrangian_set_SetFilterKernel
procedure, public :: SetFilterSize => lagrangian_set_SetFilterSize
procedure, public :: SetMPIDataTypeParams => marker_set_SetMPIDataTypeParams
procedure, public :: SetObjectType => marker_set_SetObjectType
procedure, public :: SetOverwrite => lagrangian_set_SetOverwrite
procedure, public :: SetReadFileName => lagrangian_set_SetReadFileName
procedure, public :: SetWriteFileName => lagrangian_set_SetWriteFileName
procedure, public :: SetupCollisionBlock => lagrangian_set_SetupCollisionBlock
procedure, public :: UpdateCount => lagrangian_set_UpdateCount
procedure, public :: UpdateGhostObjects => lagrangian_set_UpdateGhostObjects
procedure, public :: UpdateNeighborList => lagrangian_set_UpdateNeighborList
procedure, public :: UpdateNormals => marker_set_UpdateNormals
procedure, public :: UpdateSDF => marker_set_UpdateSDF
procedure, public :: Write => marker_set_Write
procedure, public :: lagrangian_set_Final
procedure, public :: lagrangian_set_Init

Functions

private function marker_set_CoM(this) result(CoM)

Find the center of mass

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers

Return Value real(kind=wp), (3)

Position of center of mass


Subroutines

private subroutine marker_obj_Info(this)

Print information about this marker

Arguments

Type IntentOptional Attributes Name
class(marker_obj), intent(inout) :: this

A marker_obj object

private subroutine marker_obj_assign(this, val)

Assignment

Arguments

Type IntentOptional Attributes Name
class(marker_obj), intent(inout) :: this

A marker_obj object
class(lagrangian_obj), intent(in) :: val

A marker_obj object

private subroutine marker_set_ComputeSolidVolFrac(this, VF, solver, MaxTol, MaxIt, RelaxType, Rel, It, intRHS)

Compute the solid volume fraction on the mesh

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers
type(eulerian_obj_r), intent(inout) :: VF

Volume fraction
character(len=*), intent(in) :: solver

Name of solver to be used
real(kind=wp), intent(in) :: MaxTol

Maximum relative tolerance
integer, intent(in) :: MaxIt

Maximum number of subiterations
integer, intent(in), optional :: RelaxType

Relaxation type
real(kind=wp), intent(out), optional :: Rel

Relative error at end of solve
integer, intent(out), optional :: It

Number of iterations performed
real(kind=wp), intent(out), optional :: intRHS

Magnitude of RHS

private subroutine marker_set_Filter(this, var, field)

Compute a filtered quantity on the eulerian grid

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers
character(len=*), intent(in) :: var

Variable to compute
type(eulerian_obj_r), intent(inout) :: field

Filtered quantity

private subroutine marker_set_GetIBForcing(this, Um, Vm, Wm, rhof, SA, ibF, dt)

Compute the IB forcing

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Collection of Resolved Particles
type(eulerian_obj_r), intent(in) :: Um

Velocity in 1-dir
type(eulerian_obj_r), intent(in) :: Vm

Velocity in 2-dir
type(eulerian_obj_r), intent(in) :: Wm

Velocity in 3-dir
real(kind=wp), intent(in) :: rhof

Fluid density
type(eulerian_obj_r), intent(in) :: SA

Surface area
type(eulerian_obj_r), intent(inout) :: ibF(3)

IB forcing
real(kind=wp), intent(in) :: dt

Timestep

private subroutine marker_set_LoadSTL(this, STL_file)

Load markers from a binary STL. This is a serial routine.

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers
character(len=*), intent(in) :: STL_file

Variable to compute

private subroutine marker_set_Prepare(this, timer, parser, operators, monitors)

Prepare marker_set for use with solvers

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Immersed boundary
type(timer_obj), intent(in), target :: timer

Timer utility
type(parser_obj), intent(in), target :: parser

Parser for input file
type(op_obj), intent(in), target :: operators

Operators object
type(monitor_set), intent(in), target :: monitors

Monitors to print to stdout and files

private subroutine marker_set_Read(this, iter, time)

Read marker data from file in parallel

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers
integer, intent(out) :: iter

Iteration at write
real(kind=wp), intent(out) :: time

Time at write

private subroutine marker_set_SetMPIDataTypeParams(this, types, lengths, displacement)

Set up parameters used when creating the MPI derived type Create the MPI structure

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers
type(MPI_Datatype), intent(out), allocatable :: types(:)

Array of types
integer, intent(out), allocatable :: lengths(:)

Array of lengths
integer(kind=MPI_ADDRESS_KIND), intent(out), allocatable :: displacement(:)

Array of displacements

private subroutine marker_set_SetObjectType(this)

Set the sample type used in allocation of polymorphic variables

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers

private subroutine marker_set_UpdateNormals(this, ibN)

Updates the Normals field

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Immersed boundary
type(eulerian_obj_r), intent(inout) :: ibN(3)

private subroutine marker_set_UpdateSDF(this, SA)

Updates the Surface Density Function

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Immersed boundary
type(eulerian_obj_r), intent(inout) :: SA

private subroutine marker_set_Write(this, iter, time)

Write marker data to file in parallel

Arguments

Type IntentOptional Attributes Name
class(marker_set), intent(inout) :: this

Set of ib markers
integer, intent(in) :: iter

Iteration at write
real(kind=wp), intent(in) :: time

Time at write