CN107860800A - 一种用于泡状流的电容层析成像方法 - Google Patents

一种用于泡状流的电容层析成像方法 Download PDF

Info

Publication number
CN107860800A
CN107860800A CN201710860148.6A CN201710860148A CN107860800A CN 107860800 A CN107860800 A CN 107860800A CN 201710860148 A CN201710860148 A CN 201710860148A CN 107860800 A CN107860800 A CN 107860800A
Authority
CN
China
Prior art keywords
mrow
msubsup
mfrac
msub
capacitance
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.)
Withdrawn
Application number
CN201710860148.6A
Other languages
English (en)
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.)
Qingdao Zhongwu Special Ship Research And Design Institute Co ltd
Oceanographic Instrumentation Research Institute Shandong Academy of Sciences
Institute of Oceanographic Instrumentation Shandong Academy of Sciences
Original Assignee
Qingdao Zhongwu Special Ship Research And Design Institute Co ltd
Oceanographic Instrumentation Research Institute Shandong Academy of Sciences
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 Qingdao Zhongwu Special Ship Research And Design Institute Co ltd, Oceanographic Instrumentation Research Institute Shandong Academy of Sciences filed Critical Qingdao Zhongwu Special Ship Research And Design Institute Co ltd
Priority to CN201710860148.6A priority Critical patent/CN107860800A/zh
Publication of CN107860800A publication Critical patent/CN107860800A/zh
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/22Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance
    • G01N27/221Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance by investigating the dielectric properties

Abstract

本发明公开了一种用于泡状流的电容层析成像方法,包括如下步骤:电容的测量,二阶有限元模型的建立,灵敏度矩阵建立,数据归一化处理,权值计算以及反演计算,最后成像显示。本发明提供的用于泡状流的电容层析成像方法使用二阶有限元模型计算灵敏度矩阵,对灵敏度矩阵进行两次归一化处理,再根据单元在电场中的位置计算成像权函数,对介质进行加权成像。本发明主要适用于电容层析成像技术对泡状流的成像显示,针对现有的成像算法对圆管中心区域成像质量低的问题,改进灵敏度矩阵的求解与处理流程,使用加权函数对介质分布进行成像,提高泡状流成像质量。

Description

