Description

Base class to set up a Chrono::FSI problem.

#include <ChFsiProblemSPH.h>

Inheritance diagram for chrono::fsi::ChFsiProblemSPH:

[legend]

Collaboration diagram for chrono::fsi::ChFsiProblemSPH:

[legend]

Classes
struct	CoordHash
	Hash function for a 3D integer grid coordinate. More...

struct	FeaMesh
	Specification of an FSI FEA mesh. More...

class	ParticlePropertiesCallback
	Interface for callback to set initial particle pressure, density, viscosity, and velocity. More...

struct	RigidBody
	Specification of an FSI rigid body. More...

Public Types
enum	WavemakerType { PISTON, FLAP }
	Wave generator types.

Public Member Functions
void	SetVerbose (bool verbose)
	Enable verbose output during construction of ChFsiProblemSPH (default: false).

ChFsiSystemSPH &	GetSystemFSI ()
	Access the underlying FSI system.

ChFluidSystemSPH &	GetFluidSystemSPH ()
	Access the underlying SPH system.

ChSystem &	GetMultibodySystem ()
	Access the underlying MBS system.

size_t	AddRigidBody (std::shared_ptr< ChBody > body, const chrono::utils::ChBodyGeometry &geometry, bool check_embedded, bool use_grid_bce=false)
	Add a rigid body to the FSI problem. More...

size_t	AddRigidBodySphere (std::shared_ptr< ChBody > body, const ChVector3d &pos, double radius, bool use_grid_bce=false)

size_t	AddRigidBodyBox (std::shared_ptr< ChBody > body, const ChFramed &pos, const ChVector3d &size)

size_t	AddRigidBodyCylinderX (std::shared_ptr< ChBody > body, const ChFramed &pos, double radius, double length, bool use_grid_bce=false)

size_t	AddRigidBodyMesh (std::shared_ptr< ChBody > body, const ChVector3d &pos, const std::string &obj_file, const ChVector3d &interior_point, double scale)

void	SetBcePattern1D (BcePatternMesh1D pattern, bool remove_center=false)
	Set the BCE marker pattern for 1D flexible solids for subsequent calls to AddFeaMesh. More...

void	SetBcePattern2D (BcePatternMesh2D pattern, bool remove_center=false)
	Set the BCE marker pattern for 2D flexible solids for subsequent calls to AddFeaMesh. More...

size_t	AddFeaMesh (std::shared_ptr< fea::ChMesh > mesh, bool check_embedded)
	Add an FEA mesh to the FSI problem. More...

void	RegisterParticlePropertiesCallback (std::shared_ptr< ParticlePropertiesCallback > callback)
	Register a callback for setting SPH particle initial properties.

void	SetComputationalDomainSize (ChAABB aabb)
	Explicitly set the computational domain limits. More...

void	Initialize ()
	Complete construction of the FSI problem and initialize the FSI system. More...

std::shared_ptr< ChBody >	GetGroundBody () const
	Get the ground body.

size_t	GetNumSPHParticles () const
	Get number of SPH particles.

size_t	GetNumBoundaryBCEMarkers () const
	Get number of boundary BCE markers.

const ChAABB &	GetComputationalDomainSize () const
	Get limits of computational domain.

const ChAABB &	GetSPHBoundingBox () const
	Get limits of SPH volume.

const ChVector3d &	GetFsiBodyForce (std::shared_ptr< ChBody > body) const
	Return the FSI applied force on the specified body (as returned by AddRigidBody). More...

const ChVector3d &	GetFsiBodyTorque (std::shared_ptr< ChBody > body) const
	Return the FSI applied torque on on the specified body (as returned by AddRigidBody). More...

void	SaveInitialMarkers (const std::string &out_dir) const
	Save the set of initial SPH and BCE grid locations to files in the specified output directory.

Protected Types
typedef std::unordered_set< ChVector3i, CoordHash >	GridPoints

typedef std::vector< ChVector3d >	RealPoints

Protected Member Functions
	ChFsiProblemSPH (ChSystem &sys, double spacing)
	Create a ChFsiProblemSPH object. More...

virtual ChVector3i	Snap2Grid (const ChVector3d &point)=0

virtual ChVector3d	Grid2Point (const ChVector3i &p)=0

void	ProcessBody (RigidBody &b)
	Prune SPH markers that are inside the solid body volume. More...

int	ProcessBodyMesh (RigidBody &b, ChTriangleMeshConnected trimesh, const ChVector3d &interior_point)
	Prune SPH markers that are inside a body mesh volume. More...

Protected Attributes
ChFluidSystemSPH	m_sysSPH
	underlying Chrono SPH system

ChFsiSystemSPH	m_sysFSI
	underlying Chrono FSI system

double	m_spacing
	particle and marker spacing

std::shared_ptr< ChBody >	m_ground
	ground body

GridPoints	m_sph
	SPH particle grid locations.

GridPoints	m_bce
	boundary BCE marker grid locations

ChVector3d	m_offset_sph
	SPH particles offset.

ChVector3d	m_offset_bce
	boundary BCE particles offset

ChAABB	m_domain_aabb
	computational domain bounding box

ChAABB	m_sph_aabb
	SPH volume bounding box.

std::vector< RigidBody >	m_bodies
	list of FSI rigid bodies

std::vector< FeaMesh >	m_meshes
	list of FSI FEA meshes

