OpenLB 1.7
Loading...
Searching...
No Matches
Public Member Functions | List of all members
olb::BlockLatticePhysWallShearStress3D< T, DESCRIPTOR > Class Template Referencefinal

functor returns pointwise phys wall shear stress acting on a boundary with a given material on local lattice More...

#include <latticePhysWallShearStress3D.h>

+ Inheritance diagram for olb::BlockLatticePhysWallShearStress3D< T, DESCRIPTOR >:
+ Collaboration diagram for olb::BlockLatticePhysWallShearStress3D< T, DESCRIPTOR >:

Public Member Functions

 BlockLatticePhysWallShearStress3D (BlockLattice< T, DESCRIPTOR > &blockLattice, BlockGeometry< T, 3 > &blockGeometry, int material, const UnitConverter< T, DESCRIPTOR > &converter, IndicatorF3D< T > &indicator)
 
bool operator() (T output[], const int input[]) override
 has to be implemented for 'every' derived class
 
- Public Member Functions inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
BlockLattice< T, DESCRIPTOR > & getBlock ()
 Copy Constructor.
 
- Public Member Functions inherited from olb::BlockF3D< T >
 ~BlockF3D () override
 virtual destructor for defined behaviour
 
virtual BlockStructureD< 3 > & getBlockStructure () const
 
BlockF3D< T > & operator- (BlockF3D< T > &rhs)
 
BlockF3D< T > & operator+ (BlockF3D< T > &rhs)
 
BlockF3D< T > & operator* (BlockF3D< T > &rhs)
 
BlockF3D< T > & operator/ (BlockF3D< T > &rhs)
 
- Public Member Functions inherited from olb::GenericF< T, int >
virtual ~GenericF ()=default
 
int getSourceDim () const
 read only access to member variable _m
 
int getTargetDim () const
 read only access to member variable _n
 
std::string & getName ()
 read and write access to name
 
std::string const & getName () const
 read only access to name
 
bool operator() (T output[])
 wrapper that call the pure virtual operator() (T output[], const S input[]) from above
 
bool operator() (T output[], int input0)
 
bool operator() (T output[], int input0, int input1)
 
bool operator() (T output[], int input0, int input1, int input2)
 
bool operator() (T output[], int input0, int input1, int input2, int input3)
 

Additional Inherited Members

- Public Types inherited from olb::GenericF< T, int >
using targetType
 
using sourceType
 
- Public Attributes inherited from olb::GenericF< T, int >
std::shared_ptr< GenericF< T, int > > _ptrCalcC
 memory management, frees resouces (calcClass)
 
- Protected Member Functions inherited from olb::BlockLatticePhysF3D< T, DESCRIPTOR >
 BlockLatticePhysF3D (BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter, int targetDim)
 
- Protected Member Functions inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
 BlockLatticeF3D (BlockLattice< T, DESCRIPTOR > &blockLattice, int targetDim)
 
- Protected Member Functions inherited from olb::BlockF3D< T >
 BlockF3D (BlockStructureD< 3 > &blockStructure, int targetDim)
 
- Protected Member Functions inherited from olb::GenericF< T, int >
 GenericF (int targetDim, int sourceDim)
 
- Protected Attributes inherited from olb::BlockLatticePhysF3D< T, DESCRIPTOR >
const UnitConverter< T, DESCRIPTOR > & _converter
 
- Protected Attributes inherited from olb::BlockLatticeF3D< T, DESCRIPTOR >
BlockLattice< T, DESCRIPTOR > & _blockLattice
 
- Protected Attributes inherited from olb::BlockF3D< T >
BlockStructureD< 3 > & _blockStructure
 

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::BlockLatticePhysWallShearStress3D< T, DESCRIPTOR >

functor returns pointwise phys wall shear stress acting on a boundary with a given material on local lattice

Definition at line 65 of file latticePhysWallShearStress3D.h.

