CN109801348A - 一种适用于三维流体模拟的三阶高精度对流插值算法 - Google Patents
一种适用于三维流体模拟的三阶高精度对流插值算法 Download PDFInfo
- Publication number
- CN109801348A CN109801348A CN201910036395.3A CN201910036395A CN109801348A CN 109801348 A CN109801348 A CN 109801348A CN 201910036395 A CN201910036395 A CN 201910036395A CN 109801348 A CN109801348 A CN 109801348A
- Authority
- CN
- China
- Prior art keywords
- point
- interpolation
- enabled
- dimensional
- precision
- 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.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种适用于三维流体模拟的三阶高精度对流插值算法。基于约束插值剖面法(CIP)，在使用半拉格朗日法求解三维流体运动控制方程中对流项时，针对回退点处物理量的计算发明了一种高精度插值方法，提升了对流精度。此外，为减少内存消耗，只存储物理场的值及其一阶导数作为计算变量，高阶导数则在保证计算精度不受损的前提下基于推导的泰勒展开式近似计算。本发明能够在时间和内存消耗较少的前提下保持三阶高精度，且具有紧模板特性。对比现有方法，该方法在视觉质量、速度、内存消耗等方面都有明显的改进，能够有效提升流体模拟的对流精度和速度。
Description
技术领域
本发明涉及计算机动画技术领域，是一种适用于三维流体模拟的三阶高精度对流插值算法，能有效降低对流数值耗散，减少内存开销，增强流体动画的真实感。
背景技术
在计算机图形学和虚拟现实领域，基于物理的计算机动画技术始终是科研工作者的研究重点，而基于物理的流体动态模拟技术则是其中的一个热点问题。流体模拟技术在多个领域也被广泛应用，工程领域如航天、航空、航海等的发展也离不开流体技术的支撑。
在流体模拟中，数值耗散增加了流体的粘度，使其比预期的更粘稠，并涂抹细节以致视觉质量受损。导致数值耗散的因素中，对流的精度对流体模拟的视觉质量有很大的影响，相关学者为开发精确的对流求解器进行了各种尝试。其中，高阶插值格式广泛应用于计算流体力学(CFD)领域，包括本质非振荡(ENO)、加权ENO(WENO)等，然而由于计算量太大不适用于图形学领域，而且这些插值算法在宽模板上进行，不适用于在非均匀网格上进行仿真。线性方法(Linear)、前后向误差补偿矫正方法(BFECC)、MacCormack法等对流方法广泛用于图形学中，但计算精度有限，最高只可达二阶精度。
约束插值剖面(CIP)方法只在单个网格单元构造插值函数，具备紧模板特性，而且具有三阶精度，但扩展CIP求解高维对流方程是一项困难的任务，需要大量的计算时间和内存消耗。另外，现有技术中基于CIP的多维对流求解器，如单调CIP(MCIP)、非分裂型CIP(USCIP)等，只能部分地解决这些问题，有些甚至以损失数值精度或造成不稳定为代价。如何发展一种高效、高精度的高维CIP方案对于提升对流方法的计算效率和精度颇有意义，但仍然是一个具有挑战性的问题。
发明内容
针对现有技术的不足，本发明提供一种适用于二维流体模拟的三阶高精度对流插值算法：三维基于泰勒展开的约束插值剖面法(3d Taylor Expansion based ConstrainedInterpolation Profile，3dTECIP)。现有基于CIP的方法应用到高维时，通常存在计算开销高、内存占用大、精度降低或者不稳定等问题。本发明方法只将物理量及其一阶导数作为未知变量，高阶导数并不作为变量保存，可以减少大量内存开销。其次，本发明使用局部泰勒展开按需计算高阶导数，计算速度较快且保持三阶精度，能够更有效保持丰富的流体细节。最后，方法并未破坏原CIP方法的紧模板特性，除规则网格外同样适用于非均匀网格。
本发明的技术方案为：一种适用于三维流体模拟的三阶高精度对流插值算法，本发明主要用于提高求解流体运动控制方程中对流项的精度。通过利用半拉格朗日法求解流体运动控制方程的对流项时，针对回退点处物理量的计算采用了高精度插值，降低了数值耗散，同时为减少内存消耗，只存储物理场的值及其一阶导数作为计算变量。
对于三维对流方程，设φ为对流物理量，则在每个网格点i上存储4个变量，即该物理量及其一阶偏导：φi、和回退点上相应值通过高精度插值求得。具体包括以下步骤：
S1)针对回退点P所处的三维网格单元，设网格单元边长为h，八个顶点分别为A、B、C、D、E、F、G、H，网格的所有边均与三维空间的坐标轴平行，网格面ABCD和面EFGH均平行于XY平面且二者沿Z轴正向的间距为h，网格面AEHD和面BFGC平行于YZ平面且且二者沿X轴正向的间距为h，网格面ABFE和面DCGH平行于XZ平面且二者沿Y轴正向的间距为h。设点M和N分别为回退点P在面ABCD和面EFGH上的投影点；
S2)、计算八个顶点处φ的二阶偏导的值，具体计算如下：
基于以下三阶精度泰勒展开式，
令y＝z，有
对于网格点A，针对以上公式，令Δx＝Δz＝h，得到：
对于网格点B，针对以上公式，令Δx＝-h，Δz＝h，得到：
对于网格点C，针对以上公式，令Δx＝-h，Δz＝h，得到：
对于网格点D，针对以上公式，令Δx＝Δz＝h，得到：
对于网格点E，针对以上公式，令Δx＝h，Δz＝-h，得到：
对于网格点F，针对以上公式，令Δx＝-h，Δz＝-h，得到：
对于网格点G，针对以上公式，令Δx＝-h，Δz＝-h，得到：
对于网格点H，针对以上公式，令Δx＝h，Δz＝-h，得到：
S3)、计算八个顶点处φ的二阶偏导的值，具体计算如下：
基于以下三阶精度泰勒展开式，
令x＝y,y＝z，有
对于网格点A，针对以上公式，令Δy＝Δz＝h，得到：
对于网格点B，针对以上公式，令Δy＝Δz＝h，得到：
对于网格点C，针对以上公式，令Δy＝-h，Δz＝h，得到：
对于网格点D，针对以上公式，令Δy＝-h，Δz＝h，得到：
对于网格点E，针对以上公式，令Δy＝h，Δz＝-h，得到：
对于网格点F，针对以上公式，令Δy＝h，Δz＝-h，得到：
对于网格点G，针对以上公式，令Δy＝-h，Δz＝-h，得到：
对于网格点H，针对以上公式，令Δy＝-h，Δz＝-h，得到：
S4)、通过二维基于泰勒展开的约束插值剖面算法计算投影点M和N处φ值及其一阶、二阶偏导：φi、其中i∈{M,N}；
S5)、基于一维约束插值剖面法(CIP)插值算法，通过步骤S4)求得的投影点M和N处的φM、φN、计算回退插值点P处的φP、通过步骤S4)求得的投影点M和N处的计算点P处通过步骤S4)求得的投影点M和N处的计算点P处
进一步的，步骤S4)中所述的二维基于泰勒展开的约束插值剖面算法，具体如下：
设φ为插值物理量，则在每个网格点i上存储3个变量，即该物理量及其一阶偏导：φi、插值点上相应值通过高精度插值求得。包括以下步骤：
S401)、针对插值点P′所处的二维网格单元，设网格单元边长为h′，四个顶点分别为A′、B′、C′、D′，网格的所有边均与二维维空间的坐标轴平行，边A′B′和边C′D′平行于X轴且二者沿Y轴正向的间距为h′，边A′D′和边B′C′平行于Y轴且二者沿X轴正向的间距为h′。设点E′和G′分别为插值目标点P′在边A′B′和A′D′上的投影点；
S402)、计算点A′、B′、D′处φ的二阶偏导的值。具体计算如下：
基于以下三阶精度泰勒展开式，
令Δx＝Δy＝h′，计算网格点A′处的二阶偏导值，即：
同理，令Δx＝-h′，Δy＝h′，计算网格点B′处的二阶偏导值，即：
同理，令Δx＝h′，Δy＝-h′，计算网格点D′处的二阶偏导值，即：
S403)、基于一维CIP插值算法，通过点A′和点B′处的φA′、φB′、计算点E′的φE′、通过点A′和点B′处的计算点E′的
通过网格点A′和点D′的φA′、φD′、计算网格点G′的φG′、通过网格点A′和点D′的计算点G′的
S404)、根据以下公式计算得到P′点的φP′、
其中，ξ′和η′分别为A′E′和A′G′的长度。
进一步的，步骤S4)中，针对面ABCD，由于点M为回退插值目标点P在平面ABCD上的投影，设定M点在AB边上的投影为O，在AD边上的投影为K，ξ和η分别为AO和AK的长度。基于上述二维基于泰勒展开的约束插值剖面算法，以A、B、C、D四点处的φ、值作为输入，将点M作为插值目标点，可得到M处的对应值：φM、具体如下：
同理，以A、B、C、D四点处的值为输入，将点M作为插值目标点，可得到M处的对应值：计算公式如下：
进一步的，步骤S4)中，针对面EFGH，由于点N为回退插值目标点P在平面EFGH上的投影，设定N点在EF边上的投影为S，在EH边上的投影为J，则ξ和η分别为ES和EJ的长度。基于前述二维基于泰勒展开的约束插值剖面算法，以E、F、G、H四点处的φ、 值作为输入，将点N作为插值目标点，可得到N处的对应值：φN、计算公式如下：
同理，以E、F、G、H四点处的值为输入，将点N作为插值目标点，可得到N处的对应值：计算公式如下：
进一步的，步骤S5)中，基于一维CIP插值算法，通过投影点M和N处的φM、φN、 计算回退插值点P处的φP、其计算式如下：
φP＝C3(zP-zM)3+C2(zP-zM)2+C1(zP-zM)+C0；
其中，C0＝φM，zP、zM分别为点P和点M的Z轴坐标，h为网格边长。
进一步的，步骤S5)中，基于一维CIP插值算法，通过投影点M和N处的 计算回退插值点P处的其计算式如下：
其中， zP、zM分别为点P和点M的Z轴坐标，h为网格边长。
进一步的，步骤S5)中，基于一维CIP插值算法，通过投影点M和N处的 计算回退插值点P处的其计算公式如下：
其中， zP、zM分别为点P和点M的Z轴方向坐标，h为网格边长。
进一步的，步骤S403)中，基于一维CIP插值算法，通过点A′和点B′处的φA′、φB′、 得到点E′的φE′、其计算公式如下：
φE′＝C3(xE′-xA′)3+C2(xE′-xA′)2+C1(xE′-xA′)+C0；
其中，C0＝φA′，xE′表示E′点横坐标，xA′表示A′点横坐标，h′为网格单元边长。
进一步的，步骤S403)中，基于一维CIP插值算法，通过点A′和点B′处的 得到点E′的其计算公式如下：
其中， xE′表示E′点横坐标，xA′表示A′点横坐标，h′为网格单元边长。
进一步的，步骤S403)中，基于一维CIP插值算法，通过网格点A′和点D′的φA′、φD′、得到网格点G′的φG′、其计算公式如下：
φG′＝C3(yG′-yA′)3+C2(yG′-yA′)2+C1(yG′-yA′)+C0；
其中，C0＝φA′，yG′表示G′点纵坐标，yA′表示A′点纵坐标，h′为网格单元边长。
进一步的，步骤S403)中，基于一维CIP插值算法，通过网格点A′和点D′的 得到点G′的其计算公式如下：
其中， yG′表示G′点纵坐标，yA′表示A点纵坐标，h′为网格单元边长。
本发明的有益效果为：
本发明能够在时间和内存消耗较少的前提下保持三阶高精度，且具有紧模板特性。对比现有方法，该方法在视觉质量、速度和内存消耗等方面都有明显的改进，能够有效提升流体模拟的对流精度和速度。
附图说明
图1为本发明方法用于三维流体模拟时的插值示意图；
图2为本发明方法中二维基于泰勒展开的约束插值剖面算法的插值示意图；
图3为本发明方法应用于三维流体模拟时的计算流程图；
图4为本发明与相关方法用于三维流体模拟的渲染结果对比，图中(a)、(b)、(c)、(d)、分别为采用Linear算法、BFECC算法、MCIP和USCIP算法的模拟效果图，(e)为采用本发明方法(3dTECIP)的模拟效果图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明：
提供一种适用于三维流体模拟的三阶高精度对流插值算法，用于计算流体控制方程的对流项。通过利用半拉格朗日法求解流体运动控制方程的对流项时，针对回退点处物理量的计算采用高精度插值，能够大大降低对流数值耗散，增强流体模拟的真实感。对于三维对流方程，设φ为对流物理量，则在每个网格点i上存储4个变量，即该物理量及其一阶偏导：φi、和回退点上相应值通过高精度插值求得。具体包括以下步骤：
S1)、针对泰勒展开式：
当Δx＝0，
当Δy＝0，
综合上述5式，并忽略三阶以上的余项，得到三阶精度的泰勒展开式：
整理后得到：
S2)、如图1所示，对于回退插值点P所处的三维网格单元，设网格单元边长为h，点M和N分别为点P在网格面ABCD和EFGH上的投影点。P点处的值φP、和由网格单元的八个网格点A、B、C、D、E、F、G、H的相应值通过发明的高精度插值求得；
S3)、基于步骤S1)得到的泰勒展开式，令y＝z，有
基于上式计算八个网格点处φ的二阶偏导的值，具体计算如下：
对于网格点A，针对以上公式，令Δx＝Δz＝h，得到：
对于网格点B，针对以上公式，令Δx＝-h，Δz＝h，得到：
对于网格点C，针对以上公式，令Δx＝-h，Δz＝h，得到：
对于网格点D，针对以上公式，令Δx＝Δz＝h，得到：
对于网格点E，针对以上公式，令Δx＝h，Δz＝-h，得到：
对于网格点F，针对以上公式，令Δx＝-h，Δz＝-h，得到：
对于网格点G，针对以上公式，令Δx＝-h，Δz＝-h，得到：
对于网格点H，针对以上公式，令Δx＝h，Δz＝-h，得到：
S4)、基于步骤S1)得到的泰勒展开式，令x＝y,y＝z，有
基于上式计算八个网格点处φ的二阶偏导的值，具体计算如下：对于网格点A，针对以上公式，令Δy＝Δz＝h，得到：
对于网格点B，针对以上公式，令Δy＝Δz＝h，得到：
对于网格点C，针对以上公式，令Δy＝-h，Δz＝h，得到：
对于网格点D，针对以上公式，令Δy＝-h，Δz＝h，得到：
对于网格点E，针对以上公式，令Δy＝h，Δz＝-h，得到：
对于网格点F，针对以上公式，令Δy＝h，Δz＝-h，得到：
对于网格点G，针对以上公式，令Δy＝-h，Δz＝-h，得到：
对于网格点H，针对以上公式，令Δy＝-h，Δz＝-h，得到：
S5)、通过二维基于泰勒展开的约束插值剖面算法计算投影点M和N处φ值及其一阶、二阶偏导：φi、其中i∈{M,N}；具体计算如下：
S501)、二维基于泰勒展开的约束插值剖面算法的具体实施步骤如下：
S50101)、设φ为插值物理量，则在每个网格点i上存储有3个变量，即该物理量及其一阶偏导：φi、如图2所示，针对插值点P′所处的二维网格单元A′B′C′D′，设网格单元边长为h′，点E′和G′分别为点P′在边A′B′和A′D′上的投影点，A′E′和A′G′的长度分别为ξ′和η′。P′点位置的值φP′、由网格单元的四个顶点A′、B′、C′、D′的相应值通过以下步骤插值求得；
S50102)、计算点A′、B′、D′处φ的二阶偏导的值。具体计算如下：
基于步骤S1)得到的泰勒展开式，
令Δx＝Δy＝h′，计算网格点A′处的二阶偏导值，即：
同理，令Δx＝-h′，Δy＝h′，计算网格点B′处的二阶偏导值，即：
同理，令Δx＝h′，Δy＝-h′，计算网格点D′处的二阶偏导值，即：
S50103)、基于一维CIP插值算法，通过点A′和点B′处的φA′、φB′、计算点E′的φE′、设xE′、xA′分别表示E′点、A′点的横坐标，其计算公式如下：
φE′＝C3(xE′-xA′)3+C2(xE′-xA′)2+C1(xE′-xA′)+C0；
其中，C0＝φA′，
通过点A′和点B′处的计算点E′的其计算公式如下：
其中，
S50104)、基于一维CIP插值算法，通过网格点A′和点D′的φA′、φD′、计算点G′的φG′、设yG′、yA′分别表示G′点、A′点的纵坐标，其计算公式如下：
φG′＝C3(yG′-yA′)3+C2(yG′-yA′)2+C1(yG′-yA′)+C0；
其中，C0＝φA′，
通过网格点A′和点D′的计算点G′的其计算公式如下：
其中，
S50105)、计算点P′处的φ，具体计算如下：
基于步骤S1)得到的泰勒展开式，
令Δx＝ξ′，Δy＝η′，可得点A′的值：
令Δx＝-ξ′，Δy＝η′，可得点E′的值：
令Δx＝ξ′，Δy＝-η′，可得点G′的值：
重新整理上面三式，得到以下公式以计算点P′处的φ，
S502)、如图1，针对网格面ABCD，设M点在AB边上的投影为O，在AD边上的投影为K，ξ和η分别为AO和AK的长度。
基于S501)所述二维基于泰勒展开的约束插值剖面算法，以A、B、C、D四点处的φ、值作为算法输入，将点M作为插值目标点，可得到M处的对应值：φM、具体如下：
基于S501)所述二维基于泰勒展开的约束插值剖面算法，以A、B、C、D四点处的 值作为算法输入，将点M作为插值目标点，可得到M处的对应值： 具体如下：
S503)、如图1针对网格面EFGH，N点在EF边上的投影为S，在EH边上的投影为J，则ES和EJ的长度分别为ξ和η。
基于S501)所述二维基于泰勒展开的约束插值剖面算法，以E、F、G、H四点处的φ、值作为输入，将点N作为插值目标点，可得到N处的对应值：φN、具体如下：
基于S501)所述二维基于泰勒展开的约束插值剖面算法，以E、F、G、H四点处的 值为算法输入，将点N作为插值目标点，可得到N处的对应值：具体如下：
S6)、基于一维CIP插值算法，通过投影点M和N处的φM、φN、计算回退插值点P处的φP、其计算式如下：
φP＝C3(zP-zM)3+C2(zP-zM)2+C1(zP-zM)+C0；
其中，C0＝φM，zP、zM分别为点P和点M的Z轴坐标，h为网格边长。
S7)、基于一维CIP插值算法，通过投影点M和N处的计算回退插值点P处的其计算式如下：
其中， zP、zM分别为点P和点M的Z轴坐标，h为网格边长。
S8)、基于一维CIP插值算法，通过投影点M和N处的计算回退插值点P处的其计算公式如下：
其中， zP、zM分别为点P和点M的Z轴方向坐标，h为网格边长。
图3为三维流体模拟的流程图，其中本发明提供的方法主要用于计算流体控制方程的对流项。图4为二维流体模拟结果渲染图对比，图4(e)为本发明方法用于三维流体模拟的结果渲染图，图4(a)、图4(b)、图4(c)、图4(d)分别为采用低阶精度的Linear和BFECC算法、以及MCIP和USCIP算法的结果渲染图，从图中可以看出，基于CIP的三种方法模拟效果优于低阶精度的Linear和BFECC算法。由于前述的MCIP和USCIP方法的缺陷，本发明方法在三维流体模拟中取得的效果优于MCIP和USCIP，模拟的烟雾细节更丰富。
上述实施例和说明书中描述的只是说明本发明的原理和最佳实施例，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。
Claims (8)
1.一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：通过利用半拉格朗日法求解流体运动控制方程的对流项时，针对回退点处物理量的计算采用了高精度插值，降低了数值耗散，同时为减少内存消耗，只存储物理场的值及其一阶导数作为计算变量，高阶导数则在保证计算精度不受损的前提下基于推导的泰勒展开式近似计算；
对于三维对流方程，设φ为对流物理量，则在每个网格点i上存储4个变量，即该物理量及其一阶偏导：φi、和回退点上相应值通过高精度插值求得；具体包括以下步骤：
S1)针对回退点P所处的三维网格单元，设网格单元边长为h，八个网格点分别为A、B、C、D、E、F、G、H，网格单元的所有边均与三维空间的坐标轴平行，网格面ABCD和EFGH均平行于XY平面且二者沿Z轴正向的间距为h，网格面AEHD和面BFGC平行于YZ平面且且二者沿X轴正向的间距为h，网格面ABFE和面DCGH平行于XZ平面且二者沿Y轴正向的间距为h；设点M和N分别为回退点P在面ABCD和面EFGH上的投影点；
S2)、计算八个顶点处φ的二阶偏导的值，具体计算如下：
基于以下三阶精度泰勒展开式，
令y＝z，有
对于网格点A，针对以上公式，令Δx＝Δz＝h，得到：
对于网格点B，针对以上公式，令Δx＝-h，Δz＝h，得到：
对于网格点C，针对以上公式，令Δx＝-h，Δz＝h，得到：
对于网格点D，针对以上公式，令Δx＝Δz＝h，得到：
对于网格点E，针对以上公式，令Δx＝h，Δz＝-h，得到：
对于网格点F，针对以上公式，令Δx＝-h，Δz＝-h，得到：
对于网格点G，针对以上公式，令Δx＝-h，Δz＝-h，得到：
对于网格点H，针对以上公式，令Δx＝h，Δz＝-h，得到：
S3)、计算八个顶点处φ的二阶偏导的值，具体计算如下：
基于以下三阶精度泰勒展开式，
令x＝y,y＝z，有
对于网格点A，针对以上公式，令Δy＝Δz＝h，得到：
对于网格点B，针对以上公式，令Δy＝Δz＝h，得到：
对于网格点C，针对以上公式，令Δy＝-h，Δz＝h，得到：
对于网格点D，针对以上公式，令Δy＝-h，Δz＝h，得到：
对于网格点E，针对以上公式，令Δy＝h，Δz＝-h，得到：
对于网格点F，针对以上公式，令Δy＝h，Δz＝-h，得到：
对于网格点G，针对以上公式，令Δy＝-h，Δz＝-h，得到：
对于网格点H，针对以上公式，令Δy＝-h，Δz＝-h，得到：
S4)、通过二维基于泰勒展开的约束插值剖面算法计算投影点M和N处的对流物理量φ值及其一阶、二阶偏导：φi、其中i∈{M,N}；
S5)、基于一维约束插值剖面法(CIP)插值算法，通过步骤S4)求得的投影点M和N处的φM、φN、计算回退插值点P处的φP、
通过步骤S4)求得的投影点M和N处的计算回退插值点P处
通过步骤S4)求得的投影点M和N处的计算回退插值点P处
2.根据权利要求1所述的一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：步骤S4)中，所述的二维基于泰勒展开的约束插值剖面算法，具体如下：
设φ为插值物理量，则在每个网格点i上存储3个变量，即该物理量及其一阶偏导：φi、插值点上相应值通过高精度插值求得；包括以下步骤：
S401)、针对插值点P′所处的二维网格单元，设网格单元边长为h′，四个顶点分别为A′、B′、C′、D′，网格单元的所有边均与二维维空间的坐标轴平行，边A′B′和边C′D′平行于X轴且二者沿Y轴正向的间距为h′，边A′D′和边B′C′平行于Y轴且二者沿X轴正向的间距为h′；设点E′和G′分别为插值目标点P′在边A′B′和A′D′上的投影点；
S402)、计算点A′、B′、D′处φ的二阶偏导的值。具体计算如下：
基于以下三阶精度泰勒展开式，
令Δx＝Δy＝h′，计算网格点A′处的二阶偏导值，即：
同理，令Δx＝-h′，Δy＝h′，计算网格点B′处的二阶偏导值，即：
同理，令Δx＝h′，Δy＝-h′，计算网格点D′处的二阶偏导值，即：
S403)、基于一维CIP插值算法，通过点A′和点B′处的φA′、φB′、计算点E′的φE′、通过点A′和点B′处的计算点E′的
通过网格点A′和点D′处的φA′、φD′、计算网格点G′的φG′、通过网格点A′和点D′的计算点G′的
S404)、根据以下公式计算得到P′点的φP′、
其中，ξ′和η′分别为A′E′和A′G′的长度。
3.根据权利要求2所述的一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：步骤S4)中，针对面ABCD，由于点M为回退插值目标点P在平面ABCD上的投影，设定M点在AB边上的投影为O，在AD边上的投影为K，ξ和η分别为AO和AK的长度；
基于上述二维基于泰勒展开的约束插值剖面算法，以A、B、C、D四点处的φ、值作为输入，将点M作为插值目标点，可得到M处的对应值：φM、具体如下：
同理，以A、B、C、D四点处的值为输入，将点M作为插值目标点，可得到M处的对应值：具体如下：
4.根据权利要求3所述的一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：步骤S4)中，针对面EFGH，由于点N为回退插值目标点P在平面EFGH上的投影，设定N点在EF边上的投影为S，在EH边上的投影为J，则ξ和η分别为ES和EJ的长度；
基于权上述二维基于泰勒展开的约束插值剖面算法，以E、F、G、H四点处的φ、 值作为输入，将点N作为插值目标点，可得到N处的对应值：φN、具体如下：
同理，以E、F、G、H四点处的值为输入，将点N作为插值目标点，可得到N处的对应值：具体如下：
5.根据权利要求4所述的一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：步骤S5)中，基于一维CIP插值算法，通过投影点M和N处的φM、φN、计算回退插值点P处的φP、其计算式如下：
φP＝C3(zP-zM)3+C2(zP-zM)2+C1(zP-zM)+C0；
其中，C0＝φM，zP、zM分别为点P和点M的Z轴坐标，h为网格边长。
6.根据权利要求4所述的一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：步骤S5)中，基于一维CIP插值算法，通过投影点M和N处的计算回退插值点P处的其计算式如下：
其中， zP、zM分别为点P和点M的Z轴坐标，h为网格边长。
7.根据权利要求4所述的一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：步骤S5)中，基于一维CIP插值算法，通过投影点M和N处的计算回退插值点P处的其计算公式如下：
其中， zP、zM分别为点P和点M的Z轴方向坐标，h为网格边长。
8.根据权利要求2所述的一种适用于三维流体模拟的三阶高精度对流插值算法，其特征在于：步骤S403)中，基于一维CIP插值算法，通过点A′和点B′处的φA′、φB′、得到点E′的φE′、其计算公式如下：
φE′＝C3(xE′-xA′)3+C2(xE′-xA′)2+C1(xE′-xA′)+C0；
其中，C0＝φA′，xE′表示E′点横坐标，xA′表示A′点横坐标，h′为网格单元边长；
基于一维CIP插值算法，通过点A′和点B′处的得到点E′的 其计算公式如下：
其中， xE′表示E′点横坐标，xA′表示A′点横坐标，h′为网格单元边长；
基于一维CIP插值算法，通过网格点A′和点D′的φA′、φD′、得到网格点G′的φG′、其计算公式如下：
φG′＝C3(yG′-yA′)3+C2(yG′-yA′)2+C1(yG′-yA′)+C0；
其中，C0＝φA′，yG′表示G′点纵坐标，yA′表示A′点纵坐标，h′为网格单元边长；
基于一维CIP插值算法，通过网格点A′和点D′的得到点G′的其计算公式如下：
其中， yG′表示G′点纵坐标，yA′表示A点纵坐标，h′为网格单元边长。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910036395.3A CN109801348B (zh) | 2019-01-15 | 2019-01-15 | 一种适用于三维流体模拟的三阶高精度对流插值算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910036395.3A CN109801348B (zh) | 2019-01-15 | 2019-01-15 | 一种适用于三维流体模拟的三阶高精度对流插值算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109801348A true CN109801348A (zh) | 2019-05-24 |
| CN109801348B CN109801348B (zh) | 2022-03-04 |
Family
ID=66559093
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910036395.3A Active CN109801348B (zh) | 2019-01-15 | 2019-01-15 | 一种适用于三维流体模拟的三阶高精度对流插值算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109801348B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427700A (zh) * | 2019-08-01 | 2019-11-08 | 中国人民解放军空军工程大学 | 针对尾流场不规则三维散点坐标的栅格吸附插值拟合方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6266071B1 (en) * | 1998-11-20 | 2001-07-24 | Silicon Graphics, Inc. | Method of producing fluid-like animations using a rapid and stable solver for the Navier-Stokes equations |
| US20090040220A1 (en) * | 2007-02-05 | 2009-02-12 | Jonathan Gibbs | Hybrid volume rendering in computer implemented animation |
| US20100250213A1 (en) * | 2009-03-31 | 2010-09-30 | Doyub Kim | Semi-lagrangian cip fluid solver without dimensional splitting |
| US20130027407A1 (en) * | 2011-07-27 | 2013-01-31 | Dreamworks Animation Llc | Fluid dynamics framework for animated special effects |
-
2019
- 2019-01-15 CN CN201910036395.3A patent/CN109801348B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6266071B1 (en) * | 1998-11-20 | 2001-07-24 | Silicon Graphics, Inc. | Method of producing fluid-like animations using a rapid and stable solver for the Navier-Stokes equations |
| US20090040220A1 (en) * | 2007-02-05 | 2009-02-12 | Jonathan Gibbs | Hybrid volume rendering in computer implemented animation |
| US20100250213A1 (en) * | 2009-03-31 | 2010-09-30 | Doyub Kim | Semi-lagrangian cip fluid solver without dimensional splitting |
| US20130027407A1 (en) * | 2011-07-27 | 2013-01-31 | Dreamworks Animation Llc | Fluid dynamics framework for animated special effects |
Non-Patent Citations (3)
| Title |
|---|
| 徐英等: "多步高阶隐式泰勒级数法暂态稳定计算", 《电力系统保护与控制》 * |
| 李小纲等: "三维泊松方程基于Richardson外推的高精度多重网格解法", 《兰州理工大学学报》 * |
| 詹毅等: "泰勒展开图像插值的一个改进算法", 《重庆文理学院学报(社会科学版)》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110427700A (zh) * | 2019-08-01 | 2019-11-08 | 中国人民解放军空军工程大学 | 针对尾流场不规则三维散点坐标的栅格吸附插值拟合方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109801348B (zh) | 2022-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Eisenschmidt et al. | Direct numerical simulations for multiphase flows: An overview of the multiphase code FS3D | |
| Wicke et al. | Dynamic local remeshing for elastoplastic simulation | |
| Thomaszewski et al. | Continuum‐based strain limiting | |
| Kineri et al. | B-spline surface fitting by iterative geometric interpolation/approximation algorithms | |
| Di et al. | Moving mesh finite element methods for the incompressible Navier--Stokes equations | |
| Diziol et al. | Robust real-time deformation of incompressible surface meshes | |
| CN112862972A (zh) | 一种表面结构网格生成方法 | |
| ES2770022T3 (es) | Dispositivo y procedimiento de medición de una magnitud física de un flujo de fluido | |
| Ding et al. | An efficient dynamic mesh generation method for complex multi-block structured grid | |
| CN110415352A (zh) | 基于闭合轮廓线的三维闭合曲面构造算法 | |
| Bender et al. | Adaptive cloth simulation using corotational finite elements | |
| Aubry et al. | A three-dimensional parametric mesher with surface boundary-layer capability | |
| CN109801348A (zh) | 一种适用于三维流体模拟的三阶高精度对流插值算法 | |
| Shi et al. | Filling n-sided regions with G 1 triangular Coons B-spline patches | |
| Xi et al. | Survey on smoothed particle hydrodynamics and the particle systems | |
| CN109741428A (zh) | 一种适用于二维流体模拟的三阶高精度对流插值算法 | |
| CN104572575A (zh) | 一种特大变形动态网格生成方法 | |
| Bender et al. | Efficient Cloth Simulation Using an Adaptive Finite Element Method. | |
| Li et al. | Predicting unsteady incompressible fluid dynamics with finite volume informed neural network | |
| Degtyarev et al. | Coordinate systems, numerical objects and algorithmic operations of computational experiment in fluid mechanics | |
| Marcotte et al. | Implementation of PDE models of cardiac dynamics on GPUs using OpenCL | |
| Huang et al. | A new grid deformation technology with high quality and robustness based on quaternion | |
| Yu et al. | Image deformation using modified moving least squares with outlines | |
| Selig | Curves of stationary acceleration in SE (3) | |
| Wang et al. | An efficient preprocessing and composition based finite-time Lyapunov exponent visualization algorithm for unsteady flow field |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |