arxiv: v1 [math.na] 23 Feb 2016

Transcription

1 EPJ Web of Conferences will be se by he publisher DOI: will be se by he publisher c Owned by he auhors, published by EDP Sciences, 16 arxiv:163.67v1 [mah.na] 3 Feb 16 Numerical Soluion of a Nonlinear Inegro-Differenial Equaion Ján Buša 1,a, Michal Hnaič,3,4,b, Juha Honkonen 5,c, and Tomáš Lučivjanský,6 1 Deparmen of Mahemaics and Theoreical Informaics, FEE&I, Technical Universiy, Košice, Slovakia Faculy of Sciences, P. J. Šafarik Universiy, Košice, Slovakia 3 Insiue of Experimenal Physics SAS, Košice, Slovakia 4 Bogoliubov Laboraory of Theoreical Physics, JINR, Dubna, Moscow Region, Russia 5 Deparmen of Miliary Technology, Naional Defence Universiy, Helsinki, Finland 6 Fakulä für Physik, Universiä Duisburg-Essen, D-4748 Duisburg, Germany 1 Inroducion Absrac. An algorihm for he numerical soluion of a nonlinear inegro-differenial equaion arising in he single-species annihilaion reacion A + A modeling is discussed. Finie difference mehod ogeher wih he linear approximaion of he unknown funcion is considered. For divergen inegrals presened in he equaion for dimension d a regularizaion is used. Some numerical resuls are presened. The irreversible annihilaion reacion A + A is a funmenal model of non-equilibrium physics. The reacing A paricles are assumed o perform chaoic moion due o diffusion or some exernal advecion field such as amospheric eddy [1]. Many reacions of his ype are observed in diverse chemical, biological or physical sysems [, 3]. In [1] he advecion of reacive scalar using random velociy field generaed by he sochasic Navier-Sokes equaion, which is used for producion a velociy field corresponding o hermal flucuaions [4, 5] and a urbulen velociy field wih he Kolmogorov scaling behavior [6] is sudied by hree of he presen auhors, and he inegro-differenial equaion for he number densiy is derived. No influence of he reacan on he velociy field iself is assumed. In his paper we presen some iniial experiences wih he numerical soluion of he inegrodifferenial equaion menioned above. Problem Formulaion In [1] an inegro-differenial equaion (7) for he mean number densiy a() of chemically acive molecules in anomalous kineics of single-species annihilaion reacion A + A () d a jan.busa@uke.sk b hnaic@saske.sk c juha.honkonen@helsinki.fi λuνµ Z 4 a () + 4λ u ν µ 4 a ( ) d [ 8πuν( ) ] d/, a() a (1)

2 EPJ Web of Conferences is derived. The inegral in Eq. (1) diverges a he upper limi in space dimensions d. We will sudy a numerical soluion o he Iniial Value Problem for he equaion: () d λda () + 4λ D which corresponds o D uν,, and Z Case d a ( ) d [ 8πD( ) ] d/, a() a, () For d he singulariy in he inegral on he righ side of Eq. () a is divergen. This is a consequence of he UV divergences in he model above he criical dimension d c, and near he criical dimension is remedied by he UV renormalizaion of he model [1]. In his paper we use anoher approach o overcome his problem. We will use he following regularizaion: We rewrie Eq. (3) o he form wih d λda + 4λ D a ( ) d 8π{D( ) + l }, a() a. (3) d a ( ) d αa + β + γ α λd, β α 8πD, γ l D. Le us consider he Eq. (4). We will use he difference mehod o solve i numerically. We will consider ime discreizaion wih he ime sep : k k, a( k ) a k, k, 1,,... (5) For he discreizaion of he lef side we will use wo formula he firs and he second order finie difference approximaions: a k a k 1 or 3a k 4a k 1 + a k. (6) d d k For he righ-side inegral approximaion we will use piecewise linear approximaion of he funcion a(u): k a ( ) d k 1 k u + γ a ( ) d k + γ (7) k 1 i i k [ (a i+1 a i )/ + (i + 1) a i i a i+1 ] d k. + γ Inegrals in he righ side of Eq. (7) we calculae analyically using in he following way (4) γ (8) [ (a i+1 a i )/ + (i + 1) a i i a i+1 ] d k + γ

3 Mahemaical Modeling and Compuaional Physics 15 [( k ) (a i+1 a i )/ + (k + ) (a i+1 a i ) + (i + 1) a i i a i+1 ] d k + [(k + i) a i+1 (k + i 1)a i (k + ) (a i+1 a i )/ ] d k + [(k + i)a i+1 (k + i 1)a i ] d k + [(k + i)a i+1 (k + i 1)a i ] (a i+1 a i ) + (a i+1 a i ) ( ) (k + ) d (k + ) (k + ) d (k + + ) [(k + i)a i+1 (k + i 1)a i ] [ ln(k + ) ] (i+1) [(k + i)a i+1 (k + i 1)a i ] (a i+1 a i ) + (a i+1 a i ) ( ) [ (k + ) ] (i+1) { [ ] k + i (k + i)ai+1 (k + i 1)a i ln k + i 1 ( k + i + 1 ) (ai+1 a i ) a i (a i+1 a i ) }, (9) i, 1,..., k 1. Puing approximaions Eq. (6), Eq. (7), and Eq. (9) ino Eq. (4) we arrive a he quadraic equaions Eq. (11) and Eq. (1) in more explici form wih respec o a k : k β i {[ (k + i)ai+1 (k + i 1)a i ] ln 3a k 4a k 1 + a k + α a k (1) k + i k + i 1 ( k + i + 1 ) (ai+1 a i ) a i (a i+1 a i ) } β {[ ] (1 + )ak a k 1 ln 1 + ( ) (ak a k 1 ) a k 1 (a k a k 1 ) } or in he sanrd form [ α β {(1 + ) ln 1 + ( 3 ] [ 3 + )} a k + β a k 1 { + 1 (1 + ) ln 1 + )}] a k 4a k 1 a k {[ ] (k+ i)(ai+1 a i )+a i ln k β i For he firs sep we have a 1 a β a k 1 {1 + αa 1 β{ [ (1 + )a1 a ] ln ln 1 + } (11) k + i k + i 1 (a i+1 a i ) [ (k+ i)(a i+1 a i )+a i + a i+1 + a i ]}. ( ) (a1 a ) a (a 1 a ) }

4 EPJ Web of Conferences 4 log(a()), < a() < log( a()-a/(+c) B ), A6.58, B.56, C Figure 1. Logarihms of he number densiy a() ogeher wih is fi (above), and he fiing error (below) or in he sanrd form [ α β {(1 + ) ln 1 + ( 3 ] [ 1 + )} a 1 + β a { + 1 (1 + ) ln 1 + )}] a 1 a β a {1 + ln 1 + }. (1) The algorihm consiss of he successive calculaion of he values a k, k 1,,..., K solving he quadraic equaion Eq. (1) a he firs sep using he value a a(), and furher solving Eqs. (11) using previously deermined values a i, i,..., k 1. 4 Numerical resuls Below he resuls for a() 5, λ.1, D.4 and l.1 are presened. All calculaions are done wih he uniform ime sep. Figure 1 shows he resuls on inerval [;1] for 1/3. Above he logarihm of a() is shown ogeher wih he logarihm of is approximaion a() 6.58 (.4389).56 (13)

5 Mahemaical Modeling and Compuaional Physics 15 Table 1. Comparison of he a() values a seleced poins for differen values Figure. Precision order p() for four cases:.1.5.5,..., on inerval [7.5;1]. Below he approximaion error less han cca. 1 5 is shown. Table 1 compares he resuls a seleced poins for he increasing precision for he decreasing values from.1 o.315. Even for he sep 1/3.315 he resuls may no be considered o have he sufficien precision. If he sep number is increasing wo imes, he calculaion ime is increasing four imes. If we suppose, ha he difference beween he numerical value a () for he sep and he exac value a() has he order p, e.g., a () a() C ( ) p, (14) where boh C and p depends on, hen using he successive approximaions for differen seps, /, and /4 we ge a () a / () p() log a / () a /4 (). (15)

6 EPJ Web of Conferences Figure shows he precision order p() calculaed for four grid sequences:.1.5.5, , , and One can see, ha for >.5 he p() behavior is smooh, and he orders are cca..8, 1., 1., and 1.4, respecively. For smaller -values he p() behavior is more complicaed. I seems o be reasonable o sudy analyically he behavior of he funcion a() for he small values of, and sar he numerical compuaion from some poin >. Also non-uniform grid could be considered. 5 Conclusions Numerical resuls presened above could give us some basic imaginaion abou he behavior of he number densiy funcion a(). However, furher improvemens of he algorihm are necessary. I will be also ineresing o ry o solve he renormalized inegro-differenial equaion () d { λ [ λuνµ a () + λuνµ a () γ + ln(uνµ ) ] } + 4π + λ uνµ π [ a ( ) a () ] d, (16) where (d )/ and γ.5771 is Euler s consan, presened in [1], and compare he resuls. Anoher possibiliy is o sudy a behavior of he soluion of he problem (1) for d. Acknowledgemens The work was suppored by VEGA Gran 1//13 of he Minisry of Educaion, Science, Research and Spor of he Slovak Republic. References [1] J. Honkonen, M. Hnaič and T. Lučivjanský, EPJ B, 86 : 14 (13) [] B. Derri, V. Hakim and V. Pasquier, Phys. Rev. Le. 75, 751 (1995) [3] R. Kroon, H. Fleuren and R. Sprik, Phys. Rev. E 47, 46 (1993) [4] D. Forser, D.R. Nelson and M.J. Sephen, Phys. Rev. Le. 36, 867 (1976) [5] D. Forser, D.R. Nelson and M.J. Sephen, Phys. Rev. A 16, 73 (1977) [6] L.T. Adzhemyan, A.N. Vasil ev and Y.M. Pis mak, Teor. Ma. Fiz. 57, 68 (1983)