随机参数
COMSOL Multiphysics
随机现象是真实世界中常见的一种现象,通常是指事前不可预言的现象,即,在相同条件下重复进行试验,每次结果未必相同,或知道事物过去的状况,但未来的发展却不能完全肯定。例如,以同样的方式抛置硬币却可能出现正面向上也可能出现反面向上;走到某十字路口时,可能正好是红灯,也可能正好是绿灯等。
COMSOL Multiphysics中预置了一个随机函数, 可以用来设定随机初始噪音、参数分布等;通过与第三方软件结合,可以模拟随机变化的形状(求解域)、边界等。
Part.01 随机变化的参数
COMSOL Multiphysics
两种互不混溶的物质,可以在一定条件下混合成为稳定、均匀单相状态,然后进行加热或冷却时,到达相图中的两相区域后,就会出现相分离。这个过程中充满了随机性,由于分子运动,两相的浓度不停地发生涨落,在加热或冷却过程中,当浓度的涨落满足一定条件时,就会出现成核、生长等一系列过程, 并最终产生相分离。过去人们常常采用Monte Carlo方法来模拟这种现象,最近十几年,采用相场模型(Phase Field Model)对此现象进行描述,并采用有限元方法进行模拟的分析方式逐渐得到应用。
与当前常见的MonteCarlo模拟方法类似,我们首先需要在指定的结构(求解域)中生成一个围绕着平均浓度随机波动的极小的浓度涨落分布,然后应用相场模型,计算在这种分布情况下两相的运动,从而描述出整个分相过程。
首先,在COMSOL Multiphysics中,定义一个随机函数,设定的随机函数名称为phi_init,需要两个参变量,采用统一分布形式, 围绕着平均值0的涨落范围为±0.05。
然后就可以把这个随机函数应用于初始值的定义中,在相场变量编辑框中输入: phi_init(x,y), 表示在不同的坐标点处求解随机函数。
经过求解,我们得到随时间演化的相图,刚开始,两相是完全混合的,相场变量忽略随机扰动,大约是φ=0。到t=ls相开始分离。在t=2s,开始形成纯相,两秒后,只存在纯相。t=4s 后,纯相开始接合,形成更大的相区域。
PS:后台私信回复“相分离”获取该案例操作
Part.02 随机变化的形状
COMSOL Multiphysics
上面讨论的是随机变化的参数或初始值等,有时候,我们也会遇到结构(求解域)或边界等随机变化的情况,例如随机起伏的地面、在求解域中随机分布的粒子或小区域等。结合COMSOL Multiphysics中的随机函数,动态网格技术,以及利用脚本代码,可以实现复杂的随机变化形状的仿真。
以一个二维传热分析为例,假定研究对象的上边界由于切割的关系,呈现出随机的形状,我们可以通过以下步骤来进行描述。
首先定义一个随机函数,只需要一个参变量,设定为围绕均值0变化的±0.05的统一分布。
在模型中需要耦合一个变形几何应用模式,用来描述边界的变化。其中,设定求解域为自由变形,上边界为指定的随机网格位移。
在指定y位移的编辑框中,输入随机函数ml(x),表示该处的y坐标以随机形式变化。
很明显,上边界通过随机函数演化出现随机的形状,描述出所需模拟的对象。
备注学校专业
邀进COMSOL讨论群
● 扫码添加小编
投稿专栏
COMSOL Multiphysics
“COMSOL仿真交流”微信公众号自创办以来得到了广大科研工作者和研究生的广泛关注和支持。为更好地服务仿真学习与研究,本公众号现因业务需要长期招聘供稿作者。您可发送简历到likeapoem@163.com,邮件主题请注明“姓名+供稿作者”。
欢迎“光声力热”类相关专业的科研爱好者加入“COMSOL仿真交流”团队。
来源:COMSOL仿真交流