找回密码
 注册
查看: 3932|回复: 11

DG中有悬挂点情形的一个问题

[复制链接]
发表于 2012-7-3 08:59:13 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
一般在做自适应时,DG网格就会出现悬挂点,那么这时要用限制器的话,一个三角单元的三个邻单元要怎么确定,希望有经验的能给个提示。
 楼主| 发表于 2012-7-3 16:19:01 | 显示全部楼层
这种网格称为nonconforming grid,看了一些文献,还没有发现这种网格下限制器的使用方法
发表于 2012-7-5 09:27:50 | 显示全部楼层

回复 2# zqb138 的帖子

所以说限制器是个很麻烦的东西 搞人工粘性吧 呵呵
发表于 2012-7-25 21:25:39 | 显示全部楼层

搞DG不要眼睛只盯DG

自适应产生悬挂点,对于DG,和谱有限体积和多矩积分约束有限体积方法中的处理方法一样,你可以参考王志坚老师的文章,当然你还可以参考肖锋老师有关多矩有限体积方面的文章。
       这一类通过单元内部增加自由度来提高精度的方法具有共性。
发表于 2012-8-2 21:59:07 | 显示全部楼层
好像很容易啊,没有那么多问题吧,至少我使用自适应笛卡尔网格时很容易,我想非结构网格也是一样的没有特别的本质区别
发表于 2012-8-2 21:59:49 | 显示全部楼层
好像很容易啊,没有那么多问题吧,至少我使用自适应笛卡尔网格时很容易,我想非结构网格也是一样的没有特别的本质区别
发表于 2012-8-3 17:17:36 | 显示全部楼层
原帖由 liujmljm 于 2012-8-2 13:59 发表
好像很容易啊,没有那么多问题吧,至少我使用自适应笛卡尔网格时很容易,我想非结构网格也是一样的没有特别的本质区别


  你搞得是低阶的吧,呵呵
  如果是DG或者谱之类的方法,可没那么容易的,自适应之后自由度发生了改变,不像普通低阶有限体积(二阶在内)的。
 楼主| 发表于 2012-8-6 12:41:02 | 显示全部楼层
对于非结构三角网格,限制器的设计是基于每个三角元有三个共边的相邻三角元,而自适应过程产生悬挂点之后,每个三角元的共边相邻三角元最少三个,最多六个,每个边有时有两个相邻三角元。这种情况不符合限制器的要求,所以我才提出这个问题的。
 楼主| 发表于 2012-8-14 19:21:43 | 显示全部楼层
不明白,自适应笛卡尔网格下DG的限制器的使用,看了矩形网格下限制器的定义,确实需要上下左右四个矩形内物理量的均值,但是一旦有悬挂点,那么这个矩形单元就会有4个以上的邻矩形单元,已经不符合使用限制器的要求了,怎么使用限制器呢??
 楼主| 发表于 2012-8-22 20:57:46 | 显示全部楼层
今天踩一脚,看看还能不能解决
发表于 2013-1-13 11:10:38 | 显示全部楼层
原帖由 zqb138 于 2012-8-14 19:21 发表
不明白,自适应笛卡尔网格下DG的限制器的使用,看了矩形网格下限制器的定义,确实需要上下左右四个矩形内物理量的均值,但是一旦有悬挂点,那么这个矩形单元就会有4个以上的邻矩形单元,已经不符合使用限制器的要求了 ...

四个邻居可以将他们用L2投影到大网格,这个大网格和你需要的邻居网格加密层次一样,用这个大网格上的L2投影的单元平均用于限制器,非结构网格比如三角形是类似的,如果邻居网格加密层次比当前网格低,同样使用这样的方法,只要将邻居大网格的多项式使用L2投影到一个和当前网格同样加密层次的网格就可。这个没什么难度,主要你需要推导系数公式,但要是你已经有了加密,粗化的L2投影公式,那这个是很容易做到的,可能我这个回复晚了,不知楼主是否解决这个问题

[ 本帖最后由 liujmljm 于 2013-1-13 11:15 编辑 ]
 楼主| 发表于 2013-1-20 12:14:55 | 显示全部楼层
投影方法,很有启发
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表