终于得到了计算结果

boling2000

从NACA0012到NACA0015， 计算状态为Mach = 0.4, 迎角4度。
采用了梯度下降方法，计算了20次，才得到了计算结果。
红线是目标压力分布和翼型，绿线为设计得到的翼型和压力分布曲线。

周华

翼型为什么这么厚？是不是y方向坐标被放大了？

boling2000

翼型确实放大了，这样才可以看出差别。
另外这是基于控制理论所作的翼型反设计，一位朋友把它当作是欧拉方程计算了。

周华

有个想法：我们可不可以收集一些精彩的算例，加上简单的计算说明（软件、边界条件设定、网格划分等等），就象您这个帖子一样，然后开辟一个专门的栏目供大家欣赏。如果算例比较多的话，我们还可以考虑发行小型电子刊物，这样交流会更充分一些。

boling2000

I saw several websites about the verification and validation in CFD, it includes a lot of classical cases for CFD. most time you could find the experimental data and grid data used in CFD, so I think I could do something like that.

willamscn

请教 boling2000  ，通风算例在那里啊？？

windhair

Nice to see this good result.
As a freshman about design. What kind of restrict parameter was used during design?
Hehe, I just think abuot it, you need not answer this. ;o)

boling2000

Actually what I did here is just one reverse design with given pressure distribution. It is just a validation case.

周华

***** 版主模式 *****
<a href=topic.cgi?forum=15&topic=1>该贴子已被管理员转移，请点击这里查看</a>

zhao11o8

不错哦！！能不能告诉我，边界条件怎不处理，非结构网格上的翼型绕流，怎么搞，请高手执教

streamline

对称翼型、低马赫数的情况容易些。给出历程最好。
梯度如何求得？差分法、伴随方程法？
设计变量：Hicks-Henne sine bump function？
控制方程：Euler、Navier-Stokes/RANS?
给目标压力分布的反设计前途不大！机翼的目标压力分布就不易定。
Jameson在Stanford的网站上，有关文章都可下载。

boling2000

赫赫，谢谢大家的关心，其实那是当时做优化设计时候的一个中间结果。
具体的计算方法，就是重复Jameson的方法
梯度计算采用的是共轭方程求解(adjoint equation)
设计变量采用网格点直接扰动
控制方程为euler方程。

streamline

据说网格点作设计变量，要做光滑处理才能得到光滑外形，用否？
Jameson学生的大多数工作用Hicks-Henne函数扰动外形，尽管他在后来的AIAAPaper中坚持认为网格点可取。
做过梯度精度比较吗？局部极值点烦不烦？
怎么不做下去？这好像是目前CFD的一大方向。

boling2000

谢谢楼上的仁兄，我还在继续做，目前在流场计算和共轭方程求解都引入多重网格，网格点作为设计变量，的确需要光顺，就像残值计算的时候的隐式残值光顺一样，效果还可以。Hicks-Henne函数我也采用了，亚音速的时候，网格点扰动比较好，收敛比较快，但是到了跨音速反倒是Hicks-Henne的效果比较好了。

streamline

  亚声速的反设计，MIT的Direla（？）的XFOIL做得很好了，随便下载（google一搜就出来）。国内某些人还用它挣钱哪。
  文献中有一种叫one-shot的方法，好像把优化循环与多重网格的循环“紧耦合”起来，国内没人做过。国内大多数的气动数值优化设计工作是把Euler/RANS分析嵌入优化循环中，这是一种“松耦合”的做法。