Constructor & Destructor Documentation

◆ BlockLatticePhysWallShearStress3D()

template<typename T , typename DESCRIPTOR >
olb::BlockLatticePhysWallShearStress3D< T, DESCRIPTOR >::BlockLatticePhysWallShearStress3D ( BlockLattice< T, DESCRIPTOR > & blockLattice,
BlockGeometry< T, 3 > & blockGeometry,
int material,
const UnitConverter< T, DESCRIPTOR > & converter,
IndicatorF3D< T > & indicator )

Definition at line 68 of file latticePhysWallShearStress3D.hh.

74 : BlockLatticePhysF3D<T,DESCRIPTOR>(blockLattice,converter,1),
75 _blockGeometry(blockGeometry),
76 _material(material),
77 _discreteNormal(blockGeometry.getNcells()),
78 _normal(blockGeometry.getNcells())
79{
80 this->getName() = "physWallShearStress";
81 const T scaling = this->_converter.getConversionFactorLength() * 0.1;
82 const T omega = 1. / this->_converter.getLatticeRelaxationTime();
83 const T dt = this->_converter.getConversionFactorTime();
84 const auto& blockGeometryStructure = const_cast<BlockGeometry<T,3>&>(_blockGeometry);
85 _physFactor = -omega * descriptors::invCs2<T,DESCRIPTOR>() / dt * this->_converter.getPhysDensity() * this->_converter.getPhysViscosity();
86 std::vector<int> discreteNormalOutwards(4, 0);
87
88 blockGeometryStructure.forSpatialLocations([&](auto iX, auto iY, auto iZ) {
89 if (blockGeometryStructure.getNeighborhoodRadius({iX,iY,iZ}) >= 1) {
90 if (_blockGeometry.get({iX, iY, iZ}) == _material) {
91 discreteNormalOutwards = blockGeometryStructure.getStatistics().getType(iX, iY, iZ);
92 auto discreteNormal = -1 * Vector<int,3>(&discreteNormalOutwards[1]);
93 _discreteNormal.getRowPointer(blockGeometryStructure.getCellId(iX,iY,iZ)) = discreteNormal;
94
95 T physR[3];
96 _blockGeometry.getPhysR(physR,{iX, iY, iZ});
97 Vector<T,3> origin(physR[0],physR[1],physR[2]);
98 Vector<T,3> direction(-discreteNormal[0] * scaling,
99 -discreteNormal[1] * scaling,
100 -discreteNormal[2] * scaling);
101 Vector<T,3> normal(0.,0.,0.);
102 origin[0] = physR[0];
103 origin[1] = physR[1];
104 origin[2] = physR[2];
105
106 indicator.normal(normal, origin, direction);
107 normal = normalize(normal);
108
109 _normal.getRowPointer(blockGeometryStructure.getCellId(iX,iY,iZ)) = Vector<T,3>(normal);
110 }
111 }
112 });
113}
olb::BlockGeometry::getPhysR
Vector< T, D > getPhysR(LatticeR< D > latticeR)
Definition blockGeometry.h:129
olb::BlockLatticePhysF3D::_converter
const UnitConverter< T, DESCRIPTOR > & _converter
Definition blockBaseF3D.h:170
olb::BlockStructureD::getNcells
std::size_t getNcells() const
Get number of cells.
Definition blockStructure.h:114
olb::ColumnVector::getRowPointer
ptr getRowPointer(std::size_t i)
Definition columnVector.h:114
olb::GenericF< T, int >::getName
std::string & getName()
read and write access to name
Definition genericF.hh:51
olb::IndicatorF3D::normal
virtual bool normal(Vector< S, 3 > &normal, const Vector< S, 3 > &origin, const Vector< S, 3 > &direction, int iC=-1)
returns true and the normal if there was one found for an given origin and direction
Definition indicatorBaseF3D.hh:129
olb::normalize
constexpr Vector< T, D > normalize(const ScalarVector< T, D, IMPL > &a, T scale=T{1})
Definition vector.h:245

