主题介绍
"十多年来,Ansys Fluent adjoint Solver技术逐渐增强,并作为一个完善的工具包的形状优化形式出现。目前的伴随求解器,结合配套工具,可以优化多种工程流体系统的性能,包括可压缩和不可压缩流体的流动、传热问题(包含共轭传热)、多孔介质和旋转流动。Fluent Adjoint Solver成功的关键是具有一个的鲁棒伴随求解器,它能解决各种大小流动问题,并支持广泛的应用场景、数值和边界条件。此外,精心定制制作的网格变形工具可与单/多目标约束优化算法结合起来起来使用。这些工具被集成到一个灵活的、直接的工作流中。 "
如有任何问题请点击以下链接进入答疑室与我们的技术专家进行交流互动
https://v.ansys.com.cn/live/61da320e
演讲人简介
Chris Hill
该演讲为Ansys Simulation World 虚拟大会视频
-
00:00:04.33 - 00:00:06.97 3
大家好
-
00:00:06.97 - 00:00:07.03 12
我是Chris Hill
-
00:00:07.08 - 00:00:09.73 11
我在Ansys的职位是
-
00:00:09.73 - 00:00:10.61 12
流体业务部的首席技术专家
-
00:00:10.66 - 00:00:13.31 3
工作中
-
00:00:13.31 - 00:00:13.72 11
与Min Xu博士搭档
-
00:00:13.78 - 00:00:14.13 18
Min Xu博士同样就职于流体业务部
-
00:00:14.13 - 00:00:16.96 9
担任高级研发工程师
-
00:00:16.96 - 00:00:19.73 6
今天我将讨论
-
00:00:19.73 - 00:00:21.15 9
优化技术中的伴随法
-
00:00:21.15 - 00:00:24.26 16
首先我将为您简要概述伴随求解方法
-
00:00:24.34 - 00:00:27.68 9
它如何应用于CFD
-
00:00:27.68 - 00:00:28.72 6
以及图形优化
-
00:00:28.79 - 00:00:32.13 9
可以对什么进行优化
-
00:00:32.13 - 00:00:32.65 7
支持哪类物理场
-
00:00:32.65 - 00:00:36.26 6
以及如何计算
-
00:00:36.26 - 00:00:37.15 3
伴随解
-
00:00:37.23 - 00:00:40.84 7
我们将一些约束
-
00:00:40.84 - 00:00:41.48 9
应用至设计中的更改
-
00:00:41.56 - 00:00:42.69 5
并应用变体
-
00:00:42.69 - 00:00:45.74 10
然后提供一些优化案例
-
00:00:45.74 - 00:00:46.28 3
最后再
-
00:00:46.35 - 00:00:49.41 9
讲述融入数据分析的
-
00:00:49.41 - 00:00:50.97 4
湍流建模
-
00:00:51.03 - 00:00:52.60 8
这是一个新兴领域
-
00:00:56.05 - 00:00:58.94 9
因此CFD相关问题
-
00:00:58.94 - 00:00:59.59 4
的伴随法
-
00:00:59.65 - 00:01:02.42 9
拥有很大的设计空间
-
00:01:02.48 - 00:01:05.38 9
如果您只有一些参数
-
00:01:05.38 - 00:01:05.58 12
那么伴随法或许并不适合您
-
00:01:05.58 - 00:01:08.82 10
因此这是一款能够提供
-
00:01:08.82 - 00:01:09.54 11
详细的敏感度数据的工具
-
00:01:09.61 - 00:01:12.85 7
此处的经典案例
-
00:01:12.85 - 00:01:13.35 6
如右上图所示
-
00:01:13.42 - 00:01:16.66 7
在此您可以看到
-
00:01:16.66 - 00:01:17.09 5
机翼的升力
-
00:01:17.09 - 00:01:20.33 9
对机翼形状的敏感度
-
00:01:20.40 - 00:01:23.64 14
表明了经典的趋势 即增加拱角
-
00:01:23.64 - 00:01:24.94 4
增加升力
-
00:01:25.01 - 00:01:28.25 4
增加拱角
-
00:01:28.25 - 00:01:29.47 6
并且增加迎角
-
00:01:29.47 - 00:01:31.49 8
从而生成更多升力
-
00:01:31.49 - 00:01:34.17 5
因此伴随法
-
00:01:34.17 - 00:01:34.77 11
能提供一些工程定量导数
-
00:01:34.82 - 00:01:37.51 2
以及
-
00:01:37.51 - 00:01:38.34 8
该问题的所有输入
-
00:01:38.40 - 00:01:41.09 2
包括
-
00:01:41.09 - 00:01:42.22 6
流动边界条件
-
00:01:42.22 - 00:01:45.20 7
几何结构 材料
-
00:01:45.20 - 00:01:46.26 10
甚至是使用的模型参数
-
00:01:46.32 - 00:01:49.04 10
用于定义计算中的建模
-
00:01:49.04 - 00:01:52.52 5
当然您可以
-
00:01:52.52 - 00:01:52.99 13
采用很多方式来应用这些导数
-
00:01:53.07 - 00:01:56.56 7
形状和拓扑修改
-
00:01:56.56 - 00:01:57.64 10
就是一个很不错的案例
-
00:01:57.72 - 00:01:58.73 4
网格调整
-
00:01:58.73 - 00:02:02.18 4
几何形状
-
00:02:02.18 - 00:02:03.95 12
或边界条件变化的假设研究
-
00:02:04.03 - 00:02:06.34 11
鲁棒性设计和鲁棒性优化
-
00:02:06.34 - 00:02:10.42 10
现在机翼问题出在哪里
-
00:02:10.42 - 00:02:11.60 9
或许更加明显 常见
-
00:02:11.69 - 00:02:13.15 2
案例
-
00:02:13.15 - 00:02:16.08 6
如右下方所示
-
00:02:16.08 - 00:02:16.73 7
换热器并不明显
-
00:02:16.80 - 00:02:19.73 2
此时
-
00:02:19.73 - 00:02:20.39 8
我们尽力减少压降
-
00:02:20.45 - 00:02:23.39 8
并且提高传热效率
-
00:02:23.39 - 00:02:23.71 10
因此我们拥有两个对象
-
00:02:23.71 - 00:02:26.65 8
它们可能互相竞争
-
00:02:26.65 - 00:02:27.11 3
但我们
-
00:02:27.17 - 00:02:30.11 8
可以通过优化程序
-
00:02:30.11 - 00:02:31.15 8
设法同时改善二者
-
00:02:31.22 - 00:02:34.02 6
而在本案例中
-
00:02:34.02 - 00:02:37.27 9
我们能够做到这一点
-
00:02:37.27 - 00:02:37.70 10
对设计的改变并不明显
-
00:02:37.77 - 00:02:40.23 8
从而实现这一目标
-
00:02:43.22 - 00:02:47.55 8
因此形状优化流程
-
00:02:47.55 - 00:02:50.15 8
涉及一个优化循环
-
00:02:50.15 - 00:02:53.46 7
而该循环的元素
-
00:02:53.46 - 00:02:53.97 6
是CFD分析
-
00:02:54.05 - 00:02:57.36 13
您可以首先用典型的传统方式
-
00:02:57.36 - 00:02:58.98 7
计算CFD的解
-
00:02:58.98 - 00:03:01.97 11
然后继续运行伴随求解器
-
00:03:01.97 - 00:03:02.57 7
计算一些定量的
-
00:03:02.64 - 00:03:05.04 4
导数信息
-
00:03:07.41 - 00:03:10.57 12
您可以在此获得形状敏感度
-
00:03:10.57 - 00:03:14.09 7
并应用设计约束
-
00:03:14.09 - 00:03:15.43 16
这些设计约束或许会限制敏感度应用
-
00:03:15.50 - 00:03:19.03 6
利用设计工具
-
00:03:19.03 - 00:03:19.97 8
可以更新几何结构
-
00:03:20.05 - 00:03:23.58 5
使用变形器
-
00:03:23.58 - 00:03:24.60 8
继续使用循环处理
-
00:03:24.60 - 00:03:28.44 10
您的几何结构会被更新
-
00:03:28.44 - 00:03:28.61 7
您这样做的时候
-
00:03:28.70 - 00:03:32.54 10
系统性能将会不断提升
-
00:03:32.54 - 00:03:34.08 5
直到您最终
-
00:03:34.17 - 00:03:35.71 9
得到一个局部最优值
-
00:03:35.71 - 00:03:40.86 14
然后您可以将修改后的几何结构
-
00:03:40.86 - 00:03:41.43 3
导出为
-
00:03:41.55 - 00:03:42.01 5
STL文件
-
00:03:42.01 - 00:03:44.81 13
所有这些都由我们基于梯度的
-
00:03:44.81 - 00:03:45.62 5
优化器控制
-
00:03:45.68 - 00:03:48.18 8
您可以在左侧看到
-
00:03:48.24 - 00:03:49.68 7
优化器的主面板
-
00:03:51.76 - 00:03:55.20 5
可以优化的
-
00:03:55.20 - 00:03:56.42 9
通常是一些定量的值
-
00:03:56.50 - 00:03:59.95 4
包括升力
-
00:03:59.95 - 00:04:00.86 5
阻力和压降
-
00:04:00.94 - 00:04:04.39 13
这些都是优化流程的常见目标
-
00:04:04.39 - 00:04:05.38 10
但也可能存在其他目标
-
00:04:05.38 - 00:04:08.83 7
例如流动均匀性
-
00:04:08.83 - 00:04:09.59 5
或许您希望
-
00:04:09.67 - 00:04:13.11 7
执行一些表面或
-
00:04:13.11 - 00:04:14.26 8
体积积分的最小化
-
00:04:14.34 - 00:04:17.78 4
或最大化
-
00:04:17.78 - 00:04:18.93 8
例如域的平均温度
-
00:04:18.93 - 00:04:22.29 6
以及其它选项
-
00:04:22.29 - 00:04:23.49 13
这些都能在伴随求解器中实现
-
00:04:23.56 - 00:04:26.92 9
您可以定义对您而言
-
00:04:26.92 - 00:04:27.90 9
最为重要的一些定量
-
00:04:27.90 - 00:04:31.30 11
这些定量很少是单个对象
-
00:04:31.30 - 00:04:32.36 8
通常您感兴趣的是
-
00:04:32.43 - 00:04:35.84 6
不同操作条件
-
00:04:35.84 - 00:04:37.80 6
下的一个对象
-
00:04:37.88 - 00:04:41.28 7
或不同操作条件
-
00:04:41.28 - 00:04:43.93 6
下的多个对象
-
00:04:43.93 - 00:04:47.90 22
所有这些在Fluent伴随求解器中都受到支持
-
00:04:47.90 - 00:04:53.37 10
例如在该多点扇叶优化
-
00:04:53.37 - 00:04:57.10 4
中显示了
-
00:04:57.10 - 00:04:57.85 7
多个操作点场景
-
00:04:57.85 - 00:05:01.14 7
我们感兴趣的是
-
00:05:01.14 - 00:05:02.60 8
在多个质量流率下
-
00:05:02.67 - 00:05:05.96 7
提高扇叶的效率
-
00:05:05.96 - 00:05:06.18 14
这不仅适用于在某一流率下优化
-
00:05:06.25 - 00:05:09.54 6
而且能够发现
-
00:05:09.54 - 00:05:10.20 9
其它流率下也在减少
-
00:05:10.20 - 00:05:14.01 2
因此
-
00:05:14.01 - 00:05:14.77 4
我们希望
-
00:05:14.86 - 00:05:16.05 7
拓宽提升的范围
-
00:05:16.05 - 00:05:19.35 6
如右下图所示
-
00:05:19.35 - 00:05:20.37 5
您可以看到
-
00:05:20.45 - 00:05:23.75 12
一个几何形状的升力最大化
-
00:05:23.75 - 00:05:24.41 7
以及阻力最小化
-
00:05:27.93 - 00:05:30.91 10
涉及受支持的物理场时
-
00:05:30.91 - 00:05:31.77 12
我们正在研究稳态单相流动
-
00:05:31.84 - 00:05:34.82 9
以及恒定属性的流体
-
00:05:34.82 - 00:05:36.21 9
和可压缩的理想气体
-
00:05:36.28 - 00:05:39.26 6
因此此处案例
-
00:05:39.26 - 00:05:39.72 4
展示的是
-
00:05:39.72 - 00:05:43.33 6
流过飞机机翼
-
00:05:43.33 - 00:05:43.41 5
的超音速流
-
00:05:43.49 - 00:05:46.38 16
这是一个典型的空气动力学设计问题
-
00:05:46.38 - 00:05:49.57 10
我们也可以解决热问题
-
00:05:49.57 - 00:05:49.86 2
包括
-
00:05:49.86 - 00:05:52.88 5
共轭热传递
-
00:05:52.88 - 00:05:54.02 14
流体 固体 多孔和MRF区域
-
00:05:54.09 - 00:05:56.98 12
因此我们可以解决风扇问题
-
00:05:57.05 - 00:06:00.07 12
您可以处理所有类型的网格
-
00:06:00.07 - 00:06:00.28 11
以及绝大多数的边界条件
-
00:06:00.28 - 00:06:03.75 10
我们也有伴随求解模型
-
00:06:03.75 - 00:06:05.29 9
即GEKO湍流模型
-
00:06:05.37 - 00:06:08.84 6
虽然这种模型
-
00:06:08.84 - 00:06:10.08 10
并不能解决所有的问题
-
00:06:10.15 - 00:06:13.63 2
但是
-
00:06:13.63 - 00:06:14.32 6
对于一些问题
-
00:06:14.32 - 00:06:17.65 10
例如求解湍流的伴随解
-
00:06:17.65 - 00:06:17.95 4
意义重大
-
00:06:20.74 - 00:06:23.73 5
涉及到计算
-
00:06:23.73 - 00:06:24.13 4
伴随解时
-
00:06:24.19 - 00:06:27.19 9
存在一定的计算成本
-
00:06:27.19 - 00:06:28.05 5
关于您的解
-
00:06:28.12 - 00:06:31.11 11
您将获得大量有趣的信息
-
00:06:31.11 - 00:06:31.18 6
这涉及到一个
-
00:06:31.18 - 00:06:34.47 7
迭代解提升流程
-
00:06:34.47 - 00:06:36.59 11
能够最小化一些伴随残差
-
00:06:36.66 - 00:06:39.95 2
通常
-
00:06:39.95 - 00:06:41.27 6
我们研究的是
-
00:06:41.34 - 00:06:44.63 8
针对单个流求解的
-
00:06:44.63 - 00:06:45.36 10
相同的计算时间和内存
-
00:06:45.36 - 00:06:48.64 10
因此与原始流求解相比
-
00:06:48.64 - 00:06:49.29 9
花费的时间多了一倍
-
00:06:49.37 - 00:06:52.64 3
但最终
-
00:06:52.64 - 00:06:53.08 4
您将获得
-
00:06:53.15 - 00:06:56.21 10
大量有趣的敏感度信息
-
00:06:58.73 - 00:07:02.11 7
伴随解的鲁棒性
-
00:07:02.11 - 00:07:04.06 13
在过去十年间一直是一大挑战
-
00:07:04.14 - 00:07:07.52 4
随着我们
-
00:07:07.52 - 00:07:08.28 9
研发出了伴随求解器
-
00:07:08.35 - 00:07:11.73 10
我们在很大程度上能够
-
00:07:11.73 - 00:07:12.26 6
克服这些难题
-
00:07:12.26 - 00:07:15.77 6
我们拥有一些
-
00:07:15.77 - 00:07:16.24 12
基于基耗散和残差最小化的
-
00:07:16.32 - 00:07:19.82 5
现代化方案
-
00:07:19.82 - 00:07:21.70 6
如今效果极好
-
00:07:21.77 - 00:07:25.28 7
残差最小化方案
-
00:07:25.28 - 00:07:26.69 10
确实需要一些额外内存
-
00:07:26.69 - 00:07:29.04 7
我们已经实现了
-
00:07:29.04 - 00:07:29.62 6
鼠标轻轻一点
-
00:07:29.67 - 00:07:32.02 6
就能成功解决
-
00:07:32.02 - 00:07:32.29 6
许多伴随问题
-
00:07:36.62 - 00:07:39.66 12
现在 设计中不可避免会有
-
00:07:39.66 - 00:07:40.88 6
一些设计约束
-
00:07:40.95 - 00:07:43.92 9
我们需要有能力解决
-
00:07:43.99 - 00:07:45.55 11
这些问题 我们可以做到
-
00:07:45.55 - 00:07:49.10 11
实际上 出现设计约束时
-
00:07:49.10 - 00:07:49.89 10
我们可以优化这些问题
-
00:07:49.97 - 00:07:53.52 7
典型的设计约束
-
00:07:53.52 - 00:07:56.60 9
涉及对称性或不变量
-
00:07:56.68 - 00:08:00.24 5
光滑度条件
-
00:08:00.24 - 00:08:00.48 9
或许也会应用于变形
-
00:08:00.48 - 00:08:03.86 7
我们通常会发现
-
00:08:03.86 - 00:08:04.38 8
人们渴望运动限制
-
00:08:04.46 - 00:08:07.84 13
因此您可以引入一个边界表面
-
00:08:07.84 - 00:08:08.81 10
来限制几何结构的运动
-
00:08:08.89 - 00:08:12.27 8
当然在设计定义中
-
00:08:12.27 - 00:08:13.25 4
通常需要
-
00:08:13.25 - 00:08:17.94 3
厚表面
-
00:08:17.94 - 00:08:20.70 11
现在您拥有了敏感度数据
-
00:08:20.70 - 00:08:21.38 4
以及约束
-
00:08:21.44 - 00:08:24.21 12
您需要能够将二者巧妙结合
-
00:08:24.27 - 00:08:27.03 6
在出现约束时
-
00:08:27.03 - 00:08:27.53 9
做出最优的设计调整
-
00:08:27.53 - 00:08:32.33 7
做出最优调整后
-
00:08:32.33 - 00:08:32.97 6
利用网格变形
-
00:08:33.08 - 00:08:37.89 7
来更改计算网格
-
00:08:37.89 - 00:08:39.81 7
而网格变形技术
-
00:08:39.92 - 00:08:44.73 11
在幻灯片上的单个项目上
-
00:08:44.73 - 00:08:46.22 8
适用性并不是很好
-
00:08:49.99 - 00:08:50.82 2
但是
-
00:08:50.91 - 00:08:54.67 13
也有数量可观的网格变形技术
-
00:08:54.67 - 00:08:55.76 6
能够发挥作用
-
00:08:55.84 - 00:08:59.61 12
成功调整了几何结构和网格
-
00:08:59.61 - 00:09:01.62 3
而不会
-
00:09:01.62 - 00:09:04.88 7
让网格质量变差
-
00:09:04.88 - 00:09:08.45 15
当然我们必须找到路径返回CAD
-
00:09:08.45 - 00:09:08.53 2
如果
-
00:09:08.61 - 00:09:12.18 7
您无法让新设计
-
00:09:12.18 - 00:09:12.26 7
返回CAD系统
-
00:09:12.34 - 00:09:15.67 7
就无法成功构建
-
00:09:15.67 - 00:09:19.22 13
因此我们需要能够返回CAD
-
00:09:19.22 - 00:09:19.54 7
将调整后的表面
-
00:09:19.62 - 00:09:23.16 8
导出为STL文件
-
00:09:23.16 - 00:09:23.64 10
是一种非常有用的方式
-
00:09:26.10 - 00:09:29.92 10
这就是第一个优化案例
-
00:09:29.92 - 00:09:31.19 8
这是一个内部流动
-
00:09:31.27 - 00:09:35.09 7
我们拥有两个流
-
00:09:35.09 - 00:09:35.60 6
它们彼此混合
-
00:09:35.69 - 00:09:39.42 9
一个热流 一个冷流
-
00:09:39.42 - 00:09:43.68 11
我们需要确保混合最大化
-
00:09:43.68 - 00:09:44.15 4
与此同时
-
00:09:44.24 - 00:09:48.50 10
四个外部流的速率变化
-
00:09:48.50 - 00:09:49.73 2
减弱
-
00:09:49.73 - 00:09:53.25 7
因此我们浏览了
-
00:09:53.25 - 00:09:53.87 7
三十个设计迭代
-
00:09:53.95 - 00:09:54.74 5
的优化周期
-
00:09:54.74 - 00:09:57.90 4
并且计算
-
00:09:57.90 - 00:09:59.74 13
速率变量和温度变量的伴随解
-
00:09:59.81 - 00:10:02.97 6
并将它们混合
-
00:10:02.97 - 00:10:03.40 5
因此您能够
-
00:10:03.47 - 00:10:05.16 10
将速率变化减少35%
-
00:10:05.16 - 00:10:09.33 9
温度变化减少80%
-
00:10:09.33 - 00:10:12.63 11
您可以在温度图左侧看到
-
00:10:12.63 - 00:10:12.85 4
原始形状
-
00:10:12.93 - 00:10:16.24 5
和优化形状
-
00:10:16.24 - 00:10:17.34 7
之间的几何结构
-
00:10:17.41 - 00:10:20.72 11
以及实际几何结构的改变
-
00:10:20.72 - 00:10:21.38 4
微乎其微
-
00:10:21.38 - 00:10:24.31 8
但其影响十分巨大
-
00:10:24.31 - 00:10:27.07 8
而且伴随解为我们
-
00:10:27.07 - 00:10:27.50 7
提供了战略指导
-
00:10:27.56 - 00:10:28.73 8
指导我们如何实现
-
00:10:30.96 - 00:10:34.36 10
这是一个内部流的案例
-
00:10:34.36 - 00:10:37.62 5
我们正试图
-
00:10:37.62 - 00:10:38.28 4
减少压降
-
00:10:38.28 - 00:10:42.54 5
质量流变量
-
00:10:42.54 - 00:10:43.86 7
或速度分布变量
-
00:10:43.95 - 00:10:45.19 4
位于出口
-
00:10:45.19 - 00:10:48.45 6
但是与此同时
-
00:10:48.45 - 00:10:48.53 5
有一个约束
-
00:10:48.60 - 00:10:51.87 6
限制几何结构
-
00:10:51.87 - 00:10:52.16 8
如绿色框左侧所示
-
00:10:52.23 - 00:10:55.50 7
我们的几何结构
-
00:10:55.50 - 00:10:55.79 4
一定不能
-
00:10:55.79 - 00:10:59.52 8
通过这些限制表面
-
00:10:59.52 - 00:11:00.27 5
因此这就像
-
00:11:00.35 - 00:11:02.02 6
一个封装限制
-
00:11:02.02 - 00:11:05.71 4
因此我们
-
00:11:05.71 - 00:11:06.28 11
可以进行一系列设计迭代
-
00:11:06.36 - 00:11:10.05 10
发现几何形状发生变化
-
00:11:10.05 - 00:11:10.87 4
压降减少
-
00:11:10.95 - 00:11:14.65 5
外流的速度
-
00:11:14.65 - 00:11:15.30 4
不断提升
-
00:11:15.30 - 00:11:19.03 4
直到最后
-
00:11:19.03 - 00:11:20.11 11
我们将压降提高了91%
-
00:11:20.11 - 00:11:24.75 9
流变量提高了94%
-
00:11:24.75 - 00:11:25.26 9
这些都是借助优化器
-
00:11:25.37 - 00:11:27.33 5
自动完成的
-
00:11:27.33 - 00:11:30.48 6
当您完成之后
-
00:11:30.48 - 00:11:30.69 6
将右侧的图片
-
00:11:30.76 - 00:11:33.90 9
与左侧的图片相比较
-
00:11:33.90 - 00:11:34.18 6
当您开始之时
-
00:11:34.25 - 00:11:37.40 9
这就是几何结构表面
-
00:11:37.40 - 00:11:37.54 4
的总压力
-
00:11:37.54 - 00:11:41.08 2
因此
-
00:11:41.08 - 00:11:41.47 5
您可以看到
-
00:11:41.55 - 00:11:42.18 8
这些红色的总压力
-
00:11:42.18 - 00:11:44.62 8
数值较高的总压力
-
00:11:44.62 - 00:11:44.94 6
在右边图片中
-
00:11:45.00 - 00:11:45.38 6
已经消失不见
-
00:11:47.95 - 00:11:51.27 7
这是一个典型的
-
00:11:51.27 - 00:11:52.83 13
车辆外部空气动力学优化问题
-
00:11:52.90 - 00:11:56.23 9
此时我们感兴趣的是
-
00:11:56.23 - 00:11:57.11 2
减少
-
00:11:57.19 - 00:11:57.71 6
车辆上的阻力
-
00:11:57.71 - 00:12:01.37 10
通过一连串的设计迭代
-
00:12:01.37 - 00:12:03.00 10
以及对车辆车体的调整
-
00:12:03.08 - 00:12:06.50 8
在左侧您可以看到
-
00:12:06.59 - 00:12:10.25 4
初始形状
-
00:12:10.25 - 00:12:10.99 10
以及左下方的优化形状
-
00:12:10.99 - 00:12:14.90 6
它在车体前部
-
00:12:14.90 - 00:12:15.68 7
引入了一个凹陷
-
00:12:15.77 - 00:12:19.33 12
车体的形状同样进行了调整
-
00:12:19.42 - 00:12:22.38 7
尤其是挡风玻璃
-
00:12:22.38 - 00:12:25.57 8
在右下方的图片中
-
00:12:25.64 - 00:12:28.90 15
您会看到车辆高敏感度区域的图片
-
00:12:28.90 - 00:12:29.62 5
图中描绘了
-
00:12:29.70 - 00:12:32.96 9
后视镜 立柱 轮胎
-
00:12:32.96 - 00:12:36.20 4
以及其它
-
00:12:36.20 - 00:12:36.78 7
有趣的车体区域
-
00:12:39.84 - 00:12:43.60 6
如果变形不足
-
00:12:43.60 - 00:12:44.02 9
那么 我们需要寻求
-
00:12:44.10 - 00:12:47.86 7
大幅的设计更改
-
00:12:47.86 - 00:12:49.87 12
而依靠变形和拓扑无法实现
-
00:12:49.95 - 00:12:53.71 9
通过优化则是可行的
-
00:12:53.71 - 00:12:55.30 5
因此伴随解
-
00:12:55.30 - 00:12:59.07 10
提供了对引入或者移除
-
00:12:59.07 - 00:13:00.32 13
新固体产生影响的敏感度指南
-
00:13:00.41 - 00:13:04.17 3
事实上
-
00:13:04.17 - 00:13:04.50 6
我们拥有一个
-
00:13:04.59 - 00:13:08.35 6
具有三个入口
-
00:13:08.35 - 00:13:08.77 8
和单一出口的案例
-
00:13:08.77 - 00:13:12.35 10
我们希望改变设计区域
-
00:13:12.35 - 00:13:12.83 4
如图所示
-
00:13:12.91 - 00:13:16.50 4
以便最终
-
00:13:16.50 - 00:13:16.81 10
流中将出现更大的漩涡
-
00:13:16.89 - 00:13:20.48 6
不知不觉之间
-
00:13:20.48 - 00:13:20.80 7
以这种固体结束
-
00:13:20.80 - 00:13:25.32 14
该固体看起来非常像一种混合器
-
00:13:25.32 - 00:13:28.90 6
因此这将导致
-
00:13:28.90 - 00:13:29.21 10
不寻常或者意料之外的
-
00:13:29.29 - 00:13:30.65 5
非直观设计
-
00:13:32.67 - 00:13:35.92 8
最后我们就能开展
-
00:13:35.92 - 00:13:37.30 11
融入数据分析的湍流建模
-
00:13:37.37 - 00:13:40.63 5
在该场景中
-
00:13:40.63 - 00:13:41.71 11
我们拥有可用的实验数据
-
00:13:41.78 - 00:13:45.04 7
而且我们注意到
-
00:13:45.04 - 00:13:46.13 13
计算结果和实验结果不尽相同
-
00:13:46.13 - 00:13:49.24 9
因此我们使用伴随解
-
00:13:49.31 - 00:13:52.57 8
计算与湍流模型参
-
00:13:52.57 - 00:13:53.59 9
数相关的公差的导数
-
00:13:56.92 - 00:13:58.22 7
然后对这些参数
-
00:13:58.22 - 00:14:02.27 8
做出局部最优调整
-
00:14:02.27 - 00:14:03.98 2
使其
-
00:14:04.07 - 00:14:04.70 8
与实验数据相匹配
-
00:14:04.70 - 00:14:07.93 11
现在这里的一个关键步骤
-
00:14:07.93 - 00:14:08.22 3
是普及
-
00:14:08.29 - 00:14:11.52 9
该流程中所学的内容
-
00:14:11.52 - 00:14:11.66 2
以便
-
00:14:11.74 - 00:14:14.97 6
湍流建模专家
-
00:14:14.97 - 00:14:16.12 9
能够充分利用其优势
-
00:14:16.12 - 00:14:20.07 10
设计出更好的湍流模型
-
00:14:20.07 - 00:14:21.21 2
此外
-
00:14:21.29 - 00:14:25.24 10
我们可以依靠机器学习
-
00:14:25.24 - 00:14:26.56 12
作指导 进行一次模型调整
-
00:14:26.56 - 00:14:29.47 11
如果现在您将该调整模型
-
00:14:29.47 - 00:14:29.66 8
应用至一系列问题
-
00:14:29.72 - 00:14:32.64 4
您会发现
-
00:14:32.64 - 00:14:33.35 10
这与原始问题截然不同
-
00:14:33.41 - 00:14:36.32 4
您会发现
-
00:14:36.32 - 00:14:37.16 5
与实验结果
-
00:14:37.16 - 00:14:37.81 5
会更加一致
-
00:14:40.82 - 00:14:43.92 14
现在 我们已经了解了伴随方法
-
00:14:43.92 - 00:14:44.34 6
您对其概览和
-
00:14:44.41 - 00:14:47.52 10
优化案例已经胸中有数
-
00:14:47.52 - 00:14:48.76 6
最后以对融入
-
00:14:48.83 - 00:14:51.94 11
数据分析的湍流建模进行
-
00:14:51.94 - 00:14:53.39 10
一场简单的讨论来结束
-
00:14:53.39 - 00:14:54.01 6
感谢您的关注
-
00:00:04.33 - 00:00:06.97 35
Hello my name is Chris Hill I‘m the
-
00:00:06.97 - 00:00:07.03 5
chief
-
00:00:07.08 - 00:00:09.73 41
technologist for the fluids business unit
-
00:00:09.73 - 00:00:10.61 16
here at Ansys my
-
00:00:10.66 - 00:00:13.31 42
collaborator in this work is doctor Min Xu
-
00:00:13.31 - 00:00:13.72 8
who is a
-
00:00:13.78 - 00:00:14.13 6
senior
-
00:00:14.13 - 00:00:16.96 40
R&D engineer also in the fluids business
-
00:00:16.96 - 00:00:19.73 48
I'm going to talk today about adjoint methods in
-
00:00:19.73 - 00:00:21.15 23
optimization technology
-
00:00:21.15 - 00:00:24.26 45
First of all I'll give you an overview of the
-
00:00:24.34 - 00:00:27.68 40
adjoint method how it applies to CFD and
-
00:00:27.68 - 00:00:28.72 18
shape optimization
-
00:00:28.79 - 00:00:32.13 41
what kind of things can be optimized what
-
00:00:32.13 - 00:00:32.65 10
physics is
-
00:00:32.65 - 00:00:36.26 36
supported how an adjoint solution is
-
00:00:36.26 - 00:00:37.15 16
computed and how
-
00:00:37.23 - 00:00:40.84 42
we apply constraints to the changes in the
-
00:00:40.84 - 00:00:41.48 10
design and
-
00:00:41.56 - 00:00:42.69 14
apply morphing
-
00:00:42.69 - 00:00:45.74 44
I move on to give some optimization examples
-
00:00:45.74 - 00:00:46.28 8
and then
-
00:00:46.35 - 00:00:49.41 41
close with a description of data informed
-
00:00:49.41 - 00:00:50.97 25
turbulence modeling which
-
00:00:51.03 - 00:00:52.60 23
is an emerging new area
-
00:00:56.05 - 00:00:58.94 38
So the adjoint method for CFD involves
-
00:00:58.94 - 00:00:59.59 13
problems with
-
00:00:59.65 - 00:01:02.42 43
a large design space if you only have a few
-
00:01:02.48 - 00:01:05.38 45
parameters then adjoints probably are not for
-
00:01:05.38 - 00:01:05.58 3
you
-
00:01:05.58 - 00:01:08.82 34
So it's a tool to provide detailed
-
00:01:08.82 - 00:01:09.54 20
sensitivity data and
-
00:01:09.61 - 00:01:12.85 40
the classic example here is shown in the
-
00:01:12.85 - 00:01:13.35 11
upper right
-
00:01:13.42 - 00:01:16.66 43
figure where you see the sensitivity of the
-
00:01:16.66 - 00:01:17.09 7
lift on
-
00:01:17.09 - 00:01:20.33 45
an airfoil to the shape of the airfoil and it
-
00:01:20.40 - 00:01:23.64 44
shows the classic trend of increasing camber
-
00:01:23.64 - 00:01:24.94 19
and increasing lift
-
00:01:25.01 - 00:01:28.25 41
increasing camber and increasing angle of
-
00:01:28.25 - 00:01:29.47 18
attack to generate
-
00:01:29.47 - 00:01:31.49 9
more lift
-
00:01:31.49 - 00:01:34.17 34
So the adjointmethod gives you the
-
00:01:34.17 - 00:01:34.77 19
derivatives of some
-
00:01:34.82 - 00:01:37.51 45
engineering quantity of interest with respect
-
00:01:37.51 - 00:01:38.34 13
to all of the
-
00:01:38.40 - 00:01:41.09 44
inputs to the problem this includes the flow
-
00:01:41.09 - 00:01:42.22 19
boundary conditions
-
00:01:42.22 - 00:01:45.20 45
the geometry the materials and even the model
-
00:01:45.20 - 00:01:46.26 15
parameters used
-
00:01:46.32 - 00:01:49.04 41
to define the modeling in the calculation
-
00:01:49.04 - 00:01:52.52 43
And of course you can use these derivatives
-
00:01:52.52 - 00:01:52.99 7
in many
-
00:01:53.07 - 00:01:56.56 42
ways shape and topology modification are a
-
00:01:56.56 - 00:01:57.64 17
very good example
-
00:01:57.72 - 00:01:58.73 13
mesh adaption
-
00:01:58.73 - 00:02:02.18 40
What if studies for geometry or boundary
-
00:02:02.18 - 00:02:03.95 27
condition variations robust
-
00:02:04.03 - 00:02:06.34 30
design and robust optimization
-
00:02:06.34 - 00:02:10.42 42
Now whereas the airfoil problem is perhaps
-
00:02:10.42 - 00:02:11.60 14
rather obvious
-
00:02:11.69 - 00:02:13.15 17
and classical the
-
00:02:13.15 - 00:02:16.08 42
example shown in the lower right over heat
-
00:02:16.08 - 00:02:16.73 12
exchanger is
-
00:02:16.80 - 00:02:19.73 42
not so obvious here we're trying to reduce
-
00:02:19.73 - 00:02:20.39 12
the pressure
-
00:02:20.45 - 00:02:23.39 44
drop and increase the heat transfer rate and
-
00:02:23.39 - 00:02:23.71 5
so we
-
00:02:23.71 - 00:02:26.65 45
have two objectives they may be competing but
-
00:02:26.65 - 00:02:27.11 6
we can
-
00:02:27.17 - 00:02:30.11 42
go through an optimization procedure in an
-
00:02:30.11 - 00:02:31.15 18
attempt to improve
-
00:02:31.22 - 00:02:34.02 43
both of them and in this case we're able to
-
00:02:34.02 - 00:02:37.27 44
do that and it produces rather a non obvious
-
00:02:37.27 - 00:02:37.70 6
change
-
00:02:37.77 - 00:02:40.23 34
to the design to achieve that goal
-
00:02:43.22 - 00:02:47.55 35
So the shape optimization procedure
-
00:02:47.55 - 00:02:50.15 29
involves an optimization loop
-
00:02:50.15 - 00:02:53.46 37
and the elements of that loop are CFD
-
00:02:53.46 - 00:02:53.97 14
analysis where
-
00:02:54.05 - 00:02:57.36 39
you first compute the CFD solution in a
-
00:02:57.36 - 00:02:58.98 28
classical traditional manner
-
00:02:58.98 - 00:03:01.97 43
And then proceed to run the adjoints solver
-
00:03:01.97 - 00:03:02.57 10
to compute
-
00:03:02.64 - 00:03:05.04 42
derivative information for some quantities
-
00:03:07.41 - 00:03:10.57 44
From there you can get the shape sensitivity
-
00:03:10.57 - 00:03:14.09 43
and apply design constraints that may limit
-
00:03:14.09 - 00:03:15.43 18
the application of
-
00:03:15.50 - 00:03:19.03 42
that sensitivity using the design tool you
-
00:03:19.03 - 00:03:19.97 14
can update the
-
00:03:20.05 - 00:03:23.58 43
geometry using a morpher and proceed around
-
00:03:23.58 - 00:03:24.60 14
the loop again
-
00:03:24.60 - 00:03:28.44 44
with your new updated geometry and as you do
-
00:03:28.44 - 00:03:28.61 2
so
-
00:03:28.70 - 00:03:32.54 44
the performance of the system improves until
-
00:03:32.54 - 00:03:34.08 18
eventually you get
-
00:03:34.17 - 00:03:35.71 18
to a local optimum
-
00:03:35.71 - 00:03:40.86 43
Then you can export the revised geometry as
-
00:03:40.86 - 00:03:41.43 6
an STL
-
00:03:41.55 - 00:03:42.01 4
file
-
00:03:42.01 - 00:03:44.81 41
All of this is controlled by our gradient
-
00:03:44.81 - 00:03:45.62 15
based optimizer
-
00:03:45.68 - 00:03:48.18 40
that you can see on the left as the main
-
00:03:48.24 - 00:03:49.68 23
panel for the optimizer
-
00:03:51.76 - 00:03:55.20 38
So what can be optimized are generally
-
00:03:55.20 - 00:03:56.42 22
quantities of interest
-
00:03:56.50 - 00:03:59.95 37
include things like lift and drag and
-
00:03:59.95 - 00:04:00.86 19
pressure drop those
-
00:04:00.94 - 00:04:04.39 44
are common goals for an optimization process
-
00:04:04.39 - 00:04:05.38 13
but there can
-
00:04:05.38 - 00:04:08.83 42
be many others for example flow uniformity
-
00:04:08.83 - 00:04:09.59 12
you may want
-
00:04:09.67 - 00:04:13.11 44
to perform a minimization or maximization of
-
00:04:13.11 - 00:04:14.26 15
some surface or
-
00:04:14.34 - 00:04:17.78 44
volume integral like the mean temperature in
-
00:04:17.78 - 00:04:18.93 15
your domain and
-
00:04:18.93 - 00:04:22.29 41
the many other options that are available
-
00:04:22.29 - 00:04:23.49 19
within the adjoints
-
00:04:23.56 - 00:04:26.92 47
solver for you to define the quantity that most
-
00:04:26.92 - 00:04:27.90 14
matters to you
-
00:04:27.90 - 00:04:31.30 44
Rarely is it a single objective that matters
-
00:04:31.30 - 00:04:32.36 14
usually you're
-
00:04:32.43 - 00:04:35.84 37
interested either in one objective at
-
00:04:35.84 - 00:04:37.80 33
different operating conditions or
-
00:04:37.88 - 00:04:41.28 45
multiple objectives or multiple objectives at
-
00:04:41.28 - 00:04:43.93 34
different operating conditions all
-
00:04:43.93 - 00:04:47.90 40
of these are supported in fluent adjoint
-
00:04:47.90 - 00:04:53.37 43
So for example the multiple operating point
-
00:04:53.37 - 00:04:57.10 41
scenario is shown in this multi point fan
-
00:04:57.10 - 00:04:57.85 12
optimization
-
00:04:57.85 - 00:05:01.14 45
We're interested in increasing the efficiency
-
00:05:01.14 - 00:05:02.60 24
of the fan with multiple
-
00:05:02.67 - 00:05:05.96 45
mass flow rates it would not simply be useful
-
00:05:05.96 - 00:05:06.18 2
to
-
00:05:06.25 - 00:05:09.54 43
optimize at one flow rate and find that the
-
00:05:09.54 - 00:05:10.20 10
efficiency
-
00:05:10.20 - 00:05:14.01 45
has deteriorated at others so we want a broad
-
00:05:14.01 - 00:05:14.77 8
spectrum
-
00:05:14.86 - 00:05:16.05 14
of improvement
-
00:05:16.05 - 00:05:19.35 42
Likewise on the lower right figure you see
-
00:05:19.35 - 00:05:20.37 16
the maximization
-
00:05:20.45 - 00:05:23.75 38
of lift and minimization of drag for a
-
00:05:23.75 - 00:05:24.41 15
geometric shape
-
00:05:27.93 - 00:05:30.91 40
When it comes to supported physics we're
-
00:05:30.91 - 00:05:31.77 17
looking at steady
-
00:05:31.84 - 00:05:34.82 41
single phase flows with constant property
-
00:05:34.82 - 00:05:36.21 24
fluid and a compressible
-
00:05:36.28 - 00:05:39.26 42
ideal gas so the example here is shown for
-
00:05:39.26 - 00:05:39.72 9
transonic
-
00:05:39.72 - 00:05:43.33 43
flow over an airfoil or wing of an aircraft
-
00:05:43.33 - 00:05:43.41 2
is
-
00:05:43.49 - 00:05:46.38 36
a classic aerodynamic design problem
-
00:05:46.38 - 00:05:49.57 45
We can also handle thermal problems including
-
00:05:49.57 - 00:05:49.86 3
and
-
00:05:49.86 - 00:05:52.88 43
conjugate heat transfer we're going to have
-
00:05:52.88 - 00:05:54.02 18
fluid solid porous
-
00:05:54.09 - 00:05:56.98 43
and MRF zones so we can do fan problems you
-
00:05:57.05 - 00:06:00.07 45
can handle all kinds of meshes and most kinds
-
00:06:00.07 - 00:06:00.28 2
of
-
00:06:00.28 - 00:06:03.75 44
boundary conditions We also have the adjoint
-
00:06:03.75 - 00:06:05.29 19
the GEKO turbulence
-
00:06:05.37 - 00:06:08.84 42
model available it's not something that is
-
00:06:08.84 - 00:06:10.08 18
necessarily a good
-
00:06:10.15 - 00:06:13.63 35
idea for all problems but there are
-
00:06:13.63 - 00:06:14.32 14
problems where
-
00:06:14.32 - 00:06:17.65 40
solving the adjoint to the turbulence is
-
00:06:17.65 - 00:06:17.95 8
valuable
-
00:06:20.74 - 00:06:23.73 37
When it comes to computing in a joint
-
00:06:23.73 - 00:06:24.13 8
solution
-
00:06:24.19 - 00:06:27.19 44
there is a computational cost you're getting
-
00:06:27.19 - 00:06:28.05 13
a large block
-
00:06:28.12 - 00:06:31.11 37
of interesting information about your
-
00:06:31.11 - 00:06:31.18 8
solution
-
00:06:31.18 - 00:06:34.47 37
And it involves an iterative solution
-
00:06:34.47 - 00:06:36.59 36
advancement procedure that minimizes
-
00:06:36.66 - 00:06:39.95 42
some adjoint residuals and typically we're
-
00:06:39.95 - 00:06:41.27 20
looking at about the
-
00:06:41.34 - 00:06:44.63 43
same computational time and memory as for a
-
00:06:44.63 - 00:06:45.36 11
single flow
-
00:06:45.36 - 00:06:48.64 42
solution so it's kind of doubles your time
-
00:06:48.64 - 00:06:49.29 11
compared to
-
00:06:49.37 - 00:06:52.64 41
the original flow solution but you end up
-
00:06:52.64 - 00:06:53.08 11
with a vast
-
00:06:53.15 - 00:06:56.21 44
array of interesting sensitivity information
-
00:06:58.73 - 00:07:02.11 35
Robustness of adjoint solutions has
-
00:07:02.11 - 00:07:04.06 35
traditionally been a challenge over
-
00:07:04.14 - 00:07:07.52 40
the last 10 years as we've developed the
-
00:07:07.52 - 00:07:08.28 14
adjoint solver
-
00:07:08.35 - 00:07:11.73 35
we have been able to overcome these
-
00:07:11.73 - 00:07:12.26 6
issues
-
00:07:12.26 - 00:07:15.77 45
to a large extent we have some modern schemes
-
00:07:15.77 - 00:07:16.24 5
based
-
00:07:16.32 - 00:07:19.82 45
on dissipation and residual minimization that
-
00:07:19.82 - 00:07:21.70 23
now work extremely well
-
00:07:21.77 - 00:07:25.28 45
the residual minimization scheme does require
-
00:07:25.28 - 00:07:26.69 17
some extra memory
-
00:07:26.69 - 00:07:29.04 42
We've reached the point where many adjoint
-
00:07:29.04 - 00:07:29.62 15
problems can be
-
00:07:29.67 - 00:07:32.02 44
solved successfully with just a single mouse
-
00:07:32.02 - 00:07:32.29 5
click
-
00:07:36.62 - 00:07:39.66 43
Now inevitably there are design constraints
-
00:07:39.66 - 00:07:40.88 19
that are applied to
-
00:07:40.95 - 00:07:43.92 44
in design problems and we need to be able to
-
00:07:43.99 - 00:07:45.55 23
handle those and we can
-
00:07:45.55 - 00:07:49.10 43
So we can actually optimize problems in the
-
00:07:49.10 - 00:07:49.89 11
presence of
-
00:07:49.97 - 00:07:53.52 38
design constraints typical constraints
-
00:07:53.52 - 00:07:56.60 45
involve symmetry or invariants and conditions
-
00:07:56.68 - 00:08:00.24 35
of smoothness may also apply on the
-
00:08:00.24 - 00:08:00.48 12
deformations
-
00:08:00.48 - 00:08:03.86 39
We often find that limits on motion are
-
00:08:03.86 - 00:08:04.38 12
desirable so
-
00:08:04.46 - 00:08:07.84 39
you can introduce a bounding surface to
-
00:08:07.84 - 00:08:08.81 18
confine the motion
-
00:08:08.89 - 00:08:12.27 45
of your geometry and of course thick surfaces
-
00:08:12.27 - 00:08:13.25 12
are commonly
-
00:08:13.25 - 00:08:17.94 33
required in the design definition
-
00:08:17.94 - 00:08:20.70 40
So now you have the sensitivity data and
-
00:08:20.70 - 00:08:21.38 15
constraints and
-
00:08:21.44 - 00:08:24.21 45
you need to be able to reconcile these two to
-
00:08:24.27 - 00:08:27.03 44
come up with an optimal design change in the
-
00:08:27.03 - 00:08:27.53 8
presence
-
00:08:27.53 - 00:08:32.33 40
of the constraints Once you have that we
-
00:08:32.33 - 00:08:32.97 10
apply mesh
-
00:08:33.08 - 00:08:37.89 45
morphing to change the computational mesh and
-
00:08:37.89 - 00:08:39.81 17
the mesh morphing
-
00:08:39.92 - 00:08:44.73 42
technology is not represented very well by
-
00:08:44.73 - 00:08:46.22 16
that single item
-
00:08:46.22 - 00:08:49.99 42
line item on the slide there but there's a
-
00:08:49.99 - 00:08:50.82 12
considerable
-
00:08:50.91 - 00:08:54.67 45
amount of mesh morphing technology that comes
-
00:08:54.67 - 00:08:55.76 12
into play to
-
00:08:55.84 - 00:08:59.61 45
successfully modify the geometry and the mesh
-
00:08:59.61 - 00:09:01.62 24
without introducing poor
-
00:09:01.62 - 00:09:04.88 21
quality into the mesh
-
00:09:04.88 - 00:09:08.45 43
We must of course find pathways back to CAD
-
00:09:08.45 - 00:09:08.53 2
if
-
00:09:08.61 - 00:09:12.18 42
you can't get the new design back into the
-
00:09:12.18 - 00:09:12.26 3
CAD
-
00:09:12.34 - 00:09:15.67 42
system you can't build it so we need to be
-
00:09:15.67 - 00:09:19.22 40
able to get back to CAD and exporting of
-
00:09:19.22 - 00:09:19.54 8
modified
-
00:09:19.62 - 00:09:23.16 43
surfaces as STL files is one way that works
-
00:09:23.16 - 00:09:23.64 10
quite well
-
00:09:26.10 - 00:09:29.92 45
So here's the first optimization example this
-
00:09:29.92 - 00:09:31.19 14
is an internal
-
00:09:31.27 - 00:09:35.09 44
flow where we have two flow streams that are
-
00:09:35.09 - 00:09:35.60 6
mixing
-
00:09:35.69 - 00:09:39.42 44
one hot one cool and we would like to ensure
-
00:09:39.42 - 00:09:43.68 37
that the mixing is maximized but also
-
00:09:43.68 - 00:09:44.15 8
that the
-
00:09:44.24 - 00:09:48.50 38
variance of the velocity of those four
-
00:09:48.50 - 00:09:49.73 19
outflows is reduced
-
00:09:49.73 - 00:09:53.25 43
So we go through this optimization cycle of
-
00:09:53.25 - 00:09:53.87 9
30 design
-
00:09:53.95 - 00:09:54.74 10
iterations
-
00:09:54.74 - 00:09:57.90 37
And compute the adjoints for velocity
-
00:09:57.90 - 00:09:59.74 33
variance and temperature variance
-
00:09:59.81 - 00:10:02.97 44
and blend them together and you were able to
-
00:10:02.97 - 00:10:03.40 6
reduce
-
00:10:03.47 - 00:10:05.16 24
velocity variance by 35%
-
00:10:05.16 - 00:10:09.33 35
and the temperature variance by 80%
-
00:10:09.33 - 00:10:12.63 45
And you see on the left the temperature plots
-
00:10:12.63 - 00:10:12.85 2
on
-
00:10:12.93 - 00:10:16.24 41
the geometry between the original and the
-
00:10:16.24 - 00:10:17.34 15
optimized shape
-
00:10:17.41 - 00:10:20.72 40
and the change to the actual geometry is
-
00:10:20.72 - 00:10:21.38 13
rather subtle
-
00:10:21.38 - 00:10:24.31 33
here but it has a dramatic effect
-
00:10:24.31 - 00:10:27.07 47
And the adjoint gives us the strategic guidance
-
00:10:27.07 - 00:10:27.50 6
on how
-
00:10:27.56 - 00:10:28.73 19
to make that happen
-
00:10:30.96 - 00:10:34.36 37
Here's an example of an internal flow
-
00:10:34.36 - 00:10:37.62 32
Where we're trying to reduce the
-
00:10:37.62 - 00:10:38.28 13
pressure drop
-
00:10:38.28 - 00:10:42.54 42
And the mass flow variance or the velocity
-
00:10:42.54 - 00:10:43.86 16
profile variance
-
00:10:43.95 - 00:10:45.19 13
at the outlet
-
00:10:45.19 - 00:10:48.45 41
But at the same time we have a constraint
-
00:10:48.45 - 00:10:48.53 4
that
-
00:10:48.60 - 00:10:51.87 44
there is a limiting geometry as shown on the
-
00:10:51.87 - 00:10:52.16 4
left
-
00:10:52.23 - 00:10:55.50 45
side by the green boxes and our geometry must
-
00:10:55.50 - 00:10:55.79 3
not
-
00:10:55.79 - 00:10:59.52 44
pass through those limiting surfaces so this
-
00:10:59.52 - 00:11:00.27 9
is like a
-
00:11:00.35 - 00:11:02.02 20
packaging limitation
-
00:11:02.02 - 00:11:05.71 41
So we can go through a sequence of design
-
00:11:05.71 - 00:11:06.28 10
iterations
-
00:11:06.36 - 00:11:10.05 45
and we find the geometry changes the pressure
-
00:11:10.05 - 00:11:10.87 9
drop goes
-
00:11:10.95 - 00:11:14.65 44
down the form of the velocity at the outflow
-
00:11:14.65 - 00:11:15.30 8
improves
-
00:11:15.30 - 00:11:19.03 37
until eventually we have improved the
-
00:11:19.03 - 00:11:20.11 20
pressure drop by 91%
-
00:11:20.11 - 00:11:24.75 40
on the flow variance by 94 this is fully
-
00:11:24.75 - 00:11:25.26 9
automated
-
00:11:25.37 - 00:11:27.33 19
using the optimizer
-
00:11:27.33 - 00:11:30.48 44
And when you're done you have the picture on
-
00:11:30.48 - 00:11:30.69 3
the
-
00:11:30.76 - 00:11:33.90 44
right as compared to the picture on the left
-
00:11:33.90 - 00:11:34.18 4
when
-
00:11:34.25 - 00:11:37.40 43
you began and this is the total pressure on
-
00:11:37.40 - 00:11:37.54 3
the
-
00:11:37.54 - 00:11:41.08 44
surface of the geometry so you see those red
-
00:11:41.08 - 00:11:41.47 5
total
-
00:11:41.55 - 00:11:42.18 8
pressure
-
00:11:42.18 - 00:11:44.62 40
high total pressures are now gone in the
-
00:11:44.62 - 00:11:44.94 10
right hand
-
00:11:45.00 - 00:11:45.38 7
picture
-
00:11:47.95 - 00:11:51.27 39
This is a classic external aerodynamics
-
00:11:51.27 - 00:11:52.83 26
optimization problem for a
-
00:11:52.90 - 00:11:56.23 45
vehicle here we're interested in reducing the
-
00:11:56.23 - 00:11:57.11 11
drag on the
-
00:11:57.19 - 00:11:57.71 7
vehicle
-
00:11:57.71 - 00:12:01.37 40
And by sequence of design iterations and
-
00:12:01.37 - 00:12:03.00 24
modifications being made
-
00:12:03.08 - 00:12:06.50 42
to the hood of the vehicle and on the left
-
00:12:06.59 - 00:12:10.25 38
you see the initial and lower left the
-
00:12:10.25 - 00:12:10.99 15
optimized forms
-
00:12:10.99 - 00:12:14.90 45
So it introduces an interesting dimple on the
-
00:12:14.90 - 00:12:15.68 8
front of
-
00:12:15.77 - 00:12:19.33 41
the hood and also the form of the hood is
-
00:12:19.42 - 00:12:22.38 34
modified up against the windshield
-
00:12:22.38 - 00:12:25.57 44
The lower right figure you see a plot of the
-
00:12:25.64 - 00:12:28.90 43
high sensitivity regions on the vehicle and
-
00:12:28.90 - 00:12:29.62 11
so it picks
-
00:12:29.70 - 00:12:32.96 45
out the wing mirror and the a pillars and the
-
00:12:32.96 - 00:12:36.20 42
tires and other interesting regions of the
-
00:12:36.20 - 00:12:36.78 7
vehicle
-
00:12:39.84 - 00:12:43.60 45
If morphing is not sufficient then we need to
-
00:12:43.60 - 00:12:44.02 4
look
-
00:12:44.10 - 00:12:47.86 42
for dramatic design changes not achievable
-
00:12:47.86 - 00:12:49.87 26
with morphing and topology
-
00:12:49.95 - 00:12:53.71 41
optimization offers such an avenue so the
-
00:12:53.71 - 00:12:55.30 25
adjoint solutions provide
-
00:12:55.30 - 00:12:59.07 44
guidance on the sensitivity to the effect of
-
00:12:59.07 - 00:13:00.32 15
introducing new
-
00:13:00.41 - 00:13:04.17 41
solids or removing them indeed so here we
-
00:13:04.17 - 00:13:04.50 7
have an
-
00:13:04.59 - 00:13:08.35 43
example of a plenum with three inlets and a
-
00:13:08.35 - 00:13:08.77 6
single
-
00:13:08.77 - 00:13:12.35 44
outlet and we look for changes to the design
-
00:13:12.35 - 00:13:12.83 6
region
-
00:13:12.91 - 00:13:16.50 43
as shown such that will end up with greater
-
00:13:16.50 - 00:13:16.81 5
swirl
-
00:13:16.89 - 00:13:20.48 44
in the flow and spontaneously we end up with
-
00:13:20.48 - 00:13:20.80 4
this
-
00:13:20.80 - 00:13:25.32 40
solid that looks remarkably like a mixer
-
00:13:25.32 - 00:13:28.90 45
So this can lead to unusual or unexpected and
-
00:13:28.90 - 00:13:29.21 3
non
-
00:13:29.29 - 00:13:30.65 17
intuitive designs
-
00:13:32.67 - 00:13:35.92 44
Finally we have the subject of data informed
-
00:13:35.92 - 00:13:37.30 19
turbulence modeling
-
00:13:37.37 - 00:13:40.63 41
and in this scenario we have experimental
-
00:13:40.63 - 00:13:41.71 18
data available and
-
00:13:41.78 - 00:13:45.04 44
we notice that the computed and experimental
-
00:13:45.04 - 00:13:46.13 15
results are not
-
00:13:46.13 - 00:13:49.24 42
the same and so we use the adjoint compute
-
00:13:49.31 - 00:13:52.57 38
the derivative of that difference with
-
00:13:52.57 - 00:13:53.59 20
respect to the model
-
00:13:53.66 - 00:13:56.92 44
parameters for the turbulence model and then
-
00:13:56.92 - 00:13:58.22 18
make local optimal
-
00:13:58.22 - 00:14:02.27 43
adjustments to those parameters so that the
-
00:14:02.27 - 00:14:03.98 20
experimental data is
-
00:14:04.07 - 00:14:04.70 7
matched
-
00:14:04.70 - 00:14:07.93 44
Now a key step here is the generalization of
-
00:14:07.93 - 00:14:08.22 4
what
-
00:14:08.29 - 00:14:11.52 41
has been learned by this process and so a
-
00:14:11.52 - 00:14:11.66 5
human
-
00:14:11.74 - 00:14:14.97 43
expert in turbulence modeling can take good
-
00:14:14.97 - 00:14:16.12 17
advantage of this
-
00:14:16.12 - 00:14:20.07 41
to come up with a better turbulence model
-
00:14:20.07 - 00:14:21.21 16
Alternatively we
-
00:14:21.29 - 00:14:25.24 41
could look to machine learning to guide a
-
00:14:25.24 - 00:14:26.56 18
model modification
-
00:14:26.56 - 00:14:29.47 41
If you now apply this modified model to a
-
00:14:29.47 - 00:14:29.66 6
family
-
00:14:29.72 - 00:14:32.64 43
of problems you find that are distinct from
-
00:14:32.64 - 00:14:33.35 12
the original
-
00:14:33.41 - 00:14:36.32 43
problems you can find that you get improved
-
00:14:36.32 - 00:14:37.16 14
agreement with
-
00:14:37.16 - 00:14:37.81 10
experiment
-
00:14:40.82 - 00:14:43.92 44
So we've looked at the adjoint method giving
-
00:14:43.92 - 00:14:44.34 6
you an
-
00:14:44.41 - 00:14:47.52 43
overview and some optimization examples and
-
00:14:47.52 - 00:14:48.76 19
closed with a brief
-
00:14:48.83 - 00:14:51.94 38
discussion on data informed turbulence
-
00:14:51.94 - 00:14:53.39 27
modeling thank you for your
-
00:14:53.39 - 00:14:54.01 9
attention