Entertainment & Media

A numerical approach for determination of sources in transport equations

A numerical approach for determination of sources in transport equations
of 12
All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.
  Pergamon Computers Math. Applic. Vol. 32, No. 5, pp. 31-42, 1996 Copyright~)1996 Elsevier Science Ltd Printed in Great Britain. All rights reserved 0898-1221/96 $15.00 + 0.00 S0898-1221(96)00133-2 A Numerical Approach for Determination of Sources in Transport Equations I. DIMOV Bulgarian Academy of Sciences Acad. G. Bonchev Str., bl. 25 A 1113, Sofia, Bulgaria U. JAEKEL AND H. VEREECKEN KFA-ICG-4, D-52428 Juelich, Germany (Received and accepted April 1996) Abstract--A new numerical approach for locating the source and quantifying its power is pre- sented and studied. The proposed method may also be helpful when dealing with other transport problems. The method uses a numerical integration technique to evaluate a linear functional arising from the solution of the adjoint problem. It is shown that using the solution of the adjoint equation, an efficient umerical method can be constructed (a method which is less ime consuming than the usual schemes applied for solving the srcinal problem). Our approach leads to a well-conditioned numerical problem. A number of numerical tests are performed for a one-dimensional problem with a linear advection part. The results are compared with results obtained when solving the original problem. The code which realizes he presented method is written in FORTRAN 77. Keywords--Numerical method, Adjoint problem, Fundamental solution, Source location. 1. INTRODUCTION In this paper, we discuss--for a simple, but nontrivial example--how to obtain information on unknown sources in a transport equation from observed data. We assume that the data J~ are of the form = f f dtp x, tlu(x,t), (1) where x and t are the space and time coordinates, respectively, u(x, t) is the solution of the transport equation, and pi(x, t) is a function or distribution specified by the observation process. This covers a large number of important special cases. Some examples are: 1. pi(x,t) := 5(x-xi)5(t-ti). This corresponds to a measurement at time ti and position xi. 2. p~(x,t) := 5(x - x~)H(t~ - t). In this case, Ji is the total mass which has passed the position x~ until a time t~. 3. pi(x, t) := r(xi - x)s(ti - t). This type occurs for measurements with only finite spatial and temporal resolution (characterized by the functions r and s). The information about the sources is reconstructed from the data using an adjoint formulation of the problem. We demonstrate its usefulness by solving two problems for a one-dimensional, not translation-invariant, and hence nontrivial, example: 1. We determine the set of points where a point source can be located such that the solution U(Xl,tl) at a fixed position and time is smaller or larger than a certain value C. As Typeset by ~4A~a-TEX 31  32 I. DIMOV et al. one of many possible applications, one might think about the environmental problem of determining regions where possible pollutant sources must not be located, in order to keep pollutant concentrations below a prescribed level in a protected region. 2. We reconstruct the location and power of a steady-state point source from two obser- vations U(xl,tl) and u(x2,t2). A possible application would be the identification and characterization of a pollutant source in an inaccessible area. Our formulation leads to a well-conditioned problem and can be generalized to a large class of more complicated cases. 2. FORMULATION Let us consider the following problem of 1D-transport: (~ ) ~)u O (axu ) _02u - L u = at Ox - l)-ff~x2 = Q(t)6(z - xo), (2) u = u0(x), for t = 0. (3) Assume that u(z, t) E W~, where W~ is a Sobolev space, u(x, t) > O, D > O, -oo < z < oo, t E [0, T], and the conditions providing the existence of the unique solution of the problem (2),(3) in a weak formulation are fulfilled. The problem presented by (2),(3) does not have an analytical solution, due to presence of the function Q(t). In addition, standard numerical methods used to solve equation (2) are not efficient because there is a 6-function in (2). The condition u(z, t) E W~ implies that loTt t u (2) (x, t) dx < oo, (4) which is a natural condition. Now formulate the problems under consideration. PROBLEM 1. Find the subdomain g of G (g c G) such that if at any point of g, a source of power Q(t) will be located, the concentrations of the pollutants in a given point x -- Xl at a given time t -- tl will be less than a given constant C, that is, find the set of points z0 for which it(X1, tl) __~ C. (5) PROBLEM 2. Find the place of location x0 and the power of the source Q(t) = Q = const if two measurement data u(xl, tl) and u(x2, t2) are available. (It is assumed, that Zl ~ x2 and tt ~t t2.) The first problem (the problem of finding the subdomain g) is not so complicated if the solution of the adjoint formulation of the srcinal problem (2),(3) will be used. 3. ADJOINT FORMULATION In fact, multiply equation (2) by a function u*, which will be defined later, and integrate on time and space: dt oo u* cgz axu) - O-~x 2 dx = Q t) dt oo u*~i z - xo) dz. (6) Applying the integration by part to (6), one can obtain (it is assumed that u*lt=T = O) IoTF.° .F loTfO t ~ ~ ~= = .~* ~11=o ~ - dt ~,--~- ~ £ Io - ~ ~(=)~;(=) d= - dt J_~ u--~- d=, (7) ~T ~.a2,,. ~T/.a,, a,,', ~T. t°° a2u dt ~ ~ ~= : k ~ ~ - ~T~) d=l~:__oo + dt J-~o -n'a-.~a= =, (S) and  Sources in Transport Equations ~oT f2~X) F tfT _ ~oT f °°~* t u* dx = xuu* dxltffi 0 dt - xu~ dx J-oo Oz o f5 ,T,UO ,T)U~ X)dx ~T /~ OU* _ - - dt xu-~z dx. oo oo Let us assume that u = u* = 0 when x --* oo or x --* -c~. Now from (7)-(9), one can obtain: 33 (9) T oo 02U, T ~0 dt /_ u (-~- + aX~--~ -D-~x2 ) dx- f_:(1-x)uo(x)u~(x)dx = ~o Q(t)u'(xo, t) dr. (10) Let the function u* satisfy the following adjoint equation: with an initial condition Ou* Ou* _ 02u ----~- + ~x-~, - v-y~ 2 = p(x, t), (11) u*(x,T) = 0, (12) where the function p(x, t) will be defined later. It is possible to write (10) in the next form: fo f[ f ff t p(x, t)u(x, t) dx = Q(t)u*(xo, t) dt + (1 - x)uo(x)u3(x ) dx = a. (13) oo Equation (13) will be a basic equation for our numerical method. One can choose the function p(x, t) in the form of a product of two 6-functions: p(x, t) -~ ~(X -- Xl)~ t -- tl). 14) Let us also assume that Q(t) contains the Heaviside function H(t), that is, Q(t) = q(t)g(t), where H(t) = 1 for t > 0 and H(t) = 0 for t <_ 0. This means that the second term of the right-hand side of (13) vanishes. Consider the following linear functional of the solution of the srcinal problem: /o /2 = dt p(x, t)u(x, t) dx. (15) In this case, we have oT f_: = dt u(x,t)~(x-xl)6(t-tl)dX = u(xl,tl). (16) Obviously, the value of the functional (15) is equal to the value of the concentration in point x = Xl at the time moment t = tl. From the representation (13), it follows that instead of solving the srcinal problem, one can solve the adjoint problem for u*: Ou* Ou* 02U * ----~- + ax--~- - D-y~-x2 = e(x - Xl)6(t - tl) u* = O, for t = T. (17) (18)  34 I. DIMOV et al. From the solution of the problem formulated by (17),(18), one can estimate the functional (13) as a function depending on a parameter x0 and solve the first problem under consideration. In fact, if the function u*(x,t, xl, tl) is the solution of the adjoint problem (17),(18) (the solution u*(x, t, Xl, tl) of the adjoint problem depends on parameters Xl and tl), then the presentation (13) defines a function of Xo--J(xo, Xl,tl). Now, solving the inequality J(xo,xl,tl) <C, (19) the set g of points x0 can be determined. That is the solution of the first problem under consid- eration. Now, consider the second problem. Let the power of the source q and its location x0 be unknown. Two measurement data points U(Xl, tl) and u(x2, t2) are assumed to be available. One can then define two functions of xo--F1 (Xo) and F2 (x0)--using the following formal presentations: J(Xo, zi,ti) = J(xo, xi, t,), i = 1,2; (20) q d(zo, z~, td F~(xo) = u(xi,ti) ' i= 1,2. (21) (It is possible to do this because q is strongly positive.) One can see that x0 is the solution of equation FI(x)=F2(x), (22) being the point of the location of the source of pollutants. The value of the power of the source q can be determined using the expression 1 q= Fi(xo) (23) Now, the problem consists in finding an efficient numerical method for evaluating the functional ~o ' q(t)u*(xo, t) dr, where u*(x0, t) is the solution of the adjoint problem (17),(18) and x0 is a parameter. 4. THE METHOD In a general case, one needs an efficient method for evaluating the linear functionals fo e` q(t)u*(xo, t) dt (24) of the solution of the adjoint problem (17),(18). It is known that statistical numerical methods allow us to find directly the unknown functional of the solution with a number of operations necessary to solve the problem in one point of the domain [1,2]. The statistical methods give statistical estimates for the functional of the solution by performing random sampling of a certain chance variable whose mathematical expectation is the desired functional. These methods have proved to be very efficient in solving multidimensional problems in composite domains [1,3-6]. Moreover, it is shown that for some problems (including one-dimensional ones) in the correspond- ing functional spaces, the statistical methods have a better convergence rate than the optimal deterministic methods in such functional spaces [7-9]. It is also very important that the statistical methods are very efficient when parallel or vector processors or computers are available, because  Sources n I1 ansport quations 35 the above-mentioned methods are inherently parallel nd have loose dependencies. They are also well vectorizable. Using power parallel omputers, it is possible o apply the Monte Carlo method or particle-tracking ethod for evaluating large-scale onregular problems which sometimes are difficult o be solved by the well-known numerical methods. The difficulties or finite ifference nd finite lement methods appear from nonregularity of the right-hand side of equation (17) and from so-called artificial or numerical scheme ) diffusion. So, when L is a general elliptic perator using a statistical umerical method, the functional (24) can be evaluated. In our case, the fundamental solution of the adjoint equation can be used. In fact, introduce a new variable t* -- T - t, where t* 6 [0, ] and present the adjoint problem in the next form: Ou* Ou* _ 02u * - 1)-E -2 . -t- OtX-~--~- ---- I~(X -- Xl)I~(T -- ~* -- l~l) , u* = 0, for t* = 0. (25) (26) The fundamental solution of problem (25) can be found using the techniques [10--12] 1 { [x-b(t)] 2} fi*(x,t) = [?ca(t)]1~ exp aM , (27) where the time-dependent functions a(t) and b(t) are defined below: a(t) ---- 29 (1--e -2at) t~ b(t) = xoe -~t. (28) (29) The solution of the adjoint problem can be expressed now in the following form: /oT// *(x,t*) = dr* oo (30) The result of integration (30) will be a function defined by (27) with new arguments { [X-xl-b(tl-t)] 2 ) U*(X -- Xl,~I -- $) = [71.a($1 _ $)]1/2 exp a(~ --~) (31) Substituting (31) into the functional (24), one can get ~0 ' 1 t)] 1/2 ( [XO- Xl - b(~l -t)] 2 a_(t.l.:Q q(Q [~a(tl - exp dt = J(xo,xl,tl). (32) The proposed numerical method is based on a numerical integration scheme for evaluating the last parametrized functional. Using a simple discretization scheme, one can get k_l } J" = ~qi, - 1 { [xo -- ~1 -- b tl - ti)] 2 ,=o t~ra(tl - ti)]l/2 exp ~(~--~)- (ti,+l -- ti) -F 0(7"), (33) where r is the maximal time-step ~ = max ti+ 1 -- t,), $ and q, = q t4. The numerical scheme (33) has an error O(7). It is possible to apply more complicated schemes of high-order accuracy (say, O(r2), or O(v4)). But from a computational point of view, when the
Similar documents
View more...
Related Search
We Need Your Support
Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks

We need your sign to support Project to invent "SMART AND CONTROLLABLE REFLECTIVE BALLOONS" to cover the Sun and Save Our Earth.

More details...

Sign Now!

We are very appreciated for your Prompt Action!