MONTE CARLO STUDY OF PARTICLE TRANSPORT PROBLEM IN AIR POLLUTION. R. J. Papancheva, T. V. Gurov, I. T. Dimov
|
|
- Ambrose Haynes
- 5 years ago
- Views:
Transcription
1
2 Pliska Stud. Math. Bulga. 14 (23), STUDIA MATHEMATICA BULGARICA MOTE CARLO STUDY OF PARTICLE TRASPORT PROBLEM I AIR POLLUTIO R. J. Papancheva, T. V. Guov, I. T. Dimov Abstact. The actual tanspot of the ai pollutants is due to the wind. This nomally called advection of the ai pollutants. Diffusion and deposition ae othe two majo physical pocesses, which take place duing the tanspot of pollutants in the atmosphee. In this pape we study two classes of gid-fee Monte Calo (MC) algoithms fo solving an elliptic bounday value poblem, whee the patial diffeential equation contains advection, diffusion and deposition pats. The gid-fee MC appoach to solve the above equation uses a local integal epesentation and leads to a stochastic pocess called a andom Walk on balls (WOB). In the fist class of algoithms, the choice of a tansition density function in the Makov chain depends on the adius of the maximal ball, lying inside the domain, in which the poblem is defined, and on the paametes of the diffeential opeato. While the choice of a tansition density function in the second class of algoithms does not depend on the deposition pat of the poblem. The computational complexity of both classes of gid-fee MC algoithms was investigated using vaied numeical tests on a PowePC (G4 w/altivec) 45 MHz unning YDL Mathematics Subject Classification: 65C5 Key wods: Monte Calo algoithms Suppoted by Cente of Excellence BIS-21 gant ICA and by the SF of Bulgaia unde Gants #I 811/98 and #MM 92/99
3 18 R. J. Papancheva, T. V. Guov, I. T. Dimov 1. Fomulation of the poblem Conside the functional (1) J(u) (g, u) = Ω g(x)u(x)dx, whee Ω is a domain in the Euclidean space R 3 and x = (x 1, x 2, x 3 ) Ω. The functions u(x) and g(x) belong to the Banach space X and to the adjoint space X, espectively, and u(x) is a unique solution of the following Fedholm integal equation: (2) u(x) = k(x, y)u(y)dy + f(x). Ω We ae inteested to estimate the functional (1), using MC appoach, whee Eq. (2) is a local integal epesentation of the solution of the following bounday value poblem: (3) Mu = Φ(x), x Ω, Ω R 3, (4) u = ψ(x), The opeato M is defined by: M = 3 ( i=1 2 x 2 i x Ω. + b i (x) x i ) + c(x). Such poblems appea in envionmental sciences. The functional (1) can be consideed as a measue fo the damaging effect of a dange pollution on the life in the natue. In this case g(x) is a given sensitivity function, and u(x) is the concentation of the pollutant, that can be descibed as a solution of the bounday value poblem (3) (4). Hee, we intoduce the following notations: b = max x Ω b(x), c = max x Ω c(x), and denote by R the adius of the maximal ball, lying inside Ω and by R(x) the adius of the maximal ball, lying inside Ω with cente in the point x. We assume that the conditions fo existing of an unique solution of the poblem (3) (4) ae satisfied [1, 11], and divb(x) =. Using the Levy s function L p (y, x) [3] as a Geen s function we obtain following integal fom: (5) u(x) = B(x) M y L p (y, x)u(y)dy + B(x) L p (y, x)φ(y)dy.
4 Monte Calo study of paticle tanspot poblem In the Eq.(5) B(x) is the ball inside Ω with cente in the point x and adius R(x). M = 3 ( ) b i (x) x i + c(x) is the adjoint opeato to M, and (6) i=1 2 x 2 i k(x, y) = M y L p(y, x) = µ p (R) p() 2 + µ p(r) 2 [ c(y) + 3 i=1 R µ p (R)c(y) b i (y) y i x i ] R p(ρ) ρ dρ p(ρ)dρ. The epesentation of u(x) in (5) is the basis fo the poposed MC appoach. Using it a biased estimato fo the solution can be obtained. 2. Monte Calo algoithms The MC estimato, whose mathematical expectation coincides with J(u) is (7) Θ[g] = g(ξ ) π(ξ ) Q j f(ξ j ), j= whee Q = 1, Q j = Q j 1 k(ξ j 1,ξ j ) p(ξ j 1,ξ j ), j = 1, 2, 3,... and ξ, ξ 1,... is a Makov chain in Ω with initial density function π(x) and tansition density p(x, y), which ae toleant to g(x) and k(x, y), espectively [2, 5, 11]. To ensue the convegence of the pocess, we intoduce an ε - stip of the bounday. The pocess stats at point ξ = x Ω, chosen with the initial density function π(x). The next andom point y is detemined by a tansition density function p(x, y). This pocess teminates when the point falls into the ε - stip of the bounday. In the fist class of algoithms (class A), the tansition density is taken to be popotional to the kenel (6), i.e. (8) k(x, y) p(x, y) = k(x, y)dy. B(x) In this case we need to know when the integal tansfomation kenel (6) is non-negative. It is poved in [4] that p(x, y), when the following inequality (9) p() (b + R4 ) R c p(ρ)dρ
5 11 R. J. Papancheva, T. V. Guov, I. T. Dimov is tue. In the second class (class B) the tansition density function in the Makov chain is taken to be the kenel, whee the deposition paamete is equal to zeo, i.e. (1) p(x, y) = k(x, y) = My L p(y, x) with c(y). It is poved in [1] that when the condition: (11) p() b R p(ρ)dρ is fulfilled, p(x, y) can be used as tansition density function in Makov chain. Hee we conside thee cases fo p(): whee case I p() = e k, case II p() = p = const, if kr 1, case III p() = e k, if kr ln 2, b + R 4 c fo class A, k = k = b fo class B. The function p(x, y) can be witten in spheical coodinates as: (12) p(, w) = sin θ 4π p() q p (R) p(w ), whee fo the class A we have: (13) p(w ) = [ ] R b(x + w) cos(b, w) + c(x + w) 1 + p(ρ)dρ p() R c(x + w)2 p(ρ) p() ρ dρ, and fo algoithms fom class B the conditional pobability is: (14) p(w ) = 1 + b(x + w) cos(b, w) p() R p(ρ)dρ.
6 Monte Calo study of paticle tanspot poblem The next andom point y in the Makov chain depends on both the jump into the maximal ball, and on the diection w. We sample the jump with a density function p()/q p (R), using an invese-tansfomation ule. The andom diection w with density function p(w ) = sin θ 4π p(w ) is computed, using an acceptance-ejection (AR) technique [4, 1]. 3. umeical tests As an example, the following bounday value poblem was solved, using two classes of MC algoithms: (15) (16) ( 3 2 u i=1 x 2 i + b i (x) u ) + c(x)u = in Ω = [, 1] 3, x i u = ψ(x), x Ω ε. Two test vaiants ae consideed fo ψ(x), b(x) and c(x): Test 1: ψ(x 1, x 2, x 3 ) = e a 1x 1 +a 2 x 2 +a 3 x 3, (x 1, x 2, x 3 ) Ω ε, b 1 (x) = a 2 a 3 (x 2 x 3 ), b 2 (x) = a 3 a 1 (x 3 x 1 ), b 3 (x) = a 1 a 2 (x 1 x 2 ), Test 2: c(x) = (a a2 2 + a2 3 ); ψ(x 1, x 2, x 3 ) = e a 1(x 1 +x 2 +x 3 ), (x 1, x 2, x 3 ) Ω ε, b 1 (x) = a 2 x 1 (x 2 x 3 ), b 2 (x) = a 2 x 2 (x 3 x 1 ), b 3 (x) = a 2 x 3 (x 1 x 2 ), c(x) = 3a 2 1. Hee a 1, a 2, a 3 ae paametes. It is easy to see that divb(x) =. This condition guaantees the possibility to use the local integal epesentation by Geen s function. Diffeent values fo the coefficients a i, i = 1, 2, 3 ae consideed. In ou tests the following ε-stips ae used: ε =.1, ε =.5, ε =.1. Figues 1 3 show the CPU times fo algoithms of classes A and B, with diffeent values fo coefficients a i. One can see that the algoithm B is faste than the algoithm A, with inceasing of b. This behavio can be explained with the less numbe of the moves in the Makov chains (see Table 1). When the maximum of the advection vecto b(x) inceases the aveage length of Makov
7 112 R. J. Papancheva, T. V. Guov, I. T. Dimov chain (El (ε,b) ) fo the algoithm B is shote than El (ε,a) fo algoithm A. Figue 3 shows that the algoithm B is about 1.75 times faste than the algoithm A in case I when a i = 3, b = The elative eos of the algoithms ae pactically the same fo both classes A and B, and they depend on b (see Figues 4 7). Figue 8 compaes CPU times of the MC algoithms fom the class B with a diffeent choice of the density function p(). One can see that the algoithm in case II is faste than the othes. The less CPU time of this case is due to the computational simplicity. On the othe hand the algoithm in case III is a little bit faste than the algoithm in a case I. This can be explained with esults fo the aveage moves in the WOB that ae pesented in Table 2. We conclude that the numeical esults show that the algoithms of class B ae moe efficient. Also we ecommend to use a case II o a case III, if the paametes of the poblem allow this. a i El (ε,a) El (ε,b) b p() = e k p() = e b Table 1: The aveage lengths of Makov chains with ε =.1. 5 A: p()=e -k B: p()=e -b 4 3 CPU Figue 1: Compaison of the CPU times among algoithms A and B in case when ε =.1, a =.25 fo Test 1.
8 Monte Calo study of paticle tanspot poblem A: p()=e -k B: p()=e-b 4 3 CPU Figue 2: Compaison of the CPU times among algoithms A and B in case when ε =.1, a = 2 fo Test A: p()=e -k B: p()=e -b CPU Figue 3: Compaison of the CPU times among algoithms A and B in case when ε =.1, a = 3 fo Test p()=e -k p()=e -b* Figue 4: Vaiation of the elative eo in case, when ε =.1, a =.25.
9 114 R. J. Papancheva, T. V. Guov, I. T. Dimov.25.2 p()=e -k p()=e -b* Figue 5: Vaiation of the elative eo in case, when ε =.1, a = p()=e -k p()=e -b* Figue 6: Vaiation of the elative eo in case, when ε =.1, a = A: p()=const B: p()=const Figue 7: Vaiation of the elative eo in case, when ε =.1, a =.25 fo Test 2.
10 Monte Calo study of paticle tanspot poblem B: p()=e -b* B: p()=const B: p()=e b* CPU Figue 8: Compaison of the CPU times among algoithms B in case, when ε =.1, a = 1 fo Test 2. Test 1 Test 2 El (ε,a) El (ε,b) El (ε,a) El (ε,b) case I ε = ε = ε = case III ε = ε = ε = Table 2: The aveage lengths of Makov chains in case when a =.25.
11 116 R. J. Papancheva, T. V. Guov, I. T. Dimov R E F E R E C E S [1] A. V. Bitzadze. Equations of the Mathematical Physics, auka, Moscow, [2] J. H. Cutiss, Monte Calo methods fo the iteation of the linea opeatos, J. Math. Phys (1954), [3] I. T. Dimov, T. V. Guov. Estimates of the computational complexity of iteative Monte Calo algoithm based on Geen s function appoach. Mathematics and Computes in Simulation (1998) [4] I. T. Dimov, R. J. Papancheva. Geen s function Monte Calo algoithms fo elliptic poblems, Jounal of Mathematics and Computes in simulations (submitted). [5] S. M. Emakov, G. A. Mikhailov. Statistical Simulation, auka, Moscow, [6] S. M. Emakov, V.V. ekutkin, A.S. Sipin. Random Pocesses fo Solving Classical Equations of the Mathematical Physics, auka, Moscow, [7] T. V. Guov, P. A. Withlock, I. T. Dimov. A gid fee Monte Calo algoithm fo solving elliptic bounday value poblems, Lectue notes in Comp. Sci. (eds. L.Vulkov et al.), Spinge-Velag (21), [8] V. P. Mikhailov. Patial Diffeential Equations, auka, Moscow, [9] C. Mianda, Equasioni alle diivate paziali di tipo ellipttico. Spinge- Velag, Belin, [1] R. J. Papancheva, I. T. Dimov, T. V. Guov. A new class gid-fee Monte Calo algoithms fo elliptic bounday value poblems, In: Poceeding of the MA 2 Confeence , Boovets, Bulgaia (submitted). [11] I. M. Sobol, Monte Calo umeical Methods, auka, Moscow, R. J. Papancheva, T. V. Guov, I. T. Dimov Depatment of Paallel algoithms Cental Laboatoy fo Paallel Pocessing Bulgaian Academy of Sciences Acad. G. Bonchev St., bl. 25A 1113 Sofia, Bulgaia guov@copen.bas.bg, ivdimov@bas.bg, umi@canto.bas.bg
Green s function Monte Carlo algorithms for elliptic problems
Mathematics and Computes in Simulation 63 (2003) 587 604 Geen s function Monte Calo algoithms fo elliptic poblems I.T. Dimov, R.Y. Papancheva Cental Laboatoy fo Paallel Pocessing, Depatment of Paallel
More informationTemporal-Difference Learning
.997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation
More informationSOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS
Fixed Point Theoy, Volume 5, No. 1, 2004, 71-80 http://www.math.ubbcluj.o/ nodeacj/sfptcj.htm SOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS G. ISAC 1 AND C. AVRAMESCU 2 1 Depatment of Mathematics Royal
More informationGoodness-of-fit for composite hypotheses.
Section 11 Goodness-of-fit fo composite hypotheses. Example. Let us conside a Matlab example. Let us geneate 50 obsevations fom N(1, 2): X=nomnd(1,2,50,1); Then, unning a chi-squaed goodness-of-fit test
More informationMODELING OF CYLINDRICAL COUETTE FLOW OF RAREFIED GAS. THE CASE OF ROTATING INNER CYLINDER *
th National Congess on Theoetical and Applied Mechanics, -5 Sept. 9, Boovets, Bulgaia MODELING OF CYLINDRICAL COUETTE FLOW OF RAREFIED GAS. THE CASE OF ROTATING INNER CYLINDER * PETER GOSPODINO Institute
More information15 Solving the Laplace equation by Fourier method
5 Solving the Laplace equation by Fouie method I aleady intoduced two o thee dimensional heat equation, when I deived it, ecall that it taes the fom u t = α 2 u + F, (5.) whee u: [0, ) D R, D R is the
More information12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,
th WSEAS Int. Conf. on APPLIED MATHEMATICS, Caio, Egypt, Decembe 9-3, 7 5 Magnetostatic Field calculations associated with thick Solenoids in the Pesence of Ion using a Powe Seies expansion and the Complete
More informationMONTE CARLO SIMULATION OF FLUID FLOW
MONTE CARLO SIMULATION OF FLUID FLOW M. Ragheb 3/7/3 INTRODUCTION We conside the situation of Fee Molecula Collisionless and Reflective Flow. Collisionless flows occu in the field of aefied gas dynamics.
More informationPearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationMathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden
Applied Mathematical Sciences, Vol. 7, 13, no. 7, 335-348 Mathematical Model of Magnetometic Resistivity Sounding fo a Conductive Host with a Bulge Ovebuden Teeasak Chaladgan Depatment of Mathematics Faculty
More informationComputational Methods of Solid Mechanics. Project report
Computational Methods of Solid Mechanics Poject epot Due on Dec. 6, 25 Pof. Allan F. Bowe Weilin Deng Simulation of adhesive contact with molecula potential Poject desciption In the poject, we will investigate
More informationCOMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS
Pogess In Electomagnetics Reseach, PIER 73, 93 105, 2007 COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS T.-X. Song, Y.-H. Liu, and J.-M. Xiong School of Mechanical Engineeing
More informationApplication of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates
Computational Ecology and Softwae 5 5(): 9-5 Aticle Application of homotopy petubation method to the Navie-Stokes equations in cylindical coodinates H. A. Wahab Anwa Jamal Saia Bhatti Muhammad Naeem Muhammad
More informationSurveillance Points in High Dimensional Spaces
Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage
More informationis the instantaneous position vector of any grid point or fluid
Absolute inetial, elative inetial and non-inetial coodinates fo a moving but non-defoming contol volume Tao Xing, Pablo Caica, and Fed Sten bjective Deive and coelate the govening equations of motion in
More informationSOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES
italian jounal of pue and applied mathematics n. 35 015 (433 44) 433 SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF OPERATOR MATRICES Watheq Bani-Domi Depatment of Mathematics
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationAbsorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere
Applied Mathematics, 06, 7, 709-70 Published Online Apil 06 in SciRes. http://www.scip.og/jounal/am http://dx.doi.og/0.46/am.06.77065 Absoption Rate into a Small Sphee fo a Diffusing Paticle Confined in
More informationSOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD
1 SOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD NGO By CO. H. TRAN, PHONG. T. Faculty of Mathematics & Infomatics, Univesity of Natual Sciences VNU-HCM Abstact : The pape pesents some
More informationGROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS
Annales Academiæ Scientiaum Fennicæ Mathematica Volumen 32, 2007, 595 599 GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS Teo Kilpeläinen, Henik Shahgholian and Xiao Zhong
More information, the tangent line is an approximation of the curve (and easier to deal with than the curve).
114 Tangent Planes and Linea Appoimations Back in-dimensions, what was the equation of the tangent line of f ( ) at point (, ) f ( )? (, ) ( )( ) = f Linea Appoimation (Tangent Line Appoimation) of f at
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationKOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS
Jounal of Applied Analysis Vol. 14, No. 1 2008), pp. 43 52 KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS L. KOCZAN and P. ZAPRAWA Received Mach 12, 2007 and, in evised fom,
More informationPerturbation to Symmetries and Adiabatic Invariants of Nonholonomic Dynamical System of Relative Motion
Commun. Theo. Phys. Beijing, China) 43 25) pp. 577 581 c Intenational Academic Publishes Vol. 43, No. 4, Apil 15, 25 Petubation to Symmeties and Adiabatic Invaiants of Nonholonomic Dynamical System of
More informationELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS
THE 9 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS R. Sbulati *, S. R. Atashipou Depatment of Civil, Chemical and Envionmental Engineeing,
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationResearch Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function
Abstact and Applied Analysis Volume 011, Aticle ID 697547, 7 pages doi:10.1155/011/697547 Reseach Aticle On Alze and Qiu s Conjectue fo Complete Elliptic Integal and Invese Hypebolic Tangent Function Yu-Ming
More informationME 210 Applied Mathematics for Mechanical Engineers
Tangent and Ac Length of a Cuve The tangent to a cuve C at a point A on it is defined as the limiting position of the staight line L though A and B, as B appoaches A along the cuve as illustated in the
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Seies UG Examination 2015 16 FLUID DYNAMICS WITH ADVANCED TOPICS MTH-MD59 Time allowed: 3 Hous Attempt QUESTIONS 1 and 2, and THREE othe questions.
More informationMeasure Estimates of Nodal Sets of Polyharmonic Functions
Chin. Ann. Math. Se. B 39(5), 08, 97 93 DOI: 0.007/s40-08-004-6 Chinese Annals of Mathematics, Seies B c The Editoial Office of CAM and Spinge-Velag Belin Heidelbeg 08 Measue Estimates of Nodal Sets of
More informationSolution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere
Applied Physics eseach; Vol. 9, No. 6; 7 ISSN 96-969 E-ISSN 96-9647 Published by Canadian Cente of Science and Education Solution of a Spheically Symmetic Static Poblem of Geneal elativity fo an Elastic
More information1 Explicit Explore or Exploit (E 3 ) Algorithm
2.997 Decision-Making in Lage-Scale Systems Mach 3 MIT, Sping 2004 Handout #2 Lectue Note 9 Explicit Exploe o Exploit (E 3 ) Algoithm Last lectue, we studied the Q-leaning algoithm: [ ] Q t+ (x t, a t
More informationLecture 10. Vertical coordinates General vertical coordinate
Lectue 10 Vetical coodinates We have exclusively used height as the vetical coodinate but thee ae altenative vetical coodinates in use in ocean models, most notably the teainfollowing coodinate models
More informationConservative Averaging Method and its Application for One Heat Conduction Problem
Poceedings of the 4th WSEAS Int. Conf. on HEAT TRANSFER THERMAL ENGINEERING and ENVIRONMENT Elounda Geece August - 6 (pp6-) Consevative Aveaging Method and its Application fo One Heat Conduction Poblem
More informationϕ direction (others needed boundary conditions And we have following boundary conditions in
Exact Analytical 3-D Steady-State Solution fo Cylindical Fin ANITA PILIKSERE, ANDRIS BUIKIS Institute of Mathematics and Compute Science Univesity of Latvia Raina ulv 9, Riga, LV459 LATVIA saulite@latnetlv,
More informationHomework # 3 Solution Key
PHYSICS 631: Geneal Relativity Homewok # 3 Solution Key 1. You e on you hono not to do this one by hand. I ealize you can use a compute o simply look it up. Please don t. In a flat space, the metic in
More informationNew problems in universal algebraic geometry illustrated by boolean equations
New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic
More informationPhysics 2212 GH Quiz #2 Solutions Spring 2016
Physics 2212 GH Quiz #2 Solutions Sping 216 I. 17 points) Thee point chages, each caying a chage Q = +6. nc, ae placed on an equilateal tiangle of side length = 3. mm. An additional point chage, caying
More informationarxiv: v1 [math.na] 8 Feb 2013
A mixed method fo Diichlet poblems with adial basis functions axiv:1302.2079v1 [math.na] 8 Feb 2013 Nobet Heue Thanh Tan Abstact We pesent a simple discetization by adial basis functions fo the Poisson
More informationON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS
ON THE TWO-BODY PROBLEM IN QUANTUM MECHANICS L. MICU Hoia Hulubei National Institute fo Physics and Nuclea Engineeing, P.O. Box MG-6, RO-0775 Buchaest-Maguele, Romania, E-mail: lmicu@theoy.nipne.o (Received
More informationSupplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in
Supplementay Figue 1. Cicula paallel lamellae gain size as a function of annealing time at 50 C. Eo bas epesent the σ uncetainty in the measued adii based on image pixilation and analysis uncetainty contibutions
More informationWeb-based Supplementary Materials for. Controlling False Discoveries in Multidimensional Directional Decisions, with
Web-based Supplementay Mateials fo Contolling False Discoveies in Multidimensional Diectional Decisions, with Applications to Gene Expession Data on Odeed Categoies Wenge Guo Biostatistics Banch, National
More information6 Matrix Concentration Bounds
6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom
More informationGeometry of the homogeneous and isotropic spaces
Geomety of the homogeneous and isotopic spaces H. Sonoda Septembe 2000; last evised Octobe 2009 Abstact We summaize the aspects of the geomety of the homogeneous and isotopic spaces which ae most elevant
More informationSupplementary material for the paper Platonic Scattering Cancellation for Bending Waves on a Thin Plate. Abstract
Supplementay mateial fo the pape Platonic Scatteing Cancellation fo Bending Waves on a Thin Plate M. Fahat, 1 P.-Y. Chen, 2 H. Bağcı, 1 S. Enoch, 3 S. Guenneau, 3 and A. Alù 2 1 Division of Compute, Electical,
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationTHE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2
THE LAPLACE EQUATION The Laplace (o potential) equation is the equation whee is the Laplace opeato = 2 x 2 u = 0. in R = 2 x 2 + 2 y 2 in R 2 = 2 x 2 + 2 y 2 + 2 z 2 in R 3 The solutions u of the Laplace
More informationDesigning a Sine-Coil for Measurement of Plasma Displacements in IR-T1 Tokamak
Designing a Sine-Coil fo Measuement of Plasma Displacements in IR-T Tokamak Pejman Khoshid, M. Razavi, M. Ghoanneviss, M. Molaii, A. TalebiTahe, R. Avin, S. Mohammadi and A. NikMohammadi Dept. of Physics,
More informationA Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation
Mathematical Modelling and Analysis Publishe: Taylo&Fancis and VGTU Volume 22 Numbe 3, May 27, 3 32 http://www.tandfonline.com/tmma https://doi.og/.3846/3926292.27.39329 ISSN: 392-6292 c Vilnius Gediminas
More informationGradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function
Intenational Confeence on Infomation echnology and Management Innovation (ICIMI 05) Gadient-based Neual Netwok fo Online Solution of Lyapunov Matix Equation with Li Activation unction Shiheng Wang, Shidong
More informationCollision Frequency of Adsorbed Particles
Bulg. J. Phys. 40 (2013) 214 218 Collision Fequency of Adsobed Paticles N.S. Peev Geogi Nadjakov Institute of Solid State Physics, Bulgaian Academy of Sciences, 72 Tzaigadsko Chaussee Blvd., 1784 Sofia,
More informationDuality between Statical and Kinematical Engineering Systems
Pape 00, Civil-Comp Ltd., Stiling, Scotland Poceedings of the Sixth Intenational Confeence on Computational Stuctues Technology, B.H.V. Topping and Z. Bittna (Editos), Civil-Comp Pess, Stiling, Scotland.
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationField emission of Electrons from Negatively Charged Cylindrical Particles with Nonlinear Screening in a Dusty Plasma
Reseach & Reviews: Jounal of Pue and Applied Physics Field emission of Electons fom Negatively Chaged Cylindical Paticles with Nonlinea Sceening in a Dusty Plasma Gyan Pakash* Amity School of Engineeing
More informationCOUPLED MODELS OF ROLLING, SLIDING AND WHIRLING FRICTION
ENOC 008 Saint Petesbug Russia June 30-July 4 008 COUPLED MODELS OF ROLLING SLIDING AND WHIRLING FRICTION Alexey Kieenkov Ins ti tu te fo P ob le ms in Me ch an ic s Ru ss ia n Ac ad em y of Sc ie nc es
More informationDonnishJournals
DonnishJounals 041-1189 Donnish Jounal of Educational Reseach and Reviews. Vol 1(1) pp. 01-017 Novembe, 014. http:///dje Copyight 014 Donnish Jounals Oiginal Reseach Pape Vecto Analysis Using MAXIMA Savaş
More informationDIRECT INTERBAND LIGHT ABSORPTION IN A SPHERICAL QUANTUM DOT WITH THE MODIFIED PÖSCHEL-TELLER POTENTIAL
Lase Physics Intenational Jounal of Moden Physics: Confeence Seies Vol. 5 () 4 Wold Scientific Publishing Company DOI:.4/S945767 DIRECT INTERBAND LIGHT ABSORPTION IN A SPHERICAL QANTM DOT WITH THE MODIFIED
More informationCOORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT
COORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT Link to: phsicspages home page. To leave a comment o epot an eo, please use the auilia blog. Refeence: d Inveno, Ra, Intoducing Einstein s Relativit
More informationNumerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.
Numeical solution of diffusion mass tansfe model in adsoption systems Pof., D.Sc. Agenda Mass Tansfe Mechanisms Diffusion Mass Tansfe Models Solving Diffusion Mass Tansfe Models Paamete Estimation 2 Mass
More informationAnnihilation of Relativistic Positrons in Single Crystal with production of One Photon
Annihilation of Relativistic Positons in Single Cystal with poduction of One Photon Kalashnikov N.P.,Mazu E.A.,Olczak A.S. National Reseach Nuclea Univesity MEPhI (Moscow Engineeing Physics Institute),
More informationCentral Force Motion
Cental Foce Motion Cental Foce Poblem Find the motion of two bodies inteacting via a cental foce. Examples: Gavitational foce (Keple poblem): m1m F 1, ( ) =! G ˆ Linea estoing foce: F 1, ( ) =! k ˆ Two
More informationChapter 13 Gravitation
Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department. Problem Set 10 Solutions. r s
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Depatment Physics 8.033 Decembe 5, 003 Poblem Set 10 Solutions Poblem 1 M s y x test paticle The figue above depicts the geomety of the poblem. The position
More informationMULTILAYER PERCEPTRONS
Last updated: Nov 26, 2012 MULTILAYER PERCEPTRONS Outline 2 Combining Linea Classifies Leaning Paametes Outline 3 Combining Linea Classifies Leaning Paametes Implementing Logical Relations 4 AND and OR
More informationImplicit Constraint Enforcement for Rigid Body Dynamic Simulation
Implicit Constaint Enfocement fo Rigid Body Dynamic Simulation Min Hong 1, Samuel Welch, John app, and Min-Hyung Choi 3 1 Division of Compute Science and Engineeing, Soonchunhyang Univesity, 646 Eupnae-i
More informationMultiple Experts with Binary Features
Multiple Expets with Binay Featues Ye Jin & Lingen Zhang Decembe 9, 2010 1 Intoduction Ou intuition fo the poect comes fom the pape Supevised Leaning fom Multiple Expets: Whom to tust when eveyone lies
More informationI. CONSTRUCTION OF THE GREEN S FUNCTION
I. CONSTRUCTION OF THE GREEN S FUNCTION The Helmohltz equation in 4 dimensions is 4 + k G 4 x, x = δ 4 x x. In this equation, G is the Geen s function and 4 efes to the dimensionality. In the vey end,
More informationNumerical Inversion of the Abel Integral Equation using Homotopy Perturbation Method
Numeical Invesion of the Abel Integal Equation using Homotopy Petubation Method Sunil Kuma and Om P Singh Depatment of Applied Mathematics Institute of Technology Banaas Hindu Univesity Vaanasi -15 India
More informationPhysics 2B Chapter 22 Notes - Magnetic Field Spring 2018
Physics B Chapte Notes - Magnetic Field Sping 018 Magnetic Field fom a Long Staight Cuent-Caying Wie In Chapte 11 we looked at Isaac Newton s Law of Gavitation, which established that a gavitational field
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationStress, Cauchy s equation and the Navier-Stokes equations
Chapte 3 Stess, Cauchy s equation and the Navie-Stokes equations 3. The concept of taction/stess Conside the volume of fluid shown in the left half of Fig. 3.. The volume of fluid is subjected to distibuted
More informationF(r) = r f (r) 4.8. Central forces The most interesting problems in classical mechanics are about central forces.
4.8. Cental foces The most inteesting poblems in classical mechanics ae about cental foces. Definition of a cental foce: (i) the diection of the foce F() is paallel o antipaallel to ; in othe wods, fo
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationPES 3950/PHYS 6950: Homework Assignment 6
PES 3950/PHYS 6950: Homewok Assignment 6 Handed out: Monday Apil 7 Due in: Wednesday May 6, at the stat of class at 3:05 pm shap Show all woking and easoning to eceive full points. Question 1 [5 points]
More informationMath 209 Assignment 9 Solutions
Math 9 Assignment 9 olutions 1. Evaluate 4y + 1 d whee is the fist octant pat of y x cut out by x + y + z 1. olution We need a paametic epesentation of the suface. (x, z). Now detemine the nomal vecto:
More information11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.
Fall 2007 Qualifie Pat II 12 minute questions 11) A thin, unifom od of mass M is suppoted by two vetical stings, as shown below. Find the tension in the emaining sting immediately afte one of the stings
More informationJ. Electrical Systems 1-3 (2005): Regular paper
K. Saii D. Rahem S. Saii A Miaoui Regula pape Coupled Analytical-Finite Element Methods fo Linea Electomagnetic Actuato Analysis JES Jounal of Electical Systems In this pape, a linea electomagnetic actuato
More informationGreen s Identities and Green s Functions
LECTURE 7 Geen s Identities and Geen s Functions Let us ecall The ivegence Theoem in n-dimensions Theoem 7 Let F : R n R n be a vecto field ove R n that is of class C on some closed, connected, simply
More informationr cos, and y r sin with the origin of coordinate system located at
Lectue 3-3 Kinematics of Rotation Duing ou peious lectues we hae consideed diffeent examples of motion in one and seeal dimensions. But in each case the moing object was consideed as a paticle-like object,
More informationA dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids
APCOM & ISCM 11-14 th Decembe, 013, Singapoe A dual-ecipocity bounday element method fo axisymmetic themoelastodynamic defomations in functionally gaded solids *W. T. Ang and B. I. Yun Division of Engineeing
More informationScattering in Three Dimensions
Scatteing in Thee Dimensions Scatteing expeiments ae an impotant souce of infomation about quantum systems, anging in enegy fom vey low enegy chemical eactions to the highest possible enegies at the LHC.
More informationOn the Poisson Approximation to the Negative Hypergeometric Distribution
BULLETIN of the Malaysian Mathematical Sciences Society http://mathusmmy/bulletin Bull Malays Math Sci Soc (2) 34(2) (2011), 331 336 On the Poisson Appoximation to the Negative Hypegeometic Distibution
More informationFREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS
Please cite this aticle as: Izabela Zamosa Fee tansvese vibations of non-unifom beams Scientific Reseach of the Institute of Mathematics and Compute Science Volume 9 Issue pages 3-9. The website: http://www.amcm.pcz.pl/
More informationMean Curvature and Shape Operator of Slant Immersions in a Sasakian Space Form
Mean Cuvatue and Shape Opeato of Slant Immesions in a Sasakian Space Fom Muck Main Tipathi, Jean-Sic Kim and Son-Be Kim Abstact Fo submanifolds, in a Sasakian space fom, which ae tangential to the stuctue
More informationON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},
ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability
More informationOn absence of solutions of a semi-linear elliptic equation with biharmonic operator in the exterior of a ball
Tansactions of NAS of Azebaijan, Issue Mathematics, 36, 63-69 016. Seies of Physical-Technical and Mathematical Sciences. On absence of solutions of a semi-linea elliptic euation with bihamonic opeato
More informationClassical Worm algorithms (WA)
Classical Wom algoithms (WA) WA was oiginally intoduced fo quantum statistical models by Pokof ev, Svistunov and Tupitsyn (997), and late genealized to classical models by Pokof ev and Svistunov (200).
More information763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012
763620SS STATISTICAL PHYSICS Solutions 2 Autumn 2012 1. Continuous Random Walk Conside a continuous one-dimensional andom walk. Let w(s i ds i be the pobability that the length of the i th displacement
More informationAsymptotically Lacunary Statistical Equivalent Sequence Spaces Defined by Ideal Convergence and an Orlicz Function
"Science Stays Tue Hee" Jounal of Mathematics and Statistical Science, 335-35 Science Signpost Publishing Asymptotically Lacunay Statistical Equivalent Sequence Spaces Defined by Ideal Convegence and an
More informationPROBLEM SET #3A. A = Ω 2r 2 2 Ω 1r 2 1 r2 2 r2 1
PROBLEM SET #3A AST242 Figue 1. Two concentic co-axial cylindes each otating at a diffeent angula otation ate. A viscous fluid lies between the two cylindes. 1. Couette Flow A viscous fluid lies in the
More informationThe Poisson bracket and magnetic monopoles
FYST420 Advanced electodynamics Olli Aleksante Koskivaaa Final poject ollikoskivaaa@gmail.com The Poisson backet and magnetic monopoles Abstact: In this wok magnetic monopoles ae studied using the Poisson
More informationA matrix method based on the Fibonacci polynomials to the generalized pantograph equations with functional arguments
A mati method based on the Fibonacci polynomials to the genealized pantogaph equations with functional aguments Ayşe Betül Koç*,a, Musa Çama b, Aydın Kunaz a * Coespondence: aysebetuloc @ selcu.edu.t a
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationEquilibrium Reconstruction with Grad-Shafranov Equation
Equilibium Reconstuction with Gad-Shafanov Equation Qian Peng (qp215) Apil 14, 211 1 Gad-Shafanov equation [1]Fo axial symmetic case, the tooidal B eld is B to = F φ. Whee 2π F (, θ) is the total poloidal
More informationHydroelastic Analysis of a 1900 TEU Container Ship Using Finite Element and Boundary Element Methods
TEAM 2007, Sept. 10-13, 2007,Yokohama, Japan Hydoelastic Analysis of a 1900 TEU Containe Ship Using Finite Element and Bounday Element Methods Ahmet Egin 1)*, Levent Kaydıhan 2) and Bahadı Uğulu 3) 1)
More informationQUASI-STATIONARY ELECTRON STATES IN SPHERICAL ANTI-DOT WITH DONOR IMPURITY * 1. INTRODUCTION
ATOMIC PHYSICS QUASI-STATIONARY ELECTRON STATES IN SPHERICAL ANTI-DOT ITH DONOR IMPURITY * V. HOLOVATSKY, O. MAKHANETS, I. FRANKIV Chenivtsi National Univesity, Chenivtsi, 581, Ukaine, E-mail: ktf@chnu.edu.ua
More informationHammerstein Model Identification Based On Instrumental Variable and Least Square Methods
Intenational Jounal of Emeging Tends & Technology in Compute Science (IJETTCS) Volume 2, Issue, Januay Febuay 23 ISSN 2278-6856 Hammestein Model Identification Based On Instumental Vaiable and Least Squae
More informationSensor and Simulation Notes. Note 525. Oct Lens Design for a Prolate-Spheroidal Impulse radiating Antenna (IRA)
Senso and Simulation Notes Note 55 Oct 7 Lens Design fo a Polate-Spheoidal Impulse adiating Antenna (IRA) Sehat Altunc, Cal E. Baum, Chistos G. Chistodoulou and Edl Schamiloglu Univesity of New Mexico
More informationPhysics 111 Lecture 5 (Walker: 3.3-6) Vectors & Vector Math Motion Vectors Sept. 11, 2009
Physics 111 Lectue 5 (Walke: 3.3-6) Vectos & Vecto Math Motion Vectos Sept. 11, 2009 Quiz Monday - Chap. 2 1 Resolving a vecto into x-component & y- component: Pola Coodinates Catesian Coodinates x y =
More information