Mike
Forum Replies Created
-
Posts
-
Here is my code, I do use discretnormal directions to get unknownIndexes.
As for setter code, it is like:
setHalfBounceBackBoundary(NSlattice, omega, supergeometry, 4);
4 means solid node, which is implemented with nodynamics.template<typename T, typename DESCRIPTOR, class MixinDynamics>
void setHalfBounceBackBoundary(BlockLatticeStructure2D<T,DESCRIPTOR>& block, T omega, BlockIndicatorF2D<T>& indicator, bool includeOuterCells)
{
OstreamManager clout(std::cout, “setHalfBounceBackBoundary2D”);
auto& blockGeometryStructure = indicator.getBlockGeometryStructure();
const int margin = includeOuterCells ? 0 : 1;
/*
*x0,x1,y0,y1 Range of cells to be traversed
**/
int x0 = margin;
int y0 = margin;
int x1 = blockGeometryStructure.getNx()-1 -margin;
int y1 = blockGeometryStructure.getNy()-1 -margin;
std::vector<int> discreteNormal(3, 0);
for (int iX = x0; iX <= x1; ++iX) {
for (int iY = y0; iY <= y1; ++iY) {
//momenta vector provisionally inside src/core/blockLatticeStructure3D.h
Momenta<T, DESCRIPTOR>* momenta = nullptr;
Dynamics<T, DESCRIPTOR>* dynamics = nullptr;
// Momenta<T, DESCRIPTOR>* momenta = new BulkMomenta<T, DESCRIPTOR>();
// Dynamics<T, DESCRIPTOR>* dynamics = new NoDynamics<T, DESCRIPTOR>();
PostProcessorGenerator2D<T, DESCRIPTOR>* postProcessor = nullptr;
if (indicator(iX, iY)) {
discreteNormal = indicator.getBlockGeometryStructure().getStatistics().getType(iX, iY);
clout << discreteNormal[0] << ” ” << discreteNormal[1] << ” ” << discreteNormal[2] << std::endl;if (discreteNormal[0] == 0) {
//set the momenta,dynamics and post processor on the indicated local velocity boundary cells
if (discreteNormal[1] == 1) {
postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 0, 1>(iX, iY);
}
else if (discreteNormal[1] == -1) {postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 0, -1>(iX, iY);
}
else if (discreteNormal[2] == 1) {postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 1, 1>(iX, iY);
}
else if (discreteNormal[2] == -1) {postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 1, -1>(iX, iY);
}
else {
clout << “here “<< “Could not setHalfBounceBackBoundary2D (” << iX
<< “, ” << iY << “)” << std::endl;
clout << discreteNormal[0] << ” ” << discreteNormal[1] << ” ” << discreteNormal[2] << std::endl;
}
}
else if (discreteNormal[0] == 1) {
//sets the momenta, dynamics and post processors on indicated localVelocityCornerBoundary cells
if (discreteNormal[1] == 1) {
if (discreteNormal[2] == 1) {
postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 1, 8>(iX, iY);
}
else if (discreteNormal[2] == -1) {
// clout << “5” << endl;postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 1, -8>(iX, iY);
}
else {
clout << “Could not setHalfBounceBackBoundary2D (” << iX
<< “, ” << iY << “)” << std::endl;
}
}
else if (discreteNormal[1] == -1) {
if (discreteNormal[2] == 1) {
postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 1, 9>(iX, iY);
}
else if (discreteNormal[2] == -1) {
postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 1, -9>(iX, iY);
}
else {
clout << “Could not setHalfBounceBackBoundary2D (” << iX
<< “, ” << iY << “)” << std::endl;
}
}
}
//sets momenta, dynamics and postProcessors on local InnerVelocityCornerBoundary Cells
else if (discreteNormal[0] == 2) {
if (discreteNormal[1] == 1) {
if (discreteNormal[2] == 1) {postProcessor = nullptr;
}
else if (discreteNormal[2] == -1) {postProcessor = nullptr;
}
else {
clout << “Could not setHalfBounceBackBoundary2D (” << iX
<< “, ” << iY << “)” << std::endl;
}
}
else if (discreteNormal[1] == -1) {
if (discreteNormal[2] == 1) {postProcessor = nullptr;
}
else if (discreteNormal[2] == -1) {
clout << “6” << endl;// postProcessor = new BoundaryStreamPostProcessorGenerator2D<T, DESCRIPTOR, 1, 8>(iX, iY);
}
else {
clout << “Could not setHalfBounceBackBoundary2D (” << iX
<< “, ” << iY << “)” << std::endl;
}
}
}setBoundary<T, DESCRIPTOR, MixinDynamics>(block, omega, iX,iY, momenta, dynamics, postProcessor);
}
}
}}
Sorry, I answered too fast and has used up my modified limit…
your idea seems different from what I thought before, yours in code is:
cell.neighbor(descriptors::c<DESCRIPTOR>(unknownIndexes[i]))[util::opposite<DESCRIPTOR>(unknownIndexes[i])] = cell[unknownIndexes[i]];am I right?
I will change code and see if the result is fine. Thank you and sorry to bother again…Thank you for your reply!
I agree with your idea.
Below is my code :
std::vector<int> unknownIndexes;
// east-top
// I change the oritation in self-defined sethalfbbboundary.hh
if(orientation == 8){unknownIndexes = {3};
}
//east bottom
else if(orientation == -8){unknownIndexes = {1};
}
// west top
else if(orientation == 9){unknownIndexes = {5};
}
// west bottom
else if(orientation == -9){
unknownIndexes = {7};}
else{
unknownIndexes = util::subIndexOutgoing<DESCRIPTOR, direction,
orientation>();
}
auto cell = blockLattice.get(x, y);for (unsigned i = 0; i < unknownIndexes.size(); ++i) {
cell[unknownIndexes[i]] = cell.neighbor(descriptors::c<DESCRIPTOR>(unknownIndexes[i]))[util::opposite<DESCRIPTOR>(unknownIndexes[i])];
}
as in main function, I define solid boundaries nodynamics.
Hello Adian,
Does the class BounceBack in olb-1.4 means full-way bounce back? I am running transient cases, and I learned that half-way bounce back is more accurate than full-way in time. so I wonder how to implement it in openlb.
since you have mentioned postprocessor previously, then I tried with that but it didn’t work well. Could you elaborate how to implement it? Thank you in advance.
Here is my plan, since half-way only change streaming process, I define walls with bulk dynamics and use a postprocessor to make modifications, which basically implement bounce-back rules according to their directions.
Yours,
MikeHello everyone,
I also wonder why pcm specific heat is set as 1, which I think it should be set as converter.getLatticeSpecificHeatCapacity(physSpecificHeatCapacity), it will be much appreciated if anyone can help solve this!
Best
Mike