随笔分类 - 多物理场仿真
摘要：1. 流程 COMSOL中将PDE转成ODE（瞬态仿真），再通过对时间项离散，最后获得稀疏矩阵方程，通过求解器求解。而稳态仿真则跳过上述时间离散的过程，其余与瞬态仿真求解一致。 流程如下： 瞬态： 稳态: 2. 隐式ODE，及其离散形式 其中，D：质量矩阵；K：刚度矩阵；L{Ω}：载荷；L{part
阅读全文
摘要：一、介绍 1.1 物理含义 雪花是人们最常见的枝晶。枝晶生长是一种生长的不稳定现象，常起因于过冷的液体，或晶体的生长速度受限于溶质原子向固体表面的扩散速度。造成以上条件的原因，可以是液相中负的温度梯度，也可以是负的浓度梯度。在这种模式下，晶体倾向于在其棱角处优先生长，从而形成突触状结构。 这篇博文会
阅读全文
摘要：1. 简介 能够通过 MATLAB与 COMSOL 之间的双向接口进行联合仿真。 1） 使用 MATLAB作为脚本接口来建立 COMSOL模型并进行求解。 2） 从 COMSOL中调用 MATLAB函数，从而在数据预处理、模型操作和后处理中充分利用 MATLAB及其工具箱的所有功能。 2. 工作流程
阅读全文
摘要：这一篇Blog是在A First course in FEM —— matlab代码实现求解传热问题（稳态） 基础上更进一步，求解瞬态传热问题。 两者的区别如下图所示： 1. 问题描述 求解下图图所示叶片的温度场在[0-1200s]时间段内的变化，初始条件：T(0)=25℃。 控制方程为： 2. 模
阅读全文
摘要：这篇文章会将FEM全流程走一遍，包括网格、矩阵组装、求解、后处理。内容是大三时的大作业，今天拿出来回顾下。 1. 问题简介 涡轮机叶片需要冷却以提高涡轮的性能和涡轮叶片的寿命。我们现在考虑一个如上图所示的叶片，叶片处在一个高温环境中，中间通有四个冷却孔。 假设为稳态，那么叶片内导热微分方程为： 内部
阅读全文
摘要：相变材料模型： 相态就是由温度判断（大于还是小于相变温度）；潜热使用了等效热容法 在案例库中，有几个连铸工艺（液态金属在通过冷却模具时发生降温凝固）的例子。
阅读全文
摘要：1. 计算开始时，报错“找不到一致的初始值，最后一个时间步不收敛” 解决方法： (1) 使用稳态研究的结果作为瞬态研究的初始值。单个研究可以包含多个步骤，且默认情况下，每个步骤的结果都会作为初始值传递到下一步骤。因此，在瞬态研究步骤之前添加一个稳态步骤，可以先求解稳态假设下的流场，从而为瞬态步骤提供
阅读全文
摘要：一、 Poisson方程 可转换为线性弱形式 其中， u: trial function 是近似解 v: test function 是Poisson方程两边同乘的函数，用于转换为弱形式。 在Fenics中解PDE时，必须进行下面几个步骤： 1、将PDE问题转换成离散变分问题：即寻找u∈V，使得 2
阅读全文
摘要：在Comsol中，相变可以用变形几何（deformed geometry）处理。切记不能用移动网格，因为移动网格是针对材料守恒的情况，而相变时材料不守恒。 参考链接：https://www.bilibili.com/video/BV1RC4y1W75w/?spm_id_from=333.999.0.
阅读全文
摘要：引子： 一、Maxwell方程 高斯定律（磁）运用于单纯磁场。 高斯定律（电）运用于单纯电场。 安培定律与法拉第定律描述电磁感应的，安培定律是电生磁，法拉第定律是磁生电。 二、电磁理论中的势 其中，B：磁通量；E：电场；H：磁场。 三、研究类型 有周期性的用频域求解，速度和收敛性比瞬态的好。
阅读全文
