56 auto comprCase(-(generateDFCylinderInf(rCase, makeAVec(0., 0., 1.), center)));
58 return normalize(comprCase, bl.
dx);
62int main(
int argc,
char* argv[])
69 appParamsManager.
load(argc, argv);
75 Param w(6.*3.14*2./60.);
77 AVec<int> size(makeAVec(33., 33., 90.)*(1./dx.v()));
78 asl::Block bl(size, dx.v(), makeAVec(-16.65, -26.15, -67.7));
85 std::cout <<
"Data initialization..." << flush;
90 asl::Block block(compressorMap->getInternalBlock());
96 std::cout <<
"Finished" << endl;
98 std::cout <<
"Numerics initialization..." << flush;
111 std::vector<asl::SPNumMethod> bc;
112 std::vector<asl::SPNumMethod> bcVis;
118 bc.push_back(generateBCVelocity(lbgk, vfCase, comprCaseMap));
120 bcVis.push_back(generateBCVelocityVel(lbgk, vfCase, comprCaseMap));
122 bc.push_back(generateBCNoSlip(lbgk, compressorMap));
123 bcVis.push_back(generateBCNoSlipVel(lbgk, compressorMap));
129 std::cout <<
"Finished" << endl;
130 std::cout <<
"Computing..." << endl;
134 writer.
addScalars(
"compressor", *compressorMap);
137 writer.
addVector(
"v", *lbgk->getVelocity());
145 for (
unsigned int i(1); i < 10001; ++i)
158 cout <<
"Finished" << endl;
160 cout <<
"Computation statistic:" << endl;
161 cout <<
"Real Time = " << timer.
realTime() <<
"; Processor Time = "
void load(int argc, char *argv[])
const V getBPosition() const
Numerical method for fluid flow.
contains different kernels for preprocessing and posprocessing of data used by LBGK
const double realTime() const
const double processorTime() const
const double processorLoad() const
Updatable value. This class stores value and its TimeStamp.
void addVector(std::string name, AbstractData &data)
void addScalars(std::string name, AbstractData &data)
asl::SPDistanceFunction generateCase(asl::Block &bl)
asl::UValue< double > Param
SPDataWrapperACLData generateDataContainerACL_SP(const Block &b, unsigned int n=1)
generates pointer to ACL Data field with n components
SPDataWithGhostNodesACLData readSurface(const string &fileName, double dx, acl::CommandQueue queue=acl::hardware.defaultQueue)
std::shared_ptr< DistanceFunction > SPDistanceFunction
SPPositionFunction generatePFRotationField(const AVec< double > &axis, const AVec< double > &c)
const VectorTemplate & d3q15()
Vector template.
SPBCond generateBCConstantPressure(SPLBGK nm, double p, const std::vector< SlicesNames > &sl)
VectorOfElements generateVEConstant(T a)
Generates VectorOfElements with 1 Element acl::Constant with value a.
std::shared_ptr< LBGKUtilities > SPLBGKUtilities
std::shared_ptr< LBGK > SPLBGK
void initData(SPAbstractData d, double a)