• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter iResearch App Facebook

Search Issue | RSS Feeds RSS
Previous Issue

Dec 1981

Volume 24, Issue 12, pp. 2133-2380

Page 1 of 2 Pages Next Page | Jump to Page

Modeling of vortex‐corner interaction using point vortices

A. T. Conlisk and D. Rockwell

Phys. Fluids 24, 2133 (1981); http://dx.doi.org/10.1063/1.863327 (10 pages) | Cited 8 times

Full Text: | Download PDF

Show Abstract
Using point vortices, the interaction of a single vortex, as well as patterns of vortices, with a corner is examined; comparisons are made with corresponding experiments. Trajectories of a vortex swept past the corner can be well‐approximated, provided that a sufficiently weak strength of the vortex is specified. Calculations show that the dimensionless amplitude of the pressure fluctuation at the corner is very sensitive to small variations of the initial position of the vortex, and relatively insensitive to variations in strength; this finding has important consequences for recently observed amplitude modulation of pressure at, and velocity near, impingement. In fact, by prescribing transversely staggered patterns of vortices upstream of the corner, determined from experimental flow visualization, the form of the time‐averaged pressure and velocity spectra can be approximated. The most critical feature of these spectra, a well‐defined low‐frequency component(s), confirms the hypothesized mechanism associated with low‐frequency modulation observed in experiments. In addition, shortcomings of this method, primarily due to distributed vorticity inherent in laboratory vortices, are pointed out.
Show PACS
47.32.Ef Rotating and swirling flows
47.55.dp Cavitation and boiling
47.10.-g General theory in fluid dynamics

Shape of shock wave produced by a concentrated impact on a surface

Gerald Nutt and Lewis Klein

Phys. Fluids 24, 2143 (1981); http://dx.doi.org/10.1063/1.863328 (7 pages) | Cited 1 time

Full Text: | Download PDF

Show Abstract
An approximate similarity solution, derived by Raizer, of a concentrated impact (or intense explosion) at the boundary of a semi‐infinite volume of a perfect gas is used to determine the propagation velocity of the shock front as a function of its position. This velocity function is then used to obtain the shape of the propagating shock wave. It is shown that dish‐shaped shock fronts are formed when the movement of the gas at the surface is into the gas region and that cup‐shaped shock fronts are formed when the movement is out of the gas region. Comparison of these results with the shapes of explosions and meteorite craters are discussed.
Show PACS
47.40.-x Compressible flows; shock waves
41.20.Jb Electromagnetic wave propagation; radiowave propagation

Explosions on a gas‐vacuum interface

Gerald Nutt, Lewis Klein, and Albert E. Ratcliffe

Phys. Fluids 24, 2150 (1981); http://dx.doi.org/10.1063/1.863329 (4 pages) | Cited 1 time

Full Text: | Download PDF

Show Abstract
A finite‐difference computer code is used to calculate the time development of an explosion on a gas‐vacuum interface. An analytic theory of the shape of the shock wave produced in the explosion is compared with the results of the computer simulation. The assumptions used in obtaining this analytic solution are verified, and the degree to which the variables describing the explosion are self‐similar is examined. Finally, certain consistency relations among the similarity exponents are tested.
Show PACS
47.40.-x Compressible flows; shock waves
41.20.Jb Electromagnetic wave propagation; radiowave propagation

Strong evaporation in half‐space: Integral transport solutions for one‐dimensional Bhatnagar–Gross–Krook model

S. K. Loyalka

Phys. Fluids 24, 2154 (1981); http://dx.doi.org/10.1063/1.863330 (5 pages) | Cited 4 times

Full Text: | Download PDF

Show Abstract
The strong evaporation problem in half‐space for the one‐dimensional Bhatnagar–Gross–Krook model is considered, and it is shown that highly accurate results in the lowest approximation can be obtained by a judicious use of two different integral transport methods. Extensions of the work to higher order kinetic models is discussed, and some comments regarding obtaining accurate results for the Boltzmann equation by use of the generalized Maxwell method are also included.
Show PACS
68.03.Fg Evaporation and condensation of liquids

Rarefaction shocks in spherical geometry

Dushan Mitrovich

Phys. Fluids 24, 2159 (1981); http://dx.doi.org/10.1063/1.863331 (4 pages) | Cited 3 times

Full Text: | Download PDF

Show Abstract
The uniform density region that develops behind a rarefaction shock in planar geometry (and attaches downstream to the outer expansion flow) is shown to be modified by geometric divergence into a region of self‐similar flow with nonvanishing gradients. A single set of equations is derived that describes both shock and similarity flow. The longest gradient scale length in the transition region is shown to vary linearly with the radius at the shock.
Show PACS
52.30.-q Plasma dynamics and flow
52.35.Tc Shock waves and discontinuities

Kinetic theory of the plasma sheath transition in a weakly ionized plasma

K.‐U. Riemann

Phys. Fluids 24, 2163 (1981); http://dx.doi.org/10.1063/1.863332 (10 pages) | Cited 98 times

Full Text: | Download PDF

Show Abstract
The boundary layer of a weakly ionized rare gas discharge is treated using kinetic theory. A self‐consistent two‐scale analysis is performed and the exact presheath and sheath solutions are constructed for the following model: A plasma consisting of Boltzmann distributed electrons and singly charged atomic ions is in contact with a negative absorbing wall. The ion kinetics is dominated by charge exchange with cold neutrals (T0T). The mean‐free‐path λ is constant and large compared with the electron Debye length λD. This model represents a collision‐dominated counterpart to the well known Tonks–Langmuir model of the collision‐free plasma.
Show PACS
52.40.Kh Plasma sheaths
52.25.Dg Plasma kinetic equations
52.65.-y Plasma simulation
52.20.Fs Electron collisions

Electron motion in a wave of slowly varying amplitude

G. Pocobelli

Phys. Fluids 24, 2173 (1981); http://dx.doi.org/10.1063/1.863333 (4 pages) | Cited 8 times

Full Text: | Download PDF

Show Abstract
A theory of electron motion in a spatially sinusoidal wave of slowly varying amplitude is presented. On the basis of a reformulation of classical conservative motion, a set of approximate analytical, adiabatic solutions for the case of a temporally varying amplitude is derived. The set can describe electron transitions from trapped to untrapped conditions, or vice versa, and provides a formally unified description of both kinds of motion through Jacobi’s real transformation.
Show PACS
52.35.Mw Nonlinear phenomena: waves, wave propagation, and other interactions (including parametric effects, mode coupling, ponderomotive effects, etc.)
52.35.Fp Electrostatic waves and oscillations (e.g., ion-acoustic waves)
52.20.Dq Particle orbits
52.25.Dg Plasma kinetic equations

Damping of an electron plasma wave with detrapping of the electrons

G. Pocobelli

Phys. Fluids 24, 2177 (1981); http://dx.doi.org/10.1063/1.863334 (6 pages) | Cited 9 times

Full Text: | Download PDF

Show Abstract
A self‐consistent calculation of the damping of an electron plasma wave is presented including, for the first time, the contribution of all electrons (untrapped, detrapped, and trapped). The effect of the detrapped electrons, which is strong at the higher damping rates, is to lower the first minimum of the amplitude after the initial stage of Landau damping and the general level of its oscillation afterward. At γLτ0≊−0.75, a qualitative change is observed from amplitude oscillation to near‐plateau behavior. The existence of an important secularity in the Vlasov–Poisson system is pointed out.
Show PACS
52.35.Fp Electrostatic waves and oscillations (e.g., ion-acoustic waves)
52.35.Mw Nonlinear phenomena: waves, wave propagation, and other interactions (including parametric effects, mode coupling, ponderomotive effects, etc.)
52.25.Dg Plasma kinetic equations

Universal drift waves in a low‐aspect‐ratio spheromak geometry

R. Marchand and P. N. Guzdar

Phys. Fluids 24, 2183 (1981); http://dx.doi.org/10.1063/1.863335 (3 pages) | Cited 2 times

Full Text: | Download PDF

Show Abstract
An eigenmode equation for drift waves in low‐aspect‐ratio spheromak geometry is derived and solved numerically. In contrast to the large‐aspect‐ratio tokamak limit, the strong equilibrium variation involved here, allows partially localized modes at more than one value of the poloidal angle.
Show PACS
52.35.Kt Drift waves

Dielectric tensor of a weakly relativistic, nonequilibrium, and magnetized plasma

S. T. Tsai, C. S. Wu, Y. D. Wang, and S. W. Kang

Phys. Fluids 24, 2186 (1981); http://dx.doi.org/10.1063/1.863324 (5 pages) | Cited 26 times

Full Text: | Download PDF

Show Abstract
The dielectric tensor of a weakly relativistic magnetized plasma is discussed for wave frequencies very close to the fundamental or higher harmonics of the electron or ion cyclotron frequency. Since in real situations the plasmas of interest are often in nonequilibrium states, the present paper generalizes Shkarofsky’s result to include three nonequilibrium features; namely, a loss‐cone distribution, anisotropic temperatures, and a field‐aligned drift. The discussion is motivated by current interest in the study of absorption and emission of radiation near the electron cyclotron harmonics.
Show PACS
52.27.Ny Relativistic plasmas
52.25.Os Emission, absorption, and scattering of electromagnetic radiation

Absorption and emission of extraordinary‐mode electromagnetic waves near cyclotron frequency in nonequilibrium plasmas

C. S. Wu, C. S. Lin, H. K. Wong, S. T. Tsai, and R. L. Zhou

Phys. Fluids 24, 2191 (1981); http://dx.doi.org/10.1063/1.863325 (6 pages) | Cited 21 times

Full Text: | Download PDF

Show Abstract
Two cases are investigated. In the first case the electron distribution is assumed to have a loss‐cone feature. It is found that the plasma can amplify electromagnetic waves near the electron cyclotron frequency, as expected. However, the growth rate is significantly larger than that used in an approximate calculation in which a slightly different physical model was considered. The present discussion is restricted to the fast extraordinary‐ mode radiation. One of the important conclusions is that, due to the combined relativistic and loss‐cone effects, the cutoff frequency can be lower than the electron cyclotron frequency, a result in contrast with that obtained with the usual ’’cold plasma’’ approximation. In the second case a drift motion along the ambient magnetic field is studied. It is found that the cyclotron damping can occur under certain conditions.
Show PACS
52.25.Os Emission, absorption, and scattering of electromagnetic radiation
52.40.Db Electromagnetic (nonlaser) radiation interactions with plasma

Electromagnetic radiation from an inhomogeneous plasma: Theory and experiment

R. W. Means, L. Muschietti, M. Q. Tran, and J. Vaclavik

Phys. Fluids 24, 2197 (1981); http://dx.doi.org/10.1063/1.863326 (11 pages) | Cited 27 times

Full Text: | Download PDF

Show Abstract
The problem of the linear conversion of a Langmuir wave to a transverse electromagnetic wave in the presence of a density gradient has been solved numerically with appropriate boundary conditions. A reciprocity principle was found, allowing the deduction of solutions of this problem from those obtained from the transverse‐Langmuir conversion. This model has been applied to study the spectrum emitted from an inhomogeneous plasma, including the effect of the antenna radiation pattern. Experiments have been performed in a large unmagnetized dc discharge plasma (ne∼5×1010 cm−3, Te = 1.3 eV, gradient scale length L = 100−1000 cm). The shape of the spectrum observed with a horn antenna agrees with the theoretical one, but the deduced level of Langmuir fluctuations is much higher than the thermal level. This enhancement is due to the presence of primary energetic (E = 60 eV) electrons.
Show PACS
52.25.Os Emission, absorption, and scattering of electromagnetic radiation
52.35.Fp Electrostatic waves and oscillations (e.g., ion-acoustic waves)
52.35.Hr Electromagnetic waves (e.g., electron-cyclotron, Whistler, Bernstein, upper hybrid, lower hybrid)
52.70.Gw Radio-frequency and microwave measurements

Tearing‐mode stability properties of a diffuse anisotropic field‐reversed ion layer at marginal stability

James Chen and R. C. Davidson

Phys. Fluids 24, 2208 (1981); http://dx.doi.org/10.1063/1.863336 (8 pages) | Cited 9 times

Full Text: | Download PDF

Show Abstract
Stability properties are investigated for purely growing (Reω = 0) tearing modes at marginal stability (Imω = 0) for a rotating, nonrelativistic cylindrically symmetric ion layer immersed in an axial magnetic field. The analysis is carried out within the framework of a Vlasov‐fluid model in which the electrons are described as a macroscopic cold fluid, and the layer ions are described by the Vlasov equation. Tearing‐mode stability properties are calculated numerically for azimuthally symmetric perturbations about an ion layer equilibrium with temperature anisotropy. The marginal stability eigenvalue equation is solved numerically for the perturbation amplitude and the normalized axial wavenumbers at marginal stability in terms of temperature anisotropy, layer radius, and magnetic field depression on the axis. It is found that the range of unstable wavenumbers decreases as T‖‖/T is increased, and numerical estimates are made of the anisotropy required for complete stabilization.
Show PACS
52.25.Fi Transport properties

Anomalous impurity ion transport due to magnetic fluctuations

J. Y. Hsu, R. W. Harvey, and S. K. Wong

Phys. Fluids 24, 2216 (1981); http://dx.doi.org/10.1063/1.863337 (3 pages) | Cited 5 times

Full Text: | Download PDF

Show Abstract
A positive ambipolar potential arising from large transport rates of electrons relative to ions in stochastic magnetic fields gives an outward force on the ions. It can reverse the neoclassical inward convection of the impurity ions in tokamaks. This process may by responsible for the observed decays of line emissions from high charge states of silicon in an impurity study experiment in Alcator.
Show PACS
52.25.Fi Transport properties
52.25.Gj Fluctuation and chaos phenomena
52.35.Ra Plasma turbulence

Magnetic field generation in the underdense plasma

Patrick Mora and René Pellat

Phys. Fluids 24, 2219 (1981); http://dx.doi.org/10.1063/1.863338 (8 pages) | Cited 41 times

Full Text: | Download PDF

Show Abstract
Magnetic field generation in the underdense plasma due to the nonlinear interaction of an electromagnetic wave and a plasma is investigated. It is shown that the collisional absorption of the electromagnetic wave induces a nonlinear anistropy in the electron pressure tensor, which is responsible for a large source of the dc magnetic field. Six different regimes of interest appear, depending on four characteristic parameters related to the collision rate, the thermal effects, hydrodynamic convection, and the duration of the experiment. The theory is supported by a detailed comparison with recent microwave plasma experiments.
Show PACS
52.40.Db Electromagnetic (nonlaser) radiation interactions with plasma
52.25.Os Emission, absorption, and scattering of electromagnetic radiation
52.35.Ra Plasma turbulence
52.50.Jm Plasma production and heating by laser beams (laser-foil, laser-cluster, etc.)

Ponderomotive effects and magnetic field generation due to short‐wavelength ion turbulence

Patrick Mora and René Pellat

Phys. Fluids 24, 2227 (1981); http://dx.doi.org/10.1063/1.863339 (8 pages) | Cited 7 times

Full Text: | Download PDF

Show Abstract
The ponderomotive effects and magnetic field generation due to the nonlinear interaction of short‐wavelength ion turbulence with an electromagnetic wave are studied using the Dawson and Oberman model, which exactly describes the electron‐ion interaction in the limit of infinite ion‐to‐electron mass ratio. The differences from the results of the heuristic Landau model are stressed, in particular for nonthermal ion correlations and for an electromagnetic wave frequency slightly exceeding the electron plasma frequency.
Show PACS
52.40.Db Electromagnetic (nonlaser) radiation interactions with plasma
52.25.Os Emission, absorption, and scattering of electromagnetic radiation
52.35.Ra Plasma turbulence
52.50.Jm Plasma production and heating by laser beams (laser-foil, laser-cluster, etc.)

Quasi‐linear saturation of coherent backward ionization modes

A. Thieltgen, G. Leclert, and J. R. Cussenot

Phys. Fluids 24, 2235 (1981); http://dx.doi.org/10.1063/1.863340 (4 pages) | Cited 1 time

Full Text: | Download PDF

Show Abstract
Quasi‐linear saturation of the amplitude of coherent backward ionization waves is studied using the slowly varying parameters or Van der Pol’s method. The result is a nonlinear Landau equation which describes the growth and saturation of the wave amplitude. The analytical expression of the most important nonlinear saturation term predicts results in accordance with earlier experimental observations.
Show PACS
52.35.Mw Nonlinear phenomena: waves, wave propagation, and other interactions (including parametric effects, mode coupling, ponderomotive effects, etc.)

Necessary stability conditions for field‐reversed theta pinches

J. R. Cary

Phys. Fluids 24, 2239 (1981); http://dx.doi.org/10.1063/1.863341 (3 pages) | Cited 14 times

Full Text: | Download PDF

Show Abstract
Toroidal systems of arbitrary cross section without toroidal magnetic field are analyzed via the double adiabatic fluid equations. Such systems are shown to be unstable if there exists one closed field line on which the average of κrB2 is positive, where κ is the curvature. A similar criterion is derived for linear system and is applied to a noncircular z pinch.
Show PACS
52.35.Bj Magnetohydrodynamic waves (e.g., Alfven waves)
52.35.Py Macroinstabilities (hydromagnetic, e.g., kink, fire-hose, mirror, ballooning, tearing, trapped-particle, flute, Rayleigh-Taylor, etc.)

Finite Larmor radius effect on thermosolutal instability of a plasma

R. C. Sharma and K. N. Sharma

Phys. Fluids 24, 2242 (1981); http://dx.doi.org/10.1063/1.863342 (3 pages) | Cited 6 times

Full Text: | Download PDF

Show Abstract
The thermosolutal instability of a plasma in the presence of a transverse magnetic field is studied to include the effects of the finiteness of the ion Larmor radius. The finite Larmor radius and stable solute gradient are found to have stabilizing effects on stationary convection. The case of overstable modes is also considered wherein necessary conditions for the existence of overstability are derived.
Show PACS
52.35.Py Macroinstabilities (hydromagnetic, e.g., kink, fire-hose, mirror, ballooning, tearing, trapped-particle, flute, Rayleigh-Taylor, etc.)

Variational structure of the Vlasov equation

H. L. Berk, R. R. Dominguez, and E. K. Maschke

Phys. Fluids 24, 2245 (1981); http://dx.doi.org/10.1063/1.863343 (8 pages) | Cited 15 times

Full Text: | Download PDF

Show Abstract
The variational structure of the Vlasov–Maxwell integral equations is derived for a plasma equilibrium having two ignorable coordinates. It is shown that the kernel of the Maxwell equations is a self‐adjoint integral operator. This operator may also be represented as a differential equation of arbitrary order. This representation is useful when the differential operator is truncated to finite order, yielding a system of intrinsically self‐adjoint differential equations.
Show PACS
52.55.Dy General theory and basic studies of plasma lifetime, particle and heat loss, energy balance, field structure, etc.

Most probable magnetohydrodynamic tokamak and reversed field pinch equilibria

John Ambrosiano and George Vahala

Phys. Fluids 24, 2253 (1981); http://dx.doi.org/10.1063/1.863344 (12 pages) | Cited 6 times

Full Text: | Download PDF

Show Abstract
The statistical theory of Montgomery, Turner, and Vahala, which determines the most probable ideal magnetohydrodynamic equilibrium compatible with given information on only a few global parameters (e.g., energy E, magnetic helicity H, flux Φ, current I, ⋅⋅⋅) is extended and investigated for both the tokamak regime (in which experimentally Φ≫μ0aI, with a being the plasma radius) and the reversed field pinch regime (Φ≪μ0aI). One obtains typical experimentally relevant profiles in the appropriate regimes. Most probable equilibria sequences are investigated as the energy/magnetic helicity ratio is decreased at fixed flux and current: In the tokamak regime (flux≫current) the diamagnetic toroidal field Bz becomes less diamagnetic and tends to a uniform field, while in the reversed field pinch regime (flux≪current), field reversal sets in Bz with the radial reversal position moving farther into the plasma and the eventual appearance of hollow pressure profiles. It appears that, in both regimes, the most probable equilibria are becoming more stable as μ0aE/H decreases. Linearized analytic force‐free states can also be constructed for certain regimes of the global parameters together with their nonlinear quasi‐force‐free numerical counterparts.
Show PACS
52.55.Dy General theory and basic studies of plasma lifetime, particle and heat loss, energy balance, field structure, etc.

Strongly‐localized ballooning modes in an axisymmetric tandem mirror

D. A. D’Ippolito and J. R. Myra

Phys. Fluids 24, 2265 (1981); http://dx.doi.org/10.1063/1.863345 (5 pages) | Cited 7 times

Full Text: | Download PDF

Show Abstract
A simple model of an axisymmetric tandem mirror without thermal barriers is employed to study the stability of ballooning modes which are confined to the central mirror cell. The central cell is treated as long‐thin with isotropic pressure, while the end plug is modeled by a line‐tying boundary condition. The scaling of βc with the mirror ratio R and the connection length Lc is computed, and the sensitivity of these results to the magnetic field profile chosen is discussed.
Show PACS
52.55.Jd Magnetic mirrors, gas dynamic traps
52.35.Py Macroinstabilities (hydromagnetic, e.g., kink, fire-hose, mirror, ballooning, tearing, trapped-particle, flute, Rayleigh-Taylor, etc.)

Finite Larmor radius stabilization of ballooning modes in an axisymmetric tandem mirror

D. A. D’Ippolito, G. L. Francis, J. R. Myra, and W. M. Tang

Phys. Fluids 24, 2270 (1981); http://dx.doi.org/10.1063/1.863346 (4 pages) | Cited 13 times

Full Text: | Download PDF

Show Abstract
The stability of an axisymmetric tandem mirror to ballooning modes is examined including the effect of finite ion gyroradius stabilization. The central cell of the tandem mirror is treated as a long‐thin system with isotropic pressure, while the end plug is modeled by a stabilizing boundary condition. For the particular equilibrium considered, the finite Larmor radius corrections to ideal magnetohydrodynamic theory indicate a strong stabilizing effect on high–n ballooning modes even for quite modest values of the gyroradius, k ρia/L≪1. The dependence of βc on k ρi is computed, and physical mechanisms associated with this stabilizing influence are discussed.
Show PACS
52.55.Jd Magnetic mirrors, gas dynamic traps
52.35.Py Macroinstabilities (hydromagnetic, e.g., kink, fire-hose, mirror, ballooning, tearing, trapped-particle, flute, Rayleigh-Taylor, etc.)
52.55.Dy General theory and basic studies of plasma lifetime, particle and heat loss, energy balance, field structure, etc.

Low‐m ballooning stability of an axisymmetric sharp‐boundary tandem mirror

D. A. D’Ippolito and B. Hafizi

Phys. Fluids 24, 2274 (1981); http://dx.doi.org/10.1063/1.863347 (6 pages) | Cited 10 times

Full Text: | Download PDF

Show Abstract
A simple model of an axisymmetric tandem mirror is employed to study the stability of global (low‐m) ballooning modes in the central cell. The presence of a perfectly conducting wall in the model gives an estimate of the effect of wall stabilization. The scaling of the critical β ( = 2p/B2) with mode number m, connection length Lc to the end plug, and wall position w is computed for several magnetic field profiles.
Show PACS
52.55.Jd Magnetic mirrors, gas dynamic traps
52.35.Py Macroinstabilities (hydromagnetic, e.g., kink, fire-hose, mirror, ballooning, tearing, trapped-particle, flute, Rayleigh-Taylor, etc.)
52.55.Dy General theory and basic studies of plasma lifetime, particle and heat loss, energy balance, field structure, etc.

Motion of a charged particle in a nearly axisymmetric magnetic field

Harold Weitzner

Phys. Fluids 24, 2280 (1981); http://dx.doi.org/10.1063/1.863348 (15 pages) | Cited 15 times

Full Text: | Download PDF

Show Abstract
The motion of a charged particle in a static magnetic field is studied by means of repeated canonical transformations of a Hamiltonian system. Adiabatic invariants are generated based on the assumption that the particle Larmor radius is small compared with the characteristic distance over which the magnetic field varies. Unlike many earlier treatments the transformations presented here preserve the axisymmetry of the dynamics when the magnetic field is axisymmetric. It is assumed that the magnetic field consists of a small nonaxisymmetric part plus the axisymmetric toroidal and poloidal parts. After the introduction of the magnetic moment adiabatic invariant, the motion of the guiding center is studied. The results depend sensitively on the ratio of the poloidal magnetic field to the total magnetic field. In some cases a second adiabatic invariant exists and direct inferences concerning long time particle drifts are possible. In one case where a second adiabatic invariant fails to exist, long term drifts are studied by conventional perturbation expansions. At some points, resonance or lack of resonance phenomena appear and determine the drift effects.
Show PACS
52.20.Dq Particle orbits
Page 1 of 2 Pages Next Page | Jump to Page
Close
Google Calendar
ADVERTISEMENT

close