<< scibench_getsettings Benchmark Toolbox scibench_matmul >>

Benchmark Toolbox >> Benchmark Toolbox > scibench_latticeboltz

scibench_latticeboltz

Performs a 2D Lattice Boltzmann simulation with Scilab.

Calling Sequence

scibench_latticeboltz ( lx , ly , maxT , tPlot )

Parameters

lx :

a 1-by-1 matrix of floating point integers, the number of cells in the X direction, must be greater than 1

ly :

a 1-by-1 matrix of floating point integers, the number of cells in the Y direction, must be greater than 1

maxT :

a 1-by-1 matrix of doubles, the maximum number of iterations, must be greater than 1

tPlot :

a 1-by-1 matrix of doubles, the number of iterations between 2 refresh of the graphics, must be greater than 1

Description

This function simulates the channel flow past a cylinderical obstacle, using a Lattice Boltzmann method.

Get the most recent version of this file on LBMethod.org: http://www.lbmethod.org/_media/numerics:cylinder.m

Original implementaion of Zou/He boundary condition by Adriano Sciacovelli (see example "cavity.m"). The script was translated to scilab language by Vincent Lejeune. The performance was improved by Michael Baudin, using vectorized statements.

Examples

lx     = 400;      // number of cells in x-direction
ly     = 100;      // number of cells in y-direction
tPlot  = 50;       // cycles
maxT   = 1000;     // total number of iterations
scf();
tic();
scibench_latticeboltz ( lx , ly , maxT , tPlot );
t = toc();
mprintf("lx=%d, ly=%d, maxT=%d, t=%.3f\n",lx,ly,maxT,t);

Bibliography

http://www.lbmethod.org/_media/numerics:cylinder.m

Authors

<< scibench_getsettings Benchmark Toolbox scibench_matmul >>