std::unordered_map< std::shared_ptr< ChBody >, size_t >	m_fsi_bodies

std::shared_ptr< ParticlePropertiesCallback >	m_props_cb
	callback for particle properties

bool	m_verbose
	if true, write information to standard output

bool	m_initialized
	set to 'true' once terrain is initialized

Constructor & Destructor Documentation

◆ ChFsiProblemSPH()

chrono::fsi::ChFsiProblemSPH::ChFsiProblemSPH	(	ChSystem &	sys,
		double	spacing
	)

protected

Create a ChFsiProblemSPH object.

No SPH parameters are set.

Member Function Documentation

◆ AddFeaMesh()

size_t chrono::fsi::ChFsiProblemSPH::AddFeaMesh	(	std::shared_ptr< fea::ChMesh >	mesh,
		bool	check_embedded
	)

Add an FEA mesh to the FSI problem.

BCE markers are created based on the type of elements and a corresponding FEA collision surface. To check for possible overlap with SPH particles, set 'check_embedded=true'. This function must be called before Initialize().

◆ AddRigidBody()

size_t chrono::fsi::ChFsiProblemSPH::AddRigidBody	(	std::shared_ptr< ChBody >	body,
		const chrono::utils::ChBodyGeometry &	geometry,
		bool	check_embedded,
		bool	use_grid_bce = `false`
	)

Add a rigid body to the FSI problem.

BCE markers are created for the provided geometry (which may or may not match the body collision geometry). By default, where applicable, BCE markers are created using polar coordinates (in layers starting from the shape surface). Generation of BCE markers on a uniform Cartesian grid can be enforced setting use_grid_bce=true. Creation of FSI bodies embedded in the fluid phase is allowed (SPH markers inside the body geometry volume are pruned). To check for possible overlap with SPH particles, set 'check_embedded=true'. This function must be called before Initialize().

◆ GetFsiBodyForce()

const ChVector3d & chrono::fsi::ChFsiProblemSPH::GetFsiBodyForce ( std::shared_ptr< ChBody > body ) const

Return the FSI applied force on the specified body (as returned by AddRigidBody).

The force is applied at the body COM and is expressed in the absolute frame. An exception is thrown if the given body was not added through AddRigidBody.

◆ GetFsiBodyTorque()

const ChVector3d & chrono::fsi::ChFsiProblemSPH::GetFsiBodyTorque ( std::shared_ptr< ChBody > body ) const

Return the FSI applied torque on on the specified body (as returned by AddRigidBody).

The torque is expressed in the absolute frame. An exception is thrown if the given body was not added through AddRigidBody.

◆ Initialize()

void chrono::fsi::ChFsiProblemSPH::Initialize ( )

Complete construction of the FSI problem and initialize the FSI system.

After this call, no additional solid bodies should be added to the FSI problem.

◆ ProcessBody()

void chrono::fsi::ChFsiProblemSPH::ProcessBody ( RigidBody & b )

protected

Prune SPH markers that are inside the solid body volume.

Treat separately primitive shapes (use explicit test for interior points) and mesh shapes (use ProcessBodyMesh).

◆ ProcessBodyMesh()

int chrono::fsi::ChFsiProblemSPH::ProcessBodyMesh	(	RigidBody &	b,
		ChTriangleMeshConnected	trimesh,
		const ChVector3d &	interior_point
	)

protected

Prune SPH markers that are inside a body mesh volume.

Voxelize the body mesh (at the scaling resolution) and identify grid nodes inside the boundary defined by the body BCEs. Note that this assumes the BCE markers form a watertight boundary.

◆ SetBcePattern1D()

void chrono::fsi::ChFsiProblemSPH::SetBcePattern1D	(	BcePatternMesh1D	pattern,
		bool	remove_center = `false`
	)

Set the BCE marker pattern for 1D flexible solids for subsequent calls to AddFeaMesh.

By default, a full set of BCE markers is used across each section, including a central marker.

Parameters

pattern	marker pattern in cross-section
remove_center	eliminate markers on center line

◆ SetBcePattern2D()

void chrono::fsi::ChFsiProblemSPH::SetBcePattern2D	(	BcePatternMesh2D	pattern,
		bool	remove_center = `false`
	)

Set the BCE marker pattern for 2D flexible solids for subsequent calls to AddFeaMesh.

By default, BCE markers are created centered on the mesh surface, with a layer of BCEs on the surface.

Parameters

pattern	pattern of marker locations along normal
remove_center	eliminate markers on surface

◆ SetComputationalDomainSize()

void chrono::fsi::ChFsiProblemSPH::SetComputationalDomainSize ( ChAABB aabb )

inline

Explicitly set the computational domain limits.

By default, this is set so that it encompasses all SPH particles and BCE markers.

The documentation for this class was generated from the following files:

/builds/uwsbel/chrono/src/chrono_fsi/sph/ChFsiProblemSPH.h
/builds/uwsbel/chrono/src/chrono_fsi/sph/ChFsiProblemSPH.cpp

Description

Classes

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Constructor & Destructor Documentation

◆ ChFsiProblemSPH()

Member Function Documentation

◆ AddFeaMesh()

◆ AddRigidBody()

◆ GetFsiBodyForce()

◆ GetFsiBodyTorque()

◆ Initialize()

◆ ProcessBody()

◆ ProcessBodyMesh()

◆ SetBcePattern1D()

◆ SetBcePattern2D()

◆ SetComputationalDomainSize()