Publication Cover
Numerical Heat Transfer, Part B: Fundamentals
An International Journal of Computation and Methodology
Volume 71, 2017 - Issue 3
91
Views
6
CrossRef citations to date
0
Altmetric
Original Articles

Recovering a general space-time-dependent heat source by the coupled boundary integral equation method

&
Pages 283-297 | Received 16 Aug 2016, Accepted 28 Nov 2016, Published online: 16 Mar 2017
 

ABSTRACT

As to recover a time-dependent heat source under an extra temperature measured at an interior point, we can reformulate it to be a three-point boundary value problem. We can develop a coupled boundary integral equation method, wherein by selecting two sets of adjoint test eigenfunctions in two sub domains and using polynomials as the trial functions of unknown heat source, the time-dependent heat source is recovered very well and quickly. Four numerical examples, including a discontinuous one, demonstrate the efficiency for the ill-posed inverse heat source problem in a large time duration and under a large noise up to 10–30%. Then, selecting three sets of adjoint test eigenfunctions in three domains: problem domain and two sub domains, and using the Pascal polynomials as trial functions, the unknown space-time-dependent heat source is recovered very fast and accurately from the solution of three coupled boundary integral equations.

Nomenclature

A=

coefficient matrix in Eq. (28)

aj=

coefficients in polynomials

b=

right-hand side in Eq. (28)

b1=

:= BTb

bj=

coefficients in polynomials

B=

the matrix obtained from A by equilibrated method

cij=

coefficients in the Pascal polynomials

c=

vector of coefficients

D=

:= BTB

e(t)=

:= T(t) − T0(t)

f(x)=

initial temperature

H(t)=

time-dependent heat source

H(x, t)=

general heat source

=

heat operator

ℋ*=

adjoint heat operator

=

length of the rod

m=

m − 1 highest order in the Pascal polynomials

m0=

number of test functions

m1=

m1 − 1 the highest order of polynomials

m2=

m2 − 1 the highest order of polynomials

n=

= m1 + m2 for time-dependent heat source

n=

m(m + 1)/2 for general heat source

nq=

= 2m0 for time-dependent heat source

nq=

= 3m0 for general heat source

p(t)=

:= Tm(t) − T0(t)

P=

post conditioner in Eq. (31)

q(t)=

:= T(t) − Tm(t)

R(i)=

random numbers

s=

level of noise

=

multiple scales

=

multiple scales

sij=

multiple scales in the Pascal polynomials

t=

time

tf=

final time

T(x, t)=

temperature

T0(t)=

left-boundary temperature

T(t)=

right-boundary temperature

Tm(t)=

measuring temperature at xm

u(x, t)=

transformed temperature

uk(x, t)=

adjoint Trefftz test functions in Ω

vk(x, t)=

adjoint Trefftz test functions in Ω1

wk(x, t)=

adjoint Trefftz test functions in Ω2

x=

space variable

xm=

measuring point

Greek symbols=
ε=

converegence criterion

Ω=

problem domain

Ω1=

sub domain

Ω2=

sub domain

Γ=

boundary of Ω

Subscripts and superscripts=
i=

index

j=

index

k=

index

m=

index

T=

transpose

Nomenclature

A=

coefficient matrix in Eq. (28)

aj=

coefficients in polynomials

b=

right-hand side in Eq. (28)

b1=

:= BTb

bj=

coefficients in polynomials

B=

the matrix obtained from A by equilibrated method

cij=

coefficients in the Pascal polynomials

c=

vector of coefficients

D=

:= BTB

e(t)=

:= T(t) − T0(t)

f(x)=

initial temperature

H(t)=

time-dependent heat source

H(x, t)=

general heat source

=

heat operator

ℋ*=

adjoint heat operator

=

length of the rod

m=

m − 1 highest order in the Pascal polynomials

m0=

number of test functions

m1=

m1 − 1 the highest order of polynomials

m2=

m2 − 1 the highest order of polynomials

n=

= m1 + m2 for time-dependent heat source

n=

m(m + 1)/2 for general heat source

nq=

= 2m0 for time-dependent heat source

nq=

= 3m0 for general heat source

p(t)=

:= Tm(t) − T0(t)

P=

post conditioner in Eq. (31)

q(t)=

:= T(t) − Tm(t)

R(i)=

random numbers

s=

level of noise

=

multiple scales

=

multiple scales

sij=

multiple scales in the Pascal polynomials

t=

time

tf=

final time

T(x, t)=

temperature

T0(t)=

left-boundary temperature

T(t)=

right-boundary temperature

Tm(t)=

measuring temperature at xm

u(x, t)=

transformed temperature

uk(x, t)=

adjoint Trefftz test functions in Ω

vk(x, t)=

adjoint Trefftz test functions in Ω1

wk(x, t)=

adjoint Trefftz test functions in Ω2

x=

space variable

xm=

measuring point

Greek symbols=
ε=

converegence criterion

Ω=

problem domain

Ω1=

sub domain

Ω2=

sub domain

Γ=

boundary of Ω

Subscripts and superscripts=
i=

index

j=

index

k=

index

m=

index

T=

transpose

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.