NS方程MATLAB

大气湍流的速度场模拟，MATLAB仿真，GUI，可适用于各种流场

function additionalexamples（num）

% 12 additional example meshes for mesh2d including some larger meshes 
% used in CFD studies.
%
% Example:
%
%   additionalexamples（n）          % Will run through the nth example


switch（num）
    
    case 1
        
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        node  = [00; 10; 11; 01];
        
        hdata.type = 2;
        hdata.h0   = 0.01;

        [pt] = mesh2d（node[]hdata）;
        
    case 2
        
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        theta = 0:pi/50:（2*pi-pi/50）;
        x     = cos（theta）/2;
        y     = sin（theta）/2;
        
        node = [ x‘  y‘
            -5  -5
            5  -5
            5   15
            -5   15 
        ];
        
        n = size（node1）-4;
        
        cnect = [（1:n-1）‘ （2:n）‘
            n        1
            n+1      n+2
            n+2      n+3
            n+3      n+4
            n+4      n+1
        ];    
        
        hdata.type = 2;
        hdata.h0   = 0.175;
        
        [pt] = mesh2d（nodecnecthdata）;
        
    case 3
        
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        theta = 0:pi/50:（2*pi-pi/50）;
        x     = cos（theta）/2;
        y     = sin（theta）/2;
        
        node = [ x‘  y‘
            -5  -10
            25  -10
            25   10
            -5   10 
        ];
        
        n = size（node1）-4;
        
        cnect = [（1:n-1）‘ （2:n）‘
            n        1
            n+1      n+2
            n+2      n+3
            n+3      n+4
            n+4      n+1
        ];    
        
        h0 = 0.1;
        
        hdata.type  = 1;
        hdata.slope = 0.025;
        hdata.h0    = h0;
        hdata.ref   = [00250];
        hdata.lim   = 5;
        
        options.dhmax = 0.2;
        
        [pt] = mesh2d（nodecnecthdataoptions）;
        
    case 4        
        
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        theta = 0:pi/36:（2*pi-pi/36）;
        x     = cos（theta）/2;
        y     = sin（theta）/2;
        
        cyL1 = [x‘ y‘+1];
        cyL2 = [x‘ y‘-1];
        box  = [-5  -10
            25  -10
            25   10
            -5   10 
        ];
        
        n1 = size（cyL11）;
        n2 = size（cyL21）;
        
        c1 = [（1:n1-1）‘（2:n1）‘; n11];
        c2 = [（1:n2-1）‘（2:n2）‘; n21];
        c3 = [12; 23; 34; 41];
        
        node  = [cyL1; cyL2; box];
        cnect = [c1; c2+n1; c3+n1+n2];  
        
        h0 = 0.0175;
        
        hdata（1）.type  = 1;
        hdata（1）.slope = 0;
        hdata（1）.h0    = 10*h0;
        hdata（1）.ref   = [-10250];
        hdata（1）.lim   = 5;
        
        [pt] = mesh2d（nodecnecthdata）;

    case 5
     

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1046  2007-06-27 19:23  Navier2d_matlab\Generate unstructured 2D triangular meshes.txt

     文件     930940  2007-06-27 21:08  Navier2d_matlab\Navier2d\Additional Examples.pdf

     文件      18354  2006-04-22 16:05  Navier2d_matlab\Navier2d\GNUlicense.txt

     文件     573658  2007-06-27 21:57  Navier2d_matlab\Navier2d\readme.pdf

     文件      17777  2006-05-03 19:35  Navier2d_matlab\Navier2d\Windows\build_operators.m

     文件       3682  2006-04-23 19:09  Navier2d_matlab\Navier2d\Windows\cavity_benchmark.m

     文件       1093  2006-04-22 15:44  Navier2d_matlab\Navier2d\Windows\copyright.m

     文件     111232  2006-03-28 10:15  Navier2d_matlab\Navier2d\Windows\cylinder.mat

     文件       5518  2006-04-22 15:58  Navier2d_matlab\Navier2d\Windows\data_structure.m

     文件       2592  2006-03-29 07:39  Navier2d_matlab\Navier2d\Windows\LUsubs.c

     文件       6656  2006-03-29 07:36  Navier2d_matlab\Navier2d\Windows\LUsubs.dll

     文件      48410  2007-09-03 21:24  Navier2d_matlab\Navier2d\Windows\Navier2d.asv

     文件      48410  2007-06-28 20:52  Navier2d_matlab\Navier2d\Windows\Navier2d.m

     文件       1819  2006-03-29 08:53  Navier2d_matlab\Navier2d\Windows\smvp.c

     文件       6144  2006-03-29 08:53  Navier2d_matlab\Navier2d\Windows\smvp.dll

     文件      58624  2006-03-28 10:17  Navier2d_matlab\Navier2d\Windows\square.mat

     文件      42330  2006-05-09 15:02  Navier2d_matlab\Navier2d\Windows\tvd_rk2.m

     文件      17777  2006-05-03 19:35  Navier2d_matlab\Navier2d\Not windows\build_operators.m

     文件       3682  2006-03-28 09:46  Navier2d_matlab\Navier2d\Not windows\cavity_benchmark.m

     文件       1093  2006-04-22 15:44  Navier2d_matlab\Navier2d\Not windows\copyright.m

     文件     111232  2006-03-28 10:15  Navier2d_matlab\Navier2d\Not windows\cylinder.mat

     文件       5518  2006-04-22 15:58  Navier2d_matlab\Navier2d\Not windows\data_structure.m

     文件      48277  2006-05-09 15:28  Navier2d_matlab\Navier2d\Not windows\Navier2d.m

     文件      58624  2006-03-28 10:17  Navier2d_matlab\Navier2d\Not windows\square.mat

     文件      42021  2006-05-09 15:34  Navier2d_matlab\Navier2d\Not windows\tvd_rk2.m

     文件      20347  2007-05-09 21:36  Navier2d_matlab\meshnewer\additionalexamples.m

     文件        935  2007-05-09 20:56  Navier2d_matlab\meshnewer\fixmesh.m

     文件       7181  2007-05-09 20:52  Navier2d_matlab\meshnewer\inpoly.m

     文件      19010  2007-05-09 21:43  Navier2d_matlab\meshnewer\mesh2d.m

     文件      18034  2007-04-08 17:44  Navier2d_matlab\meshnewer\meshdemo.m

............此处省略15个文件信息