Caractérisation dans le plan de matériaux isolants par thermographie infrarouge et méthode convolutive

Similar documents
Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

CHAPTER X PHASE-CHANGE PROBLEMS

Solution for singularly perturbed problems via cubic spline in tension

2 Finite difference basics

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

Interconnect Modeling

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Finite Element Method: A Short Introduction

Electrostatic Potential from Transmembrane Currents

Formal solvers of the RT equation

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Numerical Transient Heat Conduction Experiment

Multivariate Ratio Estimator of the Population Total under Stratified Random Sampling

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Implementation of the Matrix Method

A new Approach for Solving Linear Ordinary Differential Equations

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Thermal tomography on the basis of an information method

The Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially

Stanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

The Analysis of Convection Experiment

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

risk and uncertainty assessment

A Spline based computational simulations for solving selfadjoint singularly perturbed two-point boundary value problems

Model and program for the prediction of the indoor air temperature and the indoor air relative humidity

ECE 107: Electromagnetism

Diffusion Mass Transfer

8.022 (E&M) Lecture 4

PubH 7405: REGRESSION ANALYSIS. SLR: INFERENCES, Part II

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline

Supplemental Material: Causal Entropic Forces

Introduction to Antennas & Arrays

Analysis of Discrete Time Queues (Section 4.6)

Implementation of the Matrix Method

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Binomial Distribution: Tossing a coin m times. p = probability of having head from a trial. y = # of having heads from n trials (y = 0, 1,..., m).

Statistics for Business and Economics

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Digital Signal Processing

16 Reflection and transmission, TE mode

6.3.4 Modified Euler s method of integration

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).

A New Recursive Method for Solving State Equations Using Taylor Series

b ), which stands for uniform distribution on the interval a x< b. = 0 elsewhere

The Feynman path integral

π e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m

Electrical double layer: revisit based on boundary conditions

Digital Modems. Lecture 2

Probability Theory. The nth coefficient of the Taylor series of f(k), expanded around k = 0, gives the nth moment of x as ( ik) n n!

Physics 181. Particle Systems

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Quasi-Static transient Thermal Stresses in a Robin's thin Rectangular plate with internal moving heat source

Scientific Research of the Institute of Mathematics and Computer Science 1(11) 2012, 23-30

Implicit Integration Henyey Method

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Simulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM

Mean Field / Variational Approximations

Adaptive Kernel Estimation of the Conditional Quantiles

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

, rst we solve te PDE's L ad L ad n g g (x) = ; = ; ; ; n () (x) = () Ten, we nd te uncton (x), te lnearzng eedbac and coordnates transormaton are gve

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

In this section is given an overview of the common elasticity models.

THE ROYAL STATISTICAL SOCIETY 2006 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE

Chapter 3. r r. Position, Velocity, and Acceleration Revisited

4.2 Chemical Driving Force

Lab session: numerical simulations of sponateous polarization

PHYS 705: Classical Mechanics. Calculus of Variations II

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

T E C O L O T E R E S E A R C H, I N C.

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

11. Dynamics in Rotating Frames of Reference

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics

Towards the modeling of microgalvanic coupling in aluminum alloys : the choice of boundary conditions

The Finite Element Method

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Global Sensitivity. Tuesday 20 th February, 2018

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

a. (All your answers should be in the letter!

Modelli Clamfim Equazioni differenziali 22 settembre 2016

On a nonlinear compactness lemma in L p (0, T ; B).

Simulation of Turbulent Flow Using FEM

Some modelling aspects for the Matlab implementation of MMA

Video Layer Extraction and Reconstruction

How Differential Equations Arise. Newton s Second Law of Motion

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

6 Supplementary Materials

Numerical Heat and Mass Transfer

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

Transcription:

Journée SFT du Jeud 3 Jun 3 Caractérsaton Température / Température Caractérsaton dans le plan de matérau solants par termogrape nfrarouge et métode convolutve B. REMY, A. DEGIOVAI & D. MAIET B.Rémy, A.Degovann and D.Mallet " Measurement of te In-Plane Termal Dffusvty of Materals by Infrared Termograpy " (DOI:.7/s765-5-45-z Internatonal Journal of Termopyscs - Vol.6, n (5, pp 493-55. EMTA - UMR 7563 (CRS Unversté de orrane ESEM avenue de la Forêt de Haye TSA 664-5458 VADOEUVRE CEDEX Emal : benjamn.remy@unv-lorrane.fr Journée SFT du Jeud 3 Jun 3

Te n-plane dffusvty measurement tecnques Soluton n te case of a sem-nfnte medum submtted to a eat step stmulaton : φ e ( t T, Assumptons: - Sem-nfnte medum - One drectonal eat transfer - Insulated T (, t φ at ep erfc λ π 4at 4at 4at a: termal dffusvty λ: termal conductvty Metod : poston and dfferent tmes (Steere, 966 & Harmaty, 964 T T (,t (, t erfc erfc at at Anoter tecnque: at e,5τ t m m m ( τ t ln( τ t m Steere, 967 Square pulse (duraton τ

Te n-plane dffusvty measurement tecnques Metod : dfferent postons at a gven tme (Katayama, 969 T πe ep erfc T 4at 4at 4at Drawbacks: - Heat Flu Dstrbuton must be perfectly known - Heatng element s deal (no capacty and n perfect contact wt te sample - Medum s assumed to be perfectly nsulated - D Heat transfer ( pt Metod : aplace transform (Kavanpour & Beck, 977 θ ( p ( T ( t T ( t ep dt θ T ( p a T t θ a p θ ( p ( p a ϕ ep λ p a, ln θ θ (, p (, p ( Te soluton s ndependent of te n-tme eat flu dstrbuton p a

Te n-plane dffusvty measurement tecnques Metod : Fn s metod (Hadsaroyo, 99 Unknown eat flu (, t T ( t T, ( e + S e. T S λ T t θ p a S λ ( T T + θ et a e l ln θ (, p (, p θ p S ( + a λ Unknown eat flu T (, t T (, t e Bot dffusvty a and eat losses are taken nto account (sem-nfnte medum s assumed l Interest of te nfrared camera: - on-necessary n space unform stmulaton - D transfer n te case of nsulatng materals

O J D D J J A J O Te Fn s metod n Transtory Regme To measure λ, te sample s stmulated n by a non necessary unform and constant eat flu ϕ(y,t or temperature step T (y,t. Te averaged temperature s ten consdered: O 6 6 T (, t T (, t T et 6 Heat Transfer Equaton : ( e + T T T e λ a t Boundares Condtons : n T λ ϕ ( t or T T ( t Tet n T λ ϕ ( t T or ( et T t T e << or y Fn s Appromaton Intal Condton : T at t

Teoretcal Model: Soluton n te aplace doman Te soluton s obtaned n te aplace doman (ntegral transform F ( pt ( p ( f ( t f ( t ep Te Quadrupole formulaton allows to lnearly lnks te aplace transforms of te nner and outer temperatures and flues: T Φ ϕ θ ( ( dt θ Φ A C θ Φ θ Φ e [ M ] [ M ] B D θ Φ e θ Φe A B C D sn λk cos ( k ( k λk ( k sn ( k p a Two reference temperature profles are cosen n and sn θ (, p θ( p F (, p + θ ( p F (, p wt: ( ( α( F, p sn( α( and F (, p ( α p a + eλ sn sn ( α( ( α(

Teoretcal Model: Soluton n te Tme doman In te tme doman, te soluton can be wrtten as a sum of two convoluton products: - ( t T ( t f (, t + T ( t f ( t wt: f (, t ( F ( p T,,, measured by nfrared camera f and f are functons of te unknown parameters: a and H..5 *//8 *//4 */3/8 *// 4-4 6-4 *//8 *//4 */3/8 *//..5 * * * * f ( f (,t,t..4.6.8. t * at/..4.6.8. t * at/

Drect Model: Test Case Test Case: square sample 4mm and emm a 5. 7 m s. λ W. m. K W. m. K Temperature step n Temperature * θ (, p In-tme Varatons: θ T θ ( p ( T T p et (, p cos( α( T p cos( α( et Insulated n * * Profle T ( lm p. θ ( p p In-space Varatons: cos cos ( α ( ( α ( T*(T-T et /(T -T et.8.6.4. * *. * *.4 *.6 *.8 *. T*(T-T et /(T -T et.8.6.4. t*. t*.78 t*.3 t*.6 t*. t * Steady-State Profle -..5..5..5.3 t * at/ -...4.6.8 * /

Inverse Model: Parameter Estmaton Ordnary east Squares Metod: T te ( t β T T te te ( t, β, ( t, β, Tte k, M M (, t, β S n t ( Tte ( t Tep ( t ( n + profles and:, β β T ep ( t T T T ep ep ep ( t, M ( k, t ( M, t a H wt: Unknown parameters k β j k n (, K k,n S β j ( t, β nt Tte te β β j ( T ( t, T ( t ( j ep X j β j T te β ( t β, j Senstvty Coeffcent Ideal case on-deal case T ep ( t T ( t, β + ε ( t te E( ε and ( ε V σ b β ˆ β + t t ( X X X ε( t ( t, β T (, t, β + ε ( t f (, t + ( t f ( t T, ep k, te k k ε ε ε k V E( ˆ β β t ( β σ ( X X ˆ b σ b Var Cov ( a Cov( a, H ( a, H Var( H σ σ Var( a σ σ Var( H a b. H b.

Optmal coce of te reference profles statc profles: movng profles: T * (T-T /(T -T et et.8.6.4 3. 3 X. a X H 3 -..5..5..5.3 t * at/ T * (T-T /(T -T et et Var ( a Var ( H Cov ( a, H ρ ( a, H.8.6.4 X a X H -..5..5..5.3 t * at/ 3 Drect Model..36.5.743 General Model..36.5.743 k +.9.397.6.56 k.35.59.8.676.88 4.389 3.65.787 k k + atural & Optmal Coce: and

Effect of a non-unform eat transfer coeffcent T*(T-T /(T -T et et.8.6.4. Sold : Dased : (.3/4*(/ /4 -..5..5..5.3 t * at/ In-space varyng eat transfer coeffcent wt: ( 3 4. 4 ( d (FEM Software - FlePDE

Effect of a non-unform eat transfer coeffcent a 4.98e-7 H 4 a 4.98e-7 H 5456 T * (T-T /(T -T et et.8.6.4. 9,96W. m. K T * (T-T /(T -T et et.8.6.4.,4w. m. K C Ste ( -..5..5..5.3 t * at/ -..5..5..5.3 t * at/ omnal values: 7 a 5. m. s and H 5 λ W. m. K W. m. K

Epermental Results Termograms (lne source stmulaton: Temperature profles and termograms: Estmaton: profle over

Conclusons An n-plane dffusvty estmaton metod by convoluton metod for low conductve materals as been presented. It s a very smple and sem-analytcal metod tat allows us usng two reference measured temperatures as boundary condtons To measure te n-plane dffusvty of low conductve materals Watever te In-Space and In-Tme stmulaton eat flu, And eat transfer coeffcent varatons are. Tank you for your attenton