LBM Midrange Repulsion

  • x4_337103
    了解作者
  • 4.3KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • VIP专享
    资源类型
  • 0
    下载次数
  • 2022-05-14 02:08
    上传日期
基于 MATLAB 的 Midrange repulsion 伪势多相流模型。
Two phase flow - midrange repulsion.zip
  • Two phase flow - midrange repulsion
  • forces.m
    1.3KB
  • visua.m
    1.2KB
  • constant.m
    2KB
  • tflow.m
    992B
  • collision.m
    841B
  • initialization.m
    888B
内容介绍
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % constant.m: consant setting % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Shan Chen Lattice Boltzmann sample in Matlab % Copyright Wei Gong % Address: Nottingham NG7 2RD, UK % E-mail: wei.gong@nottingham.ac.uk % D2Q9 lattice constants lx = single(512); % x length of computational domain ly = single(512); % y length of computational domain lxy = lx*ly; tau = single(1); % dimensionless relaxation time cc = single(1); % speeds c_squ = cc*cc/3; % square of sound speed visc = c_squ*(tau - 0.5); % viscosity t_k = single([4/9 1/9 1/9 1/9 1/9 1/36 1/36 1/36 1/36]); % weighting factors ps = single([0 4/63 4/63 4/63 4/63 4/135 4/135 4/135 4/135]); pm = single([1/180 2/945 1/15120 2/945 1/180 2/945 1/15120 2/945 ... 1/180 2/945 1/15120 2/945 1/180 2/945 1/15120 2/945]); ee = single([0 1 0 -1 0 1 -1 -1 1; % velocity model 0 0 1 0 -1 1 1 -1 -1]); em = single([2 2 2 1 0 -1 -2 -2 -2 -2 -2 -1 0 1 2 2; 0 1 2 2 2 2 2 1 0 -1 -2 -2 -2 -2 -2 -1]); % General flow constants rho = 0.3 + rand(lx, ly); p = single(zeros(lx, ly)); % pressure psx = single(zeros(lx, ly)); % mean field potential u = single(zeros(2, lx, ly)); % velocity ue = single(zeros(2, lx, ly)); % equilibrium velocity up = single(zeros(2, lx, ly)); % real fluid velocity F = single(zeros(2, lx, ly)); % for the interaction between fluid nodes S = single(zeros(2, lx, ly)); % interaction components between the dluid nodes and solid nodes ff = single(zeros(9, lx, ly)); % distribution function psx_re = single(zeros(9, lx, ly)); % psx(x+ei*dt,t) % parameters in C-S EOS % G1 = single(-4.9); % original Shan-Chen % G2 = single(0); % G1 = single(-20); % G2 = single(15.1); G1 = single(-15); G2 = single(10.1); rho0 = single(1);
评论
    相关推荐
    • matlab 课件
      matlab 功能很强大,本课件详细的介绍了matlab在科学计算中的应用。
    • Matlab教程
      MatlaB经典ppt,主要讲解了MatlaB的入门知识,图形编程以及仿真。
    • matlab入门
      matlab的入门教程,很值得仔细看看,对提高matlab能力有帮助
    • matlab编译器
      matlab编译为exe文件所需要的文件; 使用时解压,将解压后文件拖入至matlab命令行即可
    • matlab
      matlab的基础知识进行介绍
    • MatLab
      “#MatLab
    • MATLAB编译器
      基于MATLAB 2018b版本介绍MATLAB编译器。介绍如何利用编译器将MATLAB代码编译为独立应用程序或组件,并在没有安装MATLAB的计算机上进行部署。
    • MATLAB 论文
      基于MATLAB的论文撰写!可以为你指明一条写作MATLAB论文的道路!
    • matlab
      matlab经典课堂讲义
    • matlab
      matlab