一种并行同步扰动随机近似的气动优化设计方法
【专利摘要】本发明公开了一种基于并行同步扰动随机近似的气动外形优化设计方法，针对翼型约束减阻、气动外形反设计等航空气动外形设计问题，充分利用多处理器并发处理能力，每个处理器所在进程均对设计变量进行独立的扰动，并实施同步扰动随机近似优化过程，更新设计变量以及对应目标，最后从所有进程中选出最好目标以及相应设计变量作为下一次设计的初值，如此循环迭代直到找到最优外形。本发明的积极效果是：相比于传统的串行优化设计方法，该技术发明易于实现，能够充分利用计算资源，极大提高气动优化设计效率。
【专利说明】一种并行同步扰动随机近似的气动优化设计方法
【技术领域】
[0001]本发明涉及并行的随机优化算法，具体是并行同步扰动随机近似优化方法，应用于气动外形优化设计，尤其是翼型减阻设计、叶栅设计及其反设计。
【背景技术】
[0002]气动外形优化设计是指在给定的某些约束条件下，寻找最优的气动外形，在可能的范围内极大化提高飞行器的气动性能和巡航效率。对于飞行器而已，机翼是主要升力面，也是影响气动性能的最主要部件，因此机翼的气动力是影响飞行器性能的主要因素。目前，针对各种飞行状态下，如何设计机翼使得其具有更优异的气动性能一直是研究的热点。传统的外形设计方法主要是根据经验和实验，每次设计一个新外形，都需要风洞实验或者计算来验证分析，进而修改外形。现代飞行器的外形以及绕流流场越来越复杂，对机翼的气动性能要求也越来越高，传统的方法无法胜任这些要求，这要求新的高效快速气动优化算法出现。
[0003]气动外形优化是一个极多控制参数，一般设计变量是几十个，甚至上百个，而且需要求解绕流流场，其控制方程是欧拉方程，或者更为复杂的NS方程，采用传统的优化方法，需要计算目标对于每一个设计变量的梯度，采用差分计算的话，需要多次计算扰动每个设计变量之后对于外形的流场，这个计算量是巨大的。目前，主要有基于伴随的气动优化方法，每次优化需要计算一次流场，然后计算一次伴随方程得到目标关于设计变量的导数，然后更新设计变量，该方法的缺点在于一旦更换目标约束，那么整个伴随方程都需要重新推导，不具有通用性；此外，粒子群法、遗传算法、模拟退火等现在优化算法也在气动优化中有了应用，同步扰 动随机近似是一种随机优化算法，适用于多参数的优化问题中，每一次迭代，对初始设计变量进行同步扰动，然后计算两次目标函数值获得一个近似梯度，从而更新设计变量，该方法易于实现，具有通用性，特别适用于复杂的多参数优化问题。
【发明内容】
[0004]为了解决同步扰动随机近似优化方法在气动优化设计中的效率问题，本发明提供了一种并行的同步扰动随机近似优化设计算法。
[0005]本发明的原理是通过并行的多进程实施同步扰动随机近似的寻优过程，同时保证每个进程的同步扰动量是独立随机的，从而使得每个进程上得到的当前最优解是互异的，充分搜寻可能的设计变量，最后比较每个进程上得到的目标值，将得到的最优解以及相应的设计变量通过消息传递到发送到其他进程，作为下一寻优的初始值，如此循环，直到满足迭代停止条件。
[0006]本发明对于流程图如图1所示，解决技术问题所采用的详细技术方案如下:
[0007]A、初始化:
[0008]给定优化设计的目标函数I，目标一般是气动参数，比如阻力系数等，根据参数化方法以及初始外形给出初始设计变量X(|，并计算初始流场U0,从而得到初始目标值q(!，启动并行环境，将\，U0, Q0通过消息传递方式广播到每个处理器上的进程；
[0009]B、在每个进程上均实施以下同步扰动随机近似优化方法:
[0010]以毫秒对应时间作为生成随机数种子生成独立地服从±1伯努利分布的随机变量Λ k，对设计变量进行同步扰动..= Ik + CkAkt Χζ = Xk - ckAk,根据扰动后的设计变量来计算近似梯度:
【权利要求】
1.一种基于并行同步扰动随机近似的气动优化设计方法，其特征在于，进行气动优化设计的步骤如下: A、初始化: 给定优化设计的目标函数I，根据参数化方法以及初始外形给出初始设计变量并计算初始流场U。，从而得到初始目标值％，将Xtl, U0, Q0广播到每个处理器上的进程； B、每个进程上都实施同步扰动随机近似优化方法: 对设计变量进行同步扰动:K = Xk + €kAk,X； =Xk- CkAk 计算近似梯度:
2.如权利要求1所述的方法，其特征在于，每个进程上的随机向量生成时采用的随机数种子必须是不同的。
3.如权利要求1所述的方法，其特征在于，小的正数ak，ck满足:
【文档编号】G06F17/50GK104036098SQ201410302688
【公开日】2014年9月10日 申请日期:2014年6月27日 优先权日:2014年6月27日
【发明者】刘铁钢, 王正, 余胜蛟, 何晓峰, 赵越 申请人:北京航空航天大学