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
123
124
125
126
127
128
129
130
131
132
133
134
|
/*=========================================================================
Program: ParaView
Module: vtkLagrangianIntegrationModelExample.h
Copyright (c) Kitware, Inc.
All rights reserved.
See Copyright.txt or http://www.paraview.org/HTML/Copyright.html for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
/**
* @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
|
