06 2023 档案
摘要：1. 简介 能够通过 MATLAB与 COMSOL 之间的双向接口进行联合仿真。 1） 使用 MATLAB作为脚本接口来建立 COMSOL模型并进行求解。 2） 从 COMSOL中调用 MATLAB函数，从而在数据预处理、模型操作和后处理中充分利用 MATLAB及其工具箱的所有功能。 2. 工作流程
阅读全文
摘要：1. 简介 We can classify an integration method by its stability and computational effort. The y-axis represents the stability and x-axis represents the c
阅读全文
摘要：这一篇Blog是在A First course in FEM —— matlab代码实现求解传热问题（稳态） 基础上更进一步，求解瞬态传热问题。 两者的区别如下图所示： 1. 问题描述 求解下图图所示叶片的温度场在[0-1200s]时间段内的变化，初始条件：T(0)=25℃。 控制方程为： 2. 模
阅读全文
摘要：这篇文章会将FEM全流程走一遍，包括网格、矩阵组装、求解、后处理。内容是大三时的大作业，今天拿出来回顾下。 1. 问题简介 涡轮机叶片需要冷却以提高涡轮的性能和涡轮叶片的寿命。我们现在考虑一个如上图所示的叶片，叶片处在一个高温环境中，中间通有四个冷却孔。 假设为稳态，那么叶片内导热微分方程为： 内部
阅读全文
摘要：一、介绍 这篇文章介绍一个ODE solver的实现思路，可用于一阶线性系统的求解。 首先，找到一组ODE方程，可以从网上找算例，或者从comsol等软件中导出。 其次，用matlab来求解，看计算结果。 最后，自己实现一个ODE solver，并将结果与matlab的结果进行对比。 二、ODE方程
阅读全文
摘要：一、应用背景 Krylov 子空间方法就是在Krylov子空间中寻找近似解。 二、Krylov子空间 n是原始矩阵的维度，m为降阶后矩阵的维度，通常m<<n。 三、Arnoldi 过程: 计算Km 的一组正交 对应上述MGS过程的python代码： import numpy as np def ar
阅读全文
摘要：相变材料模型： 相态就是由温度判断（大于还是小于相变温度）；潜热使用了等效热容法 在案例库中，有几个连铸工艺（液态金属在通过冷却模具时发生降温凝固）的例子。
阅读全文
摘要：1. 计算开始时，报错“找不到一致的初始值，最后一个时间步不收敛” 解决方法： (1) 使用稳态研究的结果作为瞬态研究的初始值。单个研究可以包含多个步骤，且默认情况下，每个步骤的结果都会作为初始值传递到下一步骤。因此，在瞬态研究步骤之前添加一个稳态步骤，可以先求解稳态假设下的流场，从而为瞬态步骤提供
阅读全文
