Struct for reference values of MPS method. More...

#include <refvalues.hpp>

Collaboration diagram for RefValues:

Public Member Functions
	RefValues ()=default

	RefValues (int dim, double particleDistance, double re, WeightType weightType)
	calculate reference values for MPS method

Public Attributes
double	n0 = 0
	reference value of number density for source term of pressure Poisson equation

double	lambda = 0
	coefficient for laplacian

Private Attributes
Weight	weight
	weight function instance

Detailed Description

Struct for reference values of MPS method.

Definition at line 9 of file refvalues.hpp.

Constructor & Destructor Documentation

◆ RefValues() [1/2]

RefValues::RefValues ( )

default

◆ RefValues() [2/2]

RefValues::RefValues	(	int	dim,
		double	particleDistance,
		double	re,
		WeightType	weightType )

calculate reference values for MPS method

Parameters

dim	dimension of the simulation
particleDistance	initial particle distance
re	effective radius
weightType	weight function type

Definition at line 9 of file refvalues.cpp.

                                                                                       {
    assert(dim == 2 || dim == 3);
    assert(particleDistance < re);
 
    this->weight = Weight(weightType);
 
    int iZ_start = -4;
    int iZ_end   = 5;
    if (dim == 2) {
        iZ_start = 0;
        iZ_end   = 1;
    }
 
    this->n0     = 0.0;
    this->lambda = 0.0;
    for (int iX = -4; iX < 5; iX++) {
        for (int iY = -4; iY < 5; iY++) {
            for (int iZ = iZ_start; iZ < iZ_end; iZ++) {
                if (((iX == 0) && (iY == 0)) && (iZ == 0))
                    continue;
 
                double xj   = particleDistance * (double) (iX);
                double yj   = particleDistance * (double) (iY);
                double zj   = particleDistance * (double) (iZ);
                double dis2 = xj * xj + yj * yj + zj * zj;
                double dis  = sqrt(dis2);
                n0 += weight.weightFn(dis, re);
                lambda += dis2 * weight.weightFn(dis, re);
            }
        }
    }
    this->lambda /= this->n0;
}

Here is the call graph for this function:

Member Data Documentation

◆ n0

double RefValues::n0 = 0

reference value of number density for source term of pressure Poisson equation

Definition at line 11 of file refvalues.hpp.

◆ lambda

double RefValues::lambda = 0

coefficient for laplacian

Definition at line 12 of file refvalues.hpp.

◆ weight

Weight RefValues::weight

private

weight function instance

Definition at line 26 of file refvalues.hpp.

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

src/refvalues.hpp
src/refvalues.cpp

Public Member Functions

Public Attributes

Private Attributes