CN107871337A - 一种超音速二维流场数据的可视化方法 - Google Patents

一种超音速二维流场数据的可视化方法 Download PDF

Info

Publication number
CN107871337A
CN107871337A CN201610850404.9A CN201610850404A CN107871337A CN 107871337 A CN107871337 A CN 107871337A CN 201610850404 A CN201610850404 A CN 201610850404A CN 107871337 A CN107871337 A CN 107871337A
Authority
CN
China
Prior art keywords
mrow
mfrac
mtd
mtr
density
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
Application number
CN201610850404.9A
Other languages
English (en)
Other versions
CN107871337B (zh
Inventor
刘俊
路波
王显圣
何彬华
梁锦敏
杨可
陈植
刘志勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Original Assignee
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center filed Critical High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority to CN201610850404.9A priority Critical patent/CN107871337B/zh
Publication of CN107871337A publication Critical patent/CN107871337A/zh
Application granted granted Critical
Publication of CN107871337B publication Critical patent/CN107871337B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/005General purpose rendering architectures

Abstract

本发明提供了一种超音速二维流场数据的可视化方法,该方案包括有以下步骤:a、建立流场密度值与RGB颜色的映射关系,通过插值获得网格单元的RGB颜色分布;b、将无量纲密度作为网格点的高度值,得到类似水流自由表面的三维密度场曲面,并计算各网格单元光照强度;c、将每个网格单元的RGB颜色和光照强度相乘,获得最终的颜色分布。该方案能同时扑捉空间变化剧烈的激波结构和变化较平缓的旋涡结构,并且呈现均匀区流场密度分布信息,实现高保真的二维超音速流场显示。

Description

一种超音速二维流场数据的可视化方法
技术领域
本发明涉及的是计算流体力学和计算机图形学的交叉领域,尤其是一种超音速二维流场数据的可视化方法。
背景技术
随着高性能计算技术的飞速发展,科学计算能力得到巨大提升,数值计算产生的非定常多维多场多变量数据呈几何级数增长。面对海量的仿真数据,如何开展分析研究,提取有用信息,已成为科学家亟需解决的现实难题。其中,数据可视化方法是进行数据分析的一种有效手段。这种方法将计算数据绘制成形象直观的图形图像,让研究人员的关注点从一个个离散的数据中解放出来,通过对整体图形图像的分析,提高对全局的把握能力,更易于发现一些新的规律和现象。
在计算流体力学领域,流场数据可视化技术已经获得广泛应用,市场上出现了Tecplot、ParaView、EnSight、FieldView等多款专业流场可视化软件。流场可视化目前已经成为流场数值模拟的重要环节之一。可视化方法多种多样,效果也不尽相同。一般情况下,根据研究的关注点,选取一种或几种合适的可视化方法来分析数据。一种好的可视化方法能够清晰呈现流场结构和物理量的分布特性,让研究人员清楚地观察到流动结构发展演化过程,从而对机理分析起到促进作用。
近年来,流场数据可视化研究取得了长足的进展,已经发展出多种方法来显示流场信息。针对标量场,常见的可视化方法有等值线、云图等;针对矢量场,常见的可视化方法有点矢量、流线、纹理等。每种方法都有各自的特点,如等值线在流场梯度较大的地方线条密集分布,能有效的刻画激波、接触间断等流动结构,但是对于旋涡等空间变化较为平缓的结构则难以捕捉;云图通过建立流场数据值与颜色之间的映射关系,通过颜色的变化来反应全场数据的变化规律,有助于帮助研究人员全面认识全场物理量的分布,但是对流动结构的分辨能力不足;点矢量在空间的每个网格点以箭头的方式反映当地矢量的大小和方向,对于网格点较密的情况,箭头或者相互重叠或者很小难以分辨,显示效果并不理想;流线是某一固定时刻的空间曲线,该曲线上任意一点的切向量与当地速度矢量重合,在定常流场中流线能较好的刻画旋涡等结构,在非定常情况下,流线则较为紊乱;纹理是近年来发展的一种新的数据可视化方法,纹理是兼具形状和颜色两种属性的颜色有序排列,形状属性可表示流场方向信息,颜色属性可以表示标量信息。纹理对旋涡有很好的显示效果,但是对激波结构的分辨率不足。
超音速流动一直以来都是流体力学研究的热点和难点问题。近年来,研究人员开发出高精度本质无振荡的WENO、WCNS等数值格式,实现了强激波条件下的高精度数值模拟。高精度的计算结果需要匹配高精度的可视化方法,才能最大限度的展现高精度算法在精细刻画流场方面的优势,帮助科研人员准确分析流场特性。超音速流场中一般兼具空间变化剧烈的激波结构和变化较平缓的旋涡结构,如何同时清晰扑捉这两种结构,以及呈现均匀区流场物理信息,目前为止,没有很好的可视化方法能够做到。
发明内容
本发明的目的,就是针对现有技术所存在的不足,提供一种超音速二维流场数据的可视化方法,该方案能同时扑捉空间变化剧烈的激波结构和变化较平缓的旋涡结构,并且呈现均匀区流场密度分布信息,实现高保真的二维超音速流场显示。
本方案是通过如下技术措施来实现的:
1.一种超音速二维流场数据的可视化方法,其特征是:包括有以下步骤:
a、建立流场密度值与RGB颜色的映射关系,通过插值获得网格单元的RGB颜色分布;
b、将无量纲密度作为网格点的高度值,得到类似水流自由表面的三维密度场曲面,并计算各网格单元光照强度;
c、将每个网格单元的RGB颜色和光照强度相乘,获得最终的颜色分布。
2.根据权利要求1所述的一种超音速二维流场数据的可视化方法,其特征是:所述步骤a中,通过线性插值的方法建立流场密度与RGB颜色分量的对应关系;选取的插值对象为Small-Rainbow颜色组合,该颜色组合依次由蓝(0,0,255)、青(0,255,255)、绿(0,255,0)、黄(255,255,0)和红(255,0,0)等五种颜色组成,详细差值公式如下:
式中,ρ为流场密度,ρmin为流场密度的最小值,ρmax为流场密度的最大值,r、g、b分别代表某一密度值ρ插值后的红、绿、蓝三种颜色分量。
3.根据权利要求1所述的一种超音速二维流场数据的可视化方法,其特征是:所述步骤b中,对于XY平面的二维超音速流场,选取无量纲密度作为网格点的高度值,将2D网格转换成2.5D网格;通过公式(2)计算得到三维密度场曲面一点处的法向量然后,根据郎伯特余弦定律计算当地的光照强度,选择入射光为Z轴方向,将入射光与表面法向量的点积作为漫反射光照强度因子(见式3),最终得到光照强度是密度梯度幅值的函数;当密度和坐标在进行无量纲化时,若选择不同参考值,会导致最终的密度梯度幅值发生变化,在光照强度计算公式(见式4)中引入调节参数k,从而计及这方面的影响;
具体公式如下:
式中,分别为密度曲面沿x和y方向的切向量,为曲面的单位法向量,k为调节参数,σ为光照强度因子。
4.根据权利要求1所述的一种超音速二维流场数据的可视化方法,其特征是:所述步骤c中,将插值后的RGB颜色分别与光照强度相乘就得到了合成后的颜色分布如式5所示,然后调节参数k,获得明暗适中的图像,具体涉及公式如下:
R=σr
G=σg
B=σb (5)
式中,R、G、B分别代表合成后颜色的红、绿、蓝分量。
本方案的有益效果可根据对上述方案的叙述得知,由于在该方案中选取无量纲密度作为网格点的高度值,得到类似水流自由表面的三维密度场曲面,从而实现了激波、旋涡等流动结构的立体直观显示。该方案与其他可视化方法对比,本方法对激波、旋涡结构的显示更为清晰、干净,弱激波、声波等流动细节更丰富。该方法充分利用了光强信息,使流动结构呈现出立体效果,也更加符合人类视觉的认知模式。
由此可见,本发明与现有技术相比,具有实质性特点和进步,其实施的有益效果也是显而易见的。
附图说明
图1为采用密度云图和等值线方法得到双马赫反射流场效果图。
图2为采用纹理法得到的超音速空腔流场纹理图。
图3为采用本发明方法得到的双马赫反射流场效果图。
图4为采用本发明方法得到的超音速空腔流场效果图。
具体实施方式
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
本说明书(包括任何附加权利要求、摘要和附图)中公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换。即,除非特别叙述,每个特征只是一系列等效或类似特征中的一个例子而已。
本方案包括有以下步骤:
a、建立流场密度值与RGB颜色的映射关系,通过插值获得网格单元的RGB颜色分布;
b、将无量纲密度作为网格点的高度值,得到类似水流自由表面的三维密度场曲面,并计算各网格单元光照强度;
c、将每个网格单元的RGB颜色和光照强度相乘,获得最终的颜色分布。
步骤a中,通过线性插值的方法建立密度与RGB颜色分量的对应关系;选取的插值对象为Small-Rainbow颜色组合,该颜色组合依次由蓝(0,0,255)、青(0,255,255)、绿(0,255,0)、黄(255,255,0)和红(255,0,0)等五种颜色组成,详细差值公式如下:
式中,ρ为流场密度,ρmin为流场密度的最小值,ρmax为流场密度的最大值,r、g、b分别代表某一密度值ρ插值后的红、绿、蓝三种颜色分量。
步骤b中,对于XY平面的二维超音速流场,选取无量纲密度作为网格点的高度值,2D网格转换成2.5D网格;通过式(2)计算得到三维密度场曲面一点处的法向量然后,根据郎伯特余弦定律计算当地的光照强度,选择入射光为Z轴方向,将入射光与表面法向量的点积作为漫反射光照强度因子(见式3),最终得到光照强度是密度梯度幅值的函数。当密度和坐标在进行无量纲化时,若选择不同参考值,会导致最终的密度梯度幅值发生变化,在光照强度计算公式(见式4)中引入调节参数k,从而计及这方面的影响。具体公式如下:
式中,分别为密度曲面沿x和y方向的切向量,为曲面的单位法向量,k为调节参数,σ为光照强度因子。
步骤c中,将插值后的RGB颜色分别与光照强度相乘就得到了合成后的颜色分布如式5所示,然后调节参数k,获得明暗适中的图像,具体涉及公式如下:
R=σr
G=σg
B=σb (5)
式中,R、G、B分别代表合成后颜色的红、绿、蓝分量。
通过与其他可视化方法的对比,可进一步展示本发明提出的方法的特点。图1为双马赫反射流场,与密度云图和等值线合成效果图比较,本方法中除了强激波、旋涡之外,弱压缩波等非主要流动结构也都能清晰的显示出来,流动细节更丰富,流动结构的立体感更强。图2为二维超音速空腔流场,与纹理法相比,本方法对激波结构的显示更为清晰,在均匀区对密度分布的显示也更为干净。
本发明提出的可视化方法具有十分明确的物理意义。该方法以浅水流动和可压缩气流的相似性为理论基础,通过选取密度为作为高度值,获得了类似水流自由表面的三维密度场曲面,使激波、旋涡等流动结构得以立体呈现。该方法充分利用了光强信息,使流动结构呈现出立体效果,也更加符合人类视觉的认知模式。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。

Claims (4)

1.一种超音速二维流场数据的可视化方法,其特征是:包括有以下步骤:
a、建立流场密度值与RGB颜色的映射关系,通过插值获得网格单元的RGB颜色分布;
b、将无量纲密度作为网格点的高度值,得到类似水流自由表面的三维密度场曲面,并计算各网格单元光照强度;
c、将每个网格单元的RGB颜色和光照强度相乘,获得最终的颜色分布。
2.根据权利要求1所述的一种超音速二维流场数据的可视化方法,其特征是:所述步骤a中,通过线性插值的方法建立流场密度与RGB颜色分量的对应关系;选取的插值对象为Small-Rainbow颜色组合,该颜色组合依次由蓝(0,0,255)、青(0,255,255)、绿(0,255,0)、黄(255,255,0)和红(255,0,0)等五种颜色组成,详细差值公式如下:
<mrow> <mi>r</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mn>0</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&lt;</mo> <mn>0.5</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>255</mn> <mo>&amp;times;</mo> <mfrac> <mrow> <mi>&amp;alpha;</mi> <mo>-</mo> <mn>0.5</mn> </mrow> <mn>0.25</mn> </mfrac> </mrow> </mtd> <mtd> <mrow> <mn>0.5</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&lt;</mo> <mn>0.75</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>255</mn> </mtd> <mtd> <mrow> <mn>0.75</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&amp;le;</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
<mrow> <mi>g</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mn>255</mn> <mo>&amp;times;</mo> <mfrac> <mi>&amp;alpha;</mi> <mn>0.25</mn> </mfrac> </mrow> </mtd> <mtd> <mrow> <mn>0</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&lt;</mo> <mn>0.25</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>255</mn> </mtd> <mtd> <mrow> <mn>0.25</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&lt;</mo> <mn>0.75</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>255</mn> <mo>&amp;times;</mo> <mfrac> <mrow> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>1</mn> </mrow> <mn>0.25</mn> </mfrac> </mrow> </mtd> <mtd> <mrow> <mn>0.75</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&amp;le;</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>b</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>255</mn> </mtd> <mtd> <mrow> <mn>0</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&lt;</mo> <mn>0.25</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>255</mn> <mo>&amp;times;</mo> <mfrac> <mrow> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <mn>0.5</mn> </mrow> <mn>0.25</mn> </mfrac> </mrow> </mtd> <mtd> <mrow> <mn>0.25</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&lt;</mo> <mn>0.5</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mn>0.5</mn> <mo>&amp;le;</mo> <mi>&amp;alpha;</mi> <mo>&amp;le;</mo> <mn>1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
<mrow> <mi>&amp;alpha;</mi> <mo>=</mo> <mfrac> <mrow> <mi>&amp;rho;</mi> <mo>-</mo> <msub> <mi>&amp;rho;</mi> <mi>min</mi> </msub> </mrow> <mrow> <msub> <mi>&amp;rho;</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mfrac> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>&amp;rsqb;</mo> </mrow>
式中,ρ为流场密度,ρmin为流场密度的最小值,ρmax为流场密度的最大值,r、g、b分别代表某一密度值ρ插值后的红、绿、蓝三种颜色分量。
3.根据权利要求1所述的一种超音速二维流场数据的可视化方法,其特征是:所述步骤b中,对于XY平面的二维超音速流场,选取无量纲密度作为网格点的高度值,将2D网格转换成2.5D网格;通过公式(2)计算得到三维密度场曲面一点处的法向量然后,根据郎伯特余弦定律计算当地的光照强度,选择入射光为Z轴方向,将入射光与表面法向量的点积作为漫反射光照强度因子(见式3),最终得到光照强度是密度梯度幅值的函数;当密度和坐标在进行无量纲化时,若选择不同参考值,会导致最终的密度梯度幅值发生变化,在光照强度计算公式(见式4)中引入调节参数k,从而计及这方面的影响;
具体公式如下:
<mrow> <mover> <mi>l</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow>
<mrow> <mover> <mi>m</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>=</mo> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mover> <mi>n</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mover> <mi>l</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;times;</mo> <mover> <mi>m</mi> <mo>&amp;RightArrow;</mo> </mover> </mrow> <mrow> <mo>|</mo> <mover> <mi>l</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>|</mo> <mo>|</mo> <mover> <mi>m</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>|</mo> </mrow> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>1</mn> </mrow> </msqrt> </mfrac> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>,</mo> <mo>-</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>&amp;sigma;</mi> <mo>=</mo> <mover> <mi>n</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>&amp;rho;</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>1</mn> </mrow> </msqrt> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <mo>|</mo> <mo>&amp;dtri;</mo> <mi>&amp;rho;</mi> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>+</mo> <mn>1</mn> </mrow> </msqrt> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mi>&amp;sigma;</mi> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <mi>k</mi> <mo>|</mo> <mo>&amp;dtri;</mo> <mi>&amp;rho;</mi> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>+</mo> <mn>1</mn> </mrow> </msqrt> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中,分别为密度曲面沿x和y方向的切向量,为曲面的单位法向量,k为调节参数,σ为光照强度因子。
4.根据权利要求1所述的一种超音速二维流场数据的可视化方法,其特征是:所述步骤c中,将插值后的RGB颜色分别与光照强度相乘就得到了合成后的颜色分布如式5所示,然后调节参数k,获得明暗适中的图像,具体涉及公式如下:
R=σr
G=σg
B=σb (5)
式中,R、G、B分别代表合成后颜色的红、绿、蓝分量。
CN201610850404.9A 2016-09-26 2016-09-26 一种超音速二维流场数据的可视化方法 Active CN107871337B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610850404.9A CN107871337B (zh) 2016-09-26 2016-09-26 一种超音速二维流场数据的可视化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610850404.9A CN107871337B (zh) 2016-09-26 2016-09-26 一种超音速二维流场数据的可视化方法

Publications (2)

Publication Number Publication Date
CN107871337A true CN107871337A (zh) 2018-04-03
CN107871337B CN107871337B (zh) 2020-12-08

Family

ID=61750962

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610850404.9A Active CN107871337B (zh) 2016-09-26 2016-09-26 一种超音速二维流场数据的可视化方法

Country Status (1)

Country Link
CN (1) CN107871337B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241572A (zh) * 2018-08-07 2019-01-18 北京空间技术研制试验中心 一种气动流动数值仿真流场结构的显示方法
CN110008387A (zh) * 2019-04-16 2019-07-12 成都四方伟业软件股份有限公司 流场可视化实现方法、装置及电子设备
CN112507282A (zh) * 2020-11-30 2021-03-16 中国航天空气动力技术研究院 一种基于速度梯度张量特性的流动显示方法
CN114119747A (zh) * 2021-11-23 2022-03-01 四川大学 一种基于pmd波前检测的三维流场流动显示方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101839920A (zh) * 2009-03-17 2010-09-22 英业达股份有限公司 流场可视化模组
CN102999936A (zh) * 2012-11-19 2013-03-27 北京中海新图科技有限公司 一种海洋流场数据的三维流线体绘制算法
CN103366045A (zh) * 2013-06-20 2013-10-23 华北水利水电大学 基于格子Boltzmann的流体可视化仿真方法
CN103413060A (zh) * 2013-08-26 2013-11-27 中国人民解放军国防科学技术大学 基于双控制体的格心网格数据三维激波特征定位方法
CN104867186A (zh) * 2015-04-29 2015-08-26 中国海洋大学 基于gpu的可交互海洋三维流场动态可视化算法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101839920A (zh) * 2009-03-17 2010-09-22 英业达股份有限公司 流场可视化模组
CN102999936A (zh) * 2012-11-19 2013-03-27 北京中海新图科技有限公司 一种海洋流场数据的三维流线体绘制算法
CN103366045A (zh) * 2013-06-20 2013-10-23 华北水利水电大学 基于格子Boltzmann的流体可视化仿真方法
CN103413060A (zh) * 2013-08-26 2013-11-27 中国人民解放军国防科学技术大学 基于双控制体的格心网格数据三维激波特征定位方法
CN104867186A (zh) * 2015-04-29 2015-08-26 中国海洋大学 基于gpu的可交互海洋三维流场动态可视化算法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
M. OTA等: "Quantitative density measurement of the lateral jet/cross-flow interaction field by colored-grid background oriented schlieren (CGBOS) technique", 《THE VISUALIZATION SOCIETY OF JAPAN 2015》 *
王盛波等: "二维流场可视化方法对比分析及综述", 《系统仿真学报》 *
马千里等: "基于两级采样的非结构化网格流场多激波特征可视化方法", 《计算机研究与发展》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241572A (zh) * 2018-08-07 2019-01-18 北京空间技术研制试验中心 一种气动流动数值仿真流场结构的显示方法
CN109241572B (zh) * 2018-08-07 2023-04-07 北京空间技术研制试验中心 一种气动流动数值仿真流场结构的显示方法
CN110008387A (zh) * 2019-04-16 2019-07-12 成都四方伟业软件股份有限公司 流场可视化实现方法、装置及电子设备
CN112507282A (zh) * 2020-11-30 2021-03-16 中国航天空气动力技术研究院 一种基于速度梯度张量特性的流动显示方法
CN112507282B (zh) * 2020-11-30 2023-07-28 中国航天空气动力技术研究院 一种基于速度梯度张量特性的流动显示方法
CN114119747A (zh) * 2021-11-23 2022-03-01 四川大学 一种基于pmd波前检测的三维流场流动显示方法

Also Published As

Publication number Publication date
CN107871337B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN107871337A (zh) 一种超音速二维流场数据的可视化方法
Born et al. Illustrative stream surfaces
Lee et al. View point evaluation and streamline filtering for flow visualization
Ma et al. Dynamic element textures
Ferstl et al. Interactive separating streak surfaces
CN102831275B (zh) 一种3d流体的仿真方法及系统
Tao et al. A deformation framework for focus+ context flow visualization
CN115495938B (zh) 海面波浪的三维动态模拟和可视化方法
Tušar et al. Visualizing 4D approximation sets of multiobjective optimizers with prosections
Bachthaler et al. Space-time Visualization of Dynamics in Lagrangian Coherent Structures of Time-dependent 2D Vector Fields.
Everts et al. Interactive illustrative line styles and line style transfer functions for flow visualization
Schroeder et al. Overview of Visualization.
CN103871096B (zh) 三维空间中真实感流体场景合成方法
Kuhn et al. Trajectory Density Projection for Vector Field Visualization.
Janke Mathematical structures for computer graphics
Everts et al. Illustrative line styles for flow visualization
Adams et al. Particle splatting: Interactive rendering of particle-based simulation data
Toledo et al. Visualizing 3D flow of black-oil reservoir models on arbitrary surfaces using projected 2D line integral convolution
Sadlo et al. Space-time flow visualization of dynamics in 2D Lagrangian coherent structures
Fu et al. Development and Application of Marine Environment Data Visualization Technology
Klein et al. Simulation, Modelling and Rendering of Incompressible Fluids in Real Time.
Präger et al. Visualizing Statistical Complexity in 3D Turbulent Flows using a Robust Entropy Calculation Method
CN116416409B (zh) 一种流体模拟粒子自适应分辨率表面重建方法及系统
Noguchi et al. Analytical computation and visualization of magnetic flux lines in 3-D space from hexahedral edge finite element results
Ikeda et al. Podiy: A design and production support system for pouch

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