一种用于泡状流的电容层析成像方法
技术领域
本发明涉及一种电容层析成像方法,特别涉及一种用于泡状流的电容层析成像方法。
背景技术
电容层析成像技术是一种计算机层析成像技术,因其响应速度快、非侵入式、安全性能 好、成本低廉,在物料分布、油气输运、气力输送等工业领域有广泛应用。其主要原理为: 测量管道表面阵列电极间的电容值;应用数值算法计算管道内部的灵敏度分布;应用迭代或 非迭代算法计算管道内部介电常数分布,获得管道内介质分布和相含率等信息,从而获得介 质分布图像。
用有限元法计算灵敏度矩阵,根据测量电容使用LBP、Landweber、Tikhonov、牛顿法等 等反演算法计算管道内部介电常数分布。目前的电容层析成像算法对圆管中心物质的成像精 度低。目前对电容层析成像技术的研究主要集中于对反演算法的改进,圆管中心成像质量低 的问题得不到有效的解决。
灵敏度矩阵直接影响成像质量,但由于电场本身的“软场”特性,靠近电极的区域灵敏 度系数高,距离电极较远的圆管中心区域,灵敏度系数较低。因此圆管中心区域成像质量低, 难以识别物质分布,进而影响了管道内介质相含率的计算。
发明内容
为解决上述技术问题,本发明提供了一种用于泡状流的电容层析成像方法,以达到提高 圆管中心区域泡状流成像质量的目的。
为达到上述目的,本发明的技术方案如下:
一种用于泡状流的电容层析成像方法,包括如下步骤:
步骤1:测量电容层析成像传感器各电极对间的电容值,获得N个独立电容值;
步骤2:根据有限元方法将电容层析成像系统离散为二阶单元,根据泊松方程(1)建立 电容层析成像系统的二阶有限元模型(2),根据二阶有限元模型(2)计算电势分布φ,
▽·(ε▽φ)=0, (1)
Kφ=Q, (2)
K为系统刚度矩阵,Q为自由电荷量;
步骤3:根据式(3)获得N×M维灵敏度矩阵S;
Sij(e)为单元e对电极对ij的灵敏度值; 分别为激励电极i、j施加电压Vi、Vj时单 元e内的电势分布;
步骤4:根据式(4)对灵敏度矩阵行归一化,得到行归一化灵敏度矩阵根据式(5) 对进行列归一化,得到归一化处理的灵敏度矩阵
步骤5:依据单元与各电场中心线的加权距离d(x),单元与圆管中心的距离D(x),计算成 像权函数P(x)=f(d(x),D(x)),
de(x)为单元e的加权距离,dn(x)为单元e与电场中心线n的距离,s(p)为电场中心线在 垂足点p处的取值,根据灵敏度矩阵S的分布,拟合得到;
步骤6:对于ECT线性化模型C=Sg,将g的求解转化为求解最优化问题方程(7),采用 非迭代或迭代优化算法对方程(7)进行求解,获得介电常数分布g,
C为N个电容值组成的N×1维向量,g为M×1维的介电常数分布向量;
步骤7:根据步骤5中计算出的加权函数P(x),步骤6计算的介电常数分布g,进行图像重建,将介电常数分布g加权后,对圆管中介质分布进行可视化显示。
通过上述技术方案,本发明提供的用于泡状流的电容层析成像方法使用二阶有限元模型 计算灵敏度矩阵,对灵敏度矩阵进行两次归一化处理,再根据单元在电场中的位置计算成像 权函数P(x),对介质进行加权成像。本发明主要适用于电容层析成像技术对泡状流的成像显 示,针对现有的成像算法对圆管中心区域成像质量低的问题,改进灵敏度矩阵的求解与处理 流程,使用加权函数对介质分布进行成像,提高泡状流成像质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术 描述中所需要使用的附图作简单地介绍。
图1为本发明实施例所公开的一种用于泡状流的电容层析成像方法流程示意图;
图2为本发明实施例所公开的泡状流圆管模型示意图;
图3为图2中单元e与电场中心线距离示意图;
图4为本发明实施例所公开的实际成像效果与现有技术中常规成像效果对比图。
图中,1、电极;2、气泡;3、液相;e、单元;p、垂足;n、电场中心线。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描 述。
本发明提供了一种用于泡状流的电容层析成像方法,如图1所示,具体实施例如下:
步骤1:测量电容层析成像传感器各电极对间的电容值,获得N个独立电容值;
步骤2:根据有限元方法将电容层析成像系统离散为二阶单元,根据泊松方程(1)建立 电容层析成像系统的二阶有限元模型(2),根据二阶有限元模型(2)计算电势分布φ,
▽·(ε▽φ)=0, (1)
Kφ=Q, (2)
K为系统刚度矩阵,Q为自由电荷量;
步骤3:根据式(3)获得N×M维灵敏度矩阵S;
Sij(e)为单元e对电极对ij的灵敏度值; 分别为激励电极i、j施加电压Vi、Vj时单 元e内的电势分布;
步骤4:根据式(4)对灵敏度矩阵行归一化,得到行归一化灵敏度矩阵根据式(5) 对进行列归一化,得到归一化处理的灵敏度矩阵
步骤5:依据单元与各电场中心线的加权距离d(x),单元与圆管中心的距离D(x),计算成 像权函数P(x)=f(d(x),D(x)),
de(x)为单元e的加权距离,dn(x)为单元e与电场中心线n的距离,s(p)为电场中心线在 垂足点p处的取值,根据灵敏度矩阵S的分布,拟合得到;泡状流圆管模型示意图见图2和 图3所示,圆管最外侧为电极。
步骤6:对于ECT线性化模型C=Sg,将g的求解转化为求解最优化问题方程(7),采用 非迭代或迭代优化算法对方程(7)进行求解,获得介电常数分布g,
C为N个电容值组成的N×1维向量,g为M×1维的介电常数分布向量;
步骤7:根据步骤5中计算出的加权函数P(x),步骤6计算的介电常数分布g,进行图像重建,将介电常数分布g加权后,对圆管中介质分布进行可视化显示。
按照上述步骤对多个泡状流进行成像,成像效果见图4中第三列,对应的现有的常规成 像效果分别见图4中第二列,对应的实际气泡分布示意图分别为图4中第一列所示。
从图中可以看出,本发明所公开的成像方法与现有的成像方法相比,明显提高了成像精 度,尤其是圆管中心处的成像精度,本发明有助于拓展电容层析成像技术在工业多相流测量 领域的应用。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些 实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理 可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被 限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的 范围。

Claims (1)

1.一种用于泡状流的电容层析成像方法,其特征在于,包括如下步骤:
步骤1:测量电容层析成像传感器各电极对间的电容值,获得N个独立电容值;
步骤2:根据有限元方法将电容层析成像系统离散为二阶单元,根据泊松方程(1)建立电容层析成像系统的二阶有限元模型(2),根据二阶有限元模型(2)计算电势分布φ,
<mrow> <mo>&amp;dtri;</mo> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>&amp;epsiv;</mi> <mo>&amp;dtri;</mo> <mi>&amp;phi;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
Kφ=Q, (2)
K为系统刚度矩阵,Q为自由电荷量;
步骤3:根据式(3)获得N×M维灵敏度矩阵S;
<mrow> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>V</mi> <mi>i</mi> </msub> <msub> <mi>V</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>&amp;Integral;</mo> <msub> <mo>&amp;Integral;</mo> <mi>e</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>&amp;phi;</mi> <mi>i</mi> <mi>e</mi> </msubsup> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>&amp;phi;</mi> <mi>j</mi> <mi>e</mi> </msubsup> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>&amp;phi;</mi> <mi>i</mi> <mi>e</mi> </msubsup> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>&amp;phi;</mi> <mi>j</mi> <mi>e</mi> </msubsup> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
Sij(e)为单元e对电极对ij的灵敏度值;分别为激励电极i、j施加电压Vi、Vj时单元e内的电势分布;
步骤4:根据式(4)对灵敏度矩阵行归一化,得到行归一化灵敏度矩阵根据式(5)对进行列归一化,得到归一化处理的灵敏度矩阵
<mrow> <msubsup> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mo>#</mo> </msubsup> <mo>=</mo> <mfrac> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <msub> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mo>*</mo> </msubsup> <mo>=</mo> <mfrac> <msubsup> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mo>#</mo> </msubsup> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msubsup> <mi>S</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mo>#</mo> </msubsup> <mrow> <mo>(</mo> <mi>e</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
步骤5:依据单元与各电场中心线的加权距离d(x),单元与圆管中心的距离D(x),计算成像权函数P(x)=f(d(x),D(x)),
<mrow> <msub> <mi>d</mi> <mi>e</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>d</mi> <mi>n</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mi>s</mi> <mrow> <mo>(</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
de(x)为单元e的加权距离,dn(x)为单元e与电场中心线n的距离,s(p)为电场中心线在垂足点p处的取值,根据灵敏度矩阵S的分布拟合得到;
步骤6:对于ECT线性化模型C=Sg,将g的求解转化为求解最优化问题方程(7),采用非迭代或迭代优化算法对方程(7)进行求解,获得介电常数分布g,
<mrow> <mi>min</mi> <mi> </mi> <mi>J</mi> <mrow> <mo>(</mo> <mi>g</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mo>|</mo> <mo>|</mo> <mi>S</mi> <mi>g</mi> <mo>-</mo> <mi>C</mi> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
C为N个电容值组成的N×1维向量,g为M×1维的介电常数分布向量;
步骤7:根据步骤5中计算出的加权函数P(x),步骤6计算的介电常数分布g,进行图像重建,将介电常数分布g加权后,对圆管中介质分布进行可视化显示。
CN201710860148.6A 2017-09-21 2017-09-21 一种用于泡状流的电容层析成像方法 Withdrawn CN107860800A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710860148.6A CN107860800A (zh) 2017-09-21 2017-09-21 一种用于泡状流的电容层析成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710860148.6A CN107860800A (zh) 2017-09-21 2017-09-21 一种用于泡状流的电容层析成像方法

Publications (1)

Publication Number Publication Date
CN107860800A true CN107860800A (zh) 2018-03-30

Family

ID=61698235

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710860148.6A Withdrawn CN107860800A (zh) 2017-09-21 2017-09-21 一种用于泡状流的电容层析成像方法

Country Status (1)

Country Link
CN (1) CN107860800A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110186962A (zh) * 2019-05-10 2019-08-30 天津大学 一种用于电容层析成像的不完整测量数据成像方法
CN111999346A (zh) * 2020-08-26 2020-11-27 山东省科学院海洋仪器仪表研究所 一种基于敏感场灵敏度矩阵的电容层析成像方法
CN113724352A (zh) * 2021-08-24 2021-11-30 沈阳工业大学 一种基于降维处理的电容层析成像图像重建方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101520478A (zh) * 2009-03-13 2009-09-02 北京航空航天大学 一种基于圆形传感器的电容层析成像的图像重建直接方法
CN102353703A (zh) * 2011-07-05 2012-02-15 东南大学 一种电容层析成像传感器布置方法
CN102620855A (zh) * 2012-03-30 2012-08-01 华南理工大学 一种基于电容层析成像的聚合物熔体温度场测量方法及系统
CN104535294A (zh) * 2015-01-16 2015-04-22 天津大学 一种基于二阶微分的修正l曲线电学层析成像重建方法
CN106153691A (zh) * 2016-06-21 2016-11-23 中国科学院工程热物理研究所 利用耦合微波成像信息的敏感场的电容层析成像方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101520478A (zh) * 2009-03-13 2009-09-02 北京航空航天大学 一种基于圆形传感器的电容层析成像的图像重建直接方法
CN102353703A (zh) * 2011-07-05 2012-02-15 东南大学 一种电容层析成像传感器布置方法
CN102620855A (zh) * 2012-03-30 2012-08-01 华南理工大学 一种基于电容层析成像的聚合物熔体温度场测量方法及系统
CN104535294A (zh) * 2015-01-16 2015-04-22 天津大学 一种基于二阶微分的修正l曲线电学层析成像重建方法
CN106153691A (zh) * 2016-06-21 2016-11-23 中国科学院工程热物理研究所 利用耦合微波成像信息的敏感场的电容层析成像方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
刘靖 等: "多相流测量用动态电容层析成像图像重建算法", 《仪器仪表学报》 *
张立峰 等: "一种新的电容层析成像电容测量值归一化模型", 《中国电机工程学报》 *
杜乙 等: "一种用于 CT偏置扫描重建的指数型加权函数", 《清华大学学报(自然科学版)》 *
王伊凡: "直接三维ECT传感器及重建算法的研究", 《中国博士学位论文全文数据库 信息科技辑》 *
谢代梁 等: "电容层析成像流型可视化系统研究", 《浙江大学学报(工学版)》 *
邵兴华: "电容层析成像图像重建算法的研究", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 *
邵晓寅: "电容层析成像技术在气液两相流检测中的应用研究", 《中国优秀博硕士学位论文全文数据库 (硕士) 信息科技辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110186962A (zh) * 2019-05-10 2019-08-30 天津大学 一种用于电容层析成像的不完整测量数据成像方法
CN111999346A (zh) * 2020-08-26 2020-11-27 山东省科学院海洋仪器仪表研究所 一种基于敏感场灵敏度矩阵的电容层析成像方法
CN113724352A (zh) * 2021-08-24 2021-11-30 沈阳工业大学 一种基于降维处理的电容层析成像图像重建方法
CN113724352B (zh) * 2021-08-24 2023-07-04 沈阳工业大学 一种基于降维处理的电容层析成像图像重建方法

Similar Documents

Publication Publication Date Title
Tehrani et al. L1 regularization method in electrical impedance tomography by using the L1-curve (Pareto frontier curve)
CN107860800A (zh) 一种用于泡状流的电容层析成像方法
CN102620855B (zh) 一种基于电容层析成像的聚合物熔体温度场测量方法及系统
Blake et al. Using the topology of large-scale structure in the WiggleZ Dark Energy Survey as a cosmological standard ruler
Afshar et al. A node enrichment adaptive refinement in Discrete Least Squares Meshless method for solution of elasticity problems
US10488236B2 (en) Velocity vector field mapping using electrical capacitance sensors
CN108376124B (zh) 一种用于电学成像的多导体体系导纳矩阵快速计算方法
Wajman et al. Tunnel-based method of sensitivity matrix calculation for 3D-ECT imaging
Ye et al. Resolution analysis using fully 3D electrical capacitive tomography
Cao et al. Direct recovery of the electrical admittivities in 2D electrical tomography by using Calderon's method and two-terminal/electrode excitation strategies
Li et al. A practical and effective regularized polynomial smoothing (RPS) method for high-gradient strain field measurement in digital image correlation
Xie et al. Preliminary evaluation of cryogenic two-phase flow imaging using electrical capacitance tomography
Kortschak et al. A FEM‐BEM approach using level‐sets in electrical capacitance tomography
CN106153691A (zh) 利用耦合微波成像信息的敏感场的电容层析成像方法
Ye et al. Evaluation of electrical capacitance tomography sensor based on the coupling of fluid field and electrostatic field
Zhou et al. Analytical integral algorithm in the BEM for orthotropic potential problems of thin bodies
Wang et al. Application of particle filtering algorithm in image reconstruction of EMT
Tian et al. Experimental imaging and algorithm optimization based on deep neural network for electrical capacitance tomography for LN2-VN2 flow
Xu et al. Iterative nonlinear Tikhonov algorithm with constraints for electromagnetic tomography
Lanying et al. A novel multiple-electrodes excitation method for electrical capacitance tomography system
Brański et al. A comparison of boundary methods based on inverse variational formulation
Guo et al. Sensitivity matrix construction for electrical capacitance tomography based on the difference model
Sun A meshless method based on the method of fundamental solution for solving the steady-state heat conduction problems
CN111158059A (zh) 基于三次b样条函数的重力反演方法
Zhang et al. Numerical simulation of forward problem for electrical capacitance tomography using element-free Galerkin method

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 266200, Qingdao, Shandong, Qingdao, Qingdao, the core of the blue Silicon Valley, blue Silicon Valley business center, phase one, building No. 1.

Applicant after: Inst. of Marine Apparatus & Instruments, Shandong Prov. Academy of Sciences

Applicant after: Qingdao Special Ship Research & Design Institute Co., Ltd.

Address before: Miao road Laoshan District 266071 of Shandong Province, Qingdao City, No. 37

Applicant before: Inst. of Marine Apparatus & Instruments, Shandong Prov. Academy of Sciences

Applicant before: Qingdao Special Ship Research & Design Institute Co., Ltd.

WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20180330