1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
// SPDX-FileCopyrightText: Copyright (c) Kitware Inc.
// SPDX-License-Identifier: BSD-3-Clause
/**
* @class vtkLagrangianIntegrationModelExample
* @brief Integration model example
*
* Similar to vtkLagrangianMatidaIntegrationModel
* this integration model demonstrates the capabilities
* of Lagrangian Integration Models
* The mains differences with Matida version are
* it uses a gravity constant from the data, instead of using G,
* and it also adds new variable, the particle diameter,
* and uses this diameter to make choice during interaction and free flight
* Please consult vtkLagrangianBasicIntegrationModel and vtkLagrangianMatidaIntegrationModel
* for more explanation on how this example works.
*/
#ifndef vtkLagrangianIntegrationModelExample_h
#define vtkLagrangianIntegrationModelExample_h
#include "LagrangianExampleModule.h" // for export macro
#include "vtkLagrangianBasicIntegrationModel.h"
class LAGRANGIANEXAMPLE_EXPORT vtkLagrangianIntegrationModelExample
: public vtkLagrangianBasicIntegrationModel
{
public:
vtkTypeMacro(vtkLagrangianIntegrationModelExample, vtkLagrangianBasicIntegrationModel);
void PrintSelf(ostream& os, vtkIndent indent) override;
static vtkLagrangianIntegrationModelExample* New();
using Superclass::FunctionValues;
/**
* Evaluate the integration model velocity field
* f at position x, causing data from cell in dataSet with index cellId
* This method is pure abstract at vtkLagrangianBasicIntegrationModel
* THIS IS THE MAIN METHOD TO BE DEFINED IN A LAGRANGIAN INTEGRATION MODEL PLUGIN
*/
int FunctionValues(vtkLagrangianParticle* particle, vtkDataSet* dataSet, vtkIdType cellId,
double* weights, double* x, double* f) override;
/**
* This method is called each time a particle created from the seeds
* It should be inherited in order to initialize variable data in user variables
* from seed data
*/
void InitializeParticle(vtkLagrangianParticle* particle) override;
/**
* This method is called when initializing output point data
* It should be inherited when there is some variables data needed to be put
* in output field data.
* Add some User Variable Data Array in provided particleData, allocate
* maxTuples tuples.
*/
void InitializeParticleData(vtkFieldData* particleData, int maxTuples = 0) override;
/**
* This method is called when inserting particle data in output point data
* It should be inherited when there is some variables data needed to be put in
* output field data.
* Insert user variables data in provided point data, user variables data array begins at
* arrayOffset. use stepEnum to identify which step ( prev, current or next ) should be inserted.
*/
void InsertParticleData(
vtkLagrangianParticle* particle, vtkFieldData* data, int stepEnum) override;
/**
* This method is called when checking if a particle should be terminated in free flight
* At vtkLagrangianBasicIntegrationModel this method does nothing
* Return true if particle is terminated, false otherwise
*/
bool CheckFreeFlightTermination(vtkLagrangianParticle* particle) override;
/**
* Methods used by ParaView surface helper to get default
* values for each leaf of each dataset of surface
* nComponents could be retrieved with arrayName but is
* given for simplification purposes.
* it is your responsibility to initialize all components of
* defaultValues[nComponent]
*/
void ComputeSurfaceDefaultValues(
const char* arrayName, vtkDataSet* dataset, int nComponents, double* defaultValues) override;
protected:
vtkLagrangianIntegrationModelExample();
~vtkLagrangianIntegrationModelExample() override = default;
/**
* This method is called each time a particle interact with a surface
* With an unrecodgnized surfaceType or SURFACE_TYPE_MODEL
* The particle next position is already positioned exactly on the surface and
* position of the particle is not supposed to be changed
* It is possible in this method to choose to terminate particle, alter it's variables including
* velocity,
* create new particle...
*/
bool InteractWithSurface(int surfaceType, vtkLagrangianParticle* particle, vtkDataSet* surface,
vtkIdType cellId, std::queue<vtkLagrangianParticle*>& particles) override;
/**
* This method is called when trying to find the intersection point between a particle
* and a surface, enabling to use your any intersection code. in this case it only call the
* superclass method
*/
bool IntersectWithLine(vtkLagrangianParticle* particle, vtkCell* cell, double p1[3], double p2[3],
double tol, double& t, double x[3]) override;
double GetRelaxationTime(const double& dynVisc, const double& diameter, const double& density);
double GetDragCoefficient(const double* flowVelocity, const double* particleVelocity,
const double& dynVisc, const double& particleDiameter, const double& flowDensity);
private:
vtkLagrangianIntegrationModelExample(
const vtkLagrangianIntegrationModelExample&) = delete; // Not implemented.
void operator=(const vtkLagrangianIntegrationModelExample&) = delete; // Not implemented.
};
#endif
|
