Commit 3a4d161a authored by Martin Bauer's avatar Martin Bauer

Log cleanup & more config options for complex geometry app

parent 45c1613c
Pipeline #15334 failed with stages
in 142 minutes and 13 seconds
......@@ -114,7 +114,9 @@ void initializeVelocity( const shared_ptr<StructuredBlockForest> & blocks, Block
{
auto pdfField = block.getData< lbm::PdfField< LatticeModel_T>>( pdfFieldId );
WALBERLA_FOR_ALL_CELLS_XYZ(pdfField,
Vector3< real_t > coord( x + 0.5, y + 0.5, z + 0.5 );
Vector3< real_t > coord( x + 0.5 , y + 0.5, z + 0.5 );
//Vector3< real_t > coord( x + 0.5 *blocks->dx() , y + 0.5 *blocks->dx(), z + 0.5 *blocks->dx());
blocks->transformBlockLocalToGlobal( coord, block );
auto zVel = u_s * std::pow( (coord[2] / zMax), 1.0 / 7.0);
pdfField->setDensityAndVelocity(x, y, z, Vector3<real_t>(zVel, 0, 0), 1.0);
......@@ -279,10 +281,8 @@ int main( int argc, char * argv[] )
mesh::readAndBroadcast( fileName, *mesh);
auto singleAABB = computeAABB( *mesh );
aabb.merge(singleAABB);
WALBERLA_LOG_DEVEL_ON_ROOT( "Adding distance info to mesh" );
WALBERLA_LOG_INFO_ON_ROOT("Adding meshfile " << fileName);
auto triDist = make_shared< mesh::TriangleDistance<mesh::TriangleMesh> >( mesh );
WALBERLA_LOG_DEVEL_ON_ROOT( "Building distance octree" );
auto distanceOctreeSingle = make_shared< mesh::DistanceOctree<mesh::TriangleMesh> >( triDist );
distanceOctree->add(distanceOctreeSingle);
}
......
......@@ -13,7 +13,7 @@ ComplexGeometry
coarseDx 4;
coarseOmega 1.9999454;
coarseTimeSteps 201;
numLevels 3;
numLevels 2;
bodyForce <0.0, 0, 0>;
blockSize <32,32,32>;
domainBlowUp <5,2,6>; // simulation domain is blow up factor times mesh size per dimension
......@@ -31,7 +31,7 @@ VTK {
fluid // identifier for this VTKOutput object
{
writeFrequency 10;
writeFrequency 100;
inclusion_filters
{
......@@ -39,12 +39,30 @@ VTK {
}
writers
{
Velocity;
Density;
}
}
fluid_coarse
{
samplingResolution 4;
writeFrequency 50;
inclusion_filters
{
AABB_filter
{
min <-100, -50, 0>;
max < 50 , 50, 40>;
}
}
writers
{
Velocity;
}
}
flags // identifier for this VTKOutput object
{
......@@ -65,11 +83,11 @@ VTK {
StabilityChecker
{
checkFrequency 1;
checkFrequency 10;
streamOutput false;
vtkOutput true;
vtkBaseFolder vtk_out; // stability_check;
vtkBinary true;
vtkLittleEndian true;
vtkMPIIO true;
}
\ No newline at end of file
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment