comsol中弱解形式的应用详解.pptxVIP

COMSOL中弱解形式的应用等效积分形式和等效积分弱形式（虚位移原理）微分方程域内边界上等效积分形式域内边界上等效积分弱形式COMSOL PDE模式可用于标量方程或系统注意：系数可能会变成更高阶算子系数形式系数对应于常见的物理参数 (例如，扩散、对流等)通式很灵活和紧凑弱形式作为PDE的基础的PDE形式积分形式提供更强大的灵活性 (例如，非标准化边界条件，边界方程耦合等)Lagrange算子显式求解与通式和系数形式相比，很少被采用基于弱解形式的方程式系统作为PDE的基础的一种PDE形式应用更加灵活（例如，非标准化边界条件，边界方程耦合等）形式更加紧凑对变量的连续性要求较低Lagrange算子显式求解适用范围更广泛适于求解非线性多物理场问题系数形式例如：Poisson方程域内边界上域内子域边界上隐含 c=f=h=1 和所有其他系数为 0。系数形式源扩散吸收源对流对流质量阻尼质量系数形式，波动方程初始/热应力弹性力惯性力 (重力)阻尼质量质量密度阻尼系数应力刚性，“弹簧常数”系数形式，输送扩散方程源扩散吸收对流对流源堆积/储存系数形式，稳态方程系数形式，频率响应波动方程扩散Helmholtz项源Helmholtz 方程：波数波长通式－更简练的公式域内边界上Poisson方程相应的通式为其他系数为 0W?W弱形式 (静态)通式乘以试函数 v 并积分左侧分部积分重排记住对于Poisson方程： ?=[-ux -uy], F=1, R=u (u 约束为 0)在“弱”编辑框中输入上面的求解域积分-test(ux)*ux-test(uy)*uy+test(u)*F 在边界上，设置约束：u瞬态弱形式，例子通式乘以试函数 v 并积分包含“del” 表达式的分部积分重排对于Poisson方程：?=[-ux -uy], F=1, R=u (u 约束为 0)在“weak”编辑框中输入上面的积分式：-test(ux)*ux - test(uy)*uy + test(u)*F – da*test(u)*ut边界上设置约束：u案例：输送和表面反应在不同维度耦合物理场 (唯一的)输送 2D + 吸附 1D (完全耦合)弱形式, 边界模式 (PDE)控制方程1D 吸附2D 输送耦合：弱形式PDE，边界file:///I:\..\1.Trainings\My%20Documents\NYC_July19\Ds*(-test(csTx)*csTx-test(csTy)*csTy)+test(cs)*(react_surf-cst)Ds*(-test(csTx)*csTx-test(csTy)*csTy)+test(cs)*(react_surf-cst)网格划分 (局部精细化)结果表面吸附率 (cs)体积浓度 c弹性静力学的弱形式PDE方程乘上试函数并积分分部积分整理得到域内边界一般性问题的弱形式PDE方程乘上试函数并积分分部积分整理得到域内边界弱约束优点精确的通量计算处理非线性约束处理包含微分的约束缺点引入了较多的未知量容易在Jacobian矩阵的主对角线上引入零值（鞍点）不连续约束导致较大的震荡乘子的物理意义使用弱约束Lagrangian乘子被作为独立变量求解精确的流量计算处理非线性约束和带有导数的约束变量名lm1，lm2….使用弱约束在PhysicsProperties中设置弱约束在每个边界条件中确定是否采用弱约束产生新变量，以lm＋数字命名，按照应用模式及其变量的顺序来编号弱约束类型完美（Ideal）和标准的逐点约束类似的边界条件，不涉及物理本质Lagrange乘子对于所有变量对称非完美（Non-Ideal）修改了模型的物理本质Lagrange乘子只应用到指定约束的变量假设模式A(u)和B(v)，其中A的约束Lagrange乘子及试函数变量及试函数弱约束的局限接触边界相同变量的逐点和弱约束不起作用只有Dirichlet边界条件才有效弱约束常导致线性系统的缩放比例问题使用弱约束时常需要减小单元的阶数，以减少冗余的自由度当使用迭代求解器时，如果在矩阵中引入了零对角元，需要采用Vanka或i-LU算法方程式系统中的应用在weak标签中出现的变量进入求解过程其余的为后处理变量相关的Weak项参数设置weak写入弱项公式dweak与时间相关的弱项bnd.weak超弱项，应用于内部不连续边界条件（不存在几何边界）参考模型库中Equation-Based_Models/transport_problemconstr约束（不是弱约束）Constraint type完美、非完美、用户自定义Thank you!可用在标量方程或系统-矢量方程注意：系数可能会变成更高阶算子系数对应于常见的物理参数（例如：扩散、对流等）这里可以手动打开GUI来显示如何找到系数形式的PDE通式非常灵