References olb::BlockLatticePhysF3D< T, DESCRIPTOR >::_converter, olb::GenericF< T, int >::getName(), olb::BlockGeometry< T, D >::getPhysR(), olb::ColumnVector< COLUMN, D >::getRowPointer(), olb::IndicatorF3D< S >::normal(), and olb::normalize().

+ Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename DESCRIPTOR >
bool olb::BlockLatticePhysWallShearStress3D< T, DESCRIPTOR >::operator() ( T output[],
const int input[] )
overridevirtual

has to be implemented for 'every' derived class

Implements olb::GenericF< T, int >.

Definition at line 116 of file latticePhysWallShearStress3D.hh.

117{
118 output[0] = T();
119 if (_blockGeometry.get({input[0],input[1],input[2]}) == _material) {
120 auto discreteNormal = _discreteNormal.getRowPointer(_blockGeometry.getCellId(LatticeR<3>(input)));
121 auto normal = _normal.getRowPointer(_blockGeometry.getCellId(LatticeR<3>(input)));
122
123 T traction[3];
124 T stress[6];
125 T rho = this->_blockLattice.get(input[0] + discreteNormal[0],
126 input[1] + discreteNormal[1],
127 input[2] + discreteNormal[2]).computeRho();
128 this->_blockLattice.get(input[0] + discreteNormal[0],
129 input[1] + discreteNormal[1],
130 input[2] + discreteNormal[2]).computeStress(stress);
131
132 traction[0] = stress[0]*_physFactor/rho*normal[0] +
133 stress[1]*_physFactor/rho*normal[1] +
134 stress[2]*_physFactor/rho*normal[2];
135 traction[1] = stress[1]*_physFactor/rho*normal[0] +
136 stress[3]*_physFactor/rho*normal[1] +
137 stress[4]*_physFactor/rho*normal[2];
138 traction[2] = stress[2]*_physFactor/rho*normal[0] +
139 stress[4]*_physFactor/rho*normal[1] +
140 stress[5]*_physFactor/rho*normal[2];
141
142 T traction_normal_SP;
143 T tractionNormalComponent[3];
144 // scalar product of traction and normal vector
145 traction_normal_SP = traction[0] * normal[0] +
146 traction[1] * normal[1] +
147 traction[2] * normal[2];
148 tractionNormalComponent[0] = traction_normal_SP * normal[0];
149 tractionNormalComponent[1] = traction_normal_SP * normal[1];
150 tractionNormalComponent[2] = traction_normal_SP * normal[2];
151
152 T WSS[3];
153 WSS[0] = traction[0] - tractionNormalComponent[0];
154 WSS[1] = traction[1] - tractionNormalComponent[1];
155 WSS[2] = traction[2] - tractionNormalComponent[2];
156 // magnitude of the wall shear stress vector
157 output[0] = util::sqrt(WSS[0]*WSS[0] + WSS[1]*WSS[1] + WSS[2]*WSS[2]);
158
159 return true;
160 }
161 else {
162 return true;
163 }
164}
olb::BlockGeometry::get
std::enable_if_t< sizeof...(L)==D, int > get(L... latticeR) const
Read-only access to a material number.
Definition blockGeometry.h:107
olb::BlockLatticeF3D::_blockLattice
BlockLattice< T, DESCRIPTOR > & _blockLattice
Definition blockBaseF3D.h:144
olb::BlockStructureD::getCellId
CellID getCellId(LatticeR< D > latticeR) const
Get 1D cell ID.
Definition blockStructure.h:125
olb::util::sqrt
cpu::simd::Pack< T > sqrt(cpu::simd::Pack< T > value)
Definition pack.h:100
The documentation for this class was generated from the following files:
Generated on Tue Mar 5 2024 09:57:59 for OpenLB by doxygen 1.10.0