3-D fractional single-phase pressure solver. More...

#include <jet/grid_fractional_single_phase_pressure_solver3.h>

Inheritance diagram for jet::GridFractionalSinglePhasePressureSolver3:

Public Member Functions
	GridFractionalSinglePhasePressureSolver3 ()

virtual	~GridFractionalSinglePhasePressureSolver3 ()

void	solve (const FaceCenteredGrid3 &input, double timeIntervalInSeconds, FaceCenteredGrid3 *output, const ScalarField3 &boundarySdf=ConstantScalarField3(kMaxD), const VectorField3 &boundaryVelocity=ConstantVectorField3({0, 0, 0}), const ScalarField3 &fluidSdf=ConstantScalarField3(-kMaxD), bool useCompressed=false) override
	Solves the pressure term and apply it to the velocity field. More...

GridBoundaryConditionSolver3Ptr	suggestedBoundaryConditionSolver () const override
	Returns the best boundary condition solver for this solver. More...

const FdmLinearSystemSolver3Ptr &	linearSystemSolver () const
	Returns the linear system solver. More...

void	setLinearSystemSolver (const FdmLinearSystemSolver3Ptr &solver)
	Sets the linear system solver. More...

const FdmVector3 &	pressure () const
	Returns the pressure field. More...

Public Member Functions inherited from jet::GridPressureSolver3
	GridPressureSolver3 ()
	Default constructor. More...

virtual	~GridPressureSolver3 ()
	Default destructor. More...

Detailed Description

3-D fractional single-phase pressure solver.

This class implements 3-D fractional (or variational) single-phase pressure solver. It is called fractional because the solver encodes the boundaries to the grid cells like anti-aliased pixels, meaning that a grid cell will record the partially overlapping boundary as a fractional number. Alternative apporach is to represent boundaries like Lego blocks which is the case for GridSinglePhasePressureSolver2. In addition, this class solves single-phase flow, solving the pressure for selected fluid region only and treat other area as an atmosphere region. Thus, the pressure outside the fluid will be set to a constant value and velocity field won't be altered. This solver also computes the fluid boundary in fractional manner, meaning that the solver tries to capture the subgrid structures. This class uses ghost fluid method for such calculation.

See also: Batty, Christopher, Florence Bertails, and Robert Bridson. "A fast variational framework for accurate solid-fluid coupling." ACM Transactions on Graphics (TOG). Vol. 26. No. 3. ACM, 2007.; Enright, Doug, et al. "Using the particle level set method and a second order accurate pressure boundary condition for free surface flows." ASME/JSME 2003 4th Joint Fluids Summer Engineering Conference. American Society of Mechanical Engineers, 2003.

Constructor & Destructor Documentation

◆ GridFractionalSinglePhasePressureSolver3()

jet::GridFractionalSinglePhasePressureSolver3::GridFractionalSinglePhasePressureSolver3 ( )

◆ ~GridFractionalSinglePhasePressureSolver3()

virtual jet::GridFractionalSinglePhasePressureSolver3::~GridFractionalSinglePhasePressureSolver3 ( )

virtual

Member Function Documentation

◆ linearSystemSolver()

const FdmLinearSystemSolver3Ptr& jet::GridFractionalSinglePhasePressureSolver3::linearSystemSolver ( ) const

Returns the linear system solver.

◆ pressure()

const FdmVector3& jet::GridFractionalSinglePhasePressureSolver3::pressure ( ) const

Returns the pressure field.

◆ setLinearSystemSolver()

void jet::GridFractionalSinglePhasePressureSolver3::setLinearSystemSolver ( const FdmLinearSystemSolver3Ptr & solver )

Sets the linear system solver.

◆ solve()

void jet::GridFractionalSinglePhasePressureSolver3::solve	(	const FaceCenteredGrid3 &	input,
		double	timeIntervalInSeconds,
		FaceCenteredGrid3 *	output,
		const ScalarField3 &	boundarySdf = `ConstantScalarField3(kMaxD)`,
		const VectorField3 &	boundaryVelocity = `ConstantVectorField3({0, 0, 0})`,
		const ScalarField3 &	fluidSdf = `ConstantScalarField3(-kMaxD)`,
		bool	useCompressed = `false`
	)

overridevirtual

Solves the pressure term and apply it to the velocity field.

This function takes input velocity field and outputs pressure-applied velocity field. It also accepts extra arguments such as boundarySdf and fluidSdf that represent signed-distance representation of the boundary and fluid area. The negative region of boundarySdf means it is occupied by solid object. Also, the positive / negative area of the fluidSdf means it is occupied by fluid / atmosphere. If not specified, constant scalar field with kMaxD will be used for boundarySdf meaning that no boundary at all. Similarly, a constant field with -kMaxD will be used for fluidSdf which means it's fully occupied with fluid without any atmosphere.

Parameters

[in]	input	The input velocity field.
[in]	timeIntervalInSeconds	The time interval for the sim.
[in,out]	output	The output velocity field.
[in]	boundarySdf	The SDF of the boundary.
[in]	fluidSdf	The SDF of the fluid/atmosphere.
[in]	useCompressed	True if it uses compressed system.

Implements jet::GridPressureSolver3.

◆ suggestedBoundaryConditionSolver()

GridBoundaryConditionSolver3Ptr jet::GridFractionalSinglePhasePressureSolver3::suggestedBoundaryConditionSolver ( ) const