Conversion of inhomogeneous robin boundary conditions into virtual sources for wave motions and heat conduction

Boe-Shong Hong1 , Po-Jen Su2

1Department of Mechanical Engineering, National Chung Cheng University, Chia-Yi 62102, Taiwan

2Medical Devices and Opto-Electronics Equipment Department, Metal Industries Research and Development Centre, Luzhu Dist, Kaohsiung 82151, Taiwan

1Corresponding author

Journal of Vibroengineering, Vol. 17, Issue 6, 2015, p. 2846-2857.
Received 29 June 2015; received in revised form 3 August 2015; accepted 11 August 2015; published 30 September 2015

Copyright © 2015 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License

In vibration engineering, the differential equations of wave motions and heat conduction are usually accompanied by inhomogeneous boundary conditions in practice. Boundary inhomogeneity makes the dynamics essentially nonlinear, which prevents Hilbert space from being applied for modal decomposition. To deal with this difficulty, this paper does not treat boundary inhomogeneity as a “condition”, but almost converts it into a virtual source in conjunction with homogeneous boundary. This conversion counts mostly on the Laplace-Galerkin transform, a functional tool developed in previous works. We also explore boundary topology of this virtual-source conversion, and find that its strategy is to zero the environment and simultaneously create a spatially impulsive source on the homogeneous boundary, yielding almost the same solution. In one-dimensional region, such a boundary source takes the form of Dirac delta function usually combined by its derivatives. In a sense, this paper catches how Nature really handles boundary conditions.

Keywords: inhomogeneous Robin boundary conditions, nD transfer function modelling, Sturm-Liouville systems, thermoacoustic dynamics.

1. Introduction

Vibration engineering encounter a great quantity of longitudinal waves, transverse waves and heat conduction, such as parabolic or hyperbolic heat-conduction dynamics, acoustic or thermoacoustic oscillations, structural vibrations, quantum mechanics, electromagnetic waves, and so on. Any of these dynamics is governed by a Laplacian operator or its higher orders in space, which is often spatially non-uniform. As the boundary condition, Dirichlet, von-Neumann, or Robin, thereof is homogeneous, its eigenfunctions constitute an admissible, real, orthogonal and complete basis in the Hilbert space of a bounded region. This basis provides modal decomposition of the dynamics for system identification, computational intelligence, model reduction, real-time processing and design purposes. However, boundary conditions are inhomogeneous in many occasions; for instance, the differential equation of heat conduction is always accompanied by inhomogeneous boundary conditions, since temperature is non-zero in nature. With boundary inhomogeneity, the dynamics is essentially nonlinear, which prevents Hilbert space from being directly introduced for modal decomposition.

To remedy such a situation in this paper, we realize the inhomogeneous boundary conditions as virtual sources in conjunction with homogeneous boundary conditions. In one-dimensional cases, such a source is found to be a Dirac Delta distribution combined by its spatial derivatives on boundary. Therein, with the Laplace-Galerkin transform [1-3], both equations governing the interior and the boundary are integrated into a single 2D transfer-function of two independent variables: one is from the time and the other is from the space. Performing the inverse Laplace-Galerkin transform of the 2D transfer-function realizes back the dynamics into homogeneous boundary conditions with virtual sources, both of which yield the identical solution in the interior. With the homogeneous boundary resulting from the virtual-source realization, the Sturm-Liouville properties in Hilbert space are thus applicable for further analysis and synthesis in the mode-frequency domain.

The conversion of boundary inhomogeneity into virtual source has ever been conceptually applied to identify thermal inertia by inputting von-Neumann source [2], and to derive the mechanical energy of thermoacoustics [3]. In these decades, von-Neumann boundary source was employed to obtain order-reduced modelling of combustion instabilities in rocket motors. With these hints on application, this paper systematically extends to Robin sources and studies boundary topology. Robin inhomogeneity is fascinating and necessary for real practice, since the Dirichlet or von-Neumann boundary can be considered as a degenerated version of Robin boundary [4-7]. Therein, the boundary is the complement of the union of the exterior and the interior of the domain under consideration, so boundary conditions rely on the interaction between the process and the environment. Although degenerative Robin simplifies numerical or experimental investigation, the actual Robin boundary should be identified in practice through measured data, as seen in [8-13] for examples. Therein, Robin boundaries are identified for the study of cancer destruction during hyperthermia treatment, of the optical path length in inhomogeneous tissue, and of axisymmetrical induction in heating processes, respectively. To be sure, the conversion of Robin inhomogeneity into virtual source can make these important kinds of identification more accurate and reliable, since the virtual-source conceptually suggests installing an active source in measurement to trig out desired data.

Conversion of boundary inhomogeneity into virtual source also helps computational intelligence, since boundary inhomogeneity is conventionally treated as “conditions” that constrain the spatial-temporal evolution. The legitimate spatial-temporal solutions, including the integral solutions [14-15], finite-element approximation [16-17], and series solutions [18-20], have to be developed toward matching the boundary conditions. In vibration engineering, modal solutions are particularly popular, since it reveals the spatiotemporal structure and results in model reduction in practice. In cases of time-invariant environments, the response is usually computed by shifting the origin of spatial coordinate to the steady-state response, upon which the dynamics with homogeneous boundary condition can be solved by separation-of-variable method, as in [21] for example. This method can also be extended to time-varying environments by stepwise sampling the temporal continuity, as in [22-23] for examples. Compared with the virtual-source solution, this solution is numerical tedious and incapability of capturing sudden changes in the environment.

More importantly, virtual-source conversion lead to an input-output modelling that makes possible real-time signal processing. The differential equation in conjunction with the boundary inhomogeneity can give designers computer-time solutions by taking the boundary inhomogeneity as constraining conditions, but it is unable to catch the real-time nature. That is, such a computer-time version is offline plugged into computer simulation or calculation to merely obtain the solution as a function of a preset time-span. However, a real-time version is an emulator of natural evolution, wherein the state at the next instant is only dependent on the state and the boundary inhomogeneity at the present instant. It is unnecessary to know the history of the state and the boundary inhomogeneity to predict the future, since Nature has no memory. Galerkin projection of the converted dynamics with virtual source and homogeneous boundary onto a proper basis, such as those from Proper Orthogonal Decomposition (POD) [24-28], generates an order-reduced state-space realization. Given a freely assigned sampling time, Euler discretization of the state-space realization becomes the real-time version, which interprets the dynamic nature into two times of matrix multiplication and one time of matrix addition within a sampling time.

Compared with “conditions” realization of boundary inhomogeneity, virtual-source realization provides the following merits in practice:

1) It transforms the interaction between two distributed dynamics adjacent to each other into feedback interconnection, such as thermal-acoustic interaction in the fields of thermoacoustic engines [29] and combustion instabilities [30]. The construction of feedback makes possible the application of modern or classical control theory to help design and analysis.

2) Modal decomposition is applicable for computational intelligence and order-reduced modelling. It results in real-time version of distributed dynamics with inhomogeneous in digital signal processing (DSP), which can be directly programmed into a microcontroller for real-time estimation of state distribution and environmental changes, toward a newly sensing technology.

3) The real-time fashion above can still replace computer-time versions as a numerical simulator programmed into a generous computer. For a long run, it can be employed for efficient management of computer memory.

4) Virtual-source realization makes possible the frequency-domain identification of boundary inhomogeneity.

5) As the control actuation is set on some boundary, the virtual-source realization generates an input-output model served for feedback synthesis of boundary control systems.

6) Even with temporally discontinuous or impulsive environments, exact solutions can be calculated offline with virtual-source realization.

7) Active sources on boundary can be installed for identification of Robin coefficients.

2. The considered class of dynamics

This paper considers the following three types of distributed dynamics with inhomogeneous Robin boundary conditions:

ρ 2 ψ t 2 - · k ψ = 0     in       Ω ,
α ψ + β ψ · n ^ = f       o n     Ω ,       Longitudinal wave ,


ρ 2 ψ t 2 + 2 k 2 ψ      in         Ω ,
α 1 ψ + β 1 ψ · n ^ = f 0 ,       α 2 k 2 ψ + β 2 k 2 ψ · n ^ = f 1     on   Ω ,       Transverse wave ,


ρ ψ t - · k ψ = 0    in       Ω ,
α ψ + β ψ · n ^ = f     on   Ω ,       Heat conduction .

Therein the spatial functions ρ, k are real and positive in the bounded region ΩR3 is a bounded region, wherein the distributed output is denoted by ψ; on the boundary Ω, (α,β)(0.0), (α1,β1)(0.0), (α2,β2)(0.0), and the boundary inhomogeneity is denoted by f or fj’s.

A kind of these three dynamics involves a spatial Laplacian operator a:

a ϕ = - 1 ρ · k ϕ     in      Ω ,
α ϕ + β ϕ · n ^ = 0     on   Ω ,       Elastic stiffness ,


a ϕ = 1 ρ 2 k 2 ϕ     in     Ω ,
α 1 ϕ + β 1 ϕ · n ^ = 0 ,       α 2 k 2 ϕ + β 2 k 2 ϕ · n ^ = 0     on   Ω ,       Bending stiffness .

With the inner-product metric:

ψ , ϕ = Ω ρ ( x ) ϕ * ( x ) ψ ( x ) d V ,

the Laplacian operators a in Eq. (4) and Eq. (5) belong to the Sturm-Liouville class aSL(Ω), that is, their eigenfunctions constitute a real, orthonormal, and complete basis of L2(Ω). A Sturm-Liouville operator is usually distinguished from its self-adjointness and the compactness of its inverse. With the Green’s second identity, it is easy to know that the elastic stiffness in Eq. (4) is a Sturm-Liouville operator [31]. As for the bending stiffness in Eq. (5), the following shows that it is belonging to Sturm-Liouville class.

Let two operators P and Q be defined by P=ρ-12 and Q=k2, then the bending stiffness a becomes their composite, i.e. a=PQ. For any ψ, ϕD(a):

P Q ψ , ϕ Ω = Q ψ , P ϕ Ω ,

since α1(Qψ)+β1(Qψ)·n^=0 and α1ϕ+β1ϕ·n^=0. Moreover:

ψ , P Q ϕ Ω = P ψ , Q ϕ Ω ,

since α2(Qϕ)+β2(Qϕ)·n^=0; α2ψ+β2ψ·n^=0. Observe that:

P ψ , Q ϕ Ω = Q ψ , P ϕ Ω ,

since both sides equal Ωk2ψ·2ϕdV. Therefore, the bending stiffness operator a is self-adjoint. Moreover, the inverse of a is a compact operator in L2(Ω), since a is fourth-order differential operator. Therefore, its eigenfunctions constitute a real, orthonormal, and complete basis of L2(Ω). Moreover, it can be shown that both elastic stiffness and bending stiffness are positive definite in this work.

3. Laplace-Galerkin transform- a functional tool

With respect to the eigenfunctions set ϕλλΛ of a Sturm-Liouville operator a, the Galerkin transform G from spatial functions to modal functions, F(λ)=G[f(x)], is defined by:

F λ Ω ρ x ϕ λ x f x d x .

Completeness and orthonormality of ϕλλΛ of countable cardinality jointly imply that the Galerkin transform G has a unique inverse G-1, f(x)=G-1[F(λ)]:

f x λ Λ F λ ϕ λ x .

Then, the Laplace-Galerkin transform H from spatial-temporal functions to modal-complex functions is defined by the composite of the Galerkin transform G and the Laplace transform L:

H = L G = G L ,


F λ , s H f x , t = 0 - Ω e - s t ρ x ϕ λ x f x , t d x d t .

Accordingly, the inverse of Laplace-Galerkin transform H-1 is the composite of the inverse of Laplace transform and that of Galerkin transform, that is:

H - 1 = G - 1 L - 1 = L - 1 G - 1 ,


f x , t H - 1 F λ , s = 1 2 π j λ Λ Γ F λ , s ϕ λ x e t s d s .

Here the domain Γ is an infinite line parallel to the imaginary axis, whereon the integral in Eq. (13) is converged.

Denote the temporal derivative /t by Dt, and let a be a Sturm-Liouville operator, aSL(Ω). For the set of spatial-temporal functions with homogeneous boundary and initial, the Laplace-Galerkin transform H is of:

H h a , D t f x , t = h λ , s · H f x , t ,

where h is a ratio of two expressions of finite or some infinite length constructed from two independent variables, one standing for space and the other for time, allowing for the operations of addition, subtraction, multiplication, integer exponents in time, and fraction-order exponents in space. For example:

H D t - a 1 2 D t + a 1 2 f x , t = s - λ s + λ F λ , s .

The Laplace-Galerkin transform and its inverse perform transformation between Sturm-Liouville dynamics in space-time domain and 2D transfer-function in mode-frequency domain [1-3]. As an example to explain 2D transfer-function, let us find the impulse response of the following longitudinal wave G^:

2 ψ t 2 - 2 ψ x 2 = u ,       0 x π ,       0 t < ,
ψ 0 , t = 0 ,       ψ π , t = 0 ,       0 t < ,
ψ x , 0 = 0 ,       ψ ˙ x , 0 = 0 ,       0 x π .

The elastic stiffness-2/x2 is of eigenvalues Λ=1, 4, 9, associated with eigenfunctions ϕλ(x)=2/πsinλx. Taking the Laplace-Galerkin transform on both sides of the differential equation with homogeneous boundary and initial yields:

s 2 + λ Ψ λ , s = U λ , s ,

that is, the 2D transfer-function of the dynamics G^ is:

G λ , s Ψ λ , s U λ , s = 1 s 2 + λ .

Correspondingly, the impulse response g=H-1G is to be:

g x , t = ω = 1 2 π s i n ω x L - 1 1 s 2 + ω 2 = ω = 1 2 ω π s i n ω x s i n ω t .

To check whether this solution is correct, let us give the dynamics G^ the 2D unit-pulse u(x,t)=δ(t)λΛϕλ(x), where Hu=1. Integration of the differential equation from t=0- to t=0+ yields the initial condition: ψ(x,0)=0 and ψ˙(x,0)=λΛ ϕλ(x). Thereby, the impulse response g is just the solution of the initial-value problem:

2 ψ t 2 - 2 ψ x 2 = 0 ,
ψ 0 , t = 0 ,       ψ ( π , t ) = 0 ,
ψ x , 0 = 0 ,       ψ ˙ x , 0 = λ Λ   ϕ λ x ,

which has the form solvable by the conventional separation-of-variable method. It can be found that two solutions are identical.

4. Virtual conversion of conditions into sources on boundary

This section demonstrates how to converts the boundary inhomogeneity into virtual source in conjunction with homogeneous boundary. This conversion is analogous to the strategy of Laplace transform dealing with non-zero initial conditions, wherein the initial inhomogeneity is realized as a virtual source comprising Dirac Delta function and its derivatives. Consider the following explanatory example- a normalized one-dimensional wave dynamics with initial inhomogeneity:

2 ψ t 2 - 2 ψ x 2 = 0 ,         ψ 0 , t = 0 ,       ψ ( π , t ) = 0 ,
ψ x , 0 = 0 ,         ψ ˙ x , 0 = f x .

Taking the Laplace transform on Eq. (23a) with the help of integration by parts yields:

Ψ λ , s = 1 s 2 + λ G f x .

Then taking the inverse Laplace transform on Eq. (24) yields:

2 ψ t 2 - 2 ψ x 2 = f x δ t ,
ψ 0 , t = 0 ,       ψ ( π , t ) = 0 ,
ψ x , 0 = 0 ,       ψ ˙ ( x , 0 ) = 0 ,

where δ is the Dirac delta distribution. The solution to Eq. (25) is almost the same as the solution to Eq. (23); both solutions are identical in the interior t>0 that is an open set, but different on the boundary t=0 that is a closed set.

Now consider the longitudinal wave dynamics G^ in Eq. (1), which involves the elastic stiffness a in Eq. (4). As shown in Section 2, the elastic stiffness a is a Sturm-Liouville operator aSL(Ω) under the inner-product of Eq. (6). Let Φ=ϕλλΛ denote the eigenfunctions set of a corresponding to the eigenvalues set Λ. On the boundary Ω, firstly, substitution αϕλ +βϕλ ·n^=0 for Eqs. (1b) ×(-ϕλ ) yields:

ψ ϕ λ   - ϕ λ   ψ · n ^ = - ϕ λ   β f ,         β 0 .

Secondly, substitution αϕ+βϕλn^=0 for Eq. (1b) ×ϕλ ·n^ yields:

