Skip to content

structure-fluid smoothing length and search radius are not the same #923

@efaulhaber

Description

@efaulhaber

In structure-fluid interaction, the kernel and smoothing length of the fluid is used:

# Use kernel from the fluid system in order to get the same force here in
# structure-fluid interaction as for fluid-structure interaction.
# TODO this will not use corrections if the fluid uses corrections.
grad_kernel = smoothing_kernel_grad(neighbor_system, pos_diff, distance, particle)

However, the compact support and hence the search radius for this interaction depends on the boundary model kernel:

@inline function compact_support(system::Union{TotalLagrangianSPHSystem,
WallBoundarySystem},
neighbor)
return compact_support(system, system.boundary_model, neighbor)
end

This causes incorrect kernel truncation if one uses a smaller smoothing length for the boundary model than for the fluid.

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions