Non-linear Cosmic Ray propagation close to the acceleration site
|
|
- Roland Chambers
- 5 years ago
- Views:
Transcription
1 Non-linear Cosmic Ray propagation close to the acceleration site The Hebrew University of Jerusalem, Jerusalem, 91904, Israel Stefano Gabici Astroparticule et Cosmologie (APC), CNRS, Université Paris 7 Denis Dierot, Paris, France stefano.gabici@apc.univ-paris7.fr Alexanre Marcowith Laboratoire Univers et particules e Montpellier, Université Montpellier/CNRS, Montpellier, France alexanre.marcowith@umontpellier.fr Giovanni Morlino INFN Gran Sasso Science Institute, viale F. Crispi 7, L Aquila, Italy giovanni.morlino@gssi.infn.it Vlaimir Ptuskin Pushkov Institute of Terrestrial Magnetism, Ionosphere an Raio Wave Propagation of the Russian Acaemy of Sciences, Troitsk, Moscow, Russia vptuskin@izmiran.ru Recent avances on gamma-ray observations from SuperNova Remnants an Molecular Clous offer the possibility to stuy in etail the properties of the propagation of escaping Cosmic Rays (CR). However, a complete theory for CR transport outsie the acceleration site has not been evelope yet. Two physical processes are thought to be relevant to regulate the transport: the growth of waves cause by streaming instability, an possible wave amping mechanisms that reuce the growth of the turbulence. Only a few attempts have been mae so far to incorporate these mechanisms in the theory of CR iffusion. In this work we present recent avances in this subject. In particular, we show results obtaine by solving the couple equations for the iffusion of CRs an the evolution of Alfvén waves. We iscuss the importance of streaming instabilities an wave amping in ifferent ISM phases. The 34th International Cosmic Ray Conference, 30 July- 6 August, 2015 The Hague, The Netherlans Speaker. c Copyright owne by the author(s) uner the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence.
2 1. Introuction We stuy the non-linear iffusion of a population of Cosmic Rays (CRs) after they escape from the acceleration site (for example, a supernova remnant shock). As a zero orer approximation the problem can be escribe by an impulsive injection of CRs of a given energy at a given location in the Galaxy. We stuy the transport of CRs taking into account the growth of Alfvén waves cause by particle-wave interactions, an their back-reaction on the CR iffusion. We also consier possible linear wave amping mechanisms that limit the growth of the waves. The problem is escribe by two ifferential equations that govern the evolution of the CR pressure an of the wave energy ensity as a function of time an istance from the acceleration site. We present numerical solutions of these two couple equations for ifferent values of the physical parameters involve in the problem. 2. The moel In our moel the transport of CRs is assume to be regulate by the resonant scattering off Alfvén waves, i.e. a CR of energy E resonates with waves of wave number k= 1/r L (E) where r L is the particle Larmor raius (see e.g. [1]). The (normalize) energy ensity I(k) of Alfvén waves is efine as: δb 2 8π = B2 0 I(k)ln k, (2.1) 8π where B 0 is the ambient magnetic fiel an δb the amplitue of the magnetic fiel fluctuations. Accoring to quasi linear theory, CRs iffuse along the magnetic fiel lines with a iffusion coefficient equal to [2]: D= 4 c r L(E) 3π I(k) = D B(E), (2.2) I(k) which can be expresse as the ratio between the Bohm iffusion coefficient D B (E) an the energy ensity of resonant waves I(k = 1/r L ). Formally, quasi linear theory is vali for δb/b 0 1. In this limit the iffusion perpenicular to the fiel lines can be neglecte, being suppresse by a factor of(δb k /B 0 ) 4 I(k) 2 with respect to the one parallel to B 0 (see e.g. [3, 4]). This implies that uner the conition that δb/b 0 remains small, the problem is one imensional. We consier the streaming of CR to be the main source of Alfvénic turbulence. For CRs that stream along the irection of B 0 (that we consier aligne along the coorinate z) the growth rate of Alfvén waves Γ CR is proportional to the prouct between the Alfvén spee V A an the graient of the pressure of resonant CRs, an can be efine as (see e.g. [3, 5]): 2Γ CR I = V A z. (2.3) The sign inicates that only waves traveling in the irection of the streaming are excite. For imensional consistency, the pressure of CRs P CR has to be normalize to the magnetic energy ensity B 2 0 /8π. On scales smaller than the magnetic fiel coherence length, a flux tube characterize by a constant magnetic fiel strength B 0 an irecte along the z axis can be consiere. Two couple equations (for the evolution of CRs an waves along the flux tube) must be solve: the iffusion 2
3 coefficient of CRs of energy E (Eq. 2.2), inee, epens on the energy ensity of resonant waves I(k), an in turn the growth rate of these waves (Eq. 2.3) epens on the graient of the pressure of resonant CRs. The two couple equations then rea: t +V A z = z ( DB I ) z (2.4) I t +V I A z = V A z 2Γ I+ Q (2.5) where the left sie in both expressions is the time erivative compute along the characteristic of excite waves: t = t +V A z. (2.6) The avective terms V A / z an V A I/ z are neglecte in the following since we verifie that they play little role in the situation uner examination. The last two terms in Eq. 2.5 account for possible mechanisms of wave amping, operating at a rate Γ, an for the injection Q of turbulence from an external source (i.e. other than CR streaming). In this work we consier only linear amping mechanisms, for which the amping rate is inepenent of space an time. The term representing the external source of turbulence can be set to Q = 2Γ I 0, so that when streaming instability is not relevant the level of the backgroun turbulence is at a constant level I = I 0. As alreay pointe out by [6], the approach escribe above ecouples the process of acceleration of particles, which operates, for example, at a SNR shock, from the particle escape from the acceleration site. Though such a separation might seem artificial, the problem efine above can still be useful to escribe the escape of particles form a ea accelerator, in which the acceleration mechanism oes not operate any more, or operates at a much reuce efficiency [6]. This situation woul probably apply to the case of ol SNRs. On the other han, [7] suggeste that Equations. 2.4 an 2.5 coul be also use to escribe an intermeiate phase of CR propagation in which the CRs have left the source but are not yet completely mixe with the CR backgroun. For the case of a SNR shock, the equations above coul thus be applie to those CRs characterize by a iffusion length large enough to ecouple them from the shock region. Typically, this happens uring the Seov phase to the highest energy particles confine at the SNR shock when the iffusion length D B /u s graually increases with time up to a value larger than some fraction χ of the SNR shock raius R s, where u s is the shock spee an χ [8, 9]. In both scenarios, the initial conitions for Eqs 2.4 an 2.5 can be set as follows: P CR = P 0 CR an I 1 for z<a (2.7) P CR = 0 an I = I 0 1 for z>a (2.8) where a is the spatial scale of the region fille by CRs at the time of their escape from the source. Following [6], we introuce the quantity Π, efine as: where: Π = V A D B Φ CR (2.9) Φ CR z P CR = PCR 0 a. (2.10) 0 3
4 To unerstan the physical meaning of Π, consier the initial setup of the problem, in which CRs are localize in a region of size a. The CR pressure within a is PCR 0, an then Φ CR = PCR 0 a. The initial iffusion coefficient outsie the region of size a is equal to D B /I 0. The problem of transport is characterise by three relevant timescales: i) the growth time τ g (V A /I 0 / z) 1 ai 0 /V A PCR 0, ii) the time it takes the CR clou to sprea ue to iffusion τ i f f a 2 /D a 2 I 0 /D B, an iii) the amping time τ = 1/2Γ. In orer to have a significant growth of waves ue to CR streaming, the timescale for growth must be shorter than the two other timescales: τ g < min(τ i f f,τ ). In terms of the parameter Π this conitions reas: Π>max(1,τ i f f /τ ). It is evient, then, that the parameter Π regulates the effectiveness of CR inuce growth of waves. In the absence of a amping term for waves, Π > 1 is a necessary conition for streaming instability to be relevant, while in the presence of efficient wave amping, a more stringent conition on Π applies. For Π max(1,τ i f f /τ ) CRs play no role in the generation of Alfvén waves, an Equation 2.4 can be solve analytically: P CR = I0 I 0 z 2 πd B t Φ CR e 4D B t (2.11) Equation 2.11 is referre to as test particle solution. When wave growth cannot be neglecte, the solution eviates from this analytic test-particle solution. In the extreme scenario Π >> max(1,τ i f f /τ ) waves grow so quickly that the the iffusive term in Eq. 2.4 becomes negligible when compare to avection term V A / z, an the avection term cannot be neglecte anymore. This escribes a situation in which CRs are "locke" to waves an move with them at a velocity equal to V A [10]. An ientical result was foun by [11] in a stuy of the escape of MeV CRs from sources, an also suggeste by [12]. 3. The metho In orer to solve Eqns. 2.4 an 2.5 it is convenient to perform the following change of coorinates [6]: s z ( ) VA τ t (3.1) a a an re-normalise the parameters as follow: ( ) VA a P CR P CR D B ( ) VA a W I D B ( ) a Γ Γ (3.2) V A In these notations (an after ropping the avection terms) the equations become: = ( ) 1 τ s W s (3.3) W τ = 2Γ (W W 0 ) s (3.4) Their solution (i.e., P an W as a function of the variables τ an s) epens only on three parameters: the initial values PCR 0, W 0, an Γ 0. Note that: i) P0 CR = Π, ii) W 0 = V A a/d 0, an iii) Γ = Γ 0 (since we consier here only a linear amping constant in both time an space). In these notations, the conition for effective growth of waves is Π>max(1,Γ W 0 ). 4
5 Figure 1: Left: ion-neutral amping rate (Γ IN, soli curves), Farmer & Golreich amping rate (ΓFG, ashe curves) an ratio between Alfvén spee an initial size of the source (V A /a, ot-ashe curves), as a function of the CR energy E. Different colours refer to two ifferent ISM phases: warm ionise meium (WIM, re), an warm neutral meium (WNM, blue). Right: parameter Π as a function of the CR energy for a WIM an WNM (soli lines). In absence of efficient wave amping, values of Π<1 correspon to a test particle solutions, since the wave growth rate is inefficient. Conversely, for Π larger that Π max 10 the growth rate is large an the iffusion coefficient becomes smaller than the Bohm iffusion coefficient D B. In presence of efficient wave amping, both these limits are shifte to higher values, epening on the value of Γ. We present examples of numerical solutions to the couple equations 3.3 an 3.4 for values of the parameters in the ranges: Π = [2 20], W 0 = [ ], an Γ = [0 100]. In orer to unerstan why we chose these values, an to which physical conitions they correspon, it is necessary to specify the values of V A,a,D 0, D B, an the amping mechanisms. These quantities in general epen on the properties of the ISM, an on the particle energy. We consier two ifferent ISM phases: a warm neutral meium (WNM) an a warm ionise meium (WIM). The typical values of ensity, temperature, ionisation fraction, an magnetic fiel strength for these two phases are taken from [13]. In these ISM phases, the two most relevant linear amping mechanisms are the ion-neutral friction (Γ IN ) an the amping by backgroun MHD turbulence suggeste by Farmer & Golreich [14] (Γ FG ). The value of ΓFG as a function of the particle energy is shown in Fig. 1 (left panel, long-ashe lines), for both ISM phases. The ion-neutral amping rate is estimate by numerically solving the equations in [15]. The results are shown in Fig. 1 (left panel, soli lines), for both ISM phases: at low energies Γ IN = ν IN/2 (where ν IN is the ion-neutral collision frequency), while at high energies Γ IN E 2. The level of coupling between ions an neutrals affects also the Alfvén spee, that is in the range cm/s. Since it is useful for our calculations, we show in Fig. 1 the value of the ratio V a /a, as a function of energy (ash-otte lines). In orer to estimate a we consier a R esc (E), where R esc (E) is the escape raius an epens on the CR energy. We moel its epenence on the energy following [16]: we foun that its value ranges between a few pc (for the highest energy particles) an 20 pc (for GeV particles). Since we are intereste in following the CR iffusion from the scale a where they are release up to a istance z 10 3 pc, this roughly correspons to solve the normalise equations in a range of s from 1 to several hunres. Moreover, we want to stuy the CR pressure at ifferent times t in the range 10 3 yr to 10 5 yr. Consiering the value of V a /a as a function of energy isplaye in Fig. 1, this 5
6 correspons to τ = In orer to ientify the relevant values of the parameters Π,W 0,Γ we procee as follow: Γ the amping rate Γ = max(γ IN,ΓFG ) ranges between 10 9 an s 1 (Fig. 1, left panel), epening on the particle energy an ISM phase. This correspons to 1<Γ Π the values of Π as a function of the CR energy, for the two ifferent ISM is shown in Fig. 1 (right panel, soli lines). W 0 for the backgroun iffusion D 0 we can refer to the average Galactic value D 0 = D ISM = (E/10GeV) 0.5 cm 2 /s. Since V A a ranges between cm 2 /s an cm 2 /s (at low an high energies, respectively), the range of relevant values for W 0 is 10 5 < W 0 < Results an Discussion We numerically solve Eqs. 3.3 an 3.4 an show the results for P(s,τ) vs. s (soli lines in Fig. 2) at three ifferent normalise times τ = 10 2,10 1,1 (corresponing to ifferent colours in the figure). The test particle solution is also shown for comparison (otte lines). Each row in the figure correspons to results obtaine by varying only one parameter at the time. The upper row shows the impact of the parameter Π for fixe values of Γ = 10 an W 0 = Since for Π < 1 the numerical solution coincies with the test particle solution, we show our results for values of Π>1. When Π=2, the solution starts to eviate from the test particle case: ue to the growth of wave by streaming instability the iffusion is slower, the CR pressure in the vicinity of the accelerator is larger, an the iffusion length is smaller. The effect is more an more evient for increasing values of Π (upper row, mile an right panels). The mile row shows the effect of wave amping, for fixe Π = 10 an W 0 = Since Π >> 1, when wave amping is negligible (Γ =0, left panel) the iffusion is suppresse as compare to the test particle case. Larger values of Γ limit the growth of Alfvén waves an the solution approaches the test particle solution. Note that Γ =10 (mile row, central panel) correspons to a amping mechanism that becomes important at τ > 0.1. For this reason the solution at τ = 0.01 an τ = 0.1 are unaffecte by the increase amping rate, an only the solution at late time τ = 1 is moifie, an approaches the test particle solution. The right panel shows the case of a even faster amping mechanism, Γ =100, whose effect starts to be relevant at times corresponing to τ > 0.01, explaining why the solution at intermeiate an late times approaches the test particle solution, while the solution at early time is unaffecte by wave amping. The bottom row shows calculations performe at Π=5 an Γ =10, with W 0 varying from 10 5 to A large W 0 (left panel) correspons to slow iffusion an small iffusion length, both in the test particle case an in the numerical solution. The two solutions iffer at early an intermeiate times because the growth rate is faster than the iffusion rate, an suppresses the iffusion. At later times instea the growth rate ecreases an the numerical solution oes not appreciably iffer from the test particle case. A change in the initial value of the backgroun turbulence has the same effect on both the analytical an numerical solution: the iffusion is faster for ecreasing values of W 0 (see bottom row, from left to right), an the CRs can reach larger istances. The eviation from a 6
7 P CR P CR P CR Figure 2: Normalize CR pressure as a function of the variable s, which represents the istance from the source z normalise to the size of the source a. In each row the results are shown by varying only one of the three parameters (whose value is reporte above each panel) an keeping the other two fixe to the values reporte in the shae area. The soli curves show the results from the numerical computation, while otte curves represent the test particle solution, where the role of streaming instabilities in enhancing the turbulence is neglecte. Different colours refer to ifferent normalise times τ = 10 2,10 1,1. 7
8 test particle solution however remains unchange, since both the growth rate an the iffusion rate scale as W 0. This stuy outlines the importance of consiering, in the physics of CR transport, both the growth of waves ue to streaming instability an the role of mechanisms that amp the turbulence. When amping is neglecte, even relatively small values of the integrate CR pressure Π = 10 (corresponing, for example, to particles of a few TeV, in a WIM - see Fig. 1) in a environment with small backgroun turbulence W 0 = 10 4 (relevant for TeV particles when the backgroun iffusion coefficient D 0 is similar to the average Galactic iffusion coefficient) lea to the conclusion that iffusion is strongly suppresse by the efficient growth of waves, an particles cannot travel large istances, even at late times. This situation is shown in Fig. 2 (mile row, left panel) where the numerical computation is performe in absence of a wave amping term: the iffusion is strongly suppresse also at τ = 1, corresponing to ol remnants (t age 10 5 yr). References [1] D. G. Wentzel, Cosmic-ray propagation in the Galaxy - Collective effects, ARA&A 12 (1974) 71. [2] A. R. Bell, The acceleration of cosmic rays in shock fronts. I, MNRAS 182 (1978) 147. [3] L. Drury, An introuction to the theory of iffusive shock acceleration of energetic particles in tenuous plasmas, Rep. Prog. Phys. 46 (1983) 973. [4] F. Casse, M. Lemoine, G. Pelletier, Transport of cosmic rays in chaotic magnetic fiels, Phys. Rev. D 65 (2002) [5] J. Skilling, Cosmic ray streaming. III - Self-consistent solutions, MNRAS 173 (1975) 255. [6] M. A. Malkov, P. H. Diamon, R. Z. Sageev, F. A. Aharonian, I. V. Moskalenko, Analytic Solution for Self-regulate Collective Escape of Cosmic Rays from Their Acceleration Sites, ApJ 768 (2013) 73. [7] V. S. Ptuskin, V. N. Zirakashvili, A. A. Plesser, Non-linear iffusion of cosmic rays, Av. Space Res. 42 (2008) 486. [8] V. S. Ptuskin, V. N. Zirakashvili, On the spectrum of high-energy cosmic rays prouce by supernova remnants in the presence of strong cosmic-ray streaming instability an wave issipation, A&A 429 (2005) 755. [9] S. Gabici, Cosmic ray escape from supernova remnants, MmSAI 82 (2011) 760. [10] J. Skilling, Cosmic Rays in the Galaxy: Convection or Diffusion?, MNRAS 170 (1971) 265. [11] C. J. Cesarsky, Interstellar propagation of low energy cosmic rays, ICRC 2 (1975) 634. [12] T. W. Hartquist, G. E. Morfill, Cosmic ray iffusion at energies of 1 MeV to 105 GeV, Astrophys. Space Sci. 216 (1994) 223. [13] P. Jean, W. Gillar, A. Marcowith, K. Ferriere, Positron transport in the interstellar meium, A&A 508 (2009) [14] A. J. Farmer, P. Golreich, Wave Damping by Magnetohyroynamic Turbulence an Its Effect on Cosmic-Ray Propagation in the Interstellar Meium, ApJ 604 (2004) 671. [15] E. G. Zweibel, J. M. Shull, Confinement of cosmic rays in molecular clous, ApJ 259 (1982) 859. [16] P. Cristofari, S. Gabici, S. Casanova, R. Terrier, E. Parizot, Acceleration of cosmic rays an gamma-ray emission from supernova remnants in the Galaxy, MNRAS 434 (2013)
Cosmic ray escape from supernova remnants
Mem. S.A.It. Vol. 82, 760 c SAIt 2011 Memorie della Cosmic ray escape from supernova remnants Stefano Gabici Astroparticule et Cosmologie (APC), CNRS, Université Paris 7 Denis Diderot, Paris, France, e-mail:
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationDiffusive shock acceleration with regular electric fields
Diffusive shock acceleration with regular electric fields V.N.Zirakashvili Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radiowave Propagation, Russian Academy of Sciences (IZMIRAN), 142190
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationPhysics 2212 GJ Quiz #4 Solutions Fall 2015
Physics 2212 GJ Quiz #4 Solutions Fall 215 I. (17 points) The magnetic fiel at point P ue to a current through the wire is 5. µt into the page. The curve portion of the wire is a semicircle of raius 2.
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationHigh energy radiation from molecular clouds (illuminated by a supernova remnant
High energy radiation from molecular clouds (illuminated by a supernova remnant A. Marcowith (L.P.T.A. Montpellier) collaboration with S. Gabici (D.I.A.S.) 1 Outlook Introduction: Scientific interests.
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationInteraction force in a vertical dust chain inside a glass box
Interaction force in a vertical ust chain insie a glass box Jie Kong, Ke Qiao, Lorin S. Matthews an Truell W. Hye Center for Astrophysics, Space Physics, an Engineering Research (CASPER) Baylor University
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationSupernova remnants in the very high energy sky: prospects for the Cherenkov Telescope Array
Supernova remnants in the very high energy sky: prospects for the Cherenkov Telescope Array Columbia University E-mail: pc2781@columbia.edu Roberta Zanin Max-Planck-Institut fur Kernphysik Stefano Gabici
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationSolution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010
NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationGamma rays from supernova remnants in clumpy environments.! Stefano Gabici APC, Paris
Gamma rays from supernova remnants in clumpy environments!! Stefano Gabici APC, Paris Overview of the talk Galactic cosmic rays Gamma rays from supernova remnants Hadronic or leptonic? The role of gas
More informationProduction of Secondary Cosmic Rays in Supernova Remnants
Production of Secondary Cosmic Rays in Supernova Remnants E. G. Berezhko, Yu. G. Shafer Institute of Cosmophysical Research and Aeronomy, 31 Lenin Ave., 677891 Yakutsk, Russia E-mail: ksenofon@ikfia.sbras.ru
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationThe spectacular stellar explosion - GRB A: synchrotron modeling in the wind and the ISM
The spectacular stellar explosion - GRB 17A: synchrotron moeling in the win an the ISM University of Johannesburg, Department of Physics, Aucklan Park 6, Johannesburg, South Africa E-mail: jessymolkt@uj.ac.za
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More informationGamma ray emission from supernova remnant/molecular cloud associations
Gamma ray emission from supernova remnant/molecular cloud associations Stefano Gabici APC, Paris stefano.gabici@apc.univ-paris7.fr The Origin of galactic Cosmic Rays Facts: the spectrum is (ALMOST) a single
More informationu t v t v t c a u t b a v t u t v t b a
Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationThe role of ionization in the shock acceleration theory
The role of ionization in the shock acceleration theory Giovanni Morlino INAF - L.go E. Fermi 5, Firenze, Italy E-mail: morlino@arcetri.astro.it We study the acceleration of heavy nuclei at SNR shocks
More information3-dimensional Evolution of an Emerging Flux Tube in the Sun. T. Magara
3-imensional Evolution of an Emerging Flux Tube in the Sun T. Magara (Montana State University) February 6, 2002 Introuction of the stuy Dynamical evolution of emerging fiel lines Physical process working
More informationResearch Article Relativistic Propagation of Linearly/Circularly Polarized Laser Radiation in Plasmas
ISRN Optics Volume 213, Article ID 642617, 8 pages http://x.oi.org/1.1155/213/642617 Research Article Relativistic Propagation of Linearly/Circularly Polarize Laser Raiation in Plasmas Sonu Sen, 1 Meenu
More information3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects
IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationarxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration
Prouction of pions, kaons an protons in Xe Xe collisions at s =. ev arxiv:09.0v [hep-ex] Sep 0, for the ALICE Collaboration Università i Bologna an INFN (Bologna) E-mail: simone.ragoni@cern.ch In late
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationOn Cosmic-Ray Production Efficiency at Realistic Supernova Remnant Shocks
On Cosmic-Ray Production Efficiency at Realistic Supernova Remnant Shocks, 1 T. Inoue 2, Y. Ohira 1, R. Yamazaki 1, A. Bamba 1 and J. Vink 3 1 Department of Physics and Mathematics, Aoyama-Gakuin University,
More informationCosmic-ray Acceleration and Current-Driven Instabilities
Cosmic-ray Acceleration and Current-Driven Instabilities B. Reville Max-Planck-Institut für Kernphysik, Heidelberg Sep 17 2009, KITP J.G. Kirk, P. Duffy, S.O Sullivan, Y. Ohira, F. Takahara Outline Analysis
More informationDusty Plasma Void Dynamics in Unmoving and Moving Flows
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationEffect of the Regular Galactic Magnetic Field on the Propagation of Galactic Cosmic Rays in the Galaxy
Effect of the Regular Galactic Magnetic Field on the Propagation of Galactic Cosmic Rays in the Galaxy Department of Electrical and Electronic Systems Engineering, National Institute of Technology, Ibaraki
More information1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity
AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationEvaporating droplets tracking by holographic high speed video in turbulent flow
Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationCosmic-ray acceleration by compressive plasma fluctuations in supernova shells
Cosmic-ray acceleration by compressive plasma fluctuations in supernova shells Ming Zhang Florida Institute of Technology E-mail: mzhang@fit.edu We suggest that the production of Galactic cosmic rays in
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationGravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann
Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)
More informationMicrowave Reflection from the Region of Electron Cyclotron Resonance Heating in the L-2M Stellarator )
Microwave Reflection from the Region of Electron Cyclotron Resonance Heating in the L-2M Stellarator German M. BATANOV, Valentin D. BORZOSEKOV, Nikolay K. KHARCHEV, Leoni V. KOLIK, Eugeny M. KONCHEKOV,
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationarxiv: v1 [nucl-th] 21 Jan 2019
Paramagnetic squeezing of a uniformly expaning quark-gluon plasma in an out of equilibrium N. Saooghi an S. M. A. Tabatabaee Department of Physics, Sharif University of Technology, P.O. Box 11155-9161,
More informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationarxiv:physics/ v4 [physics.class-ph] 9 Jul 1999
AIAA-99-2144 PROPULSION THROUGH ELECTROMAGNETIC SELF-SUSTAINED ACCELERATION arxiv:physics/9906059v4 [physics.class-ph] 9 Jul 1999 Abstract As is known the repulsion of the volume elements of an uniformly
More informationMATHEMATICS BONUS FILES for faculty and students
MATHMATI BONU FIL for faculty an stuents http://www.onu.eu/~mcaragiu1/bonus_files.html RIVD: May 15, 9 PUBLIHD: May 5, 9 toffel 1 Maxwell s quations through the Major Vector Theorems Joshua toffel Department
More informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationChapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies
hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation
More informationA Model of Electron-Positron Pair Formation
Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More informationQuantum optics of a Bose-Einstein condensate coupled to a quantized light field
PHYSICAL REVIEW A VOLUME 60, NUMBER 2 AUGUST 1999 Quantum optics of a Bose-Einstein conensate couple to a quantize light fiel M. G. Moore, O. Zobay, an P. Meystre Optical Sciences Center an Department
More informationLecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationPhysics 2212 K Quiz #2 Solutions Summer 2016
Physics 1 K Quiz # Solutions Summer 016 I. (18 points) A positron has the same mass as an electron, but has opposite charge. Consier a positron an an electron at rest, separate by a istance = 1.0 nm. What
More informationPCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue
PAPER View Article Online View Journal View Issue Cite this: Phys. Chem. Chem. Phys., 2014, 16, 22321 Receive 4th June 2014, Accepte 3r September 2014 DOI: 10.1039/c4cp02467b www.rsc.org/pccp 1 Introuction
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationarxiv: v1 [astro-ph.he] 15 Jan 2019
Preprint 16 January 2019 Compiled using MNRAS LATEX style file v3.0 Can a cosmic ray carrot explain the ionization level in diffuse molecular clouds? S. Recchia, 1 V. H. M. Phan, 1 S. Biswas 2 S. Gabici
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationVertical shear plus horizontal stretching as a route to mixing
Vertical shear plus horizontal stretching as a route to mixing Peter H. Haynes Department of Applie Mathematics an Theoretical Physics, University of Cambrige, UK Abstract. The combine effect of vertical
More informationMagnetic helicity evolution in a periodic domain with imposed field
PHYSICAL REVIEW E 69, 056407 (2004) Magnetic helicity evolution in a perioic omain with impose fiel Axel Branenburg* Norita, Blegamsvej 17, DK-2100 Copenhagen Ø, Denmark William H. Matthaeus University
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationChapter 11: Feedback and PID Control Theory
Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by
More informationNuclear Physics and Astrophysics
Nuclear Physics an Astrophysics PHY-302 Dr. E. Rizvi Lecture 2 - Introuction Notation Nuclies A Nuclie is a particular an is esignate by the following notation: A CN = Atomic Number (no. of Protons) A
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationRFSS: Lecture 4 Alpha Decay
RFSS: Lecture 4 Alpha Decay Reaings Nuclear an Raiochemistry: Chapter 3 Moern Nuclear Chemistry: Chapter 7 Energetics of Alpha Decay Geiger Nuttall base theory Theory of Alpha Decay Hinrance Factors Different
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationAn Approach for Design of Multi-element USBL Systems
An Approach for Design of Multi-element USBL Systems MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological University of the Mixteca Carretera a Acatlima Km. 2.5 Huajuapan e Leon Oaxaca 69000
More informationThe Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis
The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into
More informationNON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS
ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More information(3-3) = (Gauss s law) (3-6)
tatic Electric Fiels Electrostatics is the stuy of the effects of electric charges at rest, an the static electric fiels, which are cause by stationary electric charges. In the euctive approach, few funamental
More informationON THE MEANING OF LORENTZ COVARIANCE
Founations of Physics Letters 17 (2004) pp. 479 496. ON THE MEANING OF LORENTZ COVARIANCE László E. Szabó Theoretical Physics Research Group of the Hungarian Acaemy of Sciences Department of History an
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationRETROGRADE WAVES IN THE COCHLEA
August 7, 28 18:2 WSPC - Proceeings Trim Size: 9.75in x 6.5in retro wave 1 RETROGRADE WAVES IN THE COCHLEA S. T. NEELY Boys Town National Research Hospital, Omaha, Nebraska 68131, USA E-mail: neely@boystown.org
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationProject 3 Convection
12.307 Project 3 Convection Loovica Illari an John Marshall March 2010 1 Backgroun The Earth is bathe in raiation from the Sun whose intensity peaks in the visible. In orer to maintain energy balance the
More informationSources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences
Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationPhysics 115C Homework 4
Physics 115C Homework 4 Problem 1 a In the Heisenberg picture, the ynamical equation is the Heisenberg equation of motion: for any operator Q H, we have Q H = 1 t i [Q H,H]+ Q H t where the partial erivative
More information