palabos-v2.0r0.zip

/* This file is part of the Palabos library.
 *
 * Copyright （C） 2011-2017 FlowKit Sarl
 * Route d‘Oron 2
 * 1010 Lausanne Switzerland
 * E-mail contact: contact@flowkit.com
 *
 * The most recent release of Palabos can be downloaded at 
 * 
 *
 * The library Palabos is free software: you can redistribute it and/or
 * modify it under the terms of the GNU Affero General Public License as
 * published by the Free Software Foundation either version 3 of the
 * License or （at your option） any later version.
 *
 * The library is distributed in the hope that it will be useful
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU Affero General Public License for more details.
 *
 * You should have received a copy of the GNU Affero General Public License
 * along with this program.  If not see .
*/

/** file
  * Flow in a lid-driven 3D cavity. Benchmark case
**/

#include “palabos3D.h“
#include “palabos3D.hh“   // include full template code
#include 

using namespace plb;
using namespace std;

typedef double T;
#define DEscriptOR descriptors::D3Q19Descriptor

void cavitySetup（ MultiBlockLattice3DriptOR>& lattice
                  IncomprFlowParam const& parameters
                  OnLatticeBoundaryCondition3DriptOR>& boundaryCondition ）
{
    const plint nx = parameters.getNx（）;
    const plint ny = parameters.getNy（）;
    const plint nz = parameters.getNz（）;
    Box3D topLid = Box3D（0 nx-1 ny-1 ny-1 0 nz-1）;
    Box3D everythingButTopLid = Box3D（0 nx-1 0 ny-2 0 nz-1）;

    // All walls implement a Dirichlet velocity condition.
    boundaryCondition.setVelocityConditionOnBlockBoundaries（lattice）;

    T u = std::sqrt（（T）2）/（T）2 * parameters.getLatticeU（）;
    initializeAtEquilibrium（lattice everythingButTopLid （T） 1. Array（（T）0.（T）0.（T）0.） ）;
    initializeAtEquilibrium（lattice topLid （T） 1. Array（u（T）0.u） ）;
    setBoundaryVelocity（lattice topLid Array（u0.u） ）;

    lattice.initialize（）;
}

int main（int argc char* argv[]） {

    plbInit（&argc &argv）;
    //defaultMultiBlockPolicy3D（）.toggleBlockingCommunication（true）;

    plint N;
    try {
        global::argv（1）.read（N）;
    }
    catch（...）
    {
        pcout << “Wrong parameters. The syntax is “ << std::endl;
        pcout << argv[0] << “ N“ << std::endl;
        pcout << “where N is the resolution. The benchmark cases published “ << std::endl;
        pcout << “on the Palabos Wiki use N=100 N=400 N=1000 or N=4000.“ << std::endl;
        exit（1）;
    }

    pcout << “Starting benchmark with “ << N+1 << “x“ << N+1 << “x“ << N+1 << “ grid points “
          << “（approx. 2 minutes on modern processors）.“ << std::endl;


    IncomprFlowParam parameters（
            （T） 1e-2  // uMax
            （T） 1.    // Re
            N         // N
            1.        // lx
            1.

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-09-15 12:01  palabos-v2.0r0
     文件         107  2017-09-15 12:01  palabos-v2.0r0.gitignore
     文件        3605  2017-09-15 12:01  palabos-v2.0r0.gitlab-ci.yml
     文件       34520  2017-09-15 12:01  palabos-v2.0r0COPYING
     文件         339  2017-09-15 12:01  palabos-v2.0r0README.md
     文件        3118  2017-09-15 12:01  palabos-v2.0r0SConstruct
     目录           0  2017-09-15 12:01  palabos-v2.0r0codeblocks
     文件        9799  2017-09-15 12:01  palabos-v2.0r0codeblocksPalabos.cbp
     文件         397  2017-09-15 12:01  palabos-v2.0r0codeblocksPalabos.layout
     目录           0  2017-09-15 12:01  palabos-v2.0r0codeblocks mp
     文件          71  2017-09-15 12:01  palabos-v2.0r0codeblocks mp.gitignore
     目录           0  2017-09-15 12:01  palabos-v2.0r0examples
     目录           0  2017-09-15 12:01  palabos-v2.0r0examplesenchmarks
     目录           0  2017-09-15 12:01  palabos-v2.0r0examplesenchmarkscavity3d
     文件        3502  2017-09-15 12:01  palabos-v2.0r0examplesenchmarkscavity3dMakefile
     文件        4579  2017-09-15 12:01  palabos-v2.0r0examplesenchmarkscavity3dcavity3d.cpp
     目录           0  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopic
     目录           0  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamics
     文件        3508  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamicsMakefile
     文件       26305  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamicsaneurysm.cpp
     文件      625121  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamicsaneurysm.stl.tgz
     文件        7643  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamicscavity3d.cpp
     文件        2477  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamicsparam_aneurysm.xml
     文件       11269  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamicspoiseuille.cpp
     目录           0  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamics mp
     文件          71  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicasinariDynamics mp.gitignore
     目录           0  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicounceBack
     文件        3504  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicounceBackMakefile
     文件        6135  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicounceBackcomputeDrag.cpp
     文件        3354  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicounceBackinstantiateCylinder.cpp
     目录           0  2017-09-15 12:01  palabos-v2.0r0examplescodesByTopicounceBack mp
............此处省略2293个文件信息