[发明专利]基于量子算法的计算流体动力学模拟方法、装置及设备有效
| 申请号: | 202010771592.2 | 申请日: | 2020-08-04 | 
| 公开(公告)号: | CN114091363B | 公开(公告)日: | 2023-08-08 | 
| 发明(设计)人: | 李叶;安宁波 | 申请(专利权)人: | 合肥本源量子计算科技有限责任公司 | 
| 主分类号: | G06F30/28 | 分类号: | G06F30/28;G06F30/23;G06F30/15;G06N10/60;G06N10/20;G06F111/10;G06F113/08;G06F119/14 | 
| 代理公司: | 北京柏杉松知识产权代理事务所(普通合伙) 11413 | 代理人: | 马敬;丁芸 | 
| 地址: | 230088 安徽省合肥市高新*** | 国省代码: | 安徽;34 | 
| 权利要求书: | 查看更多 | 说明书: | 查看更多 | 
| 摘要: | |||
| 搜索关键词: | 基于 量子 算法 计算 流体动力学 模拟 方法 装置 设备 | ||
1.一种基于量子算法的计算流体动力学模拟方法,其特征在于,所述方法包括:
在利用有限体积法解析计算流体动力学的过程中,针对流体运动的离散化数值网格中各网格单元,构建表示所述网格单元坐标信息的第一量子线路、表示所述网格单元的状态参数的第二量子线路,其中,所述网格单元的状态参数存储在量子随机存取存储器中,所述量子随机存取存储器可操作处于量子叠加态的地址和数据;
基于所述第一量子线路、所述第二量子线路、所述量子随机存取存储器构建表示所述网格单元流体状态变化的线性系统方程AnΔUn=n的参数An和bn的第三量子线路;其中,n表示记录流体状态的第n个时刻,An表示第n个时刻的系数矩阵,bn表示第n个时刻的残差量,ΔUn表示第n+1个时刻网格单元的流体状态Un+1与表示第n个时刻网格单元的流体状态Un之差;
针对所有的网格单元,基于所述第三量子线路求解所述网格单元的线性系统方程,当所述网格单元的流体状态趋于稳定时,获得所述网格单元线性系统方程的状态参数表示的流体状态作为所述网格单元的目标状态。
2.根据权利要求1所述的方法,其特征在于,所述第一量子线路包括第一Oracle及第二Oracle;
所述第一Oracle用于提取指定网格单元的相邻网格数,所述第二Oracle用于提取网格单元的坐标信息。
3.根据权利要求1所述的方法,其特征在于,所述第二量子线路,包括:
通过量子态的振幅提取所述网格单元的状态参数的第二模拟态量子线路,和通过量子态的本征态提取所述网格单元的状态参数的第二数字态量子线路。
4.根据权利要求1所述的方法,其特征在于,所述量子随机存取存储器内存储有的地址和数据的预设映射关系,所述地址用第一叠加态编码,所述第一叠加态中的每个本征态对应一个地址信息,所述数据用第二叠加态编码,所述第二叠加态中的每个本征态对应一个数据信息。
5.根据权利要求2所述的方法,其特征在于,
其中:Un={Ui},i=0,1,2,...,N-1,Ui表示网格单元i的流体状态;且表示An为与Un,相关的函数,表示bn为与Un,相关的函数,N为网格单元的总数量,为所述网格单元的坐标信息;
所述线性系统方程的参数包括当前时刻的系数矩阵An和当前时刻的残差量bn;
所述第三量子线路包括对应所述系数矩阵An的Oracle OA,和对应所述残差量bn的Oracle Ob;
所述Oracle OA用于提取当前n时刻系数矩阵An的元素Ai′j′,An={Ai′j′},i′=4i+i1,j′=4j+j1,i,j=0,1,2,...,N-1;i1,j1是0到3的整数,对网格单元i的第j个相邻网格单元,所述Oracle OA的作用为OA| i|j=|i|j|Ai′j′;
所述Oracle Ob用于提取当前n时刻的残差量bn的元素bi,其中,bi表示网格单元i的残差量,bn={bi},i=0,1,2,...,N-1;对网格单元i,所述Oracle Ob的作用为Ob|0=∑ibi|i。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于合肥本源量子计算科技有限责任公司,未经合肥本源量子计算科技有限责任公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010771592.2/1.html,转载请声明来源钻瓜专利网。