ψ ϕ λ   - ϕ λ   ψ · n ^ = ϕ λ   · n ^ / α f       for       α 0 .

In general, the sum of Eq. (26) ×β/(α+β) and Eq. (27) ×α/(α+β) is to be:

ψ ϕ λ   - ϕ λ   ψ · n ^ = - ϕ λ   + ϕ λ   · n ^ α + β f .

Moreover, based on the Green’s second identity:

Ω ϕ λ · k ψ d V = Ω ψ · k ϕ λ d V + Ω k ψ ϕ λ - ϕ λ ψ · n ^ d S .

With Eqs. (26)-(29), performing Laplace-Galerkin transform Hon Eq. (1a) with inhomogeneous boundary conditions in Eq. (1b) yields:

s 2 + λ Ψ λ , s = Ω k ϕ λ ψ - ψ ϕ λ · n ^ d S Ω k x B λ x f ^ x , s d S Q λ , s ,

where Ψλ,sH[ψ(x,t)]f^(x,s)L[f(x,t)], and Bλ is to be:

B λ = ϕ λ   β ,   β 0 , - ϕ λ   · n ^ α , α 0 , ϕ λ   - ϕ λ   · n ^ α + β , in general.

Therefore, in the sense of virtual source, the 2D transfer-function G of the dynamics G^ in Eq. (1) is:

G λ , s Ψ λ , s Q λ , s = 1 s 2 + λ .

With the 2D transfer-function of Eq. (32), the exact solution of Eq. (1) can be obtained even when the boundary inhomogeneity f in Eq. (30) or Eq. (1b) is temporally impulsive or discontinuous.

Performing the inverse Laplace-Galerkin transform H-1 on the Eq. (32) yields:

2 ψ t 2 - 1 ρ · k ψ = q     in     Ω ,
α ψ + β ψ · n ^ = 0       on     Ω ,

where q(x,t)H-1[Q(λ,s)]. In Eq. (33), the boundary inhomogeneity f in Eq. (1b) has been realized as the virtual source q in Eq. (33a) in conjunction with homogeneous boundary of Eq. (33b). In the interior of the domain Ω, the response governed by Eq. (1) is identical to that governed by Eq. (33), since both have the same 2D transfer-function.

Next, consider the transverse wave dynamics G^ in Eq. (2). For simple explanation, let us take its on-dimensional version:

2 ψ t 2 - 1 ρ 2 x 2 k 2 ψ t 2 = 0 ,         in     0 , l ,
α 1 k ψ ' ' + β 1 k ψ ' ' ' = f 1     and       α 2 ψ + β 2 ψ ' = f 2       at      x = 0 ,
ψ = 0     a n d     ψ ' ' = 0     at    x = l .

This dynamics involves the bending stiffness a in Eq. (5):

a ϕ = 1 ρ 2 x 2 k 2 ϕ x 2     in     0 , l ,
α 1 ϕ + β 1 ϕ ' = 0     and       α 2 k ϕ ' ' + β 2 k ϕ ' ' ' = 0     at     x = 0 ,
ϕ = 0     and     ϕ ' ' = 0     at     x = l ,

which is a Sturm-Liouville operator, aSL([0,l]), under the inner-product of Eq. (6). Denote its eigenvalues set by Λ and eigenfunctions set by Φ=ϕλλΛ.

With integration by parts (one-dimensional Green’s second identity), we have:

0 l ϕ λ ( k ψ ' ' ) ' ' d x = 0 l ψ ( k ϕ λ ' ' ) ' ' d x
            + ϕ λ ' ( k ψ ' ' ) - ϕ λ ( k ψ ' ' ) ' + ψ ( k ϕ λ ' ' ) ' - ψ ' ( k ϕ λ ' ' )   x = 0 .

With Eq. (34b) and Eq. (35b), Eq. (36) can be rephrased to be:

0 l ϕ λ ( k ψ ' ' ) ' ' d x = 0 l ψ ( k ϕ λ ' ' ) ' '   d x - ϕ λ 0 β 1 f 1 - k 0 ϕ λ ' ' 0 β 2 f 2 .

With Eq. (37), performing the Laplace-Galerkin transform H on Eq. (34a) yields:

s 2 + λ Ψ λ , s = ϕ λ 0 β 1 F 1 s + k 0 ϕ λ ' ' 0 β 2 F 2 s .

Then, performing the inverse Laplace-Galerkin transform H-1 on Eq. (38) yields:

2 ψ t 2 + 1 ρ 2 x 2 k 2 ψ x 2 = 1 β 1   ρ 0 f 1 t δ x + k 0 β 2 f 2 t δ ' ' x ,
α 1 k ψ ' ' + β 1 k ψ ' ' ' = 0     and     α 2 ψ + β 2 ψ ' = 0     at     x = 0 ,
ψ = 0       and     ψ ' ' = 0     at       x = l .

The dynamics in Eq. (34) has been converted into the dynamics in Eq. (39), wherein the virtual source in Eq. (39a) comprises a Delta function δ and the second derivative of Delta function δ''ϕ,δ''=ϕ''0 distributed at x=0 in conjunction with homogeneous boundary in Eqs. (39b) and (39c). Both dynamics are identical in (0,l], since they have the identical 2D transfer-function.

As for the virtual-source realization of the dynamics of heat conduction in Eq. (3), it is similar to that of the dynamics of longitudinal wave as shown above. The above shows how the wave and heat-conduction dynamics with inhomogeneous boundary “conditions” can be almost converted to delta sources on boundary in conjunction with homogeneous boundary. Input-output modelling results therefrom and modal decomposition in Hilbert space becomes possible.

5. Investigation on boundary topology

In this section, the virtual-source solution and the separation-of-variable solution of a normalized, one-dimensional heat condition with one-side Dirichlet inhomogeneity are computed and then visualized in figures. It is expected to visualize between them the identical parts in the interior but the topological difference on the inhomogeneous boundary. Accordingly, consider the following parabolic dynamics:

y t - 2 y x 2 = 0     for    x 0 , π ,       t 0 , ,
y 0 , t = f t ,       y π , t = 0     for      t 0 , ,
y x , 0 = 0     for     x [ 0 , π ] .

Two cases of environments are considered: one is time-invariant in which f(t)=1 (Case I), and the other is time-varying in which f(t)=sinωt (Case II).

The virtual-source conversion, as shown in Section 4, is employed to solve the Case I and Case II. Based on the Green’s second identity:

y ' ' , ϕ n = - y ϕ ' n π 0 + y , ϕ n ' ' .

With Eqs. (40b) and (40c), taking the Laplace-Galerkin transform on Eq. (40a) yields:

Y n , s = ϕ n ' 0 · F s s + n 2 .

As f(t)=1 in Case I, F(s)=1/s. Then taking the inverse Laplace-Galerkin transform on Eq. (42) yields:

y x , t = 2 π n = 1 N s i n n x 1 n 1 - e - n 2 t .

In case II, let f(t) be sinωt, that is, F(s)=ω/(s2+ω2). Taking the Laplace-Galerkin transform again on Eq. (40a) yields:

y x , t = 2 π n = 1 N n n 4 + ω 2 s i n n x ( ω   e - n 2 t + n 2 s i n ω t - ω c o s ω t ) .

The separation-of-variable method is improper to solve the solution of Case II, wherein the environmental temperature is time-varying. In Case I, the separation-of-variable solution is to be:

y x , t = - 1 π x + 1 + 2 π n = 1 N b n n s i n n x e - n 2 t ,       b n ( - 1 ) n + 1 + ( - 1 ) n - 1 .

Fig. 1 shows the order-reduced solutions of Case I with 500 modes being considered. At t=0.5 the separation-of-variable solution in Eq. (45) is in juxtaposition with the virtual-source solution of Eq. (43) for comparison. It is found that both solutions are identical in the interior of the domain. However, close view on the inhomogeneous boundary x=0 as shown in Fig. 2 reveals that the exact solution (N) by the virtual-source conversion is discontinuous at x=0. With the inverse Laplace-Galerkin transform, Eq. (42) is almost equivalent to:

y t - 2 y x 2 = - δ ' x f t ,

simultaneously with homogeneous boundary. Eq. (46) implies that the virtual source is actually impulsive on the homogeneous boundary. The virtual-source conversion realizes boundary inhomogeneity as delta sources combined by their derivatives.

Figs. 3 and 4 show the virtual-source solution for time-varying environment: f(t)=sint. At t=0.5, two order-reduced solutions according to N=200 and N=1000, respectively, are juxtaposed with each other. The convergence on the boundary appears, which implies that the solution will reach discontinuity at x=0 as N. This verifies again that the virtual-source conversion is of the strategy to zero the environment and simultaneously give impulsive sources on the homogeneous boundary, yielding almost the same solution.

Fig. 1. Topological comparison of the virtual-source solution with the separation-of-variable solution under time-invariant environment

 Topological comparison of the virtual-source solution with the separation-of-variable solution under time-invariant environment

Fig. 2. Close view on the inhomogeneous boundary of Fig. 1

 Close view on the inhomogeneous  boundary of Fig. 1

Fig. 3. Topological convergence of the virtual-source solution under time-varying environment

 Topological convergence of the virtual-source solution under time-varying environment

Fig. 4. Close view on the inhomogeneous boundary of Fig. 3

 Close view on the inhomogeneous  boundary of Fig. 3

6. Conclusion

For wave and heat-conduction dynamics, we can realize inhomogeneous boundary “conditions” as virtual-sources in conjunction with homogeneous boundary. The strategy such a virtual-source conversion is to zero the environment and simultaneously to introduce an impulsive source onto the homogeneous boundary, yielding the identical response in the interior of operation domain. Therein the boundary impulsive takes the form of Dirac delta function combined by its derivatives. In practice, this virtual-source conversion helps modal decomposition, real-time processing, computational intelligence, and system identification of distributed dynamics.


  1. Hong B.-S. Construction of 2D isomorphism for 2D H-control of Sturm-Liouville systems. Asian Journal of Control, Vol. 12, Issue 2, 2010, p. 187-199. [Search CrossRef]
  2. Hong B.-S., Chou C.-Y. Realization of thermal inertia in frequency domain. Entropy, Vol. 16, 2014, p. 1101-1121. [Search CrossRef]
  3. Hong B.-S., Chou C.-Y. Energy transfer modelling of active thermoacoustic engines via Lagrangian thermoacoustic dynamics. Energy Conversion and Management, Vol. 84, 2014, p. 73-79. [Search CrossRef]
  4. Gustafson K., Abe T. (Victor) Gustave Robin: 1855-1897. Mathematical Intelligencer, Vol. 20, 1988, p. 47-53. [Search CrossRef]
  5. Gustafson K., Abe T. The third boundary condition – was it Robin’s? Mathematical Intelligencer, Vol. 20, Issue 1, 1988, p. 63-71. [Search CrossRef]
  6. Romeo A., Saharian A. A. Casimir effect for scalar fields under Robin boundary conditions on plates. Journal of Physics A: Mathematical and General, Vol. 35, 2002, p. 1297. [Search CrossRef]
  7. Mintz B., Farina C., Maia Neto P. A., Rodrigues R. B. Particle creation by a moving boundary with a Robin boundary condition. Journal of Physics A: Mathematical and General, Vol. 39, 2006, p. 11325-11333. [Search CrossRef]
  8. Okada E., Schweiger M., Arridge S. R., Firbank M., Delpy D. T. Experimental validation of Monte Carlo and finite-element methods for the estimation of the optical path length in inhomogeneous tissue. Applied Optics, Vol. 35, Issue 9, 1996, p. 3362-3371. [Search CrossRef]
  9. Bay F., Labbe V., Favennec Y., Chenot J. L. A numerical model for induction heating processes coupling electromagnetism and thermomechanics. International Journal for Numerical Methods in Engineering, Vol. 58, 2003, p. 839-867. [Search CrossRef]
  10. Jin B. Conjugate gradient method for the Robin inverse problem, associated with the Laplace equation. International Journal for Numerical Methods in Engineering, Vol. 71, 2007, p. 433-453. [Search CrossRef]
  11. Xiong X. T., Liu X. H., Yan Y. M., Guo H. B. A numerical method for identifying heat transfer coefficient. Applied Mathematical Modelling, Vol. 34, 2010, p. 1930-1938. [Search CrossRef]
  12. Majchrzak E., Paruch M. Identification of electromagnetic field parameters assuring the cancer destruction during hyperthermia treatment. Inverse Problems in Science and Engineering, Vol. 19, Issue 1, 2011, p. 45-58. [Search CrossRef]
  13. Jin B., Liu X. Numerical identification of a Robin coefficient in parabolic problems. Mathematics of Computation, Vol. 81, 2012, p. 1369-1398. [Search CrossRef]
  14. Cicalese M., Trombetti C. Asymptotic behaviour of solutions to p-Laplacian equation. Asymptotic Analysis, Vol. 35, 2003, p. 27-40. [Search CrossRef]
  15. Andreu F., Mazon J. M., Moll J. S. The total variation flow with nonlinear boundary conditions. Asymptotic Analysis, Vol. 43, Issues 1-2, 2005, p. 9-46. [Search CrossRef]
  16. Bartels S., Carstensen C., Dolzmann G. Inhomogeneous Dirichlet conditions in a priori and a posteriori finite element error analysis. Numerische Mathematik, Vol. 99, 2004, p. 1-24. [Search CrossRef]
  17. Gudi T. A new error analysis for discontinuous finite element methods for linear elliptic problems. Mathematics of Computation, Vol. 79, 2010, p. 2169-2189. [Search CrossRef]
  18. Bareket M. On an isoperimetric inequality for the first eigenvalue of a boundary value problem. Journal of Mathematical Analysis, Vol. 8, 1977, p. 280-287. [Search CrossRef]
  19. Daners D. A Faber-Krahn inequality for Robin problems in higher dimensions. Mathematische Annalen, Vol. 333, 2006, p. 767-785. [Search CrossRef]
  20. Auchmuty G. Bases and comparison results for linear elliptic eigenproblems. Journal of Mathematical Analysis and Applications, Vol. 390, 2012, p. 394-406. [Search CrossRef]
  21. Moosaie A. Axisymmetric non-Fourier temperature field in a hollow sphere. Archive of Applied Mechanics, Vol. 79, 2009, p. 679-694. [Search CrossRef]
  22. Abdel-Hamid B. Modelling non-Fourier heat conduction with periodic thermal oscillation using the finite integral transform. Applied Mathematical Modelling, Vol. 23, 1999, p. 899-914. [Search CrossRef]
  23. Heidarinejad G., Shirmohammadi R., Maerefat M. Heat wave phenomena in solids subjected to time dependent surface heat flux. Heat and Mass Transfer, Vol. 44, 2008, p. 381-392. [Search CrossRef]
  24. Joshi Y. Reduced order thermal models of multiscale microsystems. Journal of Heat Transfer, Vol. 134, 2012, p. 031008. [Search CrossRef]
  25. Agharkar P., Subramanian P., Kaisare N. S., Sujith R. I. Thermoacoustic instabilities in a ducted premixed flame: reduced-order models and control. Combustion Science and Technology, Vol. 185, Issue 6, 2013, p. 920-942. [Search CrossRef]
  26. Selimefendigil F., Oztop H. F. POD-based reduced order model of a thermoacoustic heat engine. European Journal of Mechanics – B/Fluids, Vol. 48, 2014, p. 135-142. [Search CrossRef]
  27. Behzad F., Helenbrook B. T., Ahmadi G. On the sensitivity and accuracy of proper orthogonal decomposition- based order models for Burgers equation. Computers and Fluids, Vol. 106, 2015, p. 19-32. [Search CrossRef]
  28. Kunisch K., Muller M. Uniform convergence of the POD method and application to optimal control. Discrete and Continuous Dynamical Systems, Vol. 35, Issue 9, 2015, p. 4477-4501. [Search CrossRef]
  29. Hong B.-S., Lin T.-Y. System identification and resonant control of thermoacoustic engines for robust solar power. Energies, Vol. 8, Issue 5, 2015, p. 4138-4159. [Search CrossRef]
  30. Zhang Z., Guan D., Zheng Y., Li G. Characterizing premixed laminar flame-acoustics nonlinear interaction. Energy Conversion and Management, Vol. 98, 2015, p. 331-339. [Search CrossRef]
  31. Young N. An Introduction to Hilbert Space. Cambridge University Press, 1988. [Search CrossRef]