CN105929193B - 一种基于流体力学连续方程的速度场快速修正方法及装置 - Google Patents
一种基于流体力学连续方程的速度场快速修正方法及装置 Download PDFInfo
- Publication number
- CN105929193B CN105929193B CN201610236214.8A CN201610236214A CN105929193B CN 105929193 B CN105929193 B CN 105929193B CN 201610236214 A CN201610236214 A CN 201610236214A CN 105929193 B CN105929193 B CN 105929193B
- Authority
- CN
- China
- Prior art keywords
- velocity field
- ijk
- continuity equation
- divergence
- velocity
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012937 correction Methods 0.000 title claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 32
- 238000009795 derivation Methods 0.000 claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims description 15
- 239000013598 vector Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000007935 neutral effect Effects 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 5
- 238000013383 initial experiment Methods 0.000 claims description 4
- 238000002715 modification method Methods 0.000 abstract description 5
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 4
- 238000000917 particle-image velocimetry Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于流体力学连续方程的速度场快速修正方法及装置，速度场快速修正装置，包括连续方程离散形式模块、速度场散度计算模块、求导算子分解模块、速度场重构模块，速度场快速修正方法通过对一维求导算子的矩阵的特征值分解，实现了对修正速度场的直接计算，得到的速度场在离散形式下处处满足连续方程，且是对原始速度场的最保守的修正，本发明能大大提高计算效率，降低了内存的消耗，使得修正方法具有更大的适用性。
Description
技术领域
本发明涉及一种基于流体力学连续方程的速度场快速修正方法及装置，属于流体力学速度测量技术领域。
背景技术
图像粒子测速技术(Particle Image Velocimetry，PIV)是实验流体力学领域最重要的定量测量手段之一。该技术通过分析示踪粒子在相机上的成像，获得测量区域内的速度场。随着层析粒子图像测速(Tomographic Particle Image Velocimetry，TPIV)等测量技术的发展，三维体空间内全部三个速度分量(Three-dimensional Three-component,3D3C)的获取越来越方便。3D3C速度场包含所有的速度梯度分量，允许使用流体力学的基本方程，对速度场测量精度进行检验。对于不可压缩流体，流动满足连续方程(设u,v,w是三个速度分量；表示速度场的散度)，因此速度场的散度应该处处为零。然而实验测量的结果由于存在误差，所以不能严格满足全场散度为零，有不容忽视的散度残差存在。这些残差不仅会影响流场显示的结果，造成统计量的偏差，还会传播到由速度场重构的压力场中。基于连续方程速度场修正的任务是消除原始速度场散度的残差，或者说从原始测量速度场中重构出一个新的满足连续方程的速度场。前人的研究结果表明，通过对速度场的修正可以减小测量误差10％-20％。
目前，对于速度场散度残差的修正目前已有诸多方法，包括求解泊松方程、无散基函数重构和数学优化方法。一种基于非线性优化问题的最小二乘法修正方法是最简便的方法，并被数值与实验测试证明能够有效地改善湍流统计量。这种最小二乘法修正方法以连续方程为约束条件，目标函数为修正量的二范数的优化修正方法，旨在寻求满足连续方程的最小修正。基于优化方法的修正方案，可以保证对原始速度场进行最小的修正，从而使修正速度场忠于原始实验数据。然而，目前求解这种优化问题的方法既耗内存，效率又低，缺乏实用性，对计算机硬件要求较高，很能难用于实际的速度场修正。其他的非优化的修正方法不能保证是对原始速度场的最小修正，因此可能造成对速度场的过度修正，从而引入额外的误差。
发明内容
本发明的目的是为了解决上述优化问题，提出一种快速、精确的速度场修正方法及装置，降低了对计算机硬件的要求，使速度场修正更容易实现。
本发明的一种基于流体力学连续方程的速度场快速修正方法，通过对一维求导算子的矩阵的特征值分解，实现了对修正速度场的直接计算，得到的速度场在离散形式下处处满足连续方程，且是对原始速度场的最保守的修正，本发明能大大提高计算效率，降低了内存的消耗，使得修正方法具有更大的适用性，具体包括以下几个步骤：
步骤101：确定连续方程的离散形式；
步骤102：计算原始速度场的散度；
步骤103：对矩阵DnxDnx T,DnyDnx T,DnzDnx T进行特征值分解；
步骤104：完成对速度场的修正。
本发明的一种基于流体力学连续方程的速度场快速修正装置，包括连续方程离散形式模块、速度场散度计算模块、求导算子分解模块、速度场重构模块；
连续方程离散形式模块通过对速度场网格坐标进行分析，得到连续方程的离散形式，将散度算子通过一维求导算子表示，同时得到一维求导算子矩阵的形式。
速度场散度计算模块通过连续方程的离散形式，计算原始速度场散度的残差。
求导算子分解模块通过对一维求导算子进行特征值分解，得到该算子的特征值与特征向量。
速度场重构模块根据原始速度场散度的残差，利用算子的特征值与特征向量重构出满足连续方程的速度场，即修正速度场。该模块可以通过计算机的并行计算实现，大大提高了计算效率。
本发明的优点在于：
(1)本发明方法能够有效的消除测量误差造成的速度场散度的残差，修正后的速度场处处满足离散形式的连续性方程；
(2)本发明方法是对原始速度场的最小修正，修正速度场高度忠于原始实验数据；
(3)本发明方法通过对散度算子的降维分解，实现对速度场的快速修正；方法简单、直接，容易程序实现；
(4)本发明方法对内存占用较低，允许对三维大空间内的速度场快速修正；
(5)本发明方法可以通过矩阵运算实现，便于更快速的并行计算；
(6)本发明提供的速度场修正方法及装置，能够实现对实验测量的三维速度场的快速修正，修正速度场忠实于原始实验数据，又满足流体力学的控制方程之一：连续方程，从而能更好地反应真实流场的特征。
附图说明
图1速度场修正流程图；
图2速度场修正效果图。
图3本发明实施例装置的结构示意图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
实施例一、
本发明的一种基于流体力学连续方程的速度场快速修正方法，流程如图1所示，具体实施步骤如下：
步骤101：确定连续方程的离散形式；
不可压缩流动的TPIV测量结果是分布在nx×ny×nz规则结构网格结点上的3D3C瞬时速度场，Δx,Δy,Δz是x,y,z三个方向的网格间距。设u,v,w是三个速度分量，流体力学连续方程的形式为该方程涉及速度场的沿着三个坐标方向的偏导数。为了方便后面的运算，首先给出连续方程的离散形式。假设u在网格结点(i,j,k)上的取值为uijk，相应的取值为在网格内部结点采用中心差分格式，在边界上采用前向或者后向差分格式，相应的偏导数计算公式得到
可以简写为
其中，∑是求和符号，i'是遍历的指标索引，Dnx ii'是下列矩阵Dnx的i行i'列元素
求导算子矩阵Dnx与速度场的偏导数计算有关，另外两项偏导数也可以类似地写出：
其中，j',k'是遍历的指标索引，分别为在网格结点(i,j,k)上的取值，Dny jj'Dnx kk'分别是下面矩阵的j(k)行j'(k')列元素
于是连续方程的离散形式为
最后指出，Dnx,Dny,Dnz与速度场的偏导数计算有关，这里称为求导算子矩阵。采用不同的离散格式，求导算子矩阵的具体形式也各异，例如采用更高阶精度的中心格式或者迎风格式等。
步骤102：计算原始速度场的散度。
速度场散度通过公式计算得到，So ijk是原始速度场散度在结点(i,j,k)处的数值，利用步骤101给出的Dnx ii'Dny jj'Dnx kk'可以得到速度场散度；
其中：uo,vo,wo分别是实验测量的原始速度场uo的三个速度分量，下标o表示原始实验测量的结果。
速度场散度的计算方式同样可以采用不同的离散形式。速度场散度的大小反应实验测量误差的大小。本发明的任务是消除散度的残差，因此首先需要对原始速度场散度的残差进行计算。
步骤103：对矩阵DnxDnx T,DnyDnx T,DnzDnx T进行特征值分解。
其中Φnx,Φny,Φnz分别是DnxDnx T,DnyDnx T,DnzDnx T的特征向量矩阵，
设相应的Φnx mn,Φny mn,Φnz mn分别表示Φnx,Φny,Φnz的m行n列元素。Λnx,Λny,Λnz是特征值矩阵(对角阵)。Λnx l,Λny l,Λnz l分别表示Λnx,Λny,Λnz的第l个对角元。
通过对DnxDnx T,DnyDnx T,DnzDnx T三个矩阵进行分解，得到了正交的特征向量矩阵。利用特征向量矩阵可以方便快速地对速度场进行修正。这三个矩阵维度分别为nx×nx，ny×ny，nz×nz，远远小于速度场的维度nx×ny×nz，因此可以很快地计算特征值与特征向量。而且，由于三个矩阵均是实对称阵，它们的所有特征值均为实数，所有特征向量正交，这为后面的计算提供了方便。
步骤104：依次按照如下公式完成对速度场的修正：
Γijk＝Λnx i+Λny j+Λnz k (12)
其中，Φnx mn,Φny mn,Φnz mn与Λnx l,Λny l,Λnz l为求导算子的特征值与特征向量，由步骤103得到。So i'j'k'是原始实验速度场散度，由步骤102得到。(i,j,k)(i',j',k')(l,m,n)均为三个方向的坐标索引，取值范围均为(1～nx,1～ny,1～nz)。Γijk与μijk是计算的中间变量，与速度场散度So ijk类似。另外，Γijk一般会出现一个零元素。为了方便后面的公式计算，令Γijk所有零元素均改为1。实际上，Γijk零元素改成任意非零数字均不影响修正结果。可以证明，对Γijk零元素修改后，按照上述公式计算的速度场是满足连续方程的最小修正的结果。uc ijk,vc ijk,wc ijk是修正速度场，即本发明的最终目标。
该步骤运用了步骤103矩阵分解的结果，将对速度场的修正过程通过明确的公式计算来实现。这是本发明得以快速实现的重要原因。另一方面，上述计算公式均是对原始速度场的线性运算，可以通过矩阵乘法运算来实现，这非常利于并行运算。
本发明得到的修正速度场uc ijk,vc ijk,wc ijk，消除了原始速度场散度的残差，对测量误差具有抑制作用，提供了一个更能反映真实流场结构的速度场，为后面对流动机理的研究提供了更好的素材。
实施例二、
下面结合具体应用场景对本发明方法进行描述。
下面速度场修正的实施是针对如图2a所示的施加高斯噪声的瞬时速度场进行的。也可以采用不同人工噪音源进行测试。
在直接数值模拟(Direct Numerical Simulation，DNS)的一个100×100×16速度场上施加高斯随机噪声。Δx＝Δy＝Δz＝1。高斯随机噪声的均值μ为0，方差σ是速度方差的0.1倍。高斯随机噪声的概率密度函数f(x)为
施加噪声的速度场用以模拟TPIV得到的原始速度场Uo,Vo,Wo。图2a展示了施加噪声后速度场Uo分量云图。下面对该速度场进行修正。
步骤一：
根据连续方程的离散形式，生成三个求导算子矩阵Dnx,Dny,Dnz。采用中心差分格式，边界上采用前向或后向差分格式。Dnx,Dny,Dnz的形式如下
步骤二：
根据步骤102的公式(10)计算速度场的散度So ijk。
其中，Dnx mn,Dny mn,Dnz mn分别是Dnx,Dny,Dnz的m行n列元素。
原始速度场散度的概率密度分布曲线如图c所示。可以看到未经过修正的速度场存在不可忽略的散度残差。本发明的任务是在对原始速度场尽可能少修正的前提下，消除这些残差。
步骤三：
对矩阵DnxDnx T,DnyDnx T,DnzDnx T进行特征值分解，如公式(11)所示。得到相应的特征值Λnx,Λny,Λnz与特征向量Φnx,Φny,Φnz。Φnx,Φny,Φnz分别为nx×nx，ny×ny，nz×nz的标准正交阵，这有利于后面的计算。Λnx,Λny,Λnz是对角阵，且其对角元均是不小于0的实数。由于DnxDnx T,DnyDnx T,DnzDnx T维度为nx×nx，ny×ny，nz×nz，远小于速度场的维度nx×ny×nz，该步骤计算较快。
步骤四：
按照步骤104的公式(11)-(14)依次计算Γijk，μijk与修正速度场uc ijk,vc ijk,wc ijk。Γijk的结果中有一个0元素。为了避免公式(12)对0元素取倒数的运算，这里将零元素改写为1后进行公式(12)的计算。这一操作不会改变修正速度场满足连续方程的性质，而且能够保证对速度场修正量最小。该步骤得到的修正后的速度场uc ijk在中间截面上的分布如图b所示。
本发明通过对速度场做最小修正，得到散度处处为零(<10-14)，即离散形式下满足流体力学的连续方程的速度场。在MATLAB环境下，利用本发明的方法，完成对测试速度场的修正仅需0.1s。经过本发明实施后，速度场噪声有所减小，对原始速度场的修正量较小，如图2a,b。但散度残差的概率密度曲线有明显变化，由图2c变为图2d，基本上消除了原始速度场散度的残差。当测试速度场扩大到的1000×1000×160数据点，完成修正的时间约为60s，计算机占用内存小于20Gb。而这样大小的速度场在传统的最小二乘法修正方法下是很难实现的，因为计算机内存消耗过大。
实施例三、
本发明实施例还提供了一种基于流体力学连续方程的速度场快速修正装置，如图3所示，该装置可位于计算机内部，包括连续方程离散形式模块301、速度场散度计算模块302、离散求导算子分解模块303、速度场重构模块304。
连续方程离散形式模块301通过对速度场网格坐标进行分析，得到连续方程的离散形式，将散度算子通过一维求导算子表示，至少得到一维求导算子矩阵的形式。
速度场散度计算模块302通过连续方程的离散形式，计算原始速度场散度的残差，原始速度场散度的残差包括所有测量结点上速度场散度的取值。
离散求导算子分解模块303通过对一维求导算子进行特征值分解，至少得到该算子的特征值与特征向量，为后面散度算子求逆提供方便。
速度场重构模块304根据原始速度场散度的残差，利用算子的特征值与特征向量重构出满足连续方程的速度场。该模块通过计算机的并行计算实现，大大提高计算的效率。
本发明实施例中，所述各计算模块均可通过计算机中的中央处理器(CentralProcessing Unit，CPU)、数字信号处理器(Digital Signal Processor，DSP)或可编程逻辑阵列(Field－Programmable Gate Array，FPGA)实现。
本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(装置)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。
Claims (3)
1.一种基于流体力学连续方程的速度场快速修正方法，具体包括以下几个步骤：
步骤101：确定连续方程的离散形式；
具体为：
设不可压缩流动的TPIV测量结果是分布在nx×ny×nz规则结构网格结点上的3D3C瞬时速度场，Δx,Δy,Δz是x,y,z三个方向的网格间距，
设u,v,w是三个速度分量，流体力学连续方程为假设u在网格结点(i,j,k)上的取值为uijk，相应的取值为在网格内部结点采用中心差分格式，在边界上采用前向或者后向差分格式，相应的偏导数计算公式为：
简写为
其中，∑是求和符号，i'是遍历的指标索引，Dnx ii'是下列矩阵Dnx的i行i'列元素
求导算子矩阵Dnx与速度场的偏导数计算有关，另外两项偏导数 也写为：
其中，j',k'是遍历的指标索引，分别为在网格结点(i,j,k)上的取值，Dny jj'Dnx kk'分别是下面矩阵的j(k)行j'(k')列元素
连续方程的离散形式为：
步骤102：计算原始速度场的散度；
具体为：
速度场散度通过公式计算得到，So ijk是原始速度场散度在结点(i,j,k)处的数值，根据Dnx ii'Dny jj'Dnx kk'得到速度场散度；
其中：uo,vo,wo分别是实验测量的原始速度场uo的三个速度分量，下标o表示原始实验测量的结果；
步骤103：对矩阵DnxDnx T,DnyDnx T,DnzDnx T进行特征值分解；
具体为：
其中：Φnx,Φny,Φnz分别是DnxDnx T,DnyDnx T,DnzDnx T的特征向量矩阵，
设相应的Φnx mn,Φny mn,Φnz mn分别表示Φnx,Φny,Φnz的m行n列元素；Λnx,Λny,Λnz是特征值矩阵；Λnx l,Λny l,Λnz l分别表示Λnx,Λny,Λnz的第l个对角元；
通过对DnxDnx T,DnyDnx T,DnzDnx T三个矩阵进行分解，得到了正交的特征向量矩阵；
步骤104：完成对速度场的修正；
具体为：
Γijk＝Λnx i+Λny j+Λnz k (12)
其中，Φnx mn,Φny mn,Φnz mn与Λnx l,Λny l,Λnz l为求导算子的特征值与特征向量，So i'j'k'是原始实验速度场散度，(i,j,k)(i',j',k')(l,m,n)均为三个方向的坐标索引，取值范围均为(1～nx,1～ny,1～nz)；Γijk与μijk是计算的中间变量，令Γijk所有零元素均改为任意非零数字；
最终得到修正速度场：uc ijk,vc ijk,wc ijk。
2.一种基于流体力学连续方程的速度场快速修正装置，包括连续方程离散形式模块、速度场散度计算模块、求导算子分解模块、速度场重构模块；
连续方程离散形式模块通过对速度场网格坐标进行分析，得到连续方程的离散形式，将散度算子通过一维求导算子表示，同时得到一维求导算子矩阵的形式；
速度场散度计算模块通过连续方程的离散形式，计算原始速度场散度的残差；
求导算子分解模块通过对一维求导算子进行特征值分解，得到该算子的特征值与特征向量；
速度场重构模块根据原始速度场散度的残差，利用算子的特征值与特征向量重构出满足连续方程的速度场，即修正速度场；
具体为：
Γijk＝Λnx i+Λny j+Λnz k (12)
其中，Φnx mn,Φny mn,Φnz mn与Λnx l,Λny l,Λnz l为求导算子的特征值与特征向量，So i'j'k'是原始实验速度场散度，(i,j,k)(i',j',k')(l,m,n)均为三个方向的坐标索引，取值范围均为(1～nx,1～ny,1～nz)；Γijk与μijk是计算的中间变量，令Γijk所有零元素均改为任意非零数字；
最终得到修正速度场：uc ijk,vc ijk,wc ijk。
3.根据权利要求2所示的一种基于流体力学连续方程的速度场快速修正装置，所述的速度场散度计算模块通过连续方程的离散形式，计算原始速度场散度的残差，原始速度场散度的残差包括所有测量结点上速度场散度的取值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610236214.8A CN105929193B (zh) | 2016-04-15 | 2016-04-15 | 一种基于流体力学连续方程的速度场快速修正方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610236214.8A CN105929193B (zh) | 2016-04-15 | 2016-04-15 | 一种基于流体力学连续方程的速度场快速修正方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105929193A CN105929193A (zh) | 2016-09-07 |
| CN105929193B true CN105929193B (zh) | 2019-01-04 |
Family
ID=56839094
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610236214.8A Expired - Fee Related CN105929193B (zh) | 2016-04-15 | 2016-04-15 | 一种基于流体力学连续方程的速度场快速修正方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105929193B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107064546B (zh) * | 2017-05-10 | 2019-09-10 | 东南大学 | 一种基于非圆信号频率估计算法的流体速度估计方法 |
| CN108805949B (zh) * | 2018-05-25 | 2022-04-12 | 杭州晟视科技有限公司 | 一种血流速度的修正方法、装置、终端和计算机存储介质 |
| CN108962392A (zh) * | 2018-05-29 | 2018-12-07 | 杭州晟视科技有限公司 | 一种速度场的修正方法、装置及存储介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101629965A (zh) * | 2009-08-18 | 2010-01-20 | 清华大学深圳研究生院 | 用于粒子图像测速中的多网格处理方法 |
| CN102982202A (zh) * | 2012-11-15 | 2013-03-20 | 中国电子科技集团公司第二十八研究所 | 基于缺损模态的结构模型修正方法 |
| CN103729564A (zh) * | 2014-01-06 | 2014-04-16 | 北京航空航天大学 | 一种基于粒子图像测速技术的压力场计算方法和装置 |
| US20140149055A1 (en) * | 2012-05-04 | 2014-05-29 | The Regents Of The University Of California | Multi-plane method for three-dimensional particle image velocimetry |
| CN104915928A (zh) * | 2014-03-12 | 2015-09-16 | 北京航空航天大学 | 一种基于本征正交分解的速度场坏矢量识别和修正方法 |
| CN105095555A (zh) * | 2014-07-15 | 2015-11-25 | 北京航空航天大学 | 一种速度场的无散平滑处理方法及装置 |
-
2016
- 2016-04-15 CN CN201610236214.8A patent/CN105929193B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101629965A (zh) * | 2009-08-18 | 2010-01-20 | 清华大学深圳研究生院 | 用于粒子图像测速中的多网格处理方法 |
| US20140149055A1 (en) * | 2012-05-04 | 2014-05-29 | The Regents Of The University Of California | Multi-plane method for three-dimensional particle image velocimetry |
| CN102982202A (zh) * | 2012-11-15 | 2013-03-20 | 中国电子科技集团公司第二十八研究所 | 基于缺损模态的结构模型修正方法 |
| CN103729564A (zh) * | 2014-01-06 | 2014-04-16 | 北京航空航天大学 | 一种基于粒子图像测速技术的压力场计算方法和装置 |
| CN104915928A (zh) * | 2014-03-12 | 2015-09-16 | 北京航空航天大学 | 一种基于本征正交分解的速度场坏矢量识别和修正方法 |
| CN105095555A (zh) * | 2014-07-15 | 2015-11-25 | 北京航空航天大学 | 一种速度场的无散平滑处理方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| "Divergence-free smoothing for 3D-3C PIV data";王成跃 等;《4TH INTERNATIONAL CONFERENCE ON EXPERIMENTAL FLUID MECHANICS》;20140831;第1-12页 |
| 王成跃 等."Divergence‑free smoothing for volumetric PIV data".《Springer》.2016,第57卷(第15期),第1-23页. |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105929193A (zh) | 2016-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Escalante et al. | Non-hydrostatic pressure shallow flows: GPU implementation using finite volume and finite difference scheme | |
| Erhel et al. | Flow simulation in three-dimensional discrete fracture networks | |
| Xiao et al. | Physical mechanism of the inverse energy cascade of two-dimensional turbulence: a numerical investigation | |
| CN103729564B (zh) | 一种基于粒子图像测速技术的压力场计算方法和装置 | |
| CN105929193B (zh) | 一种基于流体力学连续方程的速度场快速修正方法及装置 | |
| CN103091395B (zh) | 油藏条件下多孔介质中流体相态三维动态探测方法及装置 | |
| Menshov et al. | Free-boundary method for the numerical solution of gas-dynamic equations in domains with varying geometry | |
| Liu et al. | An interior penalty discontinuous Galerkin approach for 3D incompressible Navier–Stokes equation for permeability estimation of porous media | |
| CN105807093B (zh) | 一种基于粒子图像测速技术的加速度测量方法及装置 | |
| CN105095555B (zh) | 一种基于速度场无散平滑处理的粒子图像测速方法及装置 | |
| Cao | Binary black hole simulation with an adaptive finite element method: Solving the Einstein constraint equations | |
| US20170116774A1 (en) | Method and apparatus generating acceleration structure | |
| Assonitis et al. | A new shock-fitting technique for 2-D structured grids | |
| JP7485076B2 (ja) | 推定装置、推定方法及び推定プログラム | |
| He et al. | Accelerating full waveform inversion using HSS solver and limited memory conjugate gradient method | |
| Nouveau et al. | An adaptive, residual based, splitting approach for the penalized navier stokes equations | |
| CN108376123A (zh) | 一种卫星重力观测数据异常值提取方法、装置及电子设备 | |
| Kramer et al. | Solving the Poisson equation on small aspect ratio domains using unstructured meshes | |
| CN104777329B (zh) | 一种用于粒子图像测速三维粒子场重构的线性规划算法 | |
| Bartana* et al. | GPU implementation of minimal dispersion recursive operators for reverse time migration | |
| Bakulin et al. | Benchmarking 3D time-and frequency-domain solvers for FWI applications for different cluster sizes and variable number of sources | |
| JE | High order mimetic finite difference operators satisfying a Gauss divergence theorem | |
| Ran et al. | A new Poisson-type equation applicable to the three-dimensional non-hydrostatic model in the framework of the discontinuous Galerkin method | |
| CN116794735B (zh) | 航空磁矢量梯度数据等效源多分量联合去噪方法及装置 | |
| Baxansky et al. | Calculating geometric properties of three-dimensional objects from the spherical harmonic representation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190104 |