1 内容介绍

基于溶质扩散和界面能的作用,考虑成分过冷,曲率过冷,界面能各 向异性和界面扰动等因素,建立了单个等轴枝晶的生长模型.采用元胞自动机(cellular automata)方法模拟了枝晶生长,界面扰动和分枝的竞争演化.对液相中的溶质传输进行了计算.通过模拟发现,凝同过程中溶质容易富集在枝晶臂之间, 同时,随着凝固时间的延长,界面前沿的溶质浓度梯度逐渐下降.

2 部分代码

%

clc;clear;close all;

rand('state',0);

set(gcf,'DoubleBuffer','on');

N=256; % 生成状态矩阵大小的控制参数

S=ones(N); % 生成状态矩阵S

% S(N/2,N/4:N*3/4)=0; % 设置状态矩阵的初始值

S(N/2:N*3/4,N/2)=0; % 设置状态矩阵的初始值

Ii=imshow(S); % 显示状态矩阵

ti=title(['time = 0',', N=0'],'Fontsize',14,'Fontname','Times New Roman'); % 显示时间

T=0; % 记录时间的参数

for k=1:30000; % 循环计算

xt=4; % 产生粒子的位置

yt=N/2; % 产生粒子的位置

Ss=0; % 控制下面循环是否终止的参数

while Ss<1; % 计算当前粒子的吸附过程

[dx,dy]=choose_D(0.25,0.25,0.25,0.25); % 计算粒子位移量

xt=xt+dx; % 计算下一时刻粒子的位置

yt=yt+dy; % 计算下一时刻粒子的位置

if xt<2.5|xt>N-1.5||yt<2.5|yt>N-1.5; % 判断粒子是否逸出边界

Ss=2; % 把Ss赋值为2,从而可以退出循环

end

if [S(xt+1,yt)+S(xt-1,yt)+S(xt,yt+1)+S(xt,yt-1)]<3.5; % 判断近邻是否有粒子

S(xt,yt)=0; % 固定当前运动的粒子

Ss=2; % 把Ss赋值为2,从而可以退出循环

end

end

set(Ii,'CData',S); % 显示当前状态矩阵

T=T+1; % 累加时间参数

set(ti,'String',['time = ',num2str(T),', N=',num2str(sum(sum(1-S))-[N/4+1])]); % 更新时间参数

pause(0.01); % 暂停一下,显示动画效果

end

3 运行结果

4 参考文献

[1]黄建峰, 杨屹, 李羽晨,等. 基于元胞自动机法的等轴枝晶生长数值模拟[J]. 热加工工艺, 2010, 39(10):5.

[2]陈瑞, 许庆彦, 柳百成. 基于元胞自动机方法的定向凝固枝晶竞争生长数值模拟[J]. 物理学报, 2014(18):12.

博主简介:擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机、雷达通信、无线传感器等多种领域的Matlab仿真,相关matlab代码问题可私信交流。

部分理论引用网络文献,若有侵权联系博主删除。

【元胞自动机】基于元胞自动机模拟晶体生长附matlab代码相关推荐

  1. 【元胞自动机】基于元胞自动机实现高速公路收费站交通流问题附matlab代码

    1 简介 近年来,我国高等级公路发展迅速,截止到2012年底,我国高速公路的通车里程已达到97355公里,比2011年增长14.6%.而我国自1984年以来实行"贷款修路.收费还贷" ...

  2. 【图像增强】基于Frangi滤波器实现血管图像增强附matlab代码

    1 简介 分析视网膜血管结构的变化是诊断和检测糖尿病,高血压等血管类相关疾病的最重要步骤.针对眼底视网膜图像对比度低,受病变区域边界干扰,很难正确提取血管细节的问题提出了一种基于Frangi滤波器的视 ...

  3. 【图像识别】基于卷积神经网络实现手写汉字识别附matlab代码

    1 内容介绍 1.1. 数据集的获取 数据集的获取来自模式识别国家重点实验室共享,这个不解释直接上网址http://www.nlpr.ia.ac.cn/databases/download/featu ...

  4. 【图像加密】基于混沌系统进行灰度图像加密附Matlab代码

    1 简介 ​ 1 基于混沌系统的图像加密解密 Logistic混沌置乱,先不说有多复杂,其实很简单. Logistic函数是源于一个人口统计的动力学系统,其系统方程形式如下: **X(k+1) = u ...

  5. 【车位检测】基于计算机视觉实现停车场空位识别附matlab代码

    1 简介 为便于汽车驾驶员在室外停车场中寻找可用空车位,基于以数据采集,图像处理和目标检测等过程的计算机视觉,开发了室外停车场车位检测实验.​ 2 部分代码 clc; close all; clear ...

  6. 【火灾检测】基于HSV特征实现火灾检测附matlab代码

    1 简介 针对传统火灾监测系统对于大空间的室内场合和开阔的室外环境易失效的问题,提出了一种结合火灾火焰特征和烟雾特征来进行判断的数字图像型火灾监测算法.火焰颜色特征是基于RGB颜色模型中的R,G,B三 ...

  7. 【智能优化算法】基于遗传算法实现城市交通信号优化附matlab代码

    1 简介 本文设计实时优化的配置方案对道路畅通的应急决策管理具有重要意义.本文在分析交通控制基本理论的基础上,根据交叉口的实际情况并考虑信号灯的转换与车辆的启动损失时间,采用四相位对称式放行方案,以车 ...

  8. 【图像分割】基于计算机视觉实现脑肿瘤分割附matlab代码

    ​1 内容介绍 肿瘤被定义为组织的异常生长.脑瘤是一种异常细胞在其中不受控制地生长和繁殖的大量组织,似乎不受控制通过控制正常细胞的机制.脑肿瘤可以是原发性或转移性,恶性或良性.转移性脑肿瘤是一种癌症已 ...

  9. 【物理应用】水下浮动风力涡轮机的尾流诱导动态模拟风场附matlab代码

    1 内容介绍 风力发电机的空气动力学性能是决定风力机安全与效率的最重要因素之一.但由于影响风力机气动性能参数众多,更加高效精确地模拟风力机的气动特性一直是风力机研究的重要发展方向.本研究采用浸入边界法 ...

最新文章

  1. Linux 内核定时器使用 一 低精度定时器
  2. ADSL(Asymmetric Digital Subscriber Loop)技术
  3. Mybatis源码分析--Mapper接口的代理生成原理
  4. Lua学习教程之 可变參数数据打包与解包
  5. PHP发送数据到指定方法,php通过header发送自定义数据方法_php技巧
  6. 机器学习知识总结系列- 知识图谱(0-0)
  7. Intel 64/x86_64/x86/IA-32处理器标志寄存器详解(4) - 32位EFLAGS - 状态标志与控制标志
  8. golang mysql scan_golang 实现 mysql 数据库备份
  9. python获取某天的属于第几周等时间问题
  10. HCIE-Security Day30:IPSec:实验(五)配置基于路由的IPSec PN(采用预共享密钥认证)
  11. 基于51单片机的keli安装方法
  12. WebSql用法(自我记录)
  13. mi自学站【2018年8月】
  14. pytorch BiLSTM+CRF模型实现NER任务
  15. php 死链查询,网站最新死链查询和处理方法
  16. 计算机电缆和屏蔽线的区别,屏蔽线
  17. https免费泛域名证书申请
  18. 声网 X 牛客网 200万场视频面试背后的实时互动技术支撑
  19. 如何去掉图片上的水印,四个步骤,方法超级容易又简单
  20. Ecma 118th GA会议手记

热门文章

  1. JavaScript脚本语言快速入门
  2. [Python 高德地图] API调用学习历程(三)
  3. 【徐禾 政治经济学概论】 读书笔记 §1.3劳动的二因素:具体劳动和抽象劳动(上.价值质的规定性)
  4. thrift 基于window的安装配置
  5. 分布式强化学习之D4PG
  6. npm ERR /usr/bin/git ls-remote -h -t git://github.com/adobe-webplatform/eve.git
  7. 支付宝企业向个人付款实现
  8. 【Python爬虫】猫眼电影榜单Top100
  9. Tomcat 8(一)Tomcat常用配置
  10. 如何根据域名查询出IP地址