CN109255171B - 一种自动判定数值模拟计算收敛的方法 - Google Patents
一种自动判定数值模拟计算收敛的方法 Download PDFInfo
- Publication number
- CN109255171B CN109255171B CN201810991384.6A CN201810991384A CN109255171B CN 109255171 B CN109255171 B CN 109255171B CN 201810991384 A CN201810991384 A CN 201810991384A CN 109255171 B CN109255171 B CN 109255171B
- Authority
- CN
- China
- Prior art keywords
- residual
- steps
- residual error
- observed
- iteration step
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Geometry (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种自动判定数值模拟计算收敛的方法,包括以下步骤:(1)参考值的确定;(2)偏差的确定;(3)变化区间和恒定区间的确定;(4)收敛性评价。
Description
技术领域
此本发明涉及一种计算机辅助工程(CAE)、计算流体力学(CFD)、计算结构力学(CSD)及数值仿真软件技术,具体涉及偏微分方程离散后成为代数方程组的迭代求解,是一种自动判断解是否收敛的方法。
背景技术
随着计算机硬件技术和数值分析技术的快速发展,以计算流体力学(Computational Fluid Dynamics,CFD)软件、计算结构力学(Computational StructureDynamics,CSD)软件为代表的计算机辅助工程软件及数值仿真技术正越来越广泛地应用于各种工程技术领域中,成为重要的产品研发与技术创新平台。例如,航空飞行器的空气动力学与结构力学设计、汽车热舒适性与结构动力学设计、发动机热力学与动力学设计、船舶水动力学设计等,都离不开大量基于CFD/CSD软件的仿真计算与数值优化工作。
数值模拟通过一定的离散技术,将代表物理规律的控制方程组经过线化处理转化为大规模代数方程组进行求解,工程实际问题离散形成的代数方程组对应的往往是大规模稀疏矩阵,一般采用迭代法进行求解。在物理场层面,也是通过迭代的方式进行,即先估计一个物理场的分布,然后将物理场代入到控制方程组中求得表征方程平衡程度的残差,并根据残差来改进物理场的分布,如此迭代下去直至离散方程的平衡程度达到满足实际问题的需求。以下我们主要讨论的是自动判断物理场收敛的方法,其基本思想也可以直接应用到矩阵的迭代求解。
这里面满足问题的要求是个定性的评价标准,实际操作中有两种定量操作:一种是当残差的下降指定的量级后认为物理问题得到收敛,一个是达到指定迭代次数后认为物理问题收敛。这两种方法在实际操作中都有一定的问题。指定残差下降的量级时可能发生两种情况:指定的残差下降量级过大导致数值模拟达到物理场收敛后一直在计算而满足不了残差要求,或者是指定的残差下降量级过小,从而用户得到的物理场实际上还没有收敛。指定迭代步数的方式类似,对于快速收敛的计算,固定迭代步数造成计算资源的浪费,而对于缓慢收敛的计算,固定迭代步数后依然没达到收敛。这些方式都造成了需要用户在数值模拟过程中不停的对求解器进行干涉,以取得一个物理场收敛的结果。
根据实际数值模拟过程中残差和观测量在迭代过程中的特点,这里提出了一种能够自动判定数值模拟计算收敛的方法。
发明内容
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供一种自动判定数值模拟计算收敛的方法。
为了解决上述技术问题,本发明公开了一种自动判定数值模拟计算收敛的方法,包括以下步骤。
(1)在数值计算的初始阶段,得到残差和观测量的初始值。初始阶段指的是流场进行适当的初始化之后的一个或若干个迭代步的状态。残差是当前物理场代表的解代入离散化方程后得到的数值,其表征当前物理场与离散方程精确解之间的差别。观测量是由当前物理场给出的单点值或者统计量,例如给定空间坐标的压力、固定固体结构上的受力、入口和出口的压差和流量差等。
(2)在迭代求解过程中,获取每个迭代步的残差和监测量数值。对残差取对数后进行均值分析,对监测量不同迭代步数的均值分析,记录均值和方差等信息。
(3)根据步骤(2)得到数据寻找对数残差和监测量是否出现了常数区域,如果残差和监测量都出现了显著的常数区域,则计算收敛,否则重复(2)。
(4)输出收敛时的残差下降量级,以及监测量的数值。
步骤(1)中,得到残差和观测量的初始值,包括以下步骤:
(1a)使用合适的物理场对求解区域进行初始化;
(1b)得到物理场对应的残差和观测量;
(1c)可以根据需要对物理场进行一步或几步迭代,并记录对应的残差和观测量;
(1d)将上述(1a)至(1c)过程得到残差和观测量的统计值作为残差和观测量的参考值。
在迭代求解过程中,获取每个迭代步的残差和监测量数值,包括以下步骤:
(2a)将每个迭代步的残差除以残差的参考值,得到每个迭代步的相对残差;
(2b)将每个迭代步的相当残差取对数,得到每个迭代步的对数残差;
(2c)对从初始迭代步到当前迭代步的对数残差进行均值分析,并记录其偏差;
(2d)对从初始迭代步到当前迭代步的观测量进行均值分析,并记录其偏差。
步骤(3)中以步骤(2c)和(2d)得到的均值和偏差为依据,判断是否出现常数区域,包括以下步骤:
(3a)对数残差数据能够使用一个两段线性的曲线进行拟合,其中第一段是随迭代步下降的直线,第二段是水平直线段;
(3b)水平直线段的迭代步数达到单调下降段迭代步数的一定比例,比如1/10,水平直线段的均方差与单调下降段的均方差同一个量级或者更小;
(3c) 观测量在舍弃初始一定数目的迭代步之后趋于常数,且与常数的偏差远小于参考值,比如偏差在参考值的千分之一或者更小;
(3d)如果满足(3a)至(3c)的条件,则数值计算收敛,且对数残差的初始值与水平直线段的值之差即残差收敛量级,而达到的观测值即数值计算的观测值;若不满足(3a)至(3c)的条件,则重复步骤(2a)至(2d)直至(3a)至(3c)得到满足。
有益效果:本发明无需用户指定事先难以确定的残差下降量级,省却了指定残差下降量级和迭代步数的传统方式中用户需要不停与求解软件交互的麻烦,使得整个计算直至收敛的过程能够自动完成。
附图说明
图1为本发明自动判断数值模拟结果收敛的流程图。
具体实施方式
如图1所示,本发明包括以下步骤。
(1)在数值计算的初始阶段,得到残差和观测量的初始值。初始阶段指的是流场进行适当的初始化之后的一个或若干个迭代步的状态。残差是当前物理场代表的解代入离散化方程后得到的数值,其表征当前物理场与离散方程精确解之间的差别。观测量是由当前物理场给出的单点值或者统计量,例如给定空间坐标的压力、固定固体结构上的受力、入口和出口的压差和流量差等。
(2)在迭代求解过程中,获取每个迭代步的残差和监测量数值。对残差取对数后进行均值分析,对监测量不同迭代步数的均值分析,记录均值和方差等信息。
(3)根据步骤(2)得到数据寻找对数残差和监测量是否出现了常数区域,如果残差和监测量都出现了显著的常数区域,则计算收敛,否则重复(2)。
(4)输出收敛时的残差下降量级,以及监测量的数值。
步骤(1)中,得到残差和观测量的初始值,包括以下步骤:
(1a)使用合适的物理场对求解区域进行初始化;
(1b)得到物理场对应的残差和观测量;
(1c)可以根据需要对物理场进行一步或几步迭代,并记录对应的残差和观测量;
(1d)将上述(1a)至(1c)过程得到残差和观测量的统计值作为残差和观测量的参考值。
在迭代求解过程中,获取每个迭代步的残差和监测量数值,包括以下步骤:
(2a)将每个迭代步的残差除以残差的参考值,得到每个迭代步的相对残差;
(2b)将每个迭代步的相当残差取对数,得到每个迭代步的对数残差;
(2c)对从初始迭代步到当前迭代步的对数残差进行均值分析,并记录其偏差;
(2d)对从初始迭代步到当前迭代步的观测量进行均值分析,并记录其偏差。
步骤(3)中以步骤(2c)和(2d)得到的均值和偏差为依据,判断是否出现常数区域,包括以下步骤:
(3a)对数残差数据能够使用一个两段线性的曲线进行拟合,其中第一段是随迭代步下降的直线,第二段是水平直线段;
(3b)水平直线段的迭代步数达到单调下降段迭代步数的一定比例,比如1/10,水平直线段的均方差与单调下降段的均方差同一个量级或者更小;
(3c) 观测量在舍弃初始一定数目的迭代步之后趋于常数,且与常数的偏差远小于参考值,比如偏差在参考值的千分之一或者更小;
(3d)如果满足(3a)至(3c)的条件,则数值计算收敛,且对数残差的初始值与水平直线段的值之差即残差收敛量级,而达到的观测值即数值计算的观测值;若不满足(3a)至(3c)的条件,则重复步骤(2a)至(2d)直至(3a)至(3c)得到满足。
本发明提供了一种自动判定数值模拟计算收敛的方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (3)
1.一种自动判定数值模拟计算收敛的方法,其特征在于,包括以下步骤:
(1)在数值计算的初始阶段,得到残差和观测量的初始值,初始阶段指的是流场进行适当的初始化之后的一个或若干个迭代步的状态,残差是当前物理场代表的解代入离散化方程后得到的数值,其表征当前物理场与离散方程精确解之间的差别,观测量是由当前物理场给出的单点值或者统计量;
(2)在迭代求解过程中,获取每个迭代步的残差和观测量数值,对残差取对数后进行均值分析,对观测量不同迭代步数的均值分析,记录均值和方差等信息;
(3)根据步骤(2)得到数据寻找对数残差和观测量是否出现了常数区域,如果残差和观测量都出现了显著的常数区域,则计算收敛,否则重复步骤(2);
(4)输出收敛时的残差下降量级,以及观测量的数值;
步骤(1)中,得到残差和观测量的初始值,包括以下步骤:
(1a)使用合适的物理场对求解区域进行初始化;
(1b)得到物理场对应的残差和观测量;
(1c)对物理场进行一步或几步迭代,并记录对应的残差和观测量;
(1d)将上述(1a)至(1c)过程得到残差和观测量的统计值作为残差和观测量的参考值。
2.如权利要求1所述的一种自动判定数值模拟计算收敛的方法,其特征在于,在迭代求解过程中,获取每个迭代步的残差和观测量数值,包括以下步骤:
(2a)将每个迭代步的残差除以残差的参考值,得到每个迭代步的相对残差;
(2b)将每个迭代步的相对残差取对数,得到每个迭代步的对数残差;
(2c)对从初始迭代步到当前迭代步的对数残差进行均值分析,并记录其偏差;
(2d)对从初始迭代步到当前迭代步的观测量进行均值分析,并记录其偏差。
3.如权利要求2所述的一种自动判定数值模拟计算收敛的方法,其特征在于,步骤(3)中以步骤(2c)和(2d)得到的均值和偏差为依据,判断是否出现常数区域,包括以下步骤:
(3a)对数残差数据能够使用一个两段线性的曲线进行拟合,其中第一段是随迭代步下降的直线,第二段是水平直线段;
(3b)水平直线段的迭代步数达到单调下降段迭代步数的1/10,水平直线段的均方差与单调下降段的均方差同一个量级或者更小;
(3c)观测量在舍弃初始一定数目的迭代步之后趋于常数,且与常数的偏差远小于参考值;
(3d)如果满足(3a)至(3c)的条件,则数值计算收敛,且对数残差的初始值与水平直线段的值之差为残差收敛量级,而达到的观测值为数值计算的观测值;若不满足(3a)至(3c)的条件,则重复步骤(2a)至(2d)直至(3a)至(3c)得到满足。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810991384.6A CN109255171B (zh) | 2018-08-29 | 2018-08-29 | 一种自动判定数值模拟计算收敛的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810991384.6A CN109255171B (zh) | 2018-08-29 | 2018-08-29 | 一种自动判定数值模拟计算收敛的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109255171A CN109255171A (zh) | 2019-01-22 |
CN109255171B true CN109255171B (zh) | 2023-09-05 |
Family
ID=65050303
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810991384.6A Active CN109255171B (zh) | 2018-08-29 | 2018-08-29 | 一种自动判定数值模拟计算收敛的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109255171B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111079235B (zh) * | 2019-12-11 | 2023-04-07 | 内蒙动力机械研究所 | 一种固体火箭发动机内流场仿真快速收敛方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006268422A (ja) * | 2005-03-24 | 2006-10-05 | Matsushita Electric Ind Co Ltd | 解析ナビゲーションシステム |
US8510091B1 (en) * | 2010-09-09 | 2013-08-13 | Sas Ip, Inc. | Domain decomposition formulations for simulating electromagnetic fields |
CN105512502A (zh) * | 2016-01-13 | 2016-04-20 | 重庆大学 | 一种基于残差归一化的权函数最小二乘状态估计方法 |
CN105552904A (zh) * | 2016-01-30 | 2016-05-04 | 清华大学 | 基于双线性化的多区域电网全分布式抗差状态估计方法 |
CN107016155A (zh) * | 2015-12-28 | 2017-08-04 | 达索系统西姆利亚公司 | 非线性pde和线性求解器的收敛估计 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1642222A2 (en) * | 2003-06-25 | 2006-04-05 | BAE Systems PLC | Design optimisation of computationally intensive design problems |
JP6803173B2 (ja) * | 2015-08-24 | 2020-12-23 | エスエーエス アイピー,インコーポレーテッドSAS IP, Inc. | 時間ドメイン分解過渡シミュレーションのためのプロセッサ実行システム及び方法 |
US10223331B2 (en) * | 2016-04-08 | 2019-03-05 | Goodrich Corporation | Warp models for registering multi-spectral imagery |
-
2018
- 2018-08-29 CN CN201810991384.6A patent/CN109255171B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006268422A (ja) * | 2005-03-24 | 2006-10-05 | Matsushita Electric Ind Co Ltd | 解析ナビゲーションシステム |
US8510091B1 (en) * | 2010-09-09 | 2013-08-13 | Sas Ip, Inc. | Domain decomposition formulations for simulating electromagnetic fields |
CN107016155A (zh) * | 2015-12-28 | 2017-08-04 | 达索系统西姆利亚公司 | 非线性pde和线性求解器的收敛估计 |
CN105512502A (zh) * | 2016-01-13 | 2016-04-20 | 重庆大学 | 一种基于残差归一化的权函数最小二乘状态估计方法 |
CN105552904A (zh) * | 2016-01-30 | 2016-05-04 | 清华大学 | 基于双线性化的多区域电网全分布式抗差状态估计方法 |
Non-Patent Citations (1)
Title |
---|
安恩科等.电站锅炉燃烧过程数值模拟的收敛性分析.同济大学学报(自然科学版).2005,第33卷(第1期),第121-125页. * |
Also Published As
Publication number | Publication date |
---|---|
CN109255171A (zh) | 2019-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10489533B2 (en) | Techniques for warm starting finite element analyses with deep neural networks | |
CN112231847A (zh) | 转捩位置确定方法、装置、电子设备及存储介质 | |
CN113360983B (zh) | 一种边坡可靠度分析与风险评估方法 | |
CN109255171B (zh) | 一种自动判定数值模拟计算收敛的方法 | |
Miro et al. | Reliability analysis of an axial compressor based on one-dimensional flow modeling and survival signature | |
CN104834795B (zh) | 包带连接结构接触摩擦非线性特性模拟方法及系统 | |
CN107247831B (zh) | 一种基于迭代逐维法的多学科区间不确定性分析方法 | |
Du et al. | Aircraft design optimization with uncertainty based on fuzzy clustering analysis | |
KR102266279B1 (ko) | 비정상상태를 구현하기 위한 차수 감축 모델 구축 방법 | |
CN109325279B (zh) | 一种离散的飞机气动载荷参数化的方法 | |
CN109598059B (zh) | 一种基于代理模型的热防护系统优化设计方法及设计系统 | |
CN109635346B (zh) | 机械连接结构的可靠性分析方法 | |
CN113343523B (zh) | 一种基于空间网格的有限元数值模拟分析方法 | |
US11983467B2 (en) | Rapid aero modeling for computational experiments (RAM-C) | |
Benay | A global method of data reduction applied to seven-hole probes | |
CN114417653A (zh) | 基于不确定性的卫星组件电路板热可靠性分析方法 | |
EP3690670A1 (en) | Data processing system | |
CN114035536A (zh) | 一种基于蒙特卡洛法的飞行控制系统鲁棒性评估方法 | |
Clarich et al. | Reliability-based design optimization applying polynomial chaos expansion: theory and Applications | |
Ugryumova et al. | Synthesis of robust optimal control program for axial flow compressor turning guide vanes | |
Prots et al. | Robust Sensitivity Analysis of Complex Simulation Models Subject to Noise | |
CN110543724A (zh) | 面向总体设计的卫星结构性能预测方法 | |
Favaregh et al. | Space Launch System Aerodynamic Database Uncertainty Quantification Methodologies | |
Wright | Validation Process for LEWICE by Use of a Naviér-Stokes Solver | |
Favaregh et al. | Quantification of the Uncertainties for the Space Launch System Liftoff/Transition and Ascent Databases |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20220412 Address after: 518000 2301, building D1, Nanshan Zhiyuan, No. 1001, Xueyuan Avenue, Changyuan community, Taoyuan Street, Nanshan District, Shenzhen, Guangdong Applicant after: Shenzhen Shifeng Technology Co.,Ltd. Address before: 518000 1902, floor 19, building B1, Nanshan Zhiyuan, No. 1001, Xueyuan Avenue, Nanshan District, Shenzhen, Guangdong Applicant before: SHENZHEN QINGFENGXI TECHNOLOGY Co.,Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |