Publication Cover
Numerical Heat Transfer, Part B: Fundamentals
An International Journal of Computation and Methodology
Volume 72, 2017 - Issue 2
244
Views
9
CrossRef citations to date
0
Altmetric
Original Articles

Nonlinear inverse heat conduction: Digitally filtered space marching with phase-plane and cross-correlation analyses

, &
Pages 109-129 | Received 30 Mar 2017, Accepted 06 Jun 2017, Published online: 25 Jul 2017
 

ABSTRACT

This paper utilizes concepts extracted from phase-plane and cross-correlation analyses for estimating the optimal regularization parameter required by a nonlinear space marching inverse heat conduction method. The regularization parameter for this approach is based on the cut-off frequency used in a low-pass Gauss digital filter for preprocessing the in-depth temperature data. The simulation uses two in-depth thermocouples which are representative of a physical experiment. Numerical results indicate a high level of robustness as the best prediction is extracted from the prediction space over the spectrum of cut-off frequencies. The feasibility and effectiveness of phase-plane and cross-correlation analyses are discussed.

Nomenclature

b1, b2=

in-depth thermocouple probe positions, m

C=

specific heat capacity, J/(kg · K)

Cn=

power spectral density

f, fc=

cut-off frequency, Hz

fc,DFT=

cut-off frequency determined by DFT, Hz

fsampling=

sampling frequency, Hz

fn=

nth chosen cut-off frequency, Hz

k=

thermal conductivity, W/(m · K)

L=

thickness of the one-dimensional geometry, m

M=

total number of spatial nodes in the space marching region x ∈ [0, b1] minus 1

N=

total number of temporal nodes past the initial condition

N*=

number of data points cut in predicted surface heat flux and heat flux rate for phase-plane and the cross-correlation investigation

P=

total number of chosen cut-off frequencies

q=

heat flux, W/m2

=

predicted surface heat flux under specific cut-off frequency f, W/m2

=

predicted surface heat flux under nth chosen cut-off frequency fn, W/m2

=

predicted surface heat flux rate under nth chosen cut-off frequency fn, W/(m2 s)

=

discrete heat flux at xi = iΔx, tj = jΔt in the space marching region x ∈ [0, b1], W/m2

=

, which is heat flux at x = b1 obtained by solving direct problem in x ∈ [b1, b2] using “thermocouple” data, W/m2

=

maximum value of the double Gaussian heat flux defined in Eq. (6), W/m2

=

surface heat flux, W/m2

=

heat flux obtained by solving direct problem in x ∈ [b1, b2] using “thermocouple” data, W/m2

=

heat flux obtained by solving direct problem in x ∈ [b1, b2] using filtered “thermocouple” data by cut-off frequency fc, DFT determined by DFT, W/m2

rj=

jth random number in the interval [−1,1]

Rq=

heat flux cross-correlation coefficient defined in Eq. (5a)

=

heat flux rate cross-correlation coefficient defined in Eq. (5b)

t=

time, s

tk=

kth temporal node, s

tmax=

maximum data collection time, s

tj=

jth temporal node, s

T=

temperature, K

Tf=

predicted surface temperature under chosen cut-off frequency f, W/m2

T0=

initial temperature, K

=

Ttc(b1, tj), which is the “thermocouple” temperature at x = b1, K

=

discrete temperature at xi = iΔx, tj = jΔt in the space marching region x ∈ [0, b1], K

Ts=

surface temperature, K

Ttc=

“thermocouple” temperature, K

=

filtered “thermocouple” temperature by cut-off frequency denoted by fc, K

=

filtered “thermocouple” temperature by cut-off frequency fc,DFT determined by DFT, K

=

filtered “thermocouple” temperature by nth chosen cut-off frequency fn determined by DFT, K

x=

spatial coordinate, m

xi=

xth spatial node, m

x0=

first spatial nodes at x = b1 in the space marching region x ∈ [0, b1]

xM=

last spatial nodes at x = 0 in the space marching region x ∈ [0, b1]

Δf=

the stepping frequency for forming a finite set of predictions, Hz

Δt=

time step, s

Δx=

space step, m

β1=

parameter in the double Gauss function defined in Eq. (6), s

β2=

parameter in the double Gauss function defined in Eq. (6), s

εT=

temperature noise factor

ρ=

density, kg/m3

σ1=

parameter in the double Gauss function defined in Eq. (6), s

σ2=

parameter in the double Gauss function defined in Eq. (6), s

ωc=

circular cut-off frequency, Hz

Subscripts=
s=

surface

tc=

thermocouple

Nomenclature

b1, b2=

in-depth thermocouple probe positions, m

C=

specific heat capacity, J/(kg · K)

Cn=

power spectral density

f, fc=

cut-off frequency, Hz

fc,DFT=

cut-off frequency determined by DFT, Hz

fsampling=

sampling frequency, Hz

fn=

nth chosen cut-off frequency, Hz

k=

thermal conductivity, W/(m · K)

L=

thickness of the one-dimensional geometry, m

M=

total number of spatial nodes in the space marching region x ∈ [0, b1] minus 1

N=

total number of temporal nodes past the initial condition

N*=

number of data points cut in predicted surface heat flux and heat flux rate for phase-plane and the cross-correlation investigation

P=

total number of chosen cut-off frequencies

q=

heat flux, W/m2

=

predicted surface heat flux under specific cut-off frequency f, W/m2

=

predicted surface heat flux under nth chosen cut-off frequency fn, W/m2

=

predicted surface heat flux rate under nth chosen cut-off frequency fn, W/(m2 s)

=

discrete heat flux at xi = iΔx, tj = jΔt in the space marching region x ∈ [0, b1], W/m2

=

, which is heat flux at x = b1 obtained by solving direct problem in x ∈ [b1, b2] using “thermocouple” data, W/m2

=

maximum value of the double Gaussian heat flux defined in Eq. (6), W/m2

=

surface heat flux, W/m2

=

heat flux obtained by solving direct problem in x ∈ [b1, b2] using “thermocouple” data, W/m2

=

heat flux obtained by solving direct problem in x ∈ [b1, b2] using filtered “thermocouple” data by cut-off frequency fc, DFT determined by DFT, W/m2

rj=

jth random number in the interval [−1,1]

Rq=

heat flux cross-correlation coefficient defined in Eq. (5a)

=

heat flux rate cross-correlation coefficient defined in Eq. (5b)

t=

time, s

tk=

kth temporal node, s

tmax=

maximum data collection time, s

tj=

jth temporal node, s

T=

temperature, K

Tf=

predicted surface temperature under chosen cut-off frequency f, W/m2

T0=

initial temperature, K

=

Ttc(b1, tj), which is the “thermocouple” temperature at x = b1, K

=

discrete temperature at xi = iΔx, tj = jΔt in the space marching region x ∈ [0, b1], K

Ts=

surface temperature, K

Ttc=

“thermocouple” temperature, K

=

filtered “thermocouple” temperature by cut-off frequency denoted by fc, K

=

filtered “thermocouple” temperature by cut-off frequency fc,DFT determined by DFT, K

=

filtered “thermocouple” temperature by nth chosen cut-off frequency fn determined by DFT, K

x=

spatial coordinate, m

xi=

xth spatial node, m

x0=

first spatial nodes at x = b1 in the space marching region x ∈ [0, b1]

xM=

last spatial nodes at x = 0 in the space marching region x ∈ [0, b1]

Δf=

the stepping frequency for forming a finite set of predictions, Hz

Δt=

time step, s

Δx=

space step, m

β1=

parameter in the double Gauss function defined in Eq. (6), s

β2=

parameter in the double Gauss function defined in Eq. (6), s

εT=

temperature noise factor

ρ=

density, kg/m3

σ1=

parameter in the double Gauss function defined in Eq. (6), s

σ2=

parameter in the double Gauss function defined in Eq. (6), s

ωc=

circular cut-off frequency, Hz

Subscripts=
s=

surface

tc=

thermocouple

Log in via your institution

Log in to Taylor & Francis Online

PDF download + Online access

  • 48 hours access to article PDF & online version
  • Article PDF can be downloaded
  • Article PDF can be printed
USD 61.00 Add to cart

Issue Purchase

  • 30 days online access to complete issue
  • Article PDFs can be downloaded
  • Article PDFs can be printed
USD 486.00 Add to cart

* Local tax will be added as applicable

Related Research

People also read lists articles that other readers of this article have read.

Recommended articles lists articles that we recommend and is powered by our AI driven recommendation engine.

Cited by lists all citing articles based on Crossref citations.
Articles with the Crossref icon will open in a new tab.