CN108051472A - A kind of method for fast measuring of material at high temperature thermal physical property parameter - Google Patents
A kind of method for fast measuring of material at high temperature thermal physical property parameter Download PDFInfo
- Publication number
- CN108051472A CN108051472A CN201711264303.4A CN201711264303A CN108051472A CN 108051472 A CN108051472 A CN 108051472A CN 201711264303 A CN201711264303 A CN 201711264303A CN 108051472 A CN108051472 A CN 108051472A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- physical property
- thermal
- mfrac
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/18—Investigating or analyzing materials by the use of thermal means by investigating thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
Abstract
The invention discloses a kind of method for fast measuring of material at high temperature thermal physical property parameter：According to material thermal physical property parameter medium temperature transonic characteristic, using ultrasonic pulse-echo method, ultrasound propagation time under the conditions of acquisition Transient Heat Transfer, it, can thermal physical property parameter quick, lossless, that non-contactly measurement material varies with temperature by the material parameter in the ultrasonic echo characteristic inverting equation of heat conduction；The method of the present invention only measures once, such as test specimen heating surface carries out being warming up to such as 400 DEG C of predetermined temperature value, the multiple materials thermal physical property parameters such as thermal conductivity factor, the specific heat perhaps thermal diffusion coefficient under room temperature to 400 DEG C of different temperatures can be obtained, there is outstanding advantages of measuring speed is fast, at low cost, versatility is good, measurement range is big.
Description
Technical field
The present invention relates to ultrasound detection field, a kind of quick measurement side of material at high temperature thermal physical property parameter is related in particular to
Method.
Background technology
Thermal physical property parameter such as thermal conductivity factor, specific heat perhaps thermal diffusion coefficient etc. is the key character parameter of material, is both
Can material be weighed adapt to the key parameter of specific thermal process need of work, and be to carry out basic research and work to specific thermal process
The basic data of journey design.The quick and precisely measurement of the hot Physical Parameters of material is in the high-new skill such as aerospace, new material, microelectronics
The industrial circles such as art field and building, chemical industry, biology, medicine, metallurgy, electric power all have apparent scientific meaning and important
Engineering application value.
Steady state method and cold store enclosure are with their own characteristics in terms of material thermal physical property parameter measurement, complement one another.However it is measuring
The thermal conductivity factor that varies with temperature, specific heat perhaps thermal diffusion coefficient when, most of measuring methods often measure the parameter of a certain temperature
It needs sample being heated to the given temperature range, there are the cycle is long, control is complicated and expends the problems such as big；It is in addition, most of
Method, which also exists, is only capable of measurement one-dimensional problem or the problems such as be only applicable to measure a certain thermal physical property parameter, the versatility of method compared with
Difference；Simultaneously as sensor usually requires to contact with test specimen, thermophysical property measurement scope can be subject to sensor high temperature resistance
Limitation.
The content of the invention
The object of the present invention is to provide a kind of method for fast measuring of material at high temperature thermal physical property parameter, using ultrasonic echo
Method, according to ultrasound propagation time under the conditions of Transient Heat Transfer, based on parameter identification method, suitable for metal under condition of different temperatures,
The nondestructive measurement of the thermal physical property parameters such as composite material material thermal conductivity factor, specific heat capacity.
To achieve the above object, the present invention adopts the following technical scheme that：
Step 1：Based on calibration experiment, the relational expression V-T of ultrasonic propagation velocity V and temperature T in test specimen is obtained.
Step 2：Test specimen is heated, by common ultrasound method, obtains testee tiMoment
Ultrasonic propagation time ti,exp。
Step 3：The measures conversion of high temperature thermal physical property parameter is asked to solve the optimization of heat conduction problem border unknown parameter
The ultrasound propagation time t that topic, wherein optimization aim obtain for measurementi,expWith the ultrasound propagation time t of numerical solutioni,calDifference
Minimum, restrained boundary are heat transfer governing equation and its relevant primary condition.
The object function of optimization problem is：
In formula, k (T), Cp(T) and α (T) is respectively measured material Equivalent Thermal Conductivities, specific heat capacity and thermal expansion to be identified
Dissipate coefficient； ti,calThe t obtained for numerical computationsiThe ultrasonic propagation time at moment, the time of measuring ordinal number that subscript i is represented, n
Represent total time of measuring points；L is the length that test specimen is tested direction.
The constraints of optimization problem is:
T(x,t)|T=0=T0
In formula, T0For initial temperature field.
Step 4：Solve heat conduction inverse problem.Correlation solve the step of be：
(1) given parameters initial value；
(2) value of temperature field T (x, t) and object function J is obtained in numerical solution state equation；
(3) numerical solution sensitivity equation draws sensitivity vectors；
(4) parameter value is optimized with adjoint equation algorithm, obtains k (T) or Cp(T) or α (T)；
(5) judge whether to restrain (taking ε≤1e-6), if convergence, stops calculating；Otherwise return to step (2) iteration,
Until reaching convergence criterion；
(6) measured material Equivalent Thermal Conductivities k (T), specific heat capacity C are obtainedp(T) or thermal diffusion coefficient α (T).
In conclusion by adopting the above-described technical solution, the beneficial effects of the invention are as follows：
1st, this method only measures once, and test specimen heating surface carries out being warming up to such as 400 DEG C of predetermined temperature value, you can obtains
Material thermal physical property parameter under room temperature to 400 DEG C of different temperatures, has that measuring speed is fast, the advantages such as at low cost, particularly measures
High temperature, hyperthermal material thermal physical property parameter when this method measurement period it is short the advantages of it is more prominent；
When the 2nd, carrying out non-cpntact measurement based on electromagnetism or laser-ultrasound, material at high temperature thermophysical property measurement is hardly by sensor
The influence of heat resistance, the advantage for having measurement range big；
3rd, method versatility is good, suitable for measuring the heat such as the thermal conductivity factor, the specific heat perhaps thermal diffusion coefficient that vary with temperature
Physical parameter is also applied for the measurement of one-dimensional, two-dimentional or three-dimensional material high temperature thermal physical property parameter.
Description of the drawings
Examples of the present invention will be described by way of reference to the accompanying drawings, wherein：
The prediction result that Fig. 1 thermal conductivity factors vary with temperature；
The prediction result that Fig. 2 thermal capacity varies with temperature；
The prediction result that Fig. 3 thermal diffusion coefficients vary with temperature.
Specific embodiment
Below in conjunction with the accompanying drawings and pass through specific embodiment that material at high temperature thermal physical property parameter measures to the present invention to further
Be described in detail, following embodiment be it is descriptive, it is non-limiting, it is impossible to which that protection scope of the present invention is limited with this.
Technique according to the invention scheme and step carry out the implementation of concrete case, as follows：
The bottom of measured material is subject to constant value Uniform Heat to heat Q=1e+5W/m2, remaining surface is all adiabatic face, ultrasound
Probe is placed in upper surface, using vertical incidence mode excitation pulse ultrasonic wave, based on the variation in measurement Bottom echo propagation time,
Indirect problem, the parameters such as thermal conductivity factor, thermal capacity or the thermal diffusion coefficient of inverting material are coupled by solving thermal acoustic.Measured material
Room temperature is to thermal capacity true value ρ c=1e+6 × [1.280+0.00526 × T- (0.363e-5) × T in the range of 410 DEG C2], heat conduction
Coefficient true value k=3.195+0.00756 × T- (0.813e-5) × T2, wherein T is temperature.Above-mentioned material thermal physical property parameter by
Experiment is obtained ahead of time, and the thermal conductivity factor of material and thermal capacity etc. can be expressed as the function of temperature, but usual in practice in engineering
There is no any priori to the thermal physical property parameter of material in advance, therefore be expressed as material thermal physical property parameter in heat transfer model
Go out the function with position and the piecewise function of time change, and by the measurement result in material internal temperature field come direct inversion.
The measurement identification knot for varying with temperature thermal conductivity factor, thermal capacity and thermal diffusion coefficient is set forth in Fig. 1 to Fig. 3
Fruit.Wherein the resultant error of thermal capacity and thermal diffusion coefficient is within 1%, and thermal conductivity factor is about 2%.
The invention is not limited in foregoing specific embodiments.The present invention, which expands to, any in the present specification to be disclosed
The step of new feature or any new combination and any new method or process disclosed or any new combination.
Claims (1)
1. a kind of method for fast measuring of material at high temperature thermal physical property parameter, it is characterised in that comprise the following steps：
Step 1：Based on calibration experiment, the relational expression V-T of ultrasonic propagation velocity V and temperature T in test specimen is obtained；
Step 2：Test specimen is heated, by common ultrasound method, obtains testee tiMoment surpasses
Acoustic transit time ti,exp；
Step 3：By the measures conversion of high temperature thermal physical property parameter to solve the optimization problem of unknown parameter in heat conduction problem,
The ultrasound propagation time t that middle optimization aim obtains for measurementi,expWith the ultrasound propagation time t of numerical solutioni,calDifference is minimum,
Restrained boundary is heat transfer governing equation and its relevant primary condition；
The object function of optimization problem is：
<mrow>
<mi>J</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
<mo>,</mo>
<mi>&alpha;</mi>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>{</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>exp</mi>
</mrow>
</msub>
<mo>}</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>{</mo>
<mn>2</mn>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>L</mi>
</msubsup>
<mfrac>
<mn>1</mn>
<mrow>
<mi>V</mi>
<mo>&lsqb;</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
<mi>d</mi>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>exp</mi>
</mrow>
</msub>
<mo>}</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
In formula, k (T), Cp(T) and α (T) is respectively measured material Equivalent Thermal Conductivities to be identified, specific heat capacity and thermal diffusion system
Number；ti,calThe t obtained for numerical computationsiThe ultrasonic propagation time at moment, the time of measuring ordinal number that subscript i is represented, n represent total
Time of measuring points；L is the length that test specimen is tested direction；
The constraints of optimization problem is:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>S</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>=</mo>
<msub>
<mi>&rho;C</mi>
<mi>p</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
T(x,t)|T=0=T0
In formula, T0For initial temperature field；
Step 4：Solve heat conduction inverse problem.Correlation solve the step of be：
(1) given parameters initial value；
(2) value of temperature field T (x, t) and object function J is obtained in numerical solution state equation；
(3) numerical solution sensitivity equation draws sensitivity vectors；
(4) parameter value is optimized with adjoint equation algorithm, obtains k (T) or Cp(T) or α (T)；
(5) judge whether to restrain (taking ε≤1e-6), if convergence, stops calculating；Otherwise (2) iteration is returned to, until reaching
Convergence criterion；
(6) measured material Equivalent Thermal Conductivities k (T), specific heat capacity C are obtainedp(T) or thermal diffusion coefficient α (T).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711264303.4A CN108051472B (en) | 2017-12-05 | 2017-12-05 | Method for rapidly measuring high-temperature thermophysical property parameters of material |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711264303.4A CN108051472B (en) | 2017-12-05 | 2017-12-05 | Method for rapidly measuring high-temperature thermophysical property parameters of material |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108051472A true CN108051472A (en) | 2018-05-18 |
| CN108051472B CN108051472B (en) | 2020-08-14 |
Family
ID=62121391
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711264303.4A Active CN108051472B (en) | 2017-12-05 | 2017-12-05 | Method for rapidly measuring high-temperature thermophysical property parameters of material |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108051472B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109324079A (en) * | 2018-10-23 | 2019-02-12 | 中国特种设备检测研究院 | A kind of measurement method of the material thermal expansion coefficient based on ultrasound |
| CN109470772A (en) * | 2018-10-23 | 2019-03-15 | 中国特种设备检测研究院 | A kind of non-destructive measuring method of the internal heat resource intensity size and location based on ultrasound |
| CN113239537A (en) * | 2021-05-11 | 2021-08-10 | 扬州扬杰电子科技股份有限公司 | Method for determining high-frequency heating power of diffusion sheet alloy block |
| CN116013443A (en) * | 2023-03-22 | 2023-04-25 | 中国空气动力研究与发展中心计算空气动力研究所 | Method, device, equipment and readable storage medium for predicting heat transfer characteristic |
| CN116013442A (en) * | 2023-03-22 | 2023-04-25 | 中国空气动力研究与发展中心计算空气动力研究所 | Method, device and equipment for designing heat-resistant material and readable storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010091015A1 (en) * | 2009-02-03 | 2010-08-12 | Northrop Grumman Guidance and Electronics Company Inc. | Systems and methods for measuring at least one thermal property of materials based on a thermal brewster angle |
| US20130058448A1 (en) * | 2011-09-06 | 2013-03-07 | Battelle Energy Alliance, Llc | Vibro-acoustic sensors for materials characterization and related methods and systems |
| CN103149233A (en) * | 2013-01-28 | 2013-06-12 | 中国科学院工程热物理研究所 | Device and method of testing material thermophysical property parameter |
| CN104597078A (en) * | 2015-01-14 | 2015-05-06 | 北京科技大学 | Method for measuring anisotropic material heat conductivity based on small-plane heat source |
| CN104596667A (en) * | 2015-01-05 | 2015-05-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for detecting sensitivity of transient non-uniform temperature field in object by using ultrasonic waves |
| CN204374119U (en) * | 2014-10-29 | 2015-06-03 | 北京航天长征飞行器研究所 | Solid dynamic heat transfer coefficient measuring device and measuring system |
| CN106053596A (en) * | 2016-05-18 | 2016-10-26 | 陕西煤业化工技术研究院有限责任公司 | Nano fluid heat and mass transfer monitoring device and method |
-
2017
- 2017-12-05 CN CN201711264303.4A patent/CN108051472B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010091015A1 (en) * | 2009-02-03 | 2010-08-12 | Northrop Grumman Guidance and Electronics Company Inc. | Systems and methods for measuring at least one thermal property of materials based on a thermal brewster angle |
| US20130058448A1 (en) * | 2011-09-06 | 2013-03-07 | Battelle Energy Alliance, Llc | Vibro-acoustic sensors for materials characterization and related methods and systems |
| CN103149233A (en) * | 2013-01-28 | 2013-06-12 | 中国科学院工程热物理研究所 | Device and method of testing material thermophysical property parameter |
| CN204374119U (en) * | 2014-10-29 | 2015-06-03 | 北京航天长征飞行器研究所 | Solid dynamic heat transfer coefficient measuring device and measuring system |
| CN104596667A (en) * | 2015-01-05 | 2015-05-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for detecting sensitivity of transient non-uniform temperature field in object by using ultrasonic waves |
| CN104597078A (en) * | 2015-01-14 | 2015-05-06 | 北京科技大学 | Method for measuring anisotropic material heat conductivity based on small-plane heat source |
| CN106053596A (en) * | 2016-05-18 | 2016-10-26 | 陕西煤业化工技术研究院有限责任公司 | Nano fluid heat and mass transfer monitoring device and method |
Non-Patent Citations (2)
| Title |
|---|
| KEN KUROSAKI 等: "Thermophysical properties of SrY2O4", 《JOURNAL OF ALLOYS AND COMPOUNDS》 * |
| 魏东 等: "结构内部非均匀温度场的重建方法研究", 《工程热物理学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109324079A (en) * | 2018-10-23 | 2019-02-12 | 中国特种设备检测研究院 | A kind of measurement method of the material thermal expansion coefficient based on ultrasound |
| CN109470772A (en) * | 2018-10-23 | 2019-03-15 | 中国特种设备检测研究院 | A kind of non-destructive measuring method of the internal heat resource intensity size and location based on ultrasound |
| CN113239537A (en) * | 2021-05-11 | 2021-08-10 | 扬州扬杰电子科技股份有限公司 | Method for determining high-frequency heating power of diffusion sheet alloy block |
| CN113239537B (en) * | 2021-05-11 | 2023-10-27 | 扬州扬杰电子科技股份有限公司 | High-frequency heating power determination method for diffusion sheet alloy block |
| CN116013443A (en) * | 2023-03-22 | 2023-04-25 | 中国空气动力研究与发展中心计算空气动力研究所 | Method, device, equipment and readable storage medium for predicting heat transfer characteristic |
| CN116013442A (en) * | 2023-03-22 | 2023-04-25 | 中国空气动力研究与发展中心计算空气动力研究所 | Method, device and equipment for designing heat-resistant material and readable storage medium |
| CN116013443B (en) * | 2023-03-22 | 2023-05-26 | 中国空气动力研究与发展中心计算空气动力研究所 | Method, device, equipment and readable storage medium for predicting heat transfer characteristic |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108051472B (en) | 2020-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108051472A (en) | A kind of method for fast measuring of material at high temperature thermal physical property parameter | |
| CN101126729B (en) | Double heat flux gauge steady state method for measuring material heat conductivity | |
| CN104792435B (en) | The method for reconstructing of inside configuration non-uniform temperature field based on transient state thermal boundary inverting | |
| CN104596667B (en) | The sensitivity method of ultrasonic listening interior of articles transient state non-uniform temperature field | |
| Lu et al. | Inverse estimation of the inner wall temperature fluctuations in a pipe elbow | |
| CN108008022B (en) | Ultrasonic wave propagation speed measuring method along with temperature change | |
| CN105628790B (en) | A kind of inside configuration temperature field measurement method based on material property Parameters variation | |
| Hubble et al. | A hybrid method for measuring heat flux | |
| CN104597078A (en) | Method for measuring anisotropic material heat conductivity based on small-plane heat source | |
| CN107748205A (en) | A kind of elastic constant measurement method varied with temperature | |
| CN108051475A (en) | A kind of method for fast measuring of convection transfer rate | |
| Haghighi et al. | Three-dimensional inverse transient heat transfer analysis of thick functionally graded plates | |
| CN109324079A (en) | A kind of measurement method of the material thermal expansion coefficient based on ultrasound | |
| CN107966472A (en) | A kind of lossless method for fast measuring of high temperature contact thermal resistance | |
| Tuladhar et al. | Thermal conductivity of whey protein films undergoing swelling: Measurement by dynamic gauging | |
| Guo et al. | Inverse heat conduction estimation of inner wall temperature fluctuations under turbulent penetration | |
| Haghighi et al. | Inverse internal pressure estimation of functionally graded cylindrical shells under thermal environment | |
| CN109470772A (en) | A kind of non-destructive measuring method of the internal heat resource intensity size and location based on ultrasound | |
| Pavlík et al. | Experimental assessment of thermal conductivity of a brick block with internal cavities using a semi-scale experiment | |
| CN105403323B (en) | A kind of inside configuration temperature field measurement method based on phase-detection | |
| CN109506806B (en) | Method for simultaneously measuring internal temperature and thickness of high-temperature structure under transient condition | |
| CN109506807B (en) | Method for simultaneously measuring internal temperature and wall thickness of high-temperature structure under steady-state condition | |
| Alifanov et al. | Complex identification of thermophysical properties of anisotropic composite material | |
| RU2613194C1 (en) | Method of measuring thermal properties of anisotropic materials by linear impulse of heat source | |
| Myers et al. | An adaptive extended kalman filter incorporating state model uncertainty for localizing a high heat flux point source using an ultrasonic sensor array |
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 |