无单元galerlin方法中周期边界条件的处理.docxVIP

无单元galerlin方法中周期边界条件的处理 1 周期边界条件 在整个解算区域中，某些微分方程的解将经历如下周期，例如，磁体问题、一些流动问题、聚合物分子的定向分布问题等。在这个问题上，可以通过周期性在周期内分解问题，从而减少问题的规模，简化和计算。如何在有限的区域指定边界条件，直接关系到微分方程的解。通常，可以指定两种边界条件，一种是自由边界条件，另一种是周期性边界条件。当解同时具有某些对称能力时，可以确定对称线或对称面的位置，当解确定对称线或对称面的位置时，解区域可以指定为两个相邻线或对称面之间的区域。当解未形成对称面时，通常只能指定轮廓的边界条件。 无单元Galerkin (element free galerkin, EFG) 方法 EFG方法使用移动最小二乘法 2 移动最小二乘拟合 EFG方法使用移动最小二乘法进行未知函数的逼近和检验函数的构造．移动最小二乘法 在区域?中，未知标量函数u (x) 的移动最小二乘逼近可以表示为 其中h为节点间的平均距离，p 常数基p 式 (1) 中a (x) ={a 其中N为节点总数，u 为保证u 这里 从 (3) 式中求得a (x) 并代入 (1) 式中，便可以得到函数u (x) 的MLS逼近为 其中 3 等效积分形式 假设L为定义在区域?上的任一微分算子，u为满足一定光滑性要求的未知函数，f为已知函数，微分方程的一般形式可以写为 由于 (5) 式在区域?内任意点都满足，因此对任意函数v都有 其中，函数v称为检验函数， 式 (6) 称为微分方程 (5) 的等效积分形式． 一般地，对于给定的定解条件，方程 (5) 很难精确求解，只能采用数值的办法近似求解．设u u 如果对任意检验函数v， (6) 式都成立，此时近似解便为精确解．实际上，不可能也不需要无穷多个检验函数，而是一般将检验函数取为一组基函数的线性组合，即 其中M≥N．将 (7) , (8) 两式代入 (6) 式，并考虑到系数c 上式的意义是通过选择合适的待定参数u 对EFG方法，? 此时 (9) 式具体为 或写成 式 (11) 是以u 计算时可以使用分部积分或Green公式降低L中所含导数的阶次，从而降低对函数u 4 gagerdon变分形式 假设微分方程 (5) 具有如下周期边界条件 其中边界Γ 式中λ为Lagrange乘子 取检验函数为未知函数u的变分δu时，微分方程对应的等效积分形式 (也称Galerkin变分形式) 可以写为 将 (13) 式以变分的形式引入到 (14) 式中得到 δu和δλ分别为u和λ对应的变分函数，且δu∈H 所以同时满足 (16) 和 (17) 式的u便是微分方程的解． 计算时，u、λ分别采用MLS逼近和Lagrange插值近似 式中N 在 (15) 式中，周期边界条件是以弱形式引入的，即周期边界条件是在平均意义下满足的．也可以令边界Γ 在计算过程中，引入 (20) 式等价于在计算 (15) 式中的边界积分时采用节点积分． 5 计算值的示例 5.1 周期边界条件下的稳态分析 对均质棒状分子聚合物稀溶液，把分子的旋转限定在平面内时，聚合物分子在微观的取向概率分布函数将满足如下的一维Fokker-Planck方程 给定周期边界条件 其中ψ (θ, t) 为取向分布概率密度函数，θ∈[0, 2π]为构型变量，κ为宏观流场速度梯度张量，De为无量纲的Deborah数，C (κ, θ) =u×κ·u且u= (cos (θ) , sin (θ) ) 在二维剪切流场和单轴拉伸流场两种情形下，κ分别取 计算中Deborah数取De=10，初值取ψ (θ, 0) =1/2π，取向空间内均匀布置101个节点，每个单元上采用4点Gauss积分，时间步长取为0.005.Fokker-Planck方程的稳态解见图1．记剪切流场和单轴拉伸流场下的稳态数值解分别为ψ 为检验计算结果的正确性，图2给出了对应迎风格式下的差分解ψ 对比图1和图2可以发现，EFG方法的计算结果与差分方法的计算结果吻合．另外从 (23) - (26) 式可以看出，使用EFG方法求解时，概率密度函数的归一性可以得到更好的满足．从而说明了本文提出的周期边界条件处理方法不仅可行而且具有较高的计算精度． 5.2 周期边界条件处理 考虑如下Poisson方程 左右边界上给定Dirichlet边界条件 上下边界上给定周期边界条件 该问题具有解析解 计算时在求解区域[0, 1]×[0, 1]内均匀地布置41×41个节点，背景网格积分采用4×4点Gauss积分，计算结果见图3. 若记数值解为u 计算所得误差为ε=0计算结果表明本文提出的周期边界条件处理方法对二维问题仍然有效． 6 基于lagrange乘子法的周期边界条件 本文针对具有周期边界