Function lattice_boltzmann_cfd_demo

examples/pde/boltzmann_cfd.cpp:449–552

}

void lattice_boltzmann_cfd_demo() {
    // Define the lattice for the simulation
    const size_t len         = 128;
    const size_t grid_width  = len;
    const size_t grid_height = len;

    // Specify the image scaling displayed
    float scale = 4.0f;

    // Forge window initialization
    int height = static_cast<int>(grid_width * scale);
    int width  = static_cast<int>(grid_height * scale);
    af::Window window(height, width, "Driven Cavity Flow");

    int frame_count       = 0;
    int max_frames        = 20000;
    int simulation_frames = 100;
    float total_time      = 0;
    float total_time2     = 0;

    // CFD fluid parameters
    const float density  = 2.7f;
    const float velocity = 0.35f;
    const float reynolds = 1e5f;

    const char* ux_image = ASSETS_DIR "/examples/images/default_ux.bmp";
    const char* uy_image = ASSETS_DIR "/examples/images/default_uy.bmp";
    const char* set_boundary_image =
        ASSETS_DIR "/examples/images/default_boundary.bmp";

    // Tesla Valve Fluid Simulation - entering from constricted side
    {
        //           ux_image = ASSETS_DIR "/examples/images/left_tesla_ux.bmp";
        //           uy_image = ASSETS_DIR "/examples/images/left_tesla_uy.bmp";
        // set_boundary_image = ASSETS_DIR
        // "/examples/images/left_tesla_boundary.bmp";
    }

    // Tesla Valve Fluid Simulation - entering from transfer side
    {
        //           ux_image = ASSETS_DIR
        //           "/examples/images/right_tesla_ux.bmp"; uy_image =
        //           ASSETS_DIR "/examples/images/right_tesla_uy.bmp";
        // set_boundary_image = ASSETS_DIR
        // "/examples/images/right_tesla_boundary.bmp";
    }

    // Reads the initial values of fluid quantites and simulation parameters
    Simulation sim =
        create_simulation(grid_width, grid_height, density, velocity, reynolds,
                          ux_image, uy_image, set_boundary_image);

    // Initializes the simulation quantites
    initialize(sim);

    while (!window.close() && frame_count != max_frames) {
        af::sync();
        auto begin = std::chrono::high_resolution_clock::now();