24#ifndef LATTICE_PSM_PHYS_FORCE_2D_HH
25#define LATTICE_PSM_PHYS_FORCE_2D_HH
44template<
typename T,
typename DESCRIPTOR>
50 this->
getName() =
"PSMPhysForce";
51 for (
int iC = 0; iC < this->
_sLattice.getLoadBalancer().size(); ++iC) {
59template <
typename T,
typename DESCRIPTOR>
65 this->
getName() =
"physPSMForce";
68template<
typename T,
typename DESCRIPTOR>
71 for (
int i = 0; i < this->getTargetDim(); ++i) {
75 T epsilon = 1. - this->_blockLattice.get(input[0], input[1]).template getField<descriptors::POROSITY>();
78 if ((epsilon > 1e-5)) {
79 T rho, u[DESCRIPTOR::d], u_s[DESCRIPTOR::d];
81 for (
int i = 0; i < DESCRIPTOR::d; i++) {
82 u_s[i] = this->_blockLattice.get(input[0], input[1]).template getFieldComponent<descriptors::VELOCITY_SOLID>(i);
84 T paramA = this->_converter.getLatticeRelaxationTime() - 0.5;
86 T paramB = (epsilon * paramA) / ((1. - epsilon) + paramA);
89 T omega = 1. / this->_converter.getLatticeRelaxationTime();
91 this->_blockLattice.get(input[0], input[1]).computeRhoU(rho, u);
93 const T uSqr_s = util::normSqr<T,DESCRIPTOR::d>(u_s);
94 T uSqr = util::normSqr<T,DESCRIPTOR::d>(u);
95 for (
int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
99 - this->_blockLattice.get(input[0], input[1])[iPop])
101 * (this->_blockLattice.get(input[0], input[1])[iPop]
113 for (
int i = 0; i < this->getTargetDim(); ++i) {
114 output[i] -= descriptors::c<DESCRIPTOR>(iPop,i) * omega_s;
118 for (
int i = 0; i < this->getTargetDim(); ++i) {
119 output[i] = this->_converter.getPhysForce(output[i] * paramB);
128template<
typename T,
typename DESCRIPTOR>
136 this->
getName() =
"PSMPhysForce";
137 for (
int iC = 0; iC < this->
_sLattice.getLoadBalancer().size(); ++iC) {
149template <
typename T,
typename DESCRIPTOR>
155 :
BlockLatticePhysF2D<T,DESCRIPTOR>(blockLattice, converter, 2), _shapeIndicator(shapeIndicator), _blockGeometry(blockGeometry)
157 this->
getName() =
"physPSMForce";
161template<
typename T,
typename DESCRIPTOR>
164 for (
int i = 0; i < this->getTargetDim(); ++i) {
168 T epsilon = 1. - this->_blockLattice.get(input[0], input[1]).template getField<descriptors::POROSITY>();
169 if ((epsilon > 1e-5)) {
171 T posBlockGeo[2] = {0.};
173 _blockGeometry.getPhysR( posBlockGeo, {input[0], input[1]});
174 _shapeIndicator( insideShape, posBlockGeo );
176 if ( insideShape[0] ==
true ) {
178 T rho, u[DESCRIPTOR::d], u_s[DESCRIPTOR::d];
180 for (
int i = 0; i < DESCRIPTOR::d; i++) {
181 u_s[i] = this->_blockLattice.get(input[0], input[1]).template getFieldComponent<descriptors::VELOCITY_SOLID>(i);
183 T paramA = this->_converter.getLatticeRelaxationTime() - 0.5;
185 T paramB = (epsilon * paramA) / ((1. - epsilon) + paramA);
188 T omega = 1. / this->_converter.getLatticeRelaxationTime();
190 this->_blockLattice.get(input[0], input[1]).computeRhoU(rho, u);
192 const T uSqr_s = util::normSqr<T,DESCRIPTOR::d>(u_s);
193 T uSqr = util::normSqr<T,DESCRIPTOR::d>(u);
194 for (
int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
198 - this->_blockLattice.get(input[0], input[1])[iPop])
200 * (this->_blockLattice.get(input[0], input[1])[iPop]
212 for (
int i = 0; i < this->getTargetDim(); ++i) {
213 output[i] -= descriptors::c<DESCRIPTOR>(iPop,i) * omega_s;
217 for (
int i = 0; i < this->getTargetDim(); ++i) {
218 output[i] = this->_converter.getPhysForce(output[i] * paramB);
Representation of a block geometry.
functor returns pointwise phys force for PSM dynamics
BlockLatticePSMPhysForce2DMod(BlockLattice< T, DESCRIPTOR > &blockLattice, BlockGeometry< T, 2 > &blockGeometry, const UnitConverter< T, DESCRIPTOR > &converter, IndicatorF2D< T > &shapeIndicator)
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
functor returns pointwise phys force for PSM dynamics
BlockLatticePSMPhysForce2D(BlockLattice< T, DESCRIPTOR > &blockLattice, const UnitConverter< T, DESCRIPTOR > &converter)
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Platform-abstracted block lattice for external access and inter-block interaction.
std::string & getName()
read and write access to name
IndicatorF2D is an application from .
std::vector< std::unique_ptr< BlockF2D< T > > > _blockF
Super functors may consist of several BlockF2D<W> derived functors.
Representation of a statistic for a parallel 2D geometry.
BlockGeometry< T, D > & getBlockGeometry(int locIC)
Read and write access to a single block geometry.
SuperLattice< T, DESCRIPTOR > & _sLattice
SuperLatticePSMPhysForce2DMod(SuperLattice< T, DESCRIPTOR > &sLattice, SuperGeometry< T, 2 > &superGeometry, const UnitConverter< T, DESCRIPTOR > &converter, IndicatorF2D< T > &shapeIndicator)
SuperLatticePSMPhysForce2D(SuperLattice< T, DESCRIPTOR > &sLattice, const UnitConverter< T, DESCRIPTOR > &converter)
represents all functors that operate on a DESCRIPTOR with output in Phys, e.g. physVelocity(),...
Super class maintaining block lattices for a cuboid decomposition.
Conversion between physical and lattice units, as well as discretization.
The description of a generic interface for all functor classes – header file.
This file contains indicator functions.
Wrapper functions that simplify the use of MPI.
Top level namespace for all of OpenLB.
Collection of common computations for LBM.
Representation of a parallel 2D geometry – header file.