arxiv:astro-ph/ v2 15 Jul 2007
|
|
- Flora McKenzie
- 5 years ago
- Views:
Transcription
1 Draft version September 27, 2018 Preprint typeset using L TEX style emulateapj v. 10/09/06 CORONL HETING, WEK MHD TURBULENCE ND SCLING LWS. F. Rappazzo 1,2, M. Velli 2,3, G. Einaudi 1 and R. B. Dahlburg 4 Draft version September 27, 2018 arxiv:astro-ph/ v2 15 Jul 2007 BSTRCT Long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry are carried out, within the framework of reduced magnetohydrodynamicsrmhd), to understand coronal heating driven by motion of field lines anchored in the photosphere. We unambiguously identify MHD anisotropic turbulence as the physical mechanism responsible for the transport of energy from the large scales, where energy is injected by photospheric motions, to the small scales, where it is dissipated. s the loop parameters vary different regimes of turbulence develop: strong turbulence is found for weak axial magnetic fields and long loops, leading to Kolmogorov-like spectra in the perpendicular direction, while weaker and weaker regimes steeper spectral slopes of total energy) are found for strong axial magnetic fields and short loops. s a consequence we predict that the scaling of the heating rate with axial magnetic field intensity B 0, which depends on the spectral index of total energy for given loop parameters, must vary from B 3/2 0 for weak fields to B0 2 for strong fields at a given aspect ratio. The predicted heating rate is within the lower range of observed active region and quiet Sun coronal energy losses. Subject headings: Sun: corona Sun: magnetic fields turbulence 1. INTRODUCTION In this letter we solve, within the frameworkof RMHD in cartesian geometry, the Parker field-line tangling coronal heating) problem Parker 1972, 1988). We do this via long simulations at high resolutions, introducing hyper-resistivity models to attain extremely large Reynolds numbers. We show how small scales form and how the coronal heating rate depends on the loop and photospheric driving parameters, and derive simple formulae which may be used in the coronal heating context for other stars. Over the years a number of numerical experiments have been carried out to investigate coronal heating, with particular emphasis on exploring how photospheric field line tangling leads to current sheet formation. Mikic et al. 1989) and Hendrix & Van Hoven 1996) first carried out simulations of a loop driven by photospheric motions using a cartesian approximation a straightened out loop bounded at each end by the photosphere) imposing a time-dependent alternate direction flow pattern at the boundaries. complex coronal magnetic field results from the photospheric field line random walk, and though the field does not, strictly speaking, evolve through a sequence of static force-free equilibrium states the original Parker hypothesis), magnetic energy nonetheless tends to dominate kinetic energy in the system. In this limit the field is structured by current sheets elongated along the axial direction, separating quasi-2d flux tubes which constantly move around and interact. Longcope & Sudan 1994) focused on the current sheet formation process within the RMHD approximation, also used in the simulations 1 Dipartimento di Fisica E. Fermi, Università di Pisa, Pisa, Italy; rappazzo@jpl.nasa.gov 2 Jet Propulsion Laboratory, Pasadena, C, Dipartimento di stronomia e Scienza dello Spazio, Università di Firenze, Firenze, Italy 4 LCP&FD, Naval Research Laboratory, Washington, DC by Dmitruk & Gómez 1999). The results from these studies agreed qualitatively among themselves, in that all simulations display the development of field aligned current sheets. However, estimates of the dissipated power and its scaling characteristics differed largely, depending on the way in which extrapolations from low to large values of the plasma conductivity of the properties such as inertial range power law indices were carried out. 2D numerical simulations of incompressible MHD with magnetic forcing Einaudi et al. 1996; Georgoulis et al. 1998; Dmitruk et al. 1998; Einaudi & Velli 1999) showed that turbulent current sheets dissipation is distributed intermittently, and that the statistics of dissipation events, in terms of total energy, peak energy and event duration displays power laws not unlike the distribution of observed emission events in optical, ultraviolet and x-ray wavelengths of the quiet solar corona. More recently full 3D sections of the solar corona with a realistic geometry have been simulated by Gudiksen & Nordlund 2005). While this approach has advantages when investigating the coronal loop dynamics within its neighboring coronal region, modeling numerically a larger part of the solar corona drastically reduces the number of points occupied by the coronal loops. Thus, these simulations have not been able to shed further light on the physical mechanism responsible for the coronal heating. In 2 we introduce the coronal loop model and the simulations we have carried out; in 3 we describe our numerical results, and in 4 we give simple scaling arguments to understand the magnetic energy spectral slopes. This will lead to a quantitative asymptotic estimate of the coronal loop heating rate, and of its scaling with the axial magnetic field, photospheric velocity amplitude and coronal loop length. 2. THE MODEL coronal loop is a closed magnetic structure threaded by a strong axial field, with the footpoints rooted in
2 2 Rappazzo et al. Fig. 1. : High-resolution simulation with /u ph = 200, 512x512x200 grid points and R 1 = 800. Magnetic E M ) and kinetic E K ) energies as a function of time τ = L/ is the axial lfvénic crossing time). the photosphere. This makes it a strongly anisotropic system, as measured by the relative magnitude of the lfvén velocity 1000 kms 1 compared to the typical photospheric velocity u ph 1 kms 1. This means that the relative amplitude of the lfvén waves that are launched into the corona is very small. The loop dynamics may be studied in a simplified geometry, neglecting any curvature effect, as a straightened out cartesian box, with an orthogonal square cross section of size l, andanaxiallengthlembeddedin anaxialhomogeneous uniform magnetic field B 0 = B 0 e z. This system may be described by the reduced MHD RMHD) equations Kadomtsev & Pogutse 1974; Strauss 1976): introducing the velocity and magnetic field potentials ϕ and ψ, u = ϕe z ), b = ψe z ), and vorticity and current, ω = 2 ϕ, j = 2 ψ the non-dimensioned RMHD system is given by ψ t = v ϕ z ω t = v j z +[ϕ,ψ]+ 1)n+1 R n +[j,ψ] [ω,ϕ]+ 1)n+1 R n 2n ψ, 1) 2n ω. 2) s characteristic quantities we use the perpendicular length of the computational box l, the typical photospheric velocity u ph, and the related crossing time t = l /u ph. Theequationshavebeenrendereddimensionless using velocity units for the magnetic field the density in the loops ρ is taken to be constant) and normalizing by u ph. Then the non-dimensioned lfvén speed in eqs. 1)-2) is given by the ratio /u ph between the dimensional velocities. The Poisson bracket of two functions g and h is defined as [g,h] = x g y h y g x h, where x, y are transverse coordinates across the loop while z is the axial coordinate along the loop. simplified diffusion model is assumed and R n is the Reynolds number, with n the hyperdiffusion index dissipativity): for n = 1 ordinary diffusion is recovered. The computational box spans 0 x,y 1 and 0 z L, with L = 10, corresponding to an aspect ratio equal to 10. s boundary conditions at the photo- Fig. 2. : Same simulation of Figure 1. The integrated Poynting flux S dynamically balances the Ohmic J) and viscous Ω) dissipation. Inset shows a magnification of total dissipation and S for 200 t/τ 300. spheric surfaces z = 0, L) we impose a velocity pattern intended to mimic photospheric motions, made up of two independent large spatial scale projected convection cell flow patterns. The wave number values k excited are all those in the range 3 k 4, and the average injection wavenumber is k in RESULTS Plots of the rms magnetic and kinetic energies as a function of time, together with the dissipation due to currents, vorticity, as well as the integrated Poynting flux, are shown in Figures 1 and 2. s a result of the photospheric forcing, energy in the magnetic field first grows with time, until it dominates overthe kinetic energy by a large factor, before oscillating, chaotically, around a stationary state. Fluctuating magnetic energy E M is 35 times bigger than kinetic energy E K. The same generic features are seen in the rms current and vorticity dissipation, where however the time dependence of the signal is more strongly oscillating. The ohmic dissipation rate J is 6.5 times viscous dissipation Ω. The Poynting flux, on average, follows the current dissipation there is no accumulation of energy in the box), however a detailed examination shows that the dissipation time-series tends to lag the Poynting flux, with notable de-correlations around significant dissipation peaks. The spatial configuration of the currents which corresponds to a snapshot at a given time is displayed in Figure 5. The currents collapse into warped, torn sheets which extend almost completely along the loop. The current peaks are embedded within the 2D sheet-like structures, corresponding to an anisotropic structure for the turbulence, in agreement with previous results. dimensional analysis of eqs. 1)-2) shows that the only free nondimensional quantity is f = l /Lu ph. We fix L/l = 10 and vary the ratio of the lfvén speed to photospheric convection speed /u ph. Both runs with standard second order dissipation n = 1) as well as hyperdiffusion n = 4) have been carried out to obtain extended inertial ranges in the resulting spectra.
3 Coronal Heating 3 density and T λ is the energy transfer time at the scale λ, which is greater than the eddy turnover time τ λ λ/δz λ because of the lfvén effect Iroshnikov 1964; Kraichnan 1965). ) 2. In the classical IK case, T λ τ τλ /τ This corresponds to the fact that wave-packets interact over an lfvén crossing time with τ λ > τ ), and the collisions follow a standard random walk in energy exchange. In terms of the number of collisions N λ that a wave packet must suffer for the perturbation to build up to order unity, for IK N λ τ λ /τ ) 2. More generally, however, as the lfvén speed is increased the interaction time becomes smaller, so that turbulence becomes weaker and the number of collisions required for efficient energy transfer scales as Fig. 3. : Time-averaged total energy spectra for simulations with /u ph = 50, 200, 400, Hyperdiffusion n = 4) has been used with R 4 = , 10 20, 10 19, respectively, and a grid with points. The power spectrum of total energy in the simulation box, once a statistically stationary state has been achieved,dependsstronglyontheratio /u ph. Thiswas first found in simulations by Dmitruk et al. 2003), devoted to understanding how anisotropic regimes of MHD turbulence depend on boundary driving strength, with whom our numerical work is in broad agreement. The total energy spectrum, for values of /u ph = 50, 200, 400, 1000 is shown in Figure 3, together with fits to the inertial range power law. s /u ph increases, the spectrum steepens visibly note that the hump at the high wave-vectorvalues for the runs with large /u ph is a feature, the bottleneck effect, which is well known and documented in spectral simulations of turbulence with the hyperdiffusion used here, e.g. Falkovich1994)), with the slopes ranging from 2 to almost 3. t the same time while total energy increases, the ratio of the mean magnetic field over the axial lfvén velocity decreases, in good accordance with the theory. This steepening, which may be interpreted both as the effect of inertial line-tying of the coronal magnetic field and the progressive weakening of non-linear interactions as the magnetic field is increased, has a strong and direct bearing on the coronal heating scaling laws. 4. DISCUSSION characteristic of anisotropic MHD turbulence is that the cascade takes place mainly in the plane orthogonal to the DC magnetic guide field Shebalin et al. 1983). Consider then the anisotropic version of the Iroshnikov-Kraichnan IK) theory Sridhar & Goldreich 1994; Goldreich & Sridhar 1997). Dimensionally the energy cascade rate may be written as ρδz 2 λ /T λ, where δz λ is the rms value of the Elsässer fields z ± = u ±b at the perpendicular scale λ, where because the system is magnetically dominated δz + λ δz λ. ρ is the average N λ = τλ τ with α > 2, 3) whereαisthe scalingindex notethat α = 1corresponds to standard hydrodynamic turbulence), so that v ) 1 α λ T λ N λ τ. 4) L δz λ Integrating over the whole volume, the energy transfer rate becomes 1 ǫ l 2 L ρ δz2 λ L l 2 δz α+2 λ T L ρ λ λ α. 5) Considering the injection scale λ l, eq. 5) becomes ǫ l 2 L ρ δz2 l ρl2 Lα δz α+2 T l l α vα 1 l. 6) On the other hand the energy injection rate is given by the Poynting flux integrated across the photospheric boundaries: ǫ in = ρ dauph b. Considering that this integral is dominated by energy at the large scales, due to the characteristics of the forcing function, we can approximate it with ǫ in ρl 2 u ph δz l, 7) where the large scale component of the magnetic field can be replaced with δz l because the system is magnetically dominated. The last two equations show that the system is selforganized because both ǫ and ǫ in depend on δz l, the rms values of the fields z ± at the scale l : the internal dynamics depends on the injection of energy and the injection of energy itself depends on the internal dynamics via the boundary forcing. nother aspect of self-organization results from our simulations: the perpendicular magnetic field develops few spatial structures along the axial direction z, and in the nonlinear stage its topology substantially departs from the mapping of the boundary velocity pattern which characterizes its evolution during the linear stage. These and other features will be discussed more in depth in Rappazzo et al. 2007). In a stationary cascade the injection rate 7) is equal to the transport rate 6). Equating the two yields for the amplitude at the scale l : δz l u ph l Lu ph 8)
4 4 Rappazzo et al. Fig. 4. : The solid line showsthe exponent α/) as a function of α. Symbols show values of α corresponding to different values of /u ph, at fixed L/l = 10. Substituting this value in 6) or 7) we obtain for the energy flux ǫ l 2 ρ u 2 l ph, 9) Lu ph where = B 0 / 4πρ. This is also the dissipation rate, and hence the coronal heating scaling. The crucial parameter here is f = l /Lu ph because the scaling index α 3), upon which the strength of the stationary turbulent regime depends, must be a function of f itself. The relative amplitude of the turbulence δzl /, is a function of f, and as f increases the effect of linetying becomes stronger, decreasing the strength of turbulent interactions wave-packet collision efficiency becomes sub-diffusive) so that α increases above 2. The ratio δzl / can also be interpreted as the rms value of the Parker angle Θ P, and is given by < Θ P > δz l l L u ph ) 1. 10) This is actually an estimate of the average inclination of the magnetic field lines, while the rms value of the shear angle between neighboring field lines is at least twice that given by eq. 10), not considering that close to a current sheet an enhancement of the orthogonal magnetic field is observed which leads to a higher value for the angle). Numerical simulations determine the remaining unknown nondimensional dependence of the scaling index α on f. The power law slopes of the total energy spectra shown in Figure 3 are used to determine α. Identifying, as usual, the eddy energy with the band-integrated Fourier spectrum δz 2 λ k E k, where k l /λ, from eq. 5) we obtain E k k 3α+2 α+2, 11) where for α = 1 the 5/3 slope for the anisotropic Kolmogorov spectrum is recovered, and for α = 2 the 2 slope for the anisotropic IK case. t higher values of α correspond steeper spectral slopes up to the asymptotic value of 3. In Figure 4 we plot the values of α determined in this way, together with the resulting power dependence α/) of the amplitude 8) and of the energy flux 9) on the parameter f. The other power dependences are easily obtained from this last one, e.g. for the energy flux 9) the power of the axial lfvén speed is given by 1+α/), so that in terms of the magnetic field B 0 it scales as B 3/2 0 for weak fields and/or long loops, to B0 2 for strong fields and short loops. Dividing eq. 9) by the surface l 2 we obtain the energy flux per unit area F = ǫ /l 2. Taking for example a coronal loop 40,000 km long, with a number density of cm 3, = 2,000 kms 1 and u ph = 1 kms 1, for these parameters we can estimate a value of α/α + 1) 0.95), which models an active region loop, we obtain F ergcm 2 s 1 and a Parker angle 10) < Θ P > 4. On the other hand, for a coronal loop typical of a quiet Sun region, with a length of 100,000 km, a number density of cm 3, = 500 kms 1 and u ph = 1 kms 1, for these parameters we can estimate a value of α/α + 1) 0.7) we obtain F ergcm 2 s 1 and < Θ P > 0.9. In summary, with this paper we have shown how coronal heating rates in the Parker scenario scale with coronal loop and photospheric driving parameters, demonstrating that field line tangling can supply the coronal heating energy requirement. We also predict that there is no universal scaling with axial magnetic field intensity, a feature which can be tested by observing weak field regions on the Sun, or the atmospheres of other stars with differing levels of magnetic activity. M.V. thanks W.H. Matthaeus for useful discussions. R.B.D. is supported by NS SPTP. REFERENCES Dmitruk, P., Gómez, D. O., & DeLuca, D. D. 1998, pj, 505, 974 Dmitruk, P., & Gómez, D. O. 1999, pj, 527, L63 Dmitruk, P., Gómez, D. O., & Matthaeus, W. H. 2003, Phys. Plasmas, 10, 3584 Einaudi, G., Velli, M., Politano, H., & Pouquet,. 1996, pj, 457, L113 Einaudi, G., & Velli, M. 1999, Phys. Plasmas, 6, 4146 Falkovich, G. 1994, Phys. Fluids, 6, 1411 Georgoulis, M. K., Velli, M., & Einaudi, G., 1998, pj, 497, 957 Goldreich, P., & Sridhar, S. 1997, pj, 485, 680 Gudiksen, B. V., & Nordlund, Å. 2005, pj, 618, 1020 Hendrix, D. L., & Van Hoven, G. 1996, pj, 467, 887 Iroshnikov, P. S. 1964, Sov. stron., 7, 566 Kadomtsev, B. B., & Pogutse, O. P. 1974, Sov. J. Plasma Phys., 1, 389 Kraichnan, R. H. 1965, Phys. Fluids, 8, 1385 Longcope, D. W., & Sudan, R. N. 1994, pj, 437, 491 Mikic, Z., Schnack, D. D., & Van Hoven, G. 1989, pj, 338, 1148 Parker, E. N. 1972, pj, 174, 499 Parker, E. N. 1988, pj, 330, 474 Rappazzo,. F., Velli, M., Einaudi, G., & Dahlburg, R. B., pj, in preparation Shebalin, J. V., Matthaeus, W. H., & Montgomery, D. 1983, J. Plasma Phys., 29, 525 Sridhar, S. & Goldreich, P. 1994, pj, 432, 612 Strauss, H. R. 1976, Phys. Fluids, 19, 134
5 Coronal Heating 5 Fig. 5. : Top: side view of two isosurfaces of the squared current at a selected time for a numerical simulation with /u ph = 200, 512x512x200grid points and a Reynolds number R 1 = 800. The isosurface at the value j 2 = is represented in partially transparent yellow, while red displays the isosurface with j 2 = , well below the value of the maximum of the squared current that at this time is j 2 = N.B.: The red isosurface is always nested inside the yellow one, and appears pink in the figure. The computational box has been rescaled for an improved viewing, but the aspect ratio of the box is 10, i.e. the axial length of the box is ten times bigger than its orthogonal length. Bottom: top view of the same two isosurfaces using the same color display. The isosurfaces are extended along the axial direction, and the corresponding filling factor is small.
Date of delivery: 29 June 2011 Journal and vol/article ref: IAU Number of pages (not including this page): 7
Date of delivery: 29 June 2011 Journal and vol/article ref: IAU 1101498 Number of pages (not including this page): 7 Author queries: Typesetter queries: Non-printed material: The Physics of Sun and Star
More informationSHEAR PHOTOSPHERIC FORCING AND THE ORIGIN OF TURBULENCE IN CORONAL LOOPS
The Astrophysical Journal, 722:65-78, 2010 October 10 Preprint typeset using L A TEX style emulateapj v. 2/16/10 doi:10.1088/0004-637x/722/1/65 SHEAR PHOTOSPHERIC FORCING AND THE ORIGIN OF TURBULENCE IN
More informationRandom Walk on the Surface of the Sun
Random Walk on the Surface of the Sun Chung-Sang Ng Geophysical Institute, University of Alaska Fairbanks UAF Physics Journal Club September 10, 2010 Collaborators/Acknowledgements Amitava Bhattacharjee,
More informationFundamentals of Turbulence
Fundamentals of Turbulence Stanislav Boldyrev (University of Wisconsin - Madison) Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas What is turbulence? No exact definition.
More informationMHD turbulence in the solar corona and solar wind
MHD turbulence in the solar corona and solar wind Pablo Dmitruk Departamento de Física, FCEN, Universidad de Buenos Aires Motivations The role of MHD turbulence in several phenomena in space and solar
More informationParticle acceleration in stressed coronal magnetic fields
To be submitted to ApJ Letters Particle acceleration in stressed coronal magnetic fields R. Turkmani 1,L.Vlahos 2, K. Galsgaard 3,P.J.Cargill 1 and H. Isliker 2 ABSTRACT This letter presents an analysis
More informationMagnetohydrodynamic Turbulence
Magnetohydrodynamic Turbulence Stanislav Boldyrev (UW-Madison) Jean Carlos Perez (U. New Hampshire), Fausto Cattaneo (U. Chicago), Joanne Mason (U. Exeter, UK) Vladimir Zhdankin (UW-Madison) Konstantinos
More informationMagnetohydrodynamic Turbulence: solar wind and numerical simulations
Magnetohydrodynamic Turbulence: solar wind and numerical simulations Stanislav Boldyrev (UW-Madison) Jean Carlos Perez (U. New Hampshire) Fausto Cattaneo (U. Chicago) Joanne Mason (U. Exeter, UK) Vladimir
More informationA Constrained Tectonics Model for Coronal Heating
A Constrained Tectonics Model for Coronal Heating C. S. NG AND A. BHATTACHARJEE Center for Magnetic Self-Organization Center for Integrated Computation and Analysis of Reconnection and Turbulence Institute
More informationB.V. Gudiksen. 1. Introduction. Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della
Mem. S.A.It. Vol. 75, 282 c SAIt 2007 Memorie della À Ø Ò Ø ËÓÐ Ö ÓÖÓÒ B.V. Gudiksen Institute of Theoretical Astrophysics, University of Oslo, Norway e-mail:boris@astro.uio.no Abstract. The heating mechanism
More informationUniversità degli Studi di Pisa School of Graduate Studies Galileo Galilei Ph.D. Course in Physics XVIII Cycle
Università degli Studi di Pisa School of Graduate Studies Galileo Galilei Ph.D. Course in Physics XVIII Cycle 2003-2005 Ph.D. Thesis NONLINEAR PROCESSES IN CORONAL HEATING AND SLOW SOLAR WIND ACCELERATION
More informationTurbulent Origins of the Sun s Hot Corona and the Solar Wind
Turbulent Origins of the Sun s Hot Corona and the Solar Wind Steven R. Cranmer Harvard-Smithsonian Center for Astrophysics Turbulent Origins of the Sun s Hot Corona and the Solar Wind Outline: 1. Solar
More informationKinetic Alfvén waves in space plasmas
Kinetic Alfvén waves in space plasmas Yuriy Voitenko Belgian Institute for Space Aeronomy, Brussels, Belgium Solar-Terrestrial Center of Excellence, Space Pole, Belgium Recent results obtained in collaboration
More informationSolar Flares and Particle Acceleration
Solar Flares and Particle Acceleration Loukas Vlahos In this project many colleagues have been involved P. Cargill, H. Isliker, F. Lepreti, M. Onofri, R. Turkmani, G. Zimbardo,, M. Gkioulidou (TOSTISP
More informationAlfvén wave turbulence: new results with applications to astrophysics. Sébastien GALTIER Université Paris-Sud & Institut Universitaire de France
Alfvén wave turbulence: new results with applications to astrophysics Sébastien GALTIER Université Paris-Sud & Institut Universitaire de France 1 Recent co-workers : - Barbara Bigot (France/USA) - Ben
More informationSchool and Conference on Analytical and Computational Astrophysics
2292-3 School and Conference on Analytical and Computational Astrophysics 14-25 November - 2011 Magnetohydrodynamics in solar and space physics Daniel Osvaldo Gomez Instituto de Astronomia y Fisica del
More informationIncompressible MHD simulations
Incompressible MHD simulations Felix Spanier 1 Lehrstuhl für Astronomie Universität Würzburg Simulation methods in astrophysics Felix Spanier (Uni Würzburg) Simulation methods in astrophysics 1 / 20 Outline
More informationarxiv: v1 [astro-ph.sr] 15 May 2009
Draft version November 17, 2018 Preprint typeset using L A TEX style emulateapj v. 08/22/09 TURBULENCE IN THE SUB-ALFVÉNIC SOLAR WIND DRIVEN BY REFLECTION OF LOW-FREQUENCY ALFVÉN WAVES A. Verdini 1, M.
More informationCoronal Heating Problem
PHY 690C Project Report Coronal Heating Problem by Mani Chandra, Arnab Dhabal and Raziman T V (Y6233) (Y7081) (Y7355) Mentor: Dr. M.K. Verma 1 Contents 1 Introduction 3 2 The Coronal Heating Problem 4
More informationNatalia Tronko S.V.Nazarenko S. Galtier
IPP Garching, ESF Exploratory Workshop Natalia Tronko University of York, York Plasma Institute In collaboration with S.V.Nazarenko University of Warwick S. Galtier University of Paris XI Outline Motivations:
More informationFluctuation dynamo amplified by intermittent shear bursts
by intermittent Thanks to my collaborators: A. Busse (U. Glasgow), W.-C. Müller (TU Berlin) Dynamics Days Europe 8-12 September 2014 Mini-symposium on Nonlinear Problems in Plasma Astrophysics Introduction
More informationShell-to-shell energy transfer in magnetohydrodynamics. I. Steady state turbulence
Shell-to-shell energy transfer in magnetohydrodynamics. I. Steady state turbulence Alexandros Alexakis,* Pablo D. Mininni, and Annick Pouquet National Center for Atmospheric Research, P.O. Box 3000, Boulder,
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationTURBULENCE IN FLUIDS AND SPACE PLASMAS. Amitava Bhattacharjee Princeton Plasma Physics Laboratory, Princeton University
TURBULENCE IN FLUIDS AND SPACE PLASMAS Amitava Bhattacharjee Princeton Plasma Physics Laboratory, Princeton University What is Turbulence? Webster s 1913 Dictionary: The quality or state of being turbulent;
More informationObservable consequences
Coronal Heating through braiding of magnetic field lines Solar eclipse, 11.8.1999, Wendy Carlos & John Kern Observable consequences 3D MHD model spectral synthesis results: Doppler shifts DEM variability
More informationReduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014
Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton
More informationTurbulence models and excitation of solar oscillation modes
Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives
More informationMHD consistent cellular automata (CA) models
A&A 377, 1068 1080 (2001) DOI: 10.1051/0004-6361:20011100 c ESO 2001 Astronomy & Astrophysics MHD consistent cellular automata (CA) models II. Applications to solar flares H. Isliker 1, A. Anastasiadis
More informationDissipation Scales & Small Scale Structure
Dissipation Scales & Small Scale Structure Ellen Zweibel zweibel@astro.wisc.edu Departments of Astronomy & Physics University of Wisconsin, Madison and Center for Magnetic Self-Organization in Laboratory
More informationarxiv: v1 [astro-ph.sr] 25 Jun 2013
THE ASTROPHYSICAL JOURNAL, 2013, IN PRESS Preprint typeset using LATEX style emulateapj v. 04/17/13 THE SPATIAL AND TEMPORAL DEPENDENCE OF CORONAL HEATING BY ALFVÉN WAVE TURBULENCE M. ASGARI-TARGHI, A.
More informationThermographic analysis of turbulent non-isothermal water boundary layer
Thermographic analysis of turbulent non-isothermal water boundary layer Irina A. Znamenskaya 1,*, Ekaterina Y. Koroteeva 1, Anastasiya M. Novinskaya 1, Vladimir I. Fomichev 2 1 Faculty of Physics, Lomonosov
More informationTurbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields
Turbulent Magnetic Helicity Transport and the Rapid Growth of Large Scale Magnetic Fields Jungyeon Cho Dmitry Shapovalov MWMF Madison, Wisconsin April 2012 The Large Scale Dynamo The accumulation of magnetic
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationAlfvénic Turbulence in the Fast Solar Wind: from cradle to grave
Alfvénic Turbulence in the Fast Solar Wind: from cradle to grave, A. A. van Ballegooijen, and the UVCS/SOHO Team Harvard-Smithsonian Center for Astrophysics Alfvénic Turbulence in the Fast Solar Wind:
More informationarxiv:astro-ph/ v2 23 Oct 2001
Kolmogorov Burgers Model for Star Forming Turbulence Stanislav Boldyrev 1 Institute for Theoretical Physics, Santa Barbara, CA 93106 arxiv:astro-ph/0108300v2 23 Oct 2001 ABSTRACT The process of star formation
More information3D Reconnection of Weakly Stochastic Magnetic Field and its Implications
3D Reconnection of Weakly Stochastic Magnetic Field and its Implications Alex Lazarian Astronomy Department and Center for Magnetic Self- Organization in Astrophysical and Laboratory Plasmas Collaboration:
More informationFrom the Heliosphere into the Sun
511 th WE-Heraeus-Seminar From the Heliosphere into the Sun SailingagainsttheWind Collection of presentations Edited by Hardi Peter (peter@mps.mpg.de) Physikzentrum Bad Honnef, Germany January 31 February
More informationWaves & Turbulence in the Solar Wind: Disputed Origins & Predictions for PSP
Waves & Turbulence in the Solar Wind: Disputed Origins & Predictions for PSP Steven R. Cranmer University of Colorado Boulder, LASP A. Schiff, S. Van Kooten, C. Gilbert, L. N. Woolsey, A. A. van Ballegooijen,
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
More informationIntermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint
Intermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint Naoya Okamoto 1, Katsunori Yoshimatsu 2, Kai Schneider 3 and Marie Farge 4 1 Center for Computational Science, Nagoya University,
More informationarxiv: v1 [astro-ph.sr] 2 May 2011
THE ASTROPHYSICAL JOURNAL, 211, IN PRESS Preprint typeset using LATEX style emulateapj v. 3/7/7 HEATING OF THE SOLAR CHROMOSPHERE AND CORONA BY ALFVÉN WAVE TURBULENCE A. A. VAN BALLEGOOIJEN 1, M. ASGARI-TARGHI
More informationMHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field
MHD Simulation of Solar Chromospheric Evaporation Jets in the Oblique Coronal Magnetic Field Y. Matsui, T. Yokoyama, H. Hotta and T. Saito Department of Earth and Planetary Science, University of Tokyo,
More informationMax Planck Institut für Plasmaphysik
ASDEX Upgrade Max Planck Institut für Plasmaphysik 2D Fluid Turbulence Florian Merz Seminar on Turbulence, 08.09.05 2D turbulence? strictly speaking, there are no two-dimensional flows in nature approximately
More informationOn the decay of two-dimensional homogeneous turbulence
On the decay of two-dimensional homogeneous turbulence J. R. Chasnov a) The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong Received 19 June 1996; accepted 3 September
More informationarxiv: v1 [astro-ph.sr] 25 Jul 2016
Draft version September 9, 8 Preprint typeset using L A TEX style emulateapj v. 5// THE DYNAMICAL GENERATION OF CURRENT SHEETS IN ASTROPHYSICAL PLASMA TURBULENCE Gregory G. Howes Department of Physics
More informationMHD turbulence in the solar corona and solar wind
MHD turbulence in the solar corona and solar wind Pablo Dmitruk Departamento de Física, FCEN, Universidad de Buenos Aires Turbulence, magnetic reconnection, particle acceleration Understand the mechanisms
More informationarxiv: v1 [astro-ph.he] 6 Feb 2013
Draft version February 7, 2013 Preprint typeset using L A TEX style emulateapj v. 5/2/11 STATISTICAL ANALYSIS OF CURRENT SHEETS IN THREE-DIMENSIONAL MAGNETOHYDRODYNAMIC TURBULENCE Vladimir Zhdankin Department
More informationarxiv: v1 [astro-ph.sr] 13 Oct 2017
Draft version October 8, 2018 Typeset using L A TEX preprint style in AASTeX61 THE HEATING OF SOLAR CORONAL LOOPS BY ALFVÉN WAVE TURBULENCE A. A. van Ballegooijen, 1 M. Asgari-Targhi, 2 and A. Voss 3,
More informationPassive Scalars in Stratified Turbulence
GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, Passive Scalars in Stratified Turbulence G. Brethouwer Linné Flow Centre, KTH Mechanics, SE-100 44 Stockholm, Sweden E. Lindborg Linné Flow Centre,
More informationAstronomy. Astrophysics. Numerical modelling of 3D reconnection. II. Comparison between rotational and spinning footpoint motions
A&A 459, 627 639 (2006) DOI: 10.1051/0004-6361:20065716 c ESO 2006 Astronomy & Astrophysics Numerical modelling of 3D reconnection II. Comparison between rotational and spinning footpoint motions I. De
More informationChapter 7. Basic Turbulence
Chapter 7 Basic Turbulence The universe is a highly turbulent place, and we must understand turbulence if we want to understand a lot of what s going on. Interstellar turbulence causes the twinkling of
More informationSUPPLEMENTARY INFORMATION. Magnetic jam in the corona of the Sun. Three-dimensional coronal models compared to observational properties
Magnetic jam in the corona of the Sun F. Chen, 1 H. Peter, 1 S. Bingert, 2 M.C.M. Cheung 3 DOI:.38/NPHYS331 1 Max Planck Institute for Solar System Research, Justus-von-Liebig-Weg 3, 37077 Göttingen, Germany
More informationPhysical modeling of coronal magnetic fields and currents
Physical modeling of coronal magnetic fields and currents Participants: E. Elkina,, B. Nikutowski,, A. Otto, J. Santos (Moscow,Lindau,, Fairbanks, São José dos Campos) Goal: Forward modeling to understand
More informationWALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION
WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION Joongnyon Kim, Kyoungyoun Kim, Hyung Jin Sung Department of Mechanical Engineering, Korea Advanced Institute of Science
More informationLecture 2. Turbulent Flow
Lecture 2. Turbulent Flow Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of this turbulent water jet. If L is the size of the largest eddies, only very small
More informationReview of electron-scale current-layer dissipation in kinetic plasma turbulence
Meeting on Solar Wind Turbulence Kennebunkport, ME, June 4-7, 2013 Review of electron-scale current-layer dissipation in kinetic plasma turbulence Minping Wan University of Delaware W. H. Matthaeus, P.
More informationTheory and Modelling of Coronal Wave Heating
Theory and Modelling of Coronal Wave Heating Ineke De Moortel School of Mathematics & Statistics University of St Andrews Overview Some recent observations of [Alfvén(ic)] waves in the chromosphere and
More informationThe ECHO code: from classical MHD to GRMHD in dynamical spacetimes
The ECHO code: from classical MHD to GRMHD in dynamical spacetimes Luca Del Zanna Dipartimento di Fisica e Astronomia Università di Firenze Main collaborators: N. Bucciantini, O. Zanotti, S. Landi 09/09/2011
More informationElectron magnetohydrodynamic turbulence
PHYSICS OF PLASMAS VOLUME 6, NUMBER 3 MARCH 1999 Electron magnetohydrodynamic turbulence D. Biskamp, E. Schwarz, and A. Zeiler Max-Planck-Institut für Plasmaphysik, D-85748 Garching, Germany A. Celani
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More information6 VORTICITY DYNAMICS 41
6 VORTICITY DYNAMICS 41 6 VORTICITY DYNAMICS As mentioned in the introduction, turbulence is rotational and characterized by large uctuations in vorticity. In this section we would like to identify some
More informationThe Physics of Fluids and Plasmas
The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationCreation and destruction of magnetic fields
HAO/NCAR July 30 2007 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)
More informationScaling relations in MHD and EMHD Turbulence
Scaling relations in MHD and EMHD Turbulence Jungyeon Cho Chungnam National University, Korea Outline MHD Non-MHD E(k) MHD turb. small-scale turb. ~1/r i k Topic 1. Strong MHD Turbulence Alfven wave Suppose
More informationCoronal Heating versus Solar Wind Acceleration
SOHO 15: Coronal Heating, 6 9 September 2004, University of St. Andrews, Scotland Coronal Heating versus Solar Wind Acceleration Steven R. Cranmer Harvard-Smithsonian Center for Astrophysics, Cambridge,
More informationUNCONDITIONAL STABILITY OF A PARTITIONED IMEX METHOD FOR MAGNETOHYDRODYNAMIC FLOWS
UNCONDITIONAL STABILITY OF A PARTITIONED IMEX METHOD FOR MAGNETOHYDRODYNAMIC FLOWS CATALIN TRENCHEA Key words. magnetohydrodynamics, partitioned methods, IMEX methods, stability, Elsässer variables. Abstract.
More informationSolar Flare. A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona)
Solar Flares Solar Flare A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona) Flares release 1027-1032 ergs energy in tens of minutes. (Note: one H-bomb: 10
More informationSolar Wind Turbulent Heating by Interstellar Pickup Protons: 2-Component Model
Solar Wind Turbulent Heating by Interstellar Pickup Protons: 2-Component Model Philip A. Isenberg a, Sean Oughton b, Charles W. Smith a and William H. Matthaeus c a Inst. for Study of Earth, Oceans and
More informationJet Stability: A computational survey
Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,
More informationStatistical studies of turbulent flows: self-similarity, intermittency, and structure visualization
Statistical studies of turbulent flows: self-similarity, intermittency, and structure visualization P.D. Mininni Departamento de Física, FCEyN, UBA and CONICET, Argentina and National Center for Atmospheric
More information1 TH/P8-43 Role of Impurity Cyclotron Damping in Ion Heating and RFP Turbulence
1 Role of Impurity Cyclotron Damping in Ion Heating and RFP Turbulence P.W. Terry, V. Tangri, J.S. Sarff, G. Fiksel, A.F. Almagri, Y. Ren, and S.C. Prager Department of Physics, University of Wisconsin-Madison,
More informationarxiv: v1 [physics.flu-dyn] 7 Jul 2015
Statistical equilibria of large scales in dissipative hydrodynamic turbulence V. Dallas, S. Fauve, and A. Alexakis Laboratoire de Physique Statistique, École Normale Supérieure, CNRS, Université Pierre
More informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More informationKonvektion und solares Magnetfeld
Vorlesung Physik des Sonnensystems Univ. Göttingen, 2. Juni 2008 Konvektion und solares Magnetfeld Manfred Schüssler Max-Planck Planck-Institut für Sonnensystemforschung Katlenburg-Lindau Convection &
More informationApplied Computational Fluid Dynamics
Lecture 9 - Kolmogorov s Theory Applied Computational Fluid Dynamics Instructor: André Bakker André Bakker (2002-2005) Fluent Inc. (2002) 1 Eddy size Kolmogorov s theory describes how energy is transferred
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationA solar surface dynamo
MPS Solar Group Seminar May 8, 2007 A solar surface dynamo Alexander Vögler (Univ. of Utrecht) & Manfred Schüssler A lot of magnetic flux in the `quiet Sun Observation: Flux replenishment rates increase
More informationEffect of current sheets on the power spectrum of the solar wind magnetic field using a cell model
Available online at www.sciencedirect.com Advances in Space Research 49 (2012) 1327 1332 www.elsevier.com/locate/asr Effect of current sheets on the power spectrum of the solar wind magnetic field using
More informationMixing in Highly Compressible Turbulence
Mixing in Highly Compressible Turbulence Liubin Pan Evan Scannapieco School of Earth and Space Exploration Arizona State University Astrophysical motivation: Heavy elements from supernova and stellar winds
More informationNonequilibrium Dynamics in Astrophysics and Material Science YITP, Kyoto
Nonequilibrium Dynamics in Astrophysics and Material Science 2011-11-02 @ YITP, Kyoto Multi-scale coherent structures and their role in the Richardson cascade of turbulence Susumu Goto (Okayama Univ.)
More information3D hybrid-kinetic turbulence and phase-space cascades
3D hybrid-kinetic turbulence and phase-space cascades ( in a β = 1 plasma ) Silvio Sergio Cerri Department of Astrophysical Sciences, Princeton University, USA 11th Plasma Kinetics Working Meeting WPI
More informationTurbulence and Reconnection
Turbulence and Reconnection Jeff Tessein July 10, 2011 NASA turbulence study at Wallops Island, Virginia Outline Turbulence (length scales, Reynolds decomposition) Navier-Stokes Equation Turbulence Spectrum
More informationHigh-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows
TH/P3-3 High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows A. Bierwage 1), S. Benkadda 2), M. Wakatani 1), S. Hamaguchi 3), Q. Yu
More informationMODELLING TWISTED FLUX TUBES PHILIP BRADSHAW (ASTROPHYSICS)
MODELLING TWISTED FLUX TUBES PHILIP BRADSHAW (ASTROPHYSICS) Abstract: Twisted flux tubes are important features in the Universe and are involved in the storage and release of magnetic energy. Therefore
More informationCreation and destruction of magnetic fields
HAO/NCAR July 20 2011 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)
More informationTurbulent dissipation in the solar wind and corona
Turbulent dissipation in the solar wind and corona W. H. Matthaeus, P. Dmitruk, S. Oughton and D. Mullan Bartol Research Institute, University of Delaware, Newark, DE 19716 USA Department of Mathematics,
More informationScaling laws of free magnetic energy stored in a solar emerging flux region
Publ. Astron. Soc. Japan 2014 66 (4), L6 (1 5) doi: 10.1093/pasj/psu049 Advance Access Publication Date: 2014 July 14 Letter L6-1 Letter Scaling laws of free magnetic energy stored in a solar emerging
More informationCascade Phenomenology in Turbulence: Navier-Stokes and MHD
Cascade Phenomenology in Turbulence: Navier-Stoes and MHD WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the
More informationSolar Wind Turbulence
Solar Wind Turbulence Presentation to the Solar and Heliospheric Survey Panel W H Matthaeus Bartol Research Institute, University of Delaware 2 June 2001 Overview Context and SH Themes Scientific status
More information(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time
Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,
More informationGyrokinetic Theory and Dynamics of the Tokamak Edge
ASDEX Upgrade Gyrokinetic Theory and Dynamics of the Tokamak Edge B. Scott Max Planck Institut für Plasmaphysik D-85748 Garching, Germany PET-15, Sep 2015 these slides: basic processes in the dynamics
More informationVII. Hydrodynamic theory of stellar winds
VII. Hydrodynamic theory of stellar winds observations winds exist everywhere in the HRD hydrodynamic theory needed to describe stellar atmospheres with winds Unified Model Atmospheres: - based on the
More informationarxiv:physics/ v1 [physics.plasm-ph] 5 Nov 2004
Ion Resonance Instability in the ELTRAP electron plasma G. Bettega, 1 F. Cavaliere, 2 M. Cavenago, 3 A. Illiberi, 1 R. Pozzoli, 1 and M. Romé 1 1 INFM Milano Università, INFN Sezione di Milano, Dipartimento
More informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
More informationMagnetic Effects Change Our View of the Heliosheath
Magnetic Effects Change Our View of the Heliosheath M. Opher Λ, P. C. Liewer Λ, M. Velli, T. I. Gombosi ΛΛ, W.Manchester ΛΛ,D. L. DeZeeuw ΛΛ,G.Toth ΛΛ and I. Sokolov ΛΛ Λ Jet Propulsion Laboratory, MS
More informationObservations and Modeling of Turbulence in the Solar Wind
Observations and Modeling of Turbulence in the Solar Wind Melvyn L. Goldstein NASA Goddard Space Flight Center, USA E-mail: melvyn.l.goldstein@nasa.gov Summary. Alfvénic fluctuations are a ubiquitous component
More informationWhat can test particles tell us about magnetic reconnection in the solar corona?
What can test particles tell us about magnetic reconnection in the solar corona? James Threlfall, T. Neukirch, C. E. Parnell, A. W. Hood jwt9@st-andrews.ac.uk @JamesWThrelfall Overview Motivation (solar
More informationPrediction of solar activity cycles by assimilating sunspot data into a dynamo model
Solar and Stellar Variability: Impact on Earth and Planets Proceedings IAU Symposium No. 264, 2009 A. G. Kosovichev, A. H. Andrei & J.-P. Rozelot, eds. c International Astronomical Union 2010 doi:10.1017/s1743921309992638
More informationPoS(Extremesky 2011)056
Stochastic acceleration and the evolution of spectral distributions in SSC sources: A self consistent modeling of blazars flares ISDC, University of Geneva, Chemin d Ecogia 16, Versoix, CH-1290, Switzerland
More information