ALBERT

Beschreibung: This investigation analyzes the effect of vortex wakes on the Lagrangian displacement of particles induced by the passage of an obstacle in a two-dimensional incompressible and inviscid fluid. In addition to the trajectories of individual particles, we also study their drift and the corresponding total drift areas in the Föppl and Kirchhoff potential flow models. Our findings, which are obtained numerically and in some regimes are also supported by asymptotic analysis, are compared to the wakeless potential flow which serves as a reference. We show that in the presence of the Föppl vortex wake, some of the particles follow more complicated trajectories featuring a second loop. The appearance of an additional stagnation point in the Föppl flow is identified as a source of this effect. It is also demonstrated that, while the total drift area increases with the size of the wake for large vortex strengths, it is actually decreased for small circulation values. On the other hand, the Kirchhoff flow model is shown to have an unbounded total drift area. By providing a systematic account of the wake effects on the drift, the results of this study will allow for more accurate modeling of hydrodynamic stirring.

Beschreibung: The slow motion of a circular cylinder in a plane Poiseuille flow in a microchannel is analyzed for a wide range of cylinder radii and positions across the channel. The cylinder translates parallel to the channel walls and rotates about its axis. The Stokes approximation is used and the problem is solved analytically using the Papkovich-Fadle eigenfunction expansion and the least-squares method. The stream function and the pressure distribution of the flow field are obtained as results. The force and moment exerted on the cylinder, and the pressure change far from the cylinder, are calculated and shown as functions of the size and location of the cylinder. The results confirm some reciprocal relations exactly. In particular, the translational and rotational velocities of the drifting cylinder in the existing Poiseuille flow are determined. The induced pressure change, when the cylinder drifts in the Poiseuille flow, is also calculated. Some typical streamline patterns, depending on the size and location of the cylinder, are shown and discussed. When the cylinder translates and/or rotates in the channel blocked at infinity, a series of Moffatt eddies appears far from the cylinder in the channel, as expected.

Beschreibung: Interactions between capillary and elastic effects are relevant to a variety of applications from micro- and nano-scale manufacturing to biological systems. In this work, we investigate capillary flows in flexible, millimeter-scale cylindrical elastic tubes. We demonstrate that surface tension can cause sufficiently flexible tubes to collapse and coalesce spontaneously through non-axisymmetric buckling, and develop criteria for the initial deformation and complete collapse of a circular tube. Experimental results for capillary rise and evaporation of a liquid in a flexible tube are presented, and several regimes are seen for the equilibrium state of a flexible tube deforming under capillary pressure. Deformations of the tube walls are measured in different regimes and compared with a shell theory model. Analysis and experimental results show that despite the complex and non-axisymmetric deformed shapes of cylindrical structures, the elastocapillary length used in previous literature for flat plates and sheets can be used to predict the behavior of flexible tubes.

Beschreibung: We present laboratory experimental results demonstrating that librational forcing of an ellipsoidal container of water can produce intense motions through the mechanism of a libration driven elliptical instability (LDEI). These libration studies are conducted using an ellipsoidal acrylic container filled with water. A particle image velocimetry method is used to measure the 2D velocity field in the equatorial plane over hundreds libration cycles for a fixed Ekman number, E = 2 × 10 −5 . In doing so, we recover the libration induced base flow and a time averaged zonal flow. Further, we show that LDEI in non-axisymmetric container geometries is capable of driving both intermittent and saturated turbulent motions in the bulk fluid. Additionally, we measure the growth rate and amplitude of the LDEI induced excited flow in a fully ellipsoidal container at more extreme parameters than previously studied [Noir et al. , “Experimental study of libration-driven flows in nonaxisymmetric containers,” Phys. Earth Planet. Inter. 204-205 , 1 (2012); Cébron et al. , Phys. Fluids 24 , 061703, “Libration driven elliptical instability,” (2012)]. Excitation of bulk filling turbulence by librational forcing provides a mechanism for transferring rotational energy into turbulent fluid motion and thus can play an important role in the thermal evolution, interior dynamics, and magneto-hydrodynamics of librating bodies, as appear to be common in solar system settings [e.g., Comstock and Bills, “A solar system survey of forced librations in longitude,” J. Geophys. Res. Planets 108 , 1 (2003)].

Beschreibung: We numerically study the displacement flow of two iso-viscous Newtonian fluids in an inclined two-dimensional channel, formed by two parallel plates. The results are complementary to our previous studies on displacement flows in pipes and channels. The heavier displacing fluid moves the lighter displaced fluid in the downward direction. Three dimensionless groups largely describe these flows: the densimetric Froude number ( Fr ), the Reynolds number ( Re ), and the duct inclination (β). As a first order approximation, we are able to classify different flow regimes phenomenologically in a two-dimensional ( Fr ; Re cosβ/ Fr )-plane and provide leading order expressions for the transitions between different regimes. The stabilizing and/or de-stabilizing effects of the imposed mean flow on buoyant exchange flows (zero imposed velocity) are described for a broad range of dimensionless parameters.

Beschreibung: Compressible granular materials are involved in many applications, some of them being related to energetic porous media. Gas permeation effects are important during their compaction stage, as well as their eventual chemical decomposition. Also, many situations involve porous media separated from pure fluids through two-phase interfaces. It is thus important to develop theoretical and numerical formulations to deal with granular materials in the presence of both two-phase interfaces and gas permeation effects. Similar topic was addressed for fluid mixtures and interfaces with the Discrete Equations Method (DEM) [R. Abgrall and R. Saurel, “Discrete equations for physical and numerical compressible multiphase mixtures,” J. Comput. Phys. 186 (2), 361-396 (2003)] but it seemed impossible to extend this approach to granular media as intergranular stress [K. K. Kuo, V. Yang, and B. B. Moore, “Intragranular stress, particle-wall friction and speed of sound in granular propellant beds,” J. Ballist. 4 (1), 697-730 (1980)] and associated configuration energy [J. B. Bdzil, R. Menikoff, S. F. Son, A. K. Kapila, and D. S. Stewart, “Two-phase modeling of deflagration-to-detonation transition in granular materials: A critical examination of modeling issues,” Phys. Fluids 11 , 378 (1999)] were present with significant effects. An approach to deal with fluid-porous media interfaces was derived in Saurel et al. [“Modelling dynamic and irreversible powder compaction,” J. Fluid Mech. 664 , 348-396 (2010)] but its validity was restricted to weak velocity disequilibrium only. Thanks to a deeper analysis, the DEM is successfully extended to granular media modelling in the present paper. It results in an enhanced version of the Baer and Nunziato [“A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials,” Int. J. Multiphase Flow 12 (6), 861-889 (1986)] model as symmetry of the formulation is now preserved. Several computational examples are shown to validate and illustrate method’s capabilities.

Beschreibung: We perform a theoretical and numerical study of the Coulomb-driven electroconvection flow of a dielectric liquid between two coaxial cylinders. The specific case, where the inner to outer diameter ratio is 0.5, is analyzed. A strong unipolar injection of ions either from the inner or outer cylinder is considered to introduce free charge carriers into the system. A finite volume method is used to solve all governing equations including Navier-Stokes equations and a simplified set of Maxwell’s equations. The flow is characterized by a subcritical bifurcation in the finite amplitude regime. A linear stability criterion and a nonlinear one that correspond to the onset and stop of the flow motion, respectively, are linked with a hysteresis loop. In addition, we also explore the behavior of the system for higher values of the stability parameter. For inner injection, we observe a transition between the patterns made of 7 and 8 cells, before an oscillatory regime is attained. Such a transition leads to a second finite amplitude stability criterion. A simple modal analysis reveals that the competition of different modes is at the origin of this behavior. The charge density, as well as velocity field distributions is provided to help understand the bifurcation behavior.

Beschreibung: Phytoplankton patchiness, namely the heterogeneous distribution of microalgae over multiple spatial scales, dramatically impacts marine ecology. A spectacular example of such heterogeneity occurs in thin phytoplankton layers (TPLs), where large numbers of photosynthetic microorganisms are found within a small depth interval. Some species of motile phytoplankton can form TPLs by gyrotactic trapping due to the interplay of their particular swimming style (directed motion biased against gravity) and the transport by a flow with shear along the direction of gravity. Here we consider gyrotactic swimmers in numerical simulations of the Kolmogorov shear flow, both in laminar and turbulent regimes. In the laminar case, we show that the swimmer motion is integrable and the formation of TPLs can be fully characterized by means of dynamical systems tools. We then study the effects of rotational Brownian motion or turbulent fluctuations (appearing when the Reynolds number is large enough) on TPLs. In both cases, we show that TPLs become transient, and we characterize their persistence.

Beschreibung: The model of gas bubble growth in high-viscous gas-saturated magmatic melt, subjected to rapid decompression, is presented in the current study. It is shown that consideration of unsteady character of the process is extremely important in a wide range of supersaturation. The analytical solution is found for the profile of dissolved gas concentration and the rate of bubble growth. The model of kinetics of overall degassing is developed. This model is based on distinguishing the so-called “forbidden” zone in the melt volume with suppressed formation of the new nucleation sites. The simple analytical dependences of the number of nucleating bubbles and typical nucleation time on the value of initial decompression were derived together with time dependence of volumetric concentration of the gas phase. Our results match the available experimental data.

Beschreibung: We investigate the effect of viscosity contrast on the stability of gravitationally unstable, diffusive layers in porous media. Our analysis helps evaluate experimental observations of various diffusive (boundary) layer models that are commonly used to study the sequestration of CO 2 in brine aquifers. We evaluate the effect of viscosity contrast for two basic models that are characterized with respect to whether or not the interface between CO 2 and brine is allowed to move. We find that diffusive layers are in general more unstable when viscosity decreases with depth within the layer compared to when viscosity increases with depth. This behavior is in contrast to the one associated with the classical displacement problem of gravitationally unstable diffusive layers that are subject to mean flow. For the classical problem, a greater instability is associated with the displacement of a more viscous, lighter fluid along the direction of gravity by a less viscous, heavier fluid. We show that the contrasting behavior highlighted in this study is a special case of the classical displacement problem that depends on the relative strength of the displacement and buoyancy velocities. We demonstrate the existence of a critical viscosity ratio that determines whether the flow is buoyancy dominated or displacement dominated. We explain the new behaviors in terms of the interaction of vorticity components related to gravitational and viscous effects.

Beschreibung: ABSTRACT With the objective to understand the generation, propagation and nonlinear evolution of ion cyclotron waves (ICWs) in the corona and solar wind, we use electromagnetic hybrid (kinetic ions, fluid electrons) simulations with a non-uniform magnetic field. ICWs are generated by the temperature anisotropy of O 5+ ions as minority species in a proton-electron plasma with uniform density. A number of magnetic field models are used including radial and spiral with field strength decreasing linearly or with the square of the radial distance. O 5+ ions with perpendicular temperature larger than parallel are initially placed in the high magnetic field regions. These ions are found to expand outward along the magnetic field. Associated with this expansion, ion cyclotron waves propagating along the magnetic field are also seen to expand outward. These waves are generated at frequencies below the local gyro-frequency of O 5+ ions propagating parallel and anti-parallel to the magnetic field. Through analysis of the simulation results we demonstrate that wave generation and absorption takes place at all radial distances. Comparing the simulation results to observations of ICWs in the solar wind shows some of the observed wave characteristics may be explained by the mechanism discussed in this paper.

Beschreibung: The characteristics of nighttime medium-scale travelling ionospheric disturbance (MSTID) features observed over Yonaguni (24.5 o N, 123.0 o E; 19.3 o N dip latitude), Japan are studied using all-sky imaging of OI 630.0 nm airglow emission. The uniqueness of these observations is that the area observed by the imager covers the transition region between low to middle latitudes in the ionosphere. Typical low latitude limit of mid-latitude type nighttime MSTIDs possessing phase front alignments along the northwest to the southeast occurs in this region. These MSTID features are rarely sighted at dip latitudes below 15 o . We selected two year period for analysis in which one year corresponded to the solar minimum conditions and another year to the solar maximum conditions. The MSTIDs were observed to extend to farther lower latitudes during the solar minimum conditions than during the solar maximum periods. Their observed range of wavelengths, phase velocities, phase front alignment and propagation directions are similar to those observed at typical mid-latitude sites. However, on many occasions the phase fronts of the observed MSTIDs did not extend over the whole field of view of the imager indicating that some process inhibits their extension to further lower latitudes. Detailed investigation suggests that the poleward propagating enhancement of airglow intensity, probably associated with the midnight pressure bulge, causes the MSTID features to disappear when they reach lower latitudes later in the night. When the MSTIDs reach lower latitudes well before midnight, they are found to be inhibited by the equatorial ionization anomaly crest region.

Beschreibung: We study a magnetosphere-ionosphere coupling at low-latitudes during a moderate (CIR/HSS-driven) geomagnetic storm on 22 July 2009. Recently, it has been shown that during major (CME-driven) storms, quasi-trapped 〉30 keV electrons largely enhance below the radiation belt in the forbidden zone and produce an additional ionization in the topside ionosphere. In this work, we examine a case of the recurrent storm when the magnetosphere-ionosphere coupling through the quasi-trapped electrons also may take place. Data from NOAA/POES and Japanese GOSAT satellites were used to identify the forbidden electron enhancement (FEE). We find a positive vertical gradient of the electron fluxes that indicates to the radiation belt as a source of FEE. Using global ionospheric maps (GIM), radiotomography reconstructions from beacon data and COSMIC/FS3 radio occultation measurements, we have observed an unusually large area in the night-time ionosphere with increased total electron content (TEC) and prominent elevation of the F-layer at low-latitudes that coincides with FEEs spatially and temporarily. Ionizing particles are considered as an addition source of ionization along with generally accepted mechanisms for storm time TEC increase (a positive ionospheric storm). We discuss relative contributions of the FEE and disturbance dynamo electric field in the TEC increases during the storm recovery phase.

Beschreibung: The development of a round liquid jet under the influence of a confined coaxial flow of an immiscible liquid of comparable density (central to annular flow density ratio of 8:10) was investigated in the vicinity of the nozzle exit. Two flow regimes were considered; one where the annular flow is faster than the central jet, so the central liquid jet is accelerated and one where the annular flow is slower, so the central liquid jet is decelerated. The central jet was visualised by high speed photography. Three modes of jet development were identified and classified in terms of the Reynolds number, Re, of the central jet which was in the range of 525 〈 Re 〈 2725, a modified definition of the Weber number, We, which allows the distinction between accelerating and deceleration flows and was in the range of −22 〈 We 〈 67 and the annular to central Momentum Ratio, MR, of the two streams which was in the range of 3.6 〈 MR 〈 91. By processing the time resolved jet images using Proper Orthogonal Decomposition (POD), it was possible to reduce the description of jet morphology to a small number of spatial modes, which isolated the most significant morphologies of the jet development. In this way, the temporal and spatial characteristics of the instabilities on the interface were clearly identified which highlights the advantages of POD over direct observation of the images. Relationships between the flow parameters and the interfacial waves were established. The wavelength of the interfacial instability was found to depend on the velocity of the fastest moving stream, which is contrary to findings for fluids with large density differences.

Beschreibung: Interaction of a vortex ring impinging on multiple permeable screens orthogonal to the ring axis was studied to experimentally investigate the persistence and decay of vortical structures inside the screen array using digital particle image velocimetry in a refractive index matched environment. The permeable screens had porosities (open area ratios) of 83.8%, 69.0%, and 55.7% and were held by a transparent frame that allowed the screen spacing to be changed. Vortex rings were generated using a piston-cylinder mechanism at nominal jet Reynolds numbers of 1000, 2000, and 3000 with piston stroke length-to-diameter ratios of 2 and 3. The interaction of vortex rings with the porous medium showed a strong dependence of the overall flow evolution on the screen porosity, with a central flow being preserved and vortex ring-like structures (with smaller diameter than the primary vortex ring) being generated near the centerline. Due to the large rod size used in the screens, immediate reformation of the transmitted vortex ring with size comparable to the primary ring (as has been observed with thin screens) was not observed in most cases. Since the screens have lower complexity and high open area ratios, centerline vortex ring-like flow structures formed with comparable size to the screen pore size and penetrated through the screens. In the case of low porosity screens (55.7%) with large screen spacing, re-emergence of large scale (large separation), weak vortical structures/pairs (analogous to a transmitted vortex ring) was observed downstream of the first screen. Additional smaller scale vortical structures were generated by the interaction of the vortex ring with subsequent screens. The size distribution of the generated vortical structures were shown to be strongly affected by porosity, with smaller vortical structures playing a stronger role as porosity decreased. Finally, porosity significantly affected the decay of total energy, but the effect of screen spacing decreased as porosity decreased.

Beschreibung: Bursty bulk flow (BBF) events, frequently observed in the magnetotail, carry significant energy and mass from the tail region at distances that are often greater than 20 R E into the near-Earth plasma sheet at ~10 R E where the flow is slowed and/or diverted. This region at ~10 R E is referred to as the BBF braking region. A number of possible channels are available for the transfer or dissipation of energy in BBF events including adiabatic heating of particles, the propagation of Alfvén waves out of the BBF braking region and into the auroral region, diverted flow out of the braking region, and energy dissipation within the braking region itself. This study investigates the generation of intense high-frequency electric field activity observed within the braking region. When present, these intense electric fields have power above the ion cyclotron frequency and almost always contain nonlinear structures such as electron phase space holes and double layers, which are often associated with field-aligned currents. A hypothesis in which the observed high-frequency electric field activity is generated by field-aligned currents resulting from turbulence in the BBF braking region is considered. Although linear Alfvén waves can generate field-aligned currents, based on theoretical calculations, the required currents are likely not the result of linear waves. Observations from the THEMIS satellites support the picture of a turbulent plasma leading to the generation of nonlinear kinetic structures. This work provides a possible mechanism for energy dissipation in turbulent plasmas.

Beschreibung: Dipolarization fronts (DFs) are often associated with the leading edge of Earthward bursty bulk flows in the magnetotail plasma sheet. Here multi-spacecraft THEMIS observations are used to show that a spatially limited region of counter-propagating ion beams, whose existence is not evident in either the plasma moments or the electric field, is observed on the low density side of DFs. The THEMIS magnetic field data are used to establish appropriate comparison cuts through a particle-in-cell (PIC) simulation of reconnection, and very good agreement is found between the observed and simulated ion distributions on both sides of the DF. Self-consistent back-tracing shows that the ion beams originate from the thermal component of the pre-existing high density plasma into which the DF is propagating; they do not originate from the inflow region in the traditional sense. Forward tracing shows that some of these ions can subsequently overtake the DF and pass back into the high density pre-existing plasma sheet with an order-of-magnitude increase in energy; this process is distinct from other ion reflection processes that occur directly at the DF. The interaction of the reconnection jet with the pre-existing plasma sheet therefore occurs over a macroscopic region, rather than simply being limited to the thin DF interface. A more general consequence of this study is the conclusion that reconnection jets are not simply fed by plasma inflow across the separatrices, but are also fed by plasma from the region into which the jet is propagating; the implications of this finding are discussed.

Beschreibung: Low energy (1-10 MeV) neutrons emanating from the Sun provide unique information about accelerated ions with steep energy spectra that may be produced in weak solar flares. However, observation of these solar neutrons can only be made in the inner heliosphere where measurement is difficult due to high background rates from neutrons produced by energetic ions interacting in the spacecraft. These ions can be from solar energetic particle events or produced in passing shocks associated with fast coronal mass ejections. Therefore, it is of the utmost importance that investigators rule out these secondary neutrons before making claims about detecting neutrons from the Sun. The MErcury Surface, Space ENvironment, GEochemistry, and Ranging ( tect MESSENGER ) Neutron Spectrometer recorded an hour-long neutron transient beginning at 15:45 UTC on 2011 June 4 for which [13] claim there is “strong evidence" that the neutrons were produced by the interaction of ions in the solar atmosphere. We studied this event in detail using data from the MESSENGER neutron spectrometer, gamma-ray spectrometer, X-ray Spectrometer, and Energetic Particle Spectrometer, and from the particle spectrometers on STEREO A . We demonstrate that the transient neutrons were secondaries produced by energetic ions, probably accelerated by a passing shock, that interacted in the spacecraft. We also identify significant faults with the authors’ arguments in favor of a solar neutron origin for the transient.

Beschreibung: The response of the D-region low latitude ionosphere has been examined for extreme space weather event of 14-16 December 2006 associated with a X1.5 solar flare and an intense geomagnetic storm ( Dst  = -146 nT) using VLF signals from NWC (19.8 kHz) and NPM (21.4 kHz) transmitters monitored at Suva (Geog. 18.10 ο S, 178.40 ο E), Fiji. Modeling of flare associated amplitude and phase enhancements of NWC (3.6 dB, 223 o ) and NPM (5 dB, 153 o ) using Long Wave Propagation Capability code shows reduction in the D-region reflection height ( H ') by 11.1 km and 9.4 km, and enhancement in ionization gradients described by increases in the exponential sharpness factor ( β ) by 0.122 and 0.126 km -1 , for the NWC and NPM paths, respectively. During the storm the daytime signal strengths of the NWC and NPM signals were reduced by 3.2 dB on 15 and 16 December (for about 46 hrs) and recovered by 17 December. Modelling for the NWC path shows that storm-time values of H ' and β were reduced by 1.2 km and 0.06 km -1 , respectively. Morlet wavelet analysis of signals amplitudes shows no clearly strong signatures of gravity wave propagation to low latitudes during the main and recovery phases. The reduction in VLF signal strength is due to increased signal attenuation and absorption by the Earth-ionosphere waveguide due to storm-induced D-region ionization changes and hence changes in D-region parameters. The long duration of the storm effect results from the slow diffusion of changed composition/ionization at D-region altitudes compared with higher altitudes in the ionosphere.

Beschreibung: We study the ion density and temperature in the pre- and post-dipolarization plasma sheets in the Earth's magnetotail, using 9 years (2001–2009) of Cluster data. For our study we selected cases when Cluster observed dipolarization fronts (DFs) with an earthward plasma flow greater than 150 km/s. We perform a statistical study of the temperature and density variations during the DF crossings. Earlier studies concluded that on average the temperature increases while the densitydecreases across the DF. Our statistical results show a more diverse picture: While ~53% of the DFs follow this pattern (category A), for ~28% the temperature decreases while the density increases across the DF (category B). We found an overall decrease in thermal pressure for category A DFs with a more pronounced decrease at the DF duskside, while DFs of the category B showed no clear pattern in the pressure change. Both categories are associated with earthward plasma flows, but with some difference: (1) Category A flows are faster than category B flows. (2) The observations indicate that category B flows are directed perpendicular to the current in the near-Earth current sheet while category A flows are tilted slightly duskward from this direction. (3) The background B z of category B is higher than that of category A. Based on these results we hypothesize that after reconnection takes place, a BBF emerges with category A characteristics, and as it travels earthward it further evolves into category B characteristics, which is in a more dipolarized region with slower plasma flow (closer to the flow braking region).

Beschreibung: Vortex cavitation forming in the leading-edge vortices of a delta wing was examined to determine how the individual cavitation bubbles incepted, grew, interacted with the underlying vortical flow and produced acoustic tones. The non-cavitating vortical flow over the delta wing was chosen to be similar to those previously reported in the literature. It was found that vortex breakdown was unaffected by the presence of incipient and developed vortex cavitation bubbles in the vortex core. While some cavitation bubbles incepted, grew, and collapsed relatively quickly, others reached an equilibrium position wherein the bubble tip was stationary in the laboratory frame at a particular location along the vortex axis. For a given attack angle, the equilibrium location moved upstream with a reduction in free stream cavitation number. It is shown that the existence of these stationary vortex bubbles is possible when there is a balance between the axial growth of the bubble along the vortex axis and the opposite motion of the axial jetting flow in the vortex core, and only a single equilibrium position is possible along the axially evolving vortex for a given free stream cavitation number. These transient and stationary vortex bubbles emit significant cavitation noise upon inception, growth, and collapse. The spectral content of the noise produced was expected to be related to the interaction of the bubble with the surrounding vortical flow in a manner similar to that reported in previous studies, where sustained tones were similar to the underlying vortex frequency. However, in the present study, the dominant frequency and higher harmonics of the tones occur at a higher frequency than that of the underlying vortex. Hence, it is likely that the highly elongated stationary bubbles have higher-order volume oscillations compared to the two-dimensional radial mode of the vortex cores of vortex cavitation bubbles with much smaller diameter-to-length ratios.

Beschreibung: Laminar flow over a periodic array of cylindrical surface roughness elements is simulated with an immersed boundary spectral method both to validate the method for subsequent studies and to examine how persistent streamwise vortices are introduced by a low Reynolds number roughness element. Direct comparisons are made with prior studies at a roughness-based Reynolds number Re k (= U ( k ) k / ν ) of 205 and a diameter to spanwise spacing ratio d / λ of 1/3. Downstream velocity contours match present and past experiments very well. The shear layer developed over the top of the roughness element produces the downstream velocity deficit. Upstream of the roughness element, the vortex topology is found to be consistent with juncture flow experiments, creating three cores along the recirculation line. Streamtraces stemming from these upstream cores, however, have unexpectedly little effect on the downstream flowfield as lateral divergence of the boundary layer quickly dissipates their vorticity. Long physical relaxation time of the recirculating wake behind the roughness remains a prominent issue for simulating this type of flowfield.

Beschreibung: Ice cores are archives of climate change and possibly large solar proton events (SPEs). Wolff et al . [2012] used a single event, a nitrate peak in the GISP2-H core, which McCracken et al . [2001a] time associated with the poorly quantified 1859 Carrington event, to discredit SPE-produced, impulsive nitrate deposition in polar ice. This is not the ideal test case. We critique the Wolff et al. analysis and demonstrate that the data they used cannot detect impulsive nitrate events because of resolution limitations. We suggest re-examination of the top of the Greenland ice sheet at key intervals over the last two millennia with attention to fine resolution and replicate sampling of multiple species. This will allow further insight into polar depositional processes on a sub-seasonal scale, including atmospheric sources, transport mechanisms to the ice sheet, post-depositional interactions, and a potential SPE association.

Beschreibung: Solutal Marangoni instability (SMI) is investigated both in 2D and 3D using a combined Cahn-Hilliard and Navier-Stokes model in a finite system. Fe-Sn is chosen as a representative alloy system since the phase diagram reveals a region with a miscibility gap, where two liquid phases, namely, the Fe-rich phase L 1 and the Sn-rich phase L 2 , are in chemical equilibrium. In 3D, considering a perturbed liquid cylinder ( L 2 phase) with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H × H (length × width × height) surrounded by the phase L 1 , we find that the perturbation induced Marangoni flow is either clockwise or anti-clockwise depending on the mean curvature difference between the convex and concave regions which is affected by the ratio of λ/ R 0 . The critical ratio of λ/ R 0 for SMI is shown to be invariant for different Marangoni numbers as well as independent of the geometrical properties of the L 1 phase. In 2D, a perturbed liquid pipe with a length of λ and a radius of R 0 embedded in the middle of a simulation box of λ × H (length × height) is taken into account. Due to different curvature constitution, the critical ratio of λ/ R 0 for SMI depends on the height of the L 1 phase.

Beschreibung: For several years, a promising Plasma Synthetic Jet actuator for high-speed flow control has been under development at ONERA. So far, its confined geometry and small space-time scales at play have prevented its full experimental characterization. Complementary accurate numerical simulations are then considered in this study in order to provide a complete aerothermodynamic description of the actuator. Two major obstacles have to be overcome with this approach: the modeling of the energy deposited by the electric arc and the accurate computation of the transient response of the cavity generating the pulsed jet. To solve the first problem, an Euler solver coupled with an electric circuit model was used to evaluate the energy deposition in the cavity. Such a coupling is performed by considering the electric field between the two electrodes. The second issue was then addressed by injecting these source terms in large Eddy simulations of the entire actuator. Aerodynamic results were finally compared with Schlieren visualizations. Using the proposed methodology, the temporal evolution of the jet front is remarkably well predicted.

Beschreibung: Severe geomagnetic storms have a strong impact on space communication and satellite navigation systems. Forecasting the appearance of geomagnetically induced disturbances in the ionosphere is one of the urgent goals of the space weather community. The challenge is that the processes governing the distribution of the crucial ionospheric parameters has a rather poor quantitative description and the models, built using the empirical parameterisations, have limited capabilities for operational purposes. On the other hand, data assimilation techniques are becoming more and more popular for nowcasting the state of the large-scale geophysical systems. We present an example of an ionospheric data assimilation system performance assessment during a strong geomagnetic event, which took place on 26 September 2011. The first-principle model has assimilated slant total electron content measurements from a dense network of ground stations, provided by the Norwegian Mapping Authority. The results have shown satisfactory agreement with independent data and demonstrate that the assimilation model is accurate to about 2–4 TEC units, and can be used for operational purposes in high-latitude regions. The operational system performance assessment is the subject of future work.

Beschreibung: Snapshot and classical proper orthogonal decomposition (POD) are used to examine the large-scale, energetic motions in fully developed turbulent pipe flow at Re D = 47,000 and 93,000. The snapshot POD modes come in pairs, representing the same azimuthal mode number but with a simple phase shift. The first 10 snapshot POD modes, associated with the very large scale motions (VLSMs), contribute 43% of the average Reynolds shear stress, and for first 80 modes u ′ and v ′ are anti-correlated so that they all contribute to positive shear stress events. The attached motions are contained in the lower order modes, and detached motions do not appear until snapshot POD mode numbers ≥15. We find that snapshot POD can introduce mode mixing, which is avoided in classical POD. Classical POD also gives frequency information, confirming that the low order modes capture well the behavior of the very large scale motions.

Beschreibung: [1]  The last solar minimum period was anomalously extended and low in EUV irradiance compared with previous solar minima. It can readily be expected that the thermosphere and ionosphere must be correspondingly affected by this low solar activity. While there have been unanimous reports on the thermospheric changes, being cooler and lower in its density as expected, the ionospheric responses to low solar activity in previous studies were not consistent with each other, probably due to the limited ionospheric observations used for them. In this study, we utilized the measurements of total electron content (TEC) from TOPEX and JASON-1 satellites during the periods of 1992 to 2010, which includes both the last two solar minimum periods, in order to investigate how the ionosphere responded to the extremely low solar activity during the last solar minimum compared with previous solar minimum. Although the global daily mean TECs show negligible differences between the two solar minimum periods, the global TEC maps reveal that there are significant systematic differences ranging from about -30% to +50% depending on local time, latitude and season. The systematic variations of the ionospheric responses seem to mainly result from the relative effects of reduced solar EUV production and reduced recombination rate due to thermospheric changes during the last solar minimum period.

Beschreibung: [1]  Some of the potentially most destructive effects of severe space weather storms are caused by the Geomagnetically Induced Currents. GICs can cause failures of electric transformers and result in wide-spread blackouts. GICs are induced by the time variability of the magnetic field, and are closely related to the time derivative of the local magnetic field perturbation. Predicting dB/dt is rather challenging, since the local magnetic perturbations and their time derivatives are both highly fluctuating quantities, especially during geomagnetic storms. The currently available first-principles based and empirical models cannot predict the detailed minute-scale or even faster time variation of the local magnetic field. On the other hand, Pulkkinen et al. [2013] demonstrated recently that several models can predict with positive skill scores whether the horizontal component of dB/dt at a given magnetometer station will exceed some threshold value in a 20-minute time interval. In this paper we investigate if one can improve the efficiency of the prediction further. We find that the Space Weather Modeling Framework, the best performing among the five models compared by Pulkkinen et al. [2013], shows significantly better skill scores in predicting the magnetic perturbation than predicting its time derivative, especially for large deviations. We also find that there is a strong correlation between the magnitude of dB/dt and the magnitude of the horizontal magnetic perturbation itself. Combining these two results one can devise an algorithm that gives better skill scores for predicting dB/dt exceeding various thresholds in 20-minute time intervals than the direct approach.

Beschreibung: We explore the instabilities developed in a fluid in which viscosity depends on temperature. In particular, we consider a dependency that models a very viscous (and thus rather rigid) lithosphere over a convecting mantle. To this end, we study a 2D convection problem in which viscosity depends on temperature by abruptly changing its value by a factor of 400 within a narrow temperature gap. We conduct a study which combines bifurcation analysis and time-dependent simulations. Solutions such as limit cycles are found that are fundamentally related to the presence of symmetry. Spontaneous plate-like behaviors that rapidly evolve towards a stagnant lid regime emerge sporadically through abrupt bursts during these cycles. The plate-like evolution alternates motions towards either the right or the left, thereby introducing temporary asymmetries on the convecting styles. Further time-dependent regimes with stagnant and plate-like lids are found and described.

Beschreibung: We study the shock wave structure in a rarefied polyatomic gas based on a simplified model of extended thermodynamics in which the dissipation is due only to the dynamic pressure. In this case the differential system is very simple because it is a variant of Euler system with a new scalar equation for the dynamic pressure [T. Arima, S. Taniguchi, T. Ruggeri, and M. Sugiyama, Phys. Lett. A376, 2799–2803 (2012)]. It is shown that this theory is able to describe the three types of the shock wave structure observed in experiments: the nearly symmetric shock wave structure (Type A, small Mach number), the asymmetric structure (Type B, moderate Mach number), and the structure composed of thin and thick layers (Type C, large Mach number).

Beschreibung: In this paper, the scaling property of the inverse energy cascade and forward enstrophy cascade of the vorticity filed ω( x , y ) in two-dimensional (2D) turbulence is analyzed. This is accomplished by applying a Hilbert-based technique, namely Hilbert-Huang transform, to a vorticity field obtained from a 8192 2 grid-points direct numerical simulation of the 2D turbulence with a forcing scale k f = 100 and an Ekman friction. The measured joint probability density function p ( C , k ) of mode C i ( x ) of the vorticity ω and instantaneous wavenumber k ( x ) is separated by the forcing scale k f into two parts, which correspond to the inverse energy cascade and the forward enstrophy cascade. It is found that all conditional probability density function p ( C | k ) at given wavenumber k has an exponential tail. In the inverse energy cascade, the shape of p ( C | k ) does collapse with each other, indicating a nonintermittent cascade. The measured scaling exponent ζ ω I ( q ) is linear with the statistical order q , i.e., ζ ω I ( q ) = − q / 3 , confirming the nonintermittent cascade process. In the forward enstrophy cascade, the core part of p ( C | k ) is changing with wavenumber k , indicating an intermittent forward cascade. The measured scaling exponent ζ ω F ( q ) is nonlinear with q and can be described very well by a log-Poisson fitting: ζ ω F ( q ) = 1 3 q + 0.45 1 − 0 . 43 q . However, the extracted vorticity scaling exponents ζ ω ( q ) for both inverse energy cascade and forward enstrophy cascade are not consistent with Kraichnan's theory prediction. New theory for the vorticity field in 2D turbulence is required to interpret the observed scaling behavior.

Beschreibung: [1]  As a weakly magnetized planet, Mars ionosphere/atmosphere interacts directly with the shocked solar wind plasma flow. Even though many numerical studies have been successful in reproducing numerous features of the interaction process, these earlier studies focused mainly on interaction under steady solar wind conditions. Recent observations suggest that plasma escape fluxes are significantly enhanced in response to solar wind dynamic pressure pulses. In this study, we focus on the response of the ionosphere to pressure enhancements in the solar wind. Through modeling of two idealized events using a magnetohydrodynamics model, we find that the upper ionosphere of Mars responds almost instantaneously to solar wind pressure enhancements, while the collision dominated lower ionosphere (below ~150 km) does not have noticeable changes in density. We also find that ionospheric perturbations in density, magnetic field and velocity can last more than an hour after the solar wind returns to the quiet conditions. The topside ionosphere forms complicated transient shapes in response, which may explain unexpected ionospheric behaviors in recent observations. We also find that ionospheric escape fluxes do not correlate directly with simultaneous solar wind dynamic pressure. Rather, their intensities also depend on the earlier solar wind conditions. It takes a few hours for the ionospheric/atmospheric system to reach a new quasi-equilibrium state.

Beschreibung: [1]  Combining THEMIS wave and particle observations and a quantitative calculation of linear wave growth rate, we demonstrate that magnetosonic (MS) waves can be locally excited by ion ring distributions in the Earth's magnetosphere when the ion ring energy is comparable to the local Alfven energy. MS waves in association with ion ring distributions were observed by THEMIS A on 24 November 2010 in the afternoon sector, both outside the plasmapause where the wave spectrum varied with f LHR and inside the plasmapause where the wave frequency band remained nearly constant. Our plasma instability analysis in three different regions shows that higher and narrow frequency band MS waves are excited locally outside the plasmapause, and lower and broad frequency band MS waves are excited in the region where the density slightly increases. However, there is no evidence for wave excitation inside the plasmapause, and wave propagation from a distant source is needed to explain their existence. The simulation of the MS wave growth rate spectra during this event agrees reasonably well with the observed wave magnetic field power spectra. We also simulated a MS wave event on 19 October 2011 in the dusk sector, and found that the ion ring distribution with an ion ring energy slightly higher than the local Alfven energy can excite the typical broad band MS waves outside the plasmapause.

Beschreibung: [1]  The second Radio Aurora Explorer (RAX-2) satellite has completed more than 30 conjunction experiments with the AMISR chain of incoherent scatter radars in Alaska, and Resolute Bay, Canada. Coherent radar echoing occurred during four of the passes: three when E region electron drifts exceeded the ion acoustic speed threshold and one during HF heating of the ionosphere by the HAARP heater. In this paper, we present the results for the first three passes associated with backscatter from natural irregularities. We analyze, in detail, the largest drift case because the plasma turbulence was the most intense and because the corresponding ground-to-space bi-static scattering geometry was the most favorable for magnetic aspect sensitivity analysis. A set of data analysis procedures including interference removal, autocorrelation analysis, and the application of a radar beam deconvolution algorithm mapped the distribution of E region backscatter with 3 km resolution in altitude and ∼ 0.1 ∘ in magnetic aspect angle. To our knowledge, these are the highest resolution altitude-resolved magnetic aspect sensitivity measurements made at UHF frequencies in the auroral region. In this paper, we show that, despite the large electron drift speed of ∼ 1500 m/s, the magnetic aspect sensitivity of sub-meter scale irregularities is much higher than previously reported. The root-mean-square of the aspect angle distribution varied monotonically between 0.5-0.1 ∘ for the altitude range 100-110 km. Findings from this single but compelling event suggest that sub-meter scale waves propagating at larger angles from the main E  ×  B flow direction (secondary waves) have parallel electric fields that are too small to contribute to E region electron heating. It is possible anomalous electron heating in the auroral electrojet can be explained by (a) the dynamics of those sub-meter scale waves propagating in the E  ×  B direction (primary waves) or (b) the dynamics of longer wavelengths.

Beschreibung: We report experimental observations of the controlled deformation of a dielectric liquid jet subjected to a local high-voltage electrostatic field in the direction normal to the jet. The jet deforms to the shape of an elliptic cylinder upon application of a normal electrostatic field. As the applied electric field strength is increased, the elliptic cylindrical jet deforms permanently into a flat sheet, and eventually breaks-up into droplets. We interpret this observation—the stretch of the jet is in the normal direction to the applied electric field—qualitatively using the Taylor-Melcher leaky dielectric theory, and develop a simple scaling model that predicts the critical electric field strength for the jet-to-sheet transition. Our model shows a good agreement with experimental results, and has a form that is consistent with the classical drop deformation criterion in the Taylor-Melcher theory. Finally, we statistically analyze the resultant droplets from sheet breakup, and find that increasing the applied electric field strength improves droplet uniformity and reduces droplet size.

Beschreibung: Normal impingement of a single droplet on a thin liquid film is investigated numerically solving the axisymmetric Navier-Stokes equations. Gravity and viscosity are taken into account whereas compressibility effects are neglected. Two phases are tracked by means of volume of fluid method and adaptive mesh refinement is used to increase accuracy of the interface. Numerical results are validated both qualitatively and quantitatively using experimental measurements. Effects of gas density, gas viscosity, and film thickness on the crown behavior are studied. Influence of droplet deviation from spherical shape on the crown behavior is investigated. It is shown that increasing the gas density leads to reduction of crown radius evolution rate, while gas viscosity does not affect the rate of crown radius evolution. Development rate of crown height decreases by increasing the gas density. Reynolds number and splashing regime can change the effect of gas viscosity on the crown height evolution. Deviation of droplet from sphere can change behavior of crown completely as result of change in droplet mass center position. Difference between numerical results and experimental ones is justified using different droplet shapes.

Beschreibung: We present an analytical study of peak mode isotachophoresis (ITP), and provide closed form solutions for sample distribution and electric field, as well as for leading-, trailing-, and counter-ion concentration profiles. Importantly, the solution we present is valid not only for the case of fully ionized species, but also for systems of weak electrolytes which better represent real buffer systems and for multivalent analytes such as proteins and DNA. The model reveals two major scales which govern the electric field and buffer distributions, and an additional length scale governing analyte distribution. Using well-controlled experiments, and numerical simulations, we verify and validate the model and highlight its key merits as well as its limitations. We demonstrate the use of the model for determining the peak concentration of focused sample based on known buffer and analyte properties, and show it differs significantly from commonly used approximations based on the interface width alone. We further apply our model for studying reactions between multiple species having different effective mobilities yet co-focused at a single ITP interface. We find a closed form expression for an effective-on rate which depends on reactants distributions, and derive the conditions for optimizing such reactions. Interestingly, the model reveals that maximum reaction rate is not necessarily obtained when the concentration profiles of the reacting species perfectly overlap. In addition to the exact solutions, we derive throughout several closed form engineering approximations which are based on elementary functions and are simple to implement, yet maintain the interplay between the important scales. Both the exact and approximate solutions provide insight into sample focusing and can be used to design and optimize ITP-based assays.

Beschreibung: A novel waveform modified from the standard-sinusoidal function is adopted to enhance the virtual aeroshaping effect of the synthetic jets positioned at the front stagnation point of a circular cylinder. The waveform is characterized by a control parameter, namely, the suction duty cycle factor k , which is the ratio of the time duration of the suction cycle to that of the blowing cycle. The strength of the synthetic jet vortex pair is enhanced by increasing the suction duty cycle factor. The periodic closed envelope forms upstream of the circular cylinder for k ≤ 1.00, while the quasi-steady open envelope forms for k ≥ 2.00, acting the virtual aeroshaping effect. As a result, both the statistical characteristics and the vortex dynamics of the near-wake flow field change with the suction duty cycle factor. The recirculation region downstream of the circular cylinder becomes smaller or even disappears, and thus, the drag coefficient over the circular cylinder is reduced by increasing the suction duty cycle factor to k ≥ 1.00. The statistical mean and fluctuating velocities show corresponding changes in the near wake with the different wake patterns. For k ≤ 0.50, the wake vortex shows the antisymmetric shedding mode which is similar with the natural case. For 1.00 ≤ k ≤ 2.00, the wake vortex shows the bistable state mode, where vortex sheds with symmetric or antisymmetric mode; the antisymmetric shedding mode dominates the global flow field for k = 1.00, while it is the symmetric shedding mode that dominates the flow field for k = 2.00. For k = 4.00, it shows the antisymmetric shedding mode with a shorter vortex formation length than the natural case. The above findings indicate that the virtual aeroshaping effect of the synthetic jets can be enhanced by increasing the suction duty cycle factor so as to increase the momentum coefficient while keeping other control parameters unchanged, providing us another way for effective flow control.

Beschreibung: The dynamics of vapor-liquid interface are important because interfacial instability determines bubble growth, detachment frequency, waiting time, shape of bubbles, and the interrelationship between bubble formation sites. In this study, a detailed numerical simulation has been performed to understand the transition in bubble release pattern and multimode bubble formation in saturated pool boiling. The interfaces drop down alternatively at the nodes and antinodes of the wavelengths dictated by Rayleigh-Taylor instability and Taylor-Helmholtz instability. Due to higher degrees of superheat, vapor jets emanate from nodes and antinodes. An attempt has been made to predict the maximum and minimum heat fluxes during saturated pool boiling.

Beschreibung: We present a combination of experiment, theory, and modelling on laminar mixing at large Péclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, “The diffusive strip method for scalar mixing in two-dimensions,” J. Fluid Mech.662, 134–172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.

Beschreibung: Janus droplets are compound droplets that consist of two adhering drops of different fluids that are suspended in a third fluid. We use the Shan-Chen lattice Boltzmann method for multicomponent mixtures to simulate Janus droplets at rest and in shear. In this simulation model, interfacial tensions are not known a priori from the model parameters and must be determined using numerical experiments. We show that interfacial tensions obtained with the Young-Laplace law are consistent with those measured from the equilibrium geometry. The regimes of adhering, separated, and engulfing droplets were explored. Two different adhesion geometries were considered for two-dimensional simulations of Janus droplets in shear. The first geometry resembles two adhering circles with small overlap. In the second geometry, the two halves are semicircular. For both geometries, the rotation rate of the droplet depends on its orientation. The width of the periodic simulation domain also affects the rotation rate of both droplet types up to an aspect ratio of 6:1 (width:height). While the droplets with the first geometry oscillated about the middle of the domain, the droplets of the second geometry did not translate while rotating. A four-pole vortex structure inside droplets of the second geometry was found. These simulations of single Janus droplets reveal complex behaviour that implies a rich range of possibilities for the rheology of Janus emulsions.

Beschreibung: Understanding the physics of water evaporation from saline porous media is important in many natural and engineering applications such as durability of building materials and preservation of monuments, water quality, and mineral-fluid interactions. We applied synchrotron x-ray micro-tomography to investigate the pore-scale dynamics of dissolved salt distribution in a three dimensional drying saline porous media using a cylindrical plastic column (15 mm in height and 8 mm in diameter) packed with sand particles saturated with CaI 2 solution (5% concentration by mass) with a spatial and temporal resolution of 12 μ m and 30 min, respectively. Every time the drying sand column was set to be imaged, two different images were recorded using distinct synchrotron x-rays energies immediately above and below the K-edge value of Iodine. Taking the difference between pixel gray values enabled us to delineate the spatial and temporal distribution of CaI 2 concentration at pore scale. Results indicate that during early stages of evaporation, air preferentially invades large pores at the surface while finer pores remain saturated and connected to the wet zone at bottom via capillary-induced liquid flow acting as evaporating spots. Consequently, the salt concentration increases preferentially in finer pores where evaporation occurs. Higher salt concentration was observed close to the evaporating surface indicating a convection-driven process. The obtained salt profiles were used to evaluate the numerical solution of the convection-diffusion equation (CDE). Results show that the macro-scale CDE could capture the overall trend of the measured salt profiles but fail to produce the exact slope of the profiles. Our results shed new insight on the physics of salt transport and its complex dynamics in drying porous media and establish synchrotron x-ray tomography as an effective tool to investigate the dynamics of salt transport in porous media at high spatial and temporal resolution.

Beschreibung: [1]  Working toward a physical understanding of how solar-wind/magnetosphere coupling works, four arguments are presented indicating that the solar-wind electric field v sw  × B sw does not control the rate of reconnection between the solar wind and the magnetosphere. Those four arguments are (1) that the derived rate of dayside reconnection is not equal to solar-wind electric field, (2) that electric-field driver functions can be improved by a simple modification that disallows their interpretation as the solar-wind electric field, (3) that the electric field in the magnetosheath is not equal to the electric field in the solar wind, and (4) that the magnetosphere can mass load and reduce the dayside reconnection rate without regard for the solar-wind electric field. The data is more consistent with a coupling function based on local control of the reconnection rate than the Axford conjecture that reconnection is controlled by boundary conditions irrespective of local parameters. Physical arguments that the solar-wind electric field controls dayside reconnection are absent; it is speculated that it is a coincidence that the electric field does so well at correlations with geomagnetic indices.

Beschreibung: [1]  This paper presents a case study from a single, six-hour observing period to illustrate the application of techniques developed for interferometric radio telescopes to the spectral analysis of observations of ionospheric fluctuations with sparse arrays. We have adapted the deconvolution methods used for making high dynamic range images of cosmic sources with radio arrays to making comparably high dynamic range maps of spectral power of wavelike ionospheric phenomena. In the example presented here, we have used observations of the total electron content (TEC) gradient derived from Very Large Array (VLA) observations of synchrotron emission from two galaxy clusters at 330 MHz as well as GPS-based TEC measurements from a sparse array of 33 receivers located within New Mexico near the VLA. We show that these techniques provide a significant improvement in signal to noise (S/N) of detected wavelike structures by correcting for both measurement inaccuracies and wavefront distortions. This is especially true for the GPS data when combining all available satellite/receiver pairs, which probe a larger physical area and likely have a wider variety of measurement errors than in the single-satellite case. In this instance, we found the peak S/N of the detected waves was improved by more than an order of magnitude. The data products generated by the deconvolution procedure also allow for a reconstruction of the fluctuations as a two-dimensional waveform/phase screen that can be used to correct for their effects.

Beschreibung: [1]  In this paper, we test whether time periods with hot proton temperature anisotropy are associated with EMIC waves, and whether the plasma conditions during the observed waves satisfy the linear theory threshold condition. We identify 865 events observed by the Composition DIstribution Function (CODIF) instrument onboard Cluster spacecraft 4 (SC4) during 1 January 2001 – 1 January 2011 that exhibit a positive temperature anisotropy ( A hp  =  T ⊥  h / T ∥  h  − 1) in the 10-40 keV protons. The events occur over an L range from 4 to 10 in all magnetic local times and at magnetic latitudes (MLAT) within ±50°. Of these Hot Proton Temperature Anisotropy (HPTA) events, only 68 events have electromagnetic ion cyclotron (EMIC) waves. In these 68 HPTA events, for those at 3.8〈 L  ≤ 5 and |MLAT| ≤ 10 ° , the EMIC waves with powers 〉1.0 nT 2 /Hz mainly appear in the region with f EMIC / f H , eq  〈 0.8. Two stop bands are present, one near the region with f EMIC / f H , eq  ≈ 0.33, the other in the region with 0.8 〈  f EMIC / f H , eq  〈 0.9. Most of the EMIC waves in the He, H, and 〉 H bands satisfy A hp /( A hp  + 1) 〉  f EMIC / f H , lo , A hp /( A hp  + 1) 〉 0.45 *  f EMIC / f H , lo , and A hp /( A hp  + 1) 〈 0.45 *  f EMIC / f H , lo .   f EMIC , f H , eq   and f H , lo are the EMIC wave frequency, the magnetic equatorial and the local proton gyrofrequencies. We also find that the EMIC waves predominantly occur with A hp  〉 0.25. By testing a threshold equation for the EMIC instability based on linear theory, we find that for EMIC waves with |MLAT| ≤ 10 ° in the He, H and 〉 H bands the percentages that satisfy the predicted conditions for wave growth by the threshold equation are 15.2%, 24.6% and 25.6%. For the EMIC waves with |MLAT| 〉 10 ° the percentages that satisfy the wave growth predicted conditions are only 2.8%, 2.6% and 0.0%. Finally, possible reasons for the low forecast accuracies of EMIC waves are suggested.

Beschreibung: The behaviour of low Reynolds number, non-Boussinesq fountains from four different nozzle geometries (one circular and three rectangular nozzles) are studied. High speed laser schlieren imaging is used to study the fountain behaviour (frequency and penetration height). Bi-orthogonal decomposition and dynamic mode decomposition (DMD) are used to understand the unsteady characteristics of fountains. The flow regimes of fountains are classified as steady, flapping, and flapping-bobbing type. The DMD technique successfully separates the bobbing oscillation from the combined flapping-bobbing oscillation of the fountain. The frequency of the flapping oscillation, and the frequency of the bobbing oscillation in the flapping-bobbing regime scales as St h Fr h = C 1 and S t h F r h 2 = C 2 , respectively, where the characteristic length scale is the smallest dimension ( h ) of the nozzle. The mean steady state penetration heights ( Z s / h ) of “forced” low Reynolds number non-Boussinesq fountains are independent of nozzle shape (circular and rectangular), and scales linearly with the Froude number.

Beschreibung: The large-scale properties of self-similar unstably stratified homogeneous (USH) turbulence are investigated using an eddy-damped quasi-normal markovianized approximation of the nonlinear term. This analysis shows that a special role is played by the wave vectors contained in the equatorial plane, i.e., the plane perpendicular to gravity. It is indeed in this plane that turbulent spectra reach their maxima and evolve linearly from their initial condition when their initial infrared exponent is smaller than 4. At other angles, this property is not satisfied and turbulent spectra eventually undergo an evolution dominated by nonlinear backscattering processes. The self-similar evolution of USH turbulence is also shown to be related to the properties of large scales. In particular, the asymptotic growth rate of the mixing length depends on the initial infrared exponent in the equatorial plane. Besides, the self-similar asymptotic values of the concentration and velocity correlations also depend on the properties of large scales. This allows to derive relations between the correlations and the growth rate parameter.

Beschreibung: Secondary flow cells are commonly observed in straight laboratory channels, where they are often associated with duct corners. Here, we present velocity measurements acquired with an acoustic Doppler current profiler in a straight reach of the Seine river (France). We show that a remarkably regular series of stationary flow cells spans across the entire channel. They are arranged in pairs of counter-rotating vortices aligned with the primary flow. Their existence away from the river banks contradicts the usual interpretation of these secondary flow structures, which invokes the influence of boundaries. Based on these measurements, we use a depth-averaged model to evaluate the momentum transfer by these structures, and find that it is comparable with the classical turbulent transfer.

Beschreibung: [1]  The simultaneous onset of the preliminary impulse (PI) of the geomagnetic sudden commencement at high latitude and dayside dip equator is explained by means of the TM 0 mode waves propagating at the speed of light in the Earth-ionosphere waveguide (EIW) [Kikuchi et al., 1978]. A couple of issues remain to be addressed in the EIW model: (1) How is the TM 0 mode wave is excited by the field-aligned currents in the polar region? (2) How are the quasi-steady ionospheric currents are achieved by the TM 0 mode waves? (3) How simultaneous or delayed are the onset and peak of the equatorial PI with respect to the high latitude PI? To address these issues, we examine the TEM (TM 0 ) mode wave propagation in the finite-length transmission lines replacing the pair of FACs (magnetosphere-ionosphere (MI) transmission line) and the Earth-ionosphere waveguide (ionosphere-ground (IG) transmission line). The issue (1) is addressed by showing that a fraction of the TEM mode wave is transmitted from the MI to IG transmission lines through the polar ionosphere. To address the issues (2) and (3), we examine the properties of the finite-length IG transmission line with finite ionospheric conductivity. It is shown that the ionospheric currents start to grow instantaneously and continue to grow gradually with time constants of 1-10 sec depending on the ionospheric conductivity. The MIG transmission line enables us to explain the instantaneous onset and delayed peak time of the equatorial PI and quick electric field response of the low latitude ionosphere and inner magnetosphere.

Beschreibung: Analysis of fluxes across the turbulent/non-turbulent interface (TNTI) of turbulent boundary layers is performed using data from two-dimensional particle image velocimetry (PIV) obtained at high Reynolds numbers. The interface is identified with an iso-surface of kinetic energy, and the rate of change of total kinetic energy ( K ) inside a control volume with the TNTI as a bounding surface is investigated. Features of the growth of the turbulent region into the non-turbulent region by molecular diffusion of K , viscous nibbling, are examined in detail, focussing on correlations between interface orientation, viscous stress tensor elements, and local fluid velocity. At the level of the ensemble (Reynolds) averaged Navier-Stokes equations (RANS), the total kinetic energy K is shown to evolve predominantly due to the turbulent advective fluxes occurring through an average surface which differs considerably from the local, corrugated, sharp interface. The analysis is generalized to a hierarchy of length-scales by spatial filtering of the data as used commonly in Large-Eddy-Simulation (LES) analysis. For the same overall entrainment rate of total kinetic energy, the theoretical analysis shows that the sum of resolved viscous and subgrid-scale advective flux must be independent of scale. Within the experimental limitations of the PIV data, the results agree with these trends, namely that as the filter scale increases, the viscous resolved fluxes decrease while the subgrid-scale advective fluxes increase and tend towards the RANS values at large filter sizes. However, a definitive conclusion can only be made with fully resolved three-dimensional data, over and beyond the large dynamic spatial range presented here. The qualitative trends from the measurement results provide evidence that large-scale transport due to the energy-containing eddies determines the overall rate of entrainment, while viscous effects at the smallest scales provide the physical mechanism ultimately responsible for entrainment. Data spanning over a decade in Reynolds number suggest that the fluxes (or the entrainment velocity) scale with the friction velocity (or equivalently the local turbulent fluctuating velocity), whereas Taylor microscale and boundary-layer thickness are the appropriate length scales at small and large filter sizes, respectively.

Beschreibung: [1]  Plasmoids and other reconnection-related signatures have been observed in Jupiter's magnetotail through analysis of magnetic field and energetic particle data. Previous studies have established the spatial distribution and recurrence period of tail reconnection events, and identified the location of a statistical x-line separating inward and outward flow. Here we present new analysis focusing specifically on 43 plasmoid signatures observed in magnetometer data in order to establish the average properties and internal structure of Jovian plasmoids. We present statistics on the observed plasmoid length scale, duration, radial position, and local time distribution. On average, the observed plasmoids have a ~3 R J radial extent and ~7 minute duration, and result in the closure of ~4-8 GWb of open flux from reconnection of open field lines in the post-plasmoid plasma sheet. We also determine the amount of mass released and the magnetic flux closed in order to understand the role of tail reconnection in the transport of mass and flux in Jupiter's magnetosphere. The observed plasmoid properties are consistent with a mass loss rate of ~0.7-120 kg/s and a flux closure rate of ~7-70 GWb/day. We conclude that tail reconnection and plasmoid release is an important method of flux transport at Jupiter but likely cannot account for the mass input from Io, suggesting that additional mass loss mechanisms may be significant. Finally, we examine the plasmoid interior structure through minimum variance analysis and find that most plasmoids lack a core field and are better described by magnetic loops rather than flux ropes.

Beschreibung: [1]  Hot flow anomalies (HFAs) represent a subset of solar wind discontinuities interacting with collisionless bow shocks. They are typically formed when the normal component of the motional (convective) electric field points toward the embedded current sheet on at least one of its sides. The core region of an HFA contains hot and highly deflected ion flows and rather low and turbulent magnetic field. In this paper, we report observations of possible HFA-like events at Mercury identified over a course of two planetary years. Using data from the orbital phase of the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission, we identify a representative ensemble of active current sheets magnetically connected to Mercury's bow shock. We show that some of these events exhibit magnetic and particle signatures of HFAs similar to those observed at other planets, and present their key physical characteristics. Our analysis suggests that Mercury's bow shock does not only mediate the flow of supersonic solar wind plasma but also provides conditions for local particle acceleration and heating as predicted by previous numerical simulations. Together with earlier observations of HFA activity at Earth, Venus, Mars, and Saturn, our results suggest that hot flow anomalies could be a common property of planetary bow shocks, and show that the characteristic size of these events is controlled by the bow shock standoff distance and/or local solar wind conditions.

Beschreibung: In a recent direct numerical simulation (DNS) study [P. K. Yeung and K. R. Sreenivasan, “ Spectrum of passive scalars of high molecular diffusivity in turbulent mixing,” J. Fluid Mech.716, R14 (2013)] with Schmidt number as low as 1/2048, we verified the essential physical content of the theory of Batchelor, Howells, and Townsend [“Small-scale variation of convected quantities like temperature in turbulent fluid. 2. The case of large conductivity,” J. Fluid Mech.5, 134 (1959)] for turbulent passive scalar fields with very strong diffusivity, decaying in the absence of any production mechanism. In particular, we confirmed the existence of the −17/3 power of the scalar spectral density in the so-called inertial-diffusive range. In the present paper, we consider the DNS of the same problem, but in the presence of a uniform mean gradient, which leads to the production of scalar fluctuations at (primarily) the large scales. For the parameters of the simulations, the presence of the mean gradient alters the physics of mixing fundamentally at low Peclet numbers. While the spectrum still follows a −17/3 power law in the inertial-diffusive range, the pre-factor is non-universal and depends on the magnitude of the mean scalar gradient. Spectral transfer is greatly reduced in comparison with those for moderately and weakly diffusive scalars, leading to several distinctive features such as the absence of dissipative anomaly and a new balance of terms in the spectral transfer equation for the scalar variance, differing from the case of zero gradient. We use the DNS results to present an alternative explanation for the observed scaling behavior, and discuss a few spectral characteristics in detail.

Beschreibung: [1]  This paper investigates the influence of two solar eclipses on the ionosphere complexity measures: Tsallis entropy, Renyi entropy, Hurst exponent, beta exponent, fractal dimension. The study used GPS TEC measured at 3 locations in Japan during the solar eclipses of 22 July 2009 and 21 May 2012.This is the first effort to compare the complexity measures by comparing TEC time series of the eclipse day with those from the day before and day after the eclipse. It was found from analysis of the TEC observations that there were no abnormal variations of the complexity parameters from their expected values for either eclipse. Model calculations also show that TEC deviations during the eclipses are small.

Beschreibung: A non-equilibrium wall-model based on unsteady 3D Reynolds-averaged Navier-Stokes (RANS) equations has been implemented in an unstructured mesh environment. The method is similar to that of the wall-model for structured mesh described by Wang and Moin [Phys. Fluids14, 2043–2051 (2002)], but is supplemented by a new dynamic eddy viscosity/conductivity model that corrects the effect of the resolved Reynolds stress (resolved turbulent heat flux) on the skin friction (wall heat flux). This correction is crucial in predicting the correct level of the skin friction. Unlike earlier models, this eddy viscosity/conductivity model does not have a stress-matching procedure or a tunable free parameter, and it shows consistent performance over a wide range of Reynolds numbers. The wall-model is validated against canonical (attached) transitional and fully turbulent flows at moderate to very high Reynolds numbers: a turbulent channel flow at Re τ = 2000, an H-type transitional boundary layer up to Re θ = 3300, and a high Reynolds number boundary layer at Re θ = 31 000. Application to a separated flow over a NACA4412 airfoil operating close to maximum lift is also considered to test the performance of the wall-model in complex non-equilibrium flows.

Beschreibung: Conventional shallow water theory successfully reproduces many key features of the Jovian atmosphere: a mixture of coherent vortices and stable, large-scale, zonal jets whose amplitude decreases with distance from the equator. However, both freely decaying and forced-dissipative simulations of the shallow water equations in Jovian parameter regimes invariably yield retrograde equatorial jets, while Jupiter itself has a strong prograde equatorial jet. Simulations by Scott and Polvani [“Equatorial superrotation in shallow atmospheres,” Geophys. Res. Lett.35, L24202 (2008)] have produced prograde equatorial jets through the addition of a model for radiative relaxation in the shallow water height equation. However, their model does not conserve mass or momentum in the active layer, and produces mid-latitude jets much weaker than the equatorial jet. We present the thermal shallow water equations as an alternative model for Jovian atmospheres. These equations permit horizontal variations in the thermodynamic properties of the fluid within the active layer. We incorporate a radiative relaxation term in the separate temperature equation, leaving the mass and momentum conservation equations untouched. Simulations of this model in the Jovian regime yield a strong prograde equatorial jet, and larger amplitude mid-latitude jets than the Scott and Polvani model. For both models, the slope of the non-zonal energy spectra is consistent with the classic Kolmogorov scaling, and the slope of the zonal energy spectra is consistent with the much steeper spectrum observed for Jupiter. We also perform simulations of the thermal shallow water equations for Neptunian parameter values, with a radiative relaxation time scale calculated for the same 25 mbar pressure level we used for Jupiter. These Neptunian simulations reproduce the broad, retrograde equatorial jet and prograde mid-latitude jets seen in observations. The much longer radiative time scale for the colder planet Neptune explains the transition from a prograde to a retrograde equatorial jet, while the broader jets are due to the deformation radius being a larger fraction of the planetary radius.

Beschreibung: The hydrodynamic interaction of two deformable vesicles in shear flow induces a net displacement, in most cases an increase of their distance in the transverse direction. The statistical average of these interactions leads to shear-induced diffusion in the suspension, both at the level of individual particles which experience a random walk made of successive interactions, and at the level of suspension where a nonlinear down-gradient diffusion takes place, an important ingredient in the structuring of suspension flows. We make an experimental and computational study of the interaction of a pair of lipid vesicles in shear flow by varying physical parameters, and investigate the decay of the net lateral displacement with the distance between the streamlines on which the vesicles are initially located. This decay and its dependency upon vesicle properties can be accounted for by a simple model based on the well established law for the lateral drift of a vesicle in the vicinity of a wall. In the semi-dilute regime, a determination of self-diffusion coefficients is presented.

Beschreibung: [1]  We have analyzed the data of the world neutron monitor network for the first ground level enhancement of solar cycle 24, the GLE on May 17, 2012. A newly computed neutron monitor yield function and an inverse method are applied to estimate the energy spectrum, anisotropy axis direction and pitch-angle distribution of the high-energy solar particles in interplanetary space. The method includes the determination of the asymptotic viewing cones of neutron monitor stations through computations of trajectories of cosmic rays in a model magnetosphere. The cosmic ray particle trajectories are determined with the MAGNETOCOSMICS code using Tsyganenko 1989 and IGRF models. Subsequent calculation of the neutron monitor responses with the model function is carried out, that represents an initial guess of the inverse problem. Derivation of the solar energetic particle characteristics is fulfilled by fitting the data of the global neutron monitor network using the Levenberg-Marquardt method over the nine-dimensional parameter space. The pitch-angle distribution and rigidity spectrum of high-energy protons are obtained as function of time in the course of the GLE. The angular distribution appears quite complicated. It comprises a focused beam along the interplanetary magnetic field line from the Sun and a loss-cone feature around the opposite direction, possibly indicative of the particle transport in interplanetary magnetic field structures associated with previous coronal mass ejections.

Beschreibung: [1]  Geomagnetic activity is strongly controlled by solar wind and Interplanetary Magnetic Field (IMF) conditions, especially the southward component of IMF (IMF Bs). We analyze the statistical properties of IMF Bs at 1 AU using in situ observations for more than a solar cycle (1995 - 2010). IMF Bs-events are defined as continuous IMF Bs intervals with varying thresholds of Bs magnitude and duration, and categorized by different solar wind structures, such as magnetic cloud (MC), interplanetary small-scale magnetic flux rope (ISMFR), interplanetary coronal mass ejection (ICME) without MC signature (ejecta), stream interacting region (SIR), and shock, as well as events unrelated with well-defined solar wind structures. The statistical properties of IMF Bs-events and their geoeffectiveness are investigated in detail based on satellite and ground measurements. We find that the integrated duration and number of Bs-events follow the sunspot number when Bz 〈 -5 nT. We also find that in extreme Bs-events (t 〉 6 hours, Bz 〈 -10 nT), a majority (53 %) are related to MC and 10 % are related with ejecta, but nearly a quarter are not associated with any well-defined solar wind structure. We find different geomagnetic responsesfor Bs-events with comparable duration and magnitude depending on what type of solar wind structures they are associated with. We also find that great Bs-events (t 〉 3 hours, Bz 〈 -10 nT) do not always trigger magnetic storms.

Beschreibung: [1]  This technique paper describes a novel method for quantitatively and routinely identifying auroral breakup following substorm onset using the Time History of Events and Macroscale Interactions During Substorms (THEMIS) all-sky imagers (ASIs). Substorm onset is characterised by a brightening of the aurora that is followed by auroral poleward expansion and auroral breakup. This breakup can be identified by a sharp increase in the auroral intensity i(t) and the time derivative of auroral intensity i'(t) . Utilising both i(t) and i'(t) we have developed an algorithm for identifying the time interval and spatial location of auroral breakup during the substorm expansion phase within the field of view of ASI data based solely on quantifiable characteristics of the optical auroral emissions. We compare the time interval determined by the algorithm to independently identified auroral onset times from three previously published studies. In each case the time interval determined by the algorithm is within error of the onset independently identified by the prior studies. We further show the utility of the algorithm by comparing the breakup intervals determined using the automated algorithm to an independent list of substorm onset times. We demonstrate that up to 50% of the breakup intervals characterised by the algorithm are within the uncertainty of the times identified in the independent list. The quantitative description and routine identification of an interval of auroral brightening during the substorm expansion phase provides a foundation for unbiased statistical analysis of the aurora to probe the physics of the auroral substorm as a new scientific tool for aiding the identification of the processes leading to auroral substorm onset.

Beschreibung: [1]  We applied the Grad-Shafranov (GS) reconstruction technique to Martian magnetic flux ropes observed by Mars Global Surveyor in order to estimate their spatial structures. This technique can provide a magnetic field map of their cross section from single spacecraft data, under the assumption that the structure is two-dimensional, magneto-hydrostatic, and time-independent. We succeeded in recovering the spatial structure for 70 events observed between April 1999 and November 2006. The reconstruction results indicate that the flux rope axes were mostly oriented horizontal to the Martian surface, and were randomly distributed with respect to the typical plasma streamline. A subset of events with duration longer than 240 sec was observed at solar zenith angles larger than 75 deg. These events all occur downstream from strong crustal magnetic field in the southern hemisphere, indicating an association between the crustal fields and the detected flux ropes. Using the shape and size of the flux ropes obtained from the GS reconstruction, we estimate lower limits on their volume that span 2–3 orders of magnitude, with larger flux ropes observed downstream from strong crustal magnetic fields. Estimated ion escape rates associated with flux ropes are of the order of 10 22 –10 23 ion/sec, being approximately 10% of previously estimated escape rates during solar minimum.

Beschreibung: [1]  From a survey of Polar plasma waves conducted over the interval 1 April, 1996 to 4 April, 1997 (during solar minimum) at and inside the plasmasphere, magnetosonic waves were detected at all local times with a slight preference of occurrence in the midnight-postmidnight sector at L = 3 to 4. The waves occurred primarily during heightened geomagnetic (AE) activity. Wave occurrence (and intensities) peaked at ~ ±5° of the magnetic equator, with half-maxima at ~ ±10°. For other wave events, magnetosonic waves were also detected as far from the equator as +20° and -60° MLAT, but at lower intensities. An extreme magnetosonic wave intensity event of amplitude B w  = ~ ± 1 nT and E w  = ~ ± 25 mV/m was detected during the survey period. The event occurred near local midnight (0022 MLT), at the magnetic equator (MLAT = -0.5°), at the plasmapause (L = 3.5), and during an intense substorm/convection event (AE = 624 nT; SYM-H = -33 nT). If more stringent requirements (| MLAT| ≤ 5° and AE 〉 300 nT) are imposed, the wave occurrence rate approaches ~50% for the 23 to 00 MLT bin at L = 3 to 4. This strong local time anisotropy in the location of magnetosonic wave occurrence rate supports the idea of generation by protons injected from the plasmasheet into the midnight sector magnetosphere by substorm electric fields. Magnetosonic waves were also detected near late morning (1031 MLT) during relative geomagnetic quiet (low AE). We mention that one possible generation mechanism is a recovering/expanding plasmasphere engulfing preexisting energetic ions, which in turn leads to an ion instability. The wave magnetic component oscillations are aligned along B 0 , the ambient magnetic field direction, and the electric component oscillations are orthogonal to B 0 , indicating linear polarization. The magnetosonic wave amplitudes decreased at locations further from the magnetic equator, while transverse whistler mode wave amplitudes increased. We argue that intense magnetosonic waves are always present somewhere in the magnetosphere during strong substorm/convection events. We thus suggest that modelers use dynamic particle tracing codes and the maximum (rather than average) wave amplitudes to simulate wave-particle interactions.

Beschreibung: The magnetohydrodynamic Richtmyer-Meshkov instability is investigated for the case where the initial magnetic field is unperturbed and aligned with the mean interface location. For this initial condition, the magnetic field lines penetrate the perturbed density interface, forbidding a tangential velocity jump and therefore the presence of a vortex sheet. Through simulation, we find that the vorticity distribution present on the interface immediately after the shock acceleration breaks up into waves traveling parallel and anti-parallel to the magnetic field, which transport the vorticity. The interference of these waves as they propagate causes the perturbation amplitude of the interface to oscillate in time. This interface behavior is accurately predicted over a broad range of parameters by an incompressible linearized model derived presently by solving the corresponding impulse driven, linearized initial value problem. Our use of an equilibrium initial condition results in interface motion produced solely by the impulsive acceleration. Nonlinear compressible simulations are used to investigate the behavior of the transverse field magnetohydrodynamic Richtmyer-Meshkov instability, and the performance of the incompressible model, over a range of shock strengths, magnetic field strengths, perturbation amplitudes and Atwood numbers.

Beschreibung: We investigate experimentally the statistical properties of bedload transport induced by a steady, uniform, and laminar flow. We focus chiefly on lateral transport. The analysis is restricted to experiments where the flow-induced shear stress is just above the threshold for sediment transport. We find that, in this regime, the concentration of moving particles is low enough to neglect interactions between themselves. We can therefore represent bedload as a thin layer of independent walkers travelling over the bed surface. In addition to their downstream motion, the particles show significant fluctuations of their cross-stream velocity, likely due to the roughness of the underlying sediment bed. This causes particles to disperse laterally. Based on thousands of individual trajectories, we show that this lateral spreading is the manifestation of a random walk. The experiments are entirely consistent with Fickian diffusion.

Beschreibung: We study the behavior of heavy inertial particles in the flow field of two like-signed vortices. In a frame co-rotating with the two vortices, we find that stable fixed points exist for these heavy inertial particles; these stable frame-fixed points exist only for particle Stokes number St 〈 St cr . We estimate St cr and compare this with direct numerical simulations, and find that the addition of viscosity increases the St cr slightly. We find that the rate at which particles fall into the fixed points increases until the fixed points disappear at St = St cr . These frame-fixed points are between fixed points and limit cycles in character.

Beschreibung: [1]  On January 9, 2002 and November 14, 2001, the São Luís 30 MHz coherent backscatter radar observed unusual day-time echoes scattered from the equatorial electrojet. The electrojet echoing layers on these days, as seen in the range time intensity (RTI) maps, exhibited quasi-periodic oscillations. Time-frequency decomposition of the magnetic field perturbations ΔH , measured simultaneously by the ground-based magnetometers, also showed evidence of short period waves. The ground-based observations were aided by measurements of the brightness temperature in the water vapor and infrared bands made by the GOES 8 satellite. The GOES 8 satellite measurements indicated evidence of deep tropospheric convection activities, which are favorable for the launch of atmospheric gravity waves (AGW) near São Luís. Our multi-technique investigation, combined with an analysis of the equatorial electric field and current density, indicates that AGW forcing could have been responsible, via coupling with E-region electric fields, for the short-period electrojet oscillations observed over São Luís.

Beschreibung: [1]  The diurnal variation of the global electric circuit is investigated using the World Wide Lightning Location Network (WWLLN), which has been shown to identify nearly all thunderstorms ( [16], using WWLLN data from 2005). To create an estimate of global electric circuit activity, a clustering algorithm is applied to the WWLLN dataset to identify global thunderstorms from 2010 – 2013. Annual, seasonal, and regional thunderstorm activity is investigated in this new WWLLN thunderstorm dataset in order to estimate the source behavior of the global electric circuit. Through the clustering algorithm, the total number of active thunderstorms are counted every 30 minutes creating a measure of the global electric circuit source function. The thunderstorm clusters are compared to precipitation radar data from the Tropical Rainfall Measurement Mission satellite and with case studies of thunderstorm evolution. [2]  The clustering algorithm reveals an average of 660 ± 70 thunderstorms active at any given time with a peak-to-peak variation of 36%. The highest number of thunderstorms occurs in November (720 ± 90) and the lowest number occurs in January (610 ± 80). Thunderstorm cluster and electrified storm cloud activity are combined with thunderstorm overflight current measurements to estimate the global electric circuit thunderstorm contribution current to be 1090 ± 70 A with a variation of 24%. By utilizing the global coverage and high time resolution of WWLLN, the total active thunderstorm count and current is shown to be less than previous estimates based on compiled climatologies.

Beschreibung: [1]  Dynamics of the dayside magnetosphere and proton radiation belt was analyzed during unusual magnetic storm on 21 January 2005. We have found that during the storm from 1712 to 2400 UT, the subsolar magnetopause was continuously located inside geosynchronous orbit due to strong compression. The compression was found to be extremely strong from 1846 to 2035 UT when the dense plasma of fast erupting filament produced the solar wind dynamic pressure Pd peaked up to 〉100 nPa and, in the first time, the upstream solar wind was observed at geosynchronous orbit during almost 2 hours. Under the extreme compression, the outer magnetosphere at L  〉 5 was pushed inward and the outer radiation belt particles with energies of several tens of keV moved earthward, became adiabatically accelerated and accumulated in the inner magnetosphere at L  〈 4 that produced the intensified ring current with an exceptionally long lifetime. The observations were compared with predictions of various empirical and first principles models. All the models failed to predict the magnetospheric dynamics under the extreme compression when the minimal magnetopause distance was estimated to be ~3 Re. The inconsistencies between the model predictions and observations might result from distortions of plasma measurements by extreme heliospheric conditions consisting in very fast solar wind streams (~1000 km/s) and intense fluxes of solar energetic particles. We speculated that anomalous dynamics of the magnetosphere could be well described by the models if the He abundance in the solar wind was assumed to be 〉20%, which is well appropriate for erupting filaments and which is in agreement with the upper 27% threshold for the He/H ratio obtained from Cluster measurements.

Beschreibung: The paper reports a new phenomenon—vortex flows in isothermal magnetic fluids in the vicinity of the localized source of magnetic field (magnetized iron sphere) induced by the drift of drop-like aggregates. Although the observed magnetic precipitation of drop-like aggregates resembles an ordinary rainfall in the Earth atmosphere, its origin and nature are quite different. In magnetic fluids this “rain” is induced by the non-uniform magnetic field and occurs at the scale of 1 mm, not at the scale of several kilometers as in the Earth atmosphere. The reason of this phenomenon is that the applied magnetic field initiates phase transition of “gas-liquid” type which is accompanied by formation of condensed phase represented by drop-like aggregates with the characteristic dimension of about tens of micrometers elongated along the field lines. Inhomogeneous spatial distribution of drop-like aggregates leads to deviation of the ponderomotive force, which is responsible for the formation of vortex flows in the fluid. The “rain” is the primary reason for the vortex flows and it lasts until all magnetic particles capable of condensing into drop-like aggregates precipitate at the surface of the condensation core (iron sphere). Thus, vortex flows induced by drop-like aggregate magnetophoresis represent one variant of “gas-liquid” phase transition. Hydrodynamic flows intensify mass transfer in vicinity of magnetic condensation core and considerably speed it up.

Beschreibung: [1]  Type III radio bursts are produced near the local electron plasma frequency and/or near its harmonic by fast electrons ejected from the solar active regions and moving through the corona and solar wind. These bursts have dynamic spectra with frequency rapidly falling with time. This paper presents two new methods developed to detect type III bursts automatically in the data from High Frequency Receiver (HFR) of the STEREO/WAVES (S/WAVES) radio instrument onboard the STEREO spacecraft. The first technique is applicable to the low frequency band (HFR-1: 125 kHz to 1.975 MHz) only. This technique can possibly be implemented in on-board satellite software aimed at preliminary detection of bursts and identification of time intervals with relatively high solar activity. In the second technique the bursts are detected in both the low frequency band and the high frequency band (HFR-2: 2.025 MHz to 16.025 MHz), with the computational burden being higher by one order of magnitude as compared with that for the first technique. Preliminary tests of the method show that the performance of the first technique is quite high, P dL  = 72 % ± 3 %. The performance of the second technique is considerably higher, P dL  +  H  = 81 % ± 1%, while the number of false alarms does not exceed 10% for one daily spectrum.

Beschreibung: [1]  A dipolarizing flux bundle (DFB) is a small magnetotail flux tube (typically 〈 ~3 R E in X GSM and Y GSM ) with a significantly more dipolar magnetic field than its background. Dipolarizing flux bundles typically propagate earthward at a high speed from the near-Earth reconnection region. Knowledge of a DFB's flux transport properties leads to better understanding of near-Earth (X = -6 to -30 R E ) magnetotail flux transport and thus conversion of magnetic energy to kinetic and thermal plasma energy following magnetic reconnection. We explore DFB properties with a statistical study using data from the Time History of Events and Macroscale Interactions during Substorms (THEMIS) mission. To establish the importance of DFB flux transport, we compare it with transport by bursty bulk flows (BBFs) that typically envelop DFBs. Because DFBs coexist with flow bursts inside BBFs, they contribute 〉65% of BBF flux transport, even though they last only ~30% as long as BBFs. The rate of DFB flux transport increases with proximity to Earth and to the pre-midnight sector, as well as with geomagnetic activity and distance from the neutral sheet. Under the latter two conditions the total flux transport by a typical DFB also increases. Dipolarizing flux bundles appear more often during increased geomagnetic activity. Since BBFs have been previously shown to be the major flux transporters in the tail, we conclude that DFBs are the dominant drivers of this transport. The occurrence rate of DFBs as a function of location and geomagnetic activity informs us about processes that shape global convection and energy conversion

Beschreibung: [1]  The substorm current wedge (SCW) is a fundamental component of geomagnetic substorms. Models tend to describe the SCW as a simple line current flowing into the ionosphere towards dawn and out of the ionosphere towards dusk, linked by a westward electrojet. We use multi-spacecraft observations from perigee passes of the Cluster 1 and 4 spacecraft during a substorm on 15 Jan 2010, in conjunction with ground-based observations, to examine the spatial structuring and temporal variability of the SCW. At this time, the spacecraft travelled east-west azimuthally above the auroral region. We show that the SCW has significant azimuthal sub-structure on scales of 100 km at altitudes of 4,000-7,000 km. We identify 26 individual current sheets in the Cluster 4 data and 34 individual current sheets in the Cluster 1 data, with Cluster 1 passing through the SCW 120-240 s after Cluster 4 at 1,300-2,000 km higher altitude. Both spacecraft observed large-scale regions of net upward and downward field-aligned current, consistent with the large-scale characteristics of the SCW, although sheets of oppositely directed currents were observed within both regions. We show that the majority of these current sheets were closely aligned to a north-south direction,in contrast to the expected east-west orientation of the pre-onset aurora. Comparing our results with observations of the field-aligned current associated with bursty bulk flows (BBFs) we conclude that significant questions remain for the explanation of SCW structuring by BBF driven “wedgelets”. Our results therefore represent constraints on future modelling and theoretical frameworks on the generation of the SCW.

Beschreibung: [1]  Large horizontal winds and wind shears have been measured in the lower thermosphere by rockets, lidars, and non-specular meteor radars. This paper describes a detailed analysis of 3 multi-hour non-specular meteor radar data sets collected at the Jicamarca Radio Observatory. This provides some of the highest resolution sustained measurements in this part of the atmosphere. These show: (1) intense wind speeds, maintaining 180 m/s for half an hour and 160 m/s for another half an hour; (2) winds structured in layers that move up or, more commonly, down in the pre-dawn hours at rates of a few km/hr; (3) intense wind shears that typically persist at around 50 m/s/km but, in one instance, sustains values approaching 100 m/s/km for a few hours.

Beschreibung: [1]  We present results from an analysis of high-latitude ionosphere-thermosphere (IT) coupling to the solar wind during a moderate magnetic storm which occurred on 5-6 August 2011. During the storm, a multi-point set of observations of the ionosphere and thermosphere was available. We make use of ionospheric measurements of electromagnetic and particle energy made by the Defense Meteorological Satellite Program (DMSP), and neutral densities measured by the Gravity Recovery and Climate Experiment (GRACE) satellite to infer: (1) the energy budget and (2) timing of the energy transfer process during the storm. We conclude that the primary location for energy input to the IT system may be the extremely high latitude region. We suggest that the total energy available to the IT system is not completely captured either by observation or empirical models.

Beschreibung: [1]  The level of solar activity varies from cycle to cycle. This variability is probably caused by a combination of nonlinear and random effects. Based on surface flux transport simulations, we show that the observed inflows into active regions and towards the activity belts provide an important nonlinearity in the framework of Babcock-Leighton model for the solar dynamo. Inclusion of these inflows also leads to a reproduction of the observed proportionality between the open heliospheric flux during activity minima and the maximum sunspot number of the following cycle. A substantial component of the random variability of the cycle strength is associated with the cross-equatorial flux plumes that occur when large, highly-tilted sunspot groups emerge close to the equator. We show that the flux transported by these events is important for the amplitude of the polar fields and open flux during activity minima. The combined action of inflows and cross-equatorial flux plumes provides an explanation for the weakness of the polar fields at the end of solar cycle 23 (and hence for the relative weakness of solar cycle 24).

Beschreibung: [1]  We present the first direct measurement of neutral oxygen in the lunar exosphere, detected by the Chandrayaan-1 Energetic Neutral Analyzer (CENA). With the lunar surface consisting of about 60% of oxygen in number, the neutral oxygen detected in CENA's energy range (11 eV – 3.3 keV) is attributed to have originated from the lunar surface, where it is released through solar wind ion sputtering. We verify this proposition by comparing the measured oxygen content in two different mass spectra groups with ion sputtering theory. One group contains mass spectra that were recorded when the solar wind consisted of almost pure hydrogen and the other group contains mass spectra that were recorded when the helium content in the solar wind was very high (〉3.5%). Since helium is a much more effective sputtering agent than hydrogen (5% of alpha particles present in the solar wind typically contribute 30% of the total sputter yield), these two groups should show clear differences in the oxygen sputter yield. Fitting of CENA's mass spectra with calibration spectra from ground and in-flight data resulted in the detection of a robust oxygen signal, with a flux of 0.2 to 0.4 times the flux of backscattered hydrogen, depending, as expected, on the solar wind helium content and particle velocity. These measurements present the first in-situ detection of oxygen in the lunar exosphere. For the two solar wind types observed, we derive sub-solar surface oxygen atom densities of N 0  = (1.1 ± 0.3) ⋅ 10 7 m − 3 and (1.4 ± 0.4) ⋅ 10 7 m − 3 , respectively, which agree well with earlier model predictions and measured upper limits. From these surface densities we derive, by modeling, column densities of N C  = (1.5 ± 0.5) ⋅ 10 13 m − 2 and (1.6 ± 0.5) ⋅ 10 13 m − 2 . [2]  In addition, in the CENA mass spectra, we identified for the first time a helium component. This helium is attributed to backscattering of solar wind helium (alpha particles) from the lunar surface as neutral energetic helium atoms, which has been observed for the first time. This identification is supported by the characteristic energy of the measured helium atoms, which is roughly four times the energy of reflected solar wind hydrogen, and the correlation with solar wind helium content.

Beschreibung: [1]  The diffuse aurora is an almost permanent feature in the Earth's upper polar atmosphere, providing the major source of ionizing energy input into the high-latitude region. Previous theoretical and observational studies have demonstrated that whistler-mode chorus scattering primarily accounts for intense nightside diffuse auroral precipitation within ~ 8 R E , but what causes the dayside diffuse aurora remains poorly understood. Using conjugate satellite wave and particle observations on 13 August 2009 from the THEMIS spacecraft and ground-based all-sky imager measurements at the South Pole on the dayside, we perform a quantitative analysis of wave driven diffusion and electron precipitation. Our results demonstrate that the dayside chorus scattering was the dominant contributor to the observed dayside diffuse auroral precipitation and that the chorus wave intensity primarily controlled its brightness, indicating that dayside chorus can be the major driver of the Earth's dayside diffuse aurora. While further investigations are required to bring closure to the origin of the dayside diffuse aurora under differing solar wind conditions and geomagnetic situations, our finding is an important complement to recent work on the formation mechanism of the diffuse aurora and provides improved understanding of the roles of resonant wave-particle interactions in diffuse auroral precipitation pattern on a global scale.

Beschreibung: [1]  Ultra Low Frequency (ULF: 0.001-5 Hz) magnetic records have recently been used in the search for short term earthquake prediction methods. The separation of local and global effects in the magnetic records is the greatest challenge in this research area. Geomagnetic indices are often used to predict global ULF magnetic behavior where it is assumed that increases in a geomagnetic index correspond with an increase in ULF power. This paper examines the relationships between geomagnetic indices and ULF power, spectral polarization ratio and the relationship between the spectral polarization ratio and solar wind parameters. The power in the ULF, Pc3-5 bands (10-600 s) shows a linear correlation coefficient of  0.2 with the Kp magnetic activity index. The correlation varies with magnetic local time (MLT) and latitude. The correlation coefficient is inversely related to the integrated power in the ULF Pc3 band (10-45 s) over MLT and magnetic latitude. The ratio of spectral powers Z ( ω )/ G ( ω ) is discussed and shown to be a promising parameter in the search for earthquake precursor signals in ULF records.

Beschreibung: [1]  Recent testing of a quantitative model describing the classical (region-1-sense, referred to as the R1 current loop) substorm current wedge (SCW) revealed systematic discrepancies between the observed and predicted amplitudes, which suggested us to include additional region-2-sense currents (R2-loop) earthward of the dipolarized region (SCW2L model). In this paper we discuss alternative circuit geometries of the 3d substorm current system and interpret simultaneous observations of the magnetic field dipolarizations by NOAA Geostationary Operational Environmental Satellite (GOES) and by NASA Time History of Events and Macroscale Interactions during Substorms (THEMIS) spacecraft, to quantitatively investigate the SCW2L model parameters. During two cases of a dipole-like magnetotail configuration, the dipolarization/injection front fortuitously stopped at r  ∼ 9 Re for the entire duration of ∼ 30 min-long SCW-related dipolarization within a unique, radially-distributed multi-spacecraft constellation, which allowed us to determine the locations and total currents of both SCW2L loops. In addition, we conducted a survey of the dipolarization amplitudes in events, simultaneously observed at 6.6 Re (GOES) and 11 Re (THEMIS) under a wide range of magnetotail conditions. We infer that the ratio I 2 / I 1 varies in the range 0.2 to 0.6 (median value 0.4) and that the equatorial part of the R2 current loop stays at the distance r  〉 6.6 Re in the case of a dipole-like field geometry ( BZ 0  〉 75 nT at 6.6 Re prior to the onset), but it is located at r  〈 6.6 Re in the case of a stretched magnetic field configuration (with BZ 0  〈 60 nT). Since the ground midlatitude perturbations are sensitive to the combined effect of the R1- and R2-sense current loops with the total current roughly equal to I 1  −  I 2 , the ratio I 2 / I 1 becomes an important issue when attempting to monitor the current disruption intensity from ground observations.

Beschreibung: Compressive sampling is well-known to be a useful tool used to resolve the energetic content of signals that admit a sparse representation. The broadband temporal spectrum acquired from point measurements in wall-bounded turbulence has precluded the prior use of compressive sampling in this kind of flow, however it is shown here that the frequency content of flow fields that have been Fourier transformed in the homogeneous spatial (wall-parallel) directions is approximately sparse, giving rise to a compact representation of the velocity field. As such, compressive sampling is an ideal tool for reducing the amount of information required to approximate the velocity field. Further, success of the compressive sampling approach provides strong evidence that this representation is both physically meaningful and indicative of special properties of wall turbulence. Another advantage of compressive sampling over periodic sampling becomes evident at high Reynolds numbers, since the number of samples required to resolve a given bandwidth with compressive sampling scales as the logarithm of the dynamically significant bandwidth instead of linearly for periodic sampling. The combination of the Fourier decomposition in the wall-parallel directions, the approximate sparsity in frequency, and empirical bounds on the convection velocity leads to a compact representation of an otherwise broadband distribution of energy in the space defined by streamwise and spanwise wavenumber, frequency, and wall-normal location. The data storage requirements for reconstruction of the full field using compressive sampling are shown to be significantly less than for periodic sampling, in which the Nyquist criterion limits the maximum frequency that can be resolved. Conversely, compressive sampling maximizes the frequency range that can be recovered if the number of samples is limited, resolving frequencies up to several times higher than the mean sampling rate. It is proposed that the approximate sparsity in frequency and the corresponding structure in the spatial domain can be exploited to design simulation schemes for canonical wall turbulence with significantly reduced computational expense compared with current techniques.

Beschreibung: This paper presents extensive acoustic measurements on jets impinging on surfaces of various surface roughness values. Besides surface roughness, the effects of nozzle-to-plate spacing distance and nozzle pressure ratio are also investigated. Turbulent mixing noise and tonal noise are explained using far-field wall-jet velocity and impingement region temperature fields. The results demonstrate that roughness of the impingement plate widens the staging region of impingement noise. In general, high speed jet impinging on a rough plate generates less noise compared to a smooth plate. When tones are removed from the spectra, it is found that acoustic power monotonically decreases with increasing surface roughness. Thermal imaging in the stagnation region indicates that whenever tones are present, the temperature at the stagnation region is high. Further, sound pressure directivity pattern of impingement noise is constructed by superimposing a wall-jet and a free jet in the appropriate orientations.

Beschreibung: We investigate the elasticity of an isolated, threefold junction of soap films (Plateau border), which displays static undulations when liquid rapidly flows into it. By analyzing the shape of the Plateau border (thickness R and transverse displacement) as a function of the liquid flow rate Q , we show experimentally and theoretically that the elasticity of the Plateau border is dominated by the bending of the soap films pulling on the Plateau border. In this asymptotic regime, the undulation wavelength obeys the scaling law ∼ Q 2   R −2 and the decay length ∼ Q 2   R −4 .

Beschreibung: We present Cluster observations of wave-particle interactions during an earthward-propagating dipolarization front (DF) and associated fast plasma bulk flows detected at the central current sheet in Earth's magnetotail. During this period, flux tubes behind the DF frequently contain more energetic or hotter ions than did the pre-existing flux tubes ahead of the DF. On the other hand, electrons within the DF flux tubes heat less, or are even colder, than were the pre-existing populations, and are often accompanied by superposed isolated beams. At the same time, electrostatic emissions are strongly enhanced over a wide range of frequencies (up to several times the electron cyclotron frequency) behind the DFs. This low-frequency electrostatic wave power is well correlated with ion energization. From linear theory, we find two wave modes: a high-frequency beam mode and a low-frequency whistler mode that are associated with the electron beam component. We attribute the generation of whistlers to electron beams that persist for a while before undergoing rapid thermalization. The existence of isolated beam components behind DFs detected during the 4-sec Cluster spin period indicates that DFs either provide a continuous source of electron beams or facilitate a physical process that maintains the beams against rapid thermalization. Our analysis suggests that the earthward motion of the DF flux tube, via Fermi acceleration as the magnetic field lines behind the DF shorten, can lead to the persistent electron beams that generate whistler mode waves, which in turn can heat ions. This scenario, by which free energy in electron beams generates waves that then heat ions, accounts for the Cluster observations of different energization behaviors between electrons and ions behind DFs.

Beschreibung: We examine a unique data set from seven Hubble Space Telescope (HST) ‘visits’ that imaged Saturn's northern dayside ultraviolet emissions exhibiting usual circumpolar ‘auroral oval’ morphologies, during which Cassini measured the interplanetary magnetic field (IMF) upstream of Saturn's bow shock over intervals of several hours. The auroras generally consist of a dawn arc extending towards noon centered near ~15º co-latitude, together with intermittent patchy forms at ~10º co-latitude and poleward thereof, located between noon and dusk. The dawn arc is a persistent feature, but exhibits variations in position, width, and intensity, which have no clear relationship with the concurrent IMF. However, the patchy post-noon auroras are found to relate to the (suitably lagged and averaged) IMF B z , being present during all four visits with positive B z and absent during all three visits with negative B z . The most continuous such forms occur in the case of strongest positive B z . These results suggest that the post-noon forms are associated with reconnection and open flux production at Saturn's magnetopause, related to the similarly-interpreted bifurcated auroral arc structures previously observed in this LT sector in Cassini UVIS data, whose details remain unresolved in these HST images. One of the intervals with negative IMF B z , however, exhibits a pre-noon patch of very high latitude emission extending poleward of the dawn arc to the magnetic/spin pole, suggestive of the occurrence of lobe reconnection. Overall, these data provide evidence of significant IMF-dependence in the morphology of Saturn's dayside auroras.

Beschreibung: A three-dimensional (3D) particle-in-cell (PIC) simulation of the whistler anisotropy instability is carried out for a collisionless, homogeneous, magnetized plasma with β e  = 0.10. This is the first 3D PIC simulation of the evolution of enchanced fluctuations from this growing mode driven by an anisotropic electron velocity distribution with T ⊥  e / T ‖ e  〉 1 where ⊥ and ∥ represent directions perpendicular and parallel to the background magnetic field B o , respectively. The early-time magnetic fluctuation spectrum grows with properties reflecting the predictions of linear theory with narrowband maxima at kc / ω e  ≃ 1 and k  ×  B o  = 0, and a wavevector anisotropy in the sense of k ⊥  〈 〈  k ∥ . Here ω e represents the electron plasma frequency. At later times the fluctuations undergo both a forward transfer to shorter wavelengths, also with k ⊥  〈 〈  k ∥ , and an inverse transfer to longer wavelengths with wavevector anisotropy k ⊥  〉 〉  k ∥ . The inverse transfer is consistent with a prediction of nonlinear three-wave coupling theory.

Beschreibung: Auroral ionospheric F-region density depletions observed by PFISR (Poker Flat Incoherent Scatter Radar) during the MICA (Magnetosphere-Ionosphere Coupling in the Alfvén Resonator) sounding rocket campaign are critically examined alongside complementary numerical simulations. Particular processes of interest include cavity formation due to intense frictional heating and Pedersen drifts, evolution in the presence of structured precipitation, and refilling due to impact ionization and downflows. Our analysis uses an ionospheric fluid model which solves conservation of mass, momentum, and energy equations for all major ionospheric species. These fluid equations are coupled to an electrostatic current continuity equation to self-consistently describe auroral electric fields. Energetic electron precipitation inputs for the model are specified by inverting optical data, and electric field boundary conditions are obtained from direct PFISR measurements. Thus, the model is driven in as realistic a manner as possible. Both ISR data and simulations indicate that the conversion of the F-region plasma to molecular ions and subsequent recombination is the dominant process contributing to the formation of the observed cavities, all of which occur in conjuction with electric fields exceeding ∼ 90 mV/m. Furthermore, the cavities often persist several minutes past the point when the frictional heating stops. Impact ionization and field-aligned plasma flows modulate the cavity depth in a significant way, but are of secondary importance to the molecular generation process. Informal comparisons of the ISR density and temperature fits to the model verify that the simulations reproduce most of the observed cavity features to a reasonable level of detail.

Beschreibung: Energetic particle injections in the near-Earth plasma sheet are critical for supplying particles and energy to the inner magnetosphere. Recent case studies have demonstrated a good correlation between injections and transient, narrow, fast flow channels as well as earthward reconnection (dipolarization) fronts in the magnetotail, but statistical observations beyond geosynchronous orbit (GEO) to verify the findings were lacking. By surveying trans-geosynchronous injections using THEMIS, we show that their likely origin is the earthward-traveling, dipolarizing flux bundles following near-Earth reconnection. The good correlation between injections and fast flows, reconnection fronts and impulsive, dawn-dusk electric field increases is not limited to within 12 R E , but extends out to 30 R E . Like near-Earth reconnection, both ion and electron injections are most probable in the pre-midnight sector. Similar to bursty bulk flows (BBFs), injection-time flow speeds are faster with increasing distance from Earth. With faster flows, injection intensity generally increases and extends to higher energy channels. With increased geomagnetic activity, injection occurrence rate increases (akin to that of BBFs) and spectral hardening occurs (κ decreases). The occurrence rate increase within the inner magnetosphere suggests that injections populate the radiation belts more effectively under enhanced geomagnetic activity. Our results are inconsistent with the classical concept of an azimuthally wide injection boundary moving earthward from ~9-12 R E to GEO under an enhanced cross-tail electric field. Rather, particle injection and transport occur along a large range of radial distances due to effects from earthward-penetrating, azimuthally localized, transient, strong electric fields of recently reconnected, dipolarizing flux bundles.

Beschreibung: A method for estimating the vector neutral wind profiles in the mesosphere and lower thermosphere (MLT) region of the upper atmosphere from Arecibo dual-beam incoherent scatter radar data is presented. The method yields continuous estimates of both the altitude-averaged F -region plasma drifts and all three components of the altitude-resolved neutral wind profiles in the MLT using data taken while the Arecibo feed system swings in azimuth. The problem is mixed determined, and its solution is not inherently unique. Second order Tikhonov regularization is used to find solutions consistent with the available data while being minimally structured, additional structure being unsupported by the data. The solution is found using the method of conjugate gradient least squares and sparse matrix mathematics. Example data acquired during an interval of midlatitude spread F are used to illustrate the method. The estimated wind profiles exhibit characteristics broadly consistent with gravity waves but are impulsive, with features that generally persist for less than one and a half wave periods.

Beschreibung: The case of a turbulent round jet impinging perpendicularly onto a rotating, heated disc is investigated, in order to understand the mechanisms at the origin of the influence of rotation on the radial wall jet and associated heat transfer. The present study is based on the complementary use of an analysis of the orders of magnitude of the terms of the mean momentum and Reynolds stress transport equations, available experiments, and dedicated Reynolds-averaged Navier–Stokes computations with refined turbulence models. The Reynolds number Re j = 14 500, the orifice-to-plate distance H = 5 D , where D is the jet-orifice diameter, and the four rotation rates were chosen to match the experiments of Minagawa and Obi [“Development of turbulent impinging jet on a rotating disk,” Int. J. Heat Fluid Flow25, 759–766 (2004)] and comparisons are made with the Nusselt number distribution measured by Popiel and Boguslawski [“Local heat transfer from a rotating disk in an impinging round jet,” J. Heat Transfer108, 357–364 (1986)], at a higher Reynolds number. The overestimation of turbulent mixing in the free-jet before the impact on the disk is detrimental to the prediction of the impingement region, in particular of the Nusselt number close to the symmetry axis, but the self-similar wall jet developing along the disk is correctly reproduced by the models. The analysis, experiments, and computations show that the rotational effect do not directly affect the outer layer, but only the inner layer of the wall jet. A noteworthy consequence is that entrainment at the outer edge of the wall jet is insensitive to rotation, which explains the dependence of the wall-jet thickness on the inverse of the non-dimensional rotation rate, observed in the experiments and the Reynolds stress model computations, but not reproduced by the eddy-viscosity models, due to the algebraic dependence to the mean flow. The analysis makes moreover possible the identification of a scenario for the appearance of rotational effects when the rotation rate is gradually increased. For weak rotation rates, the rotation-induced boundary layer appears but does not break the self-similar solution observed for the case without rotation. For intermediate rotation rates, the production of the azimuthal Reynolds stress becomes much stronger than other components, leading to a complete modification of the turbulence anisotropy which is reproduced only by Reynolds stress models. For strong rotation rates, centrifugal effects dominate, leading to an acceleration and thinning of the jet, and consequently an increase of turbulent production and heat transfer, reproduced by all the turbulence models.

Beschreibung: Microfluidic channels are powerful means of control of minute volumes such as droplets. These droplets are usually conveyed at will in an externally imposed flow which follows the geometry of the micro-channel. It has recently been pointed out by Dangla et al. [“Trapping microfluidic drops in wells of surface energy,” Phys. Rev. Lett.107(12), 124501 (2011)] that the motion of transported droplets may also be stopped in the flow, when they are anchored to grooves which are etched in the channels top wall. This feature of the channel geometry explores a direction that is usually uniform in microfluidics. Herein, this anchoring effect exploiting the three spatial directions is studied combining a depth averaged fluid description and a geometrical model that accounts for the shape of the droplet in the anchor. First, the presented method is shown to enable the capture and release droplets in numerical simulations. Second, this tool is used in a numerical investigation of the physical mechanisms at play in the capture of the droplet: a localized reduced Laplace pressure jump is found on its interface when the droplet penetrates the groove. This modified boundary condition helps the droplet cope with the linear pressure drop in the surrounding fluid. Held on the anchor the droplet deforms and stretches in the flow. The combination of these ingredients leads to recover the scaling law for the critical capillary number at which the droplets exit the anchors C a ★ ∝ h 2 / R 2 where h is the channel height and R the droplet undeformed radius.

Beschreibung: We report the peculiar interaction of two type III bursts observed in the solar wind. As electrons beams propagating on the same magnetic field lines cross, a spectacular depletion of the type III radio emission is observed. We combine observations from the WAVES experiment on board the STEREO mission together with kinetic plasma simulations to study the extinction of type III radio emission resulting from the interaction between two electron beams. The remote observations enable to follow the electron beams in the interplanetary medium and show that the level of radiated radio waves is recovered after the beam crossing. The in situ observations of beam-driven Langmuir waves give evidence for Langmuir decay. The density fluctuations are extracted from in situ observations. The velocity of the beams is independently evaluated from in situ observations of decaying Langmuir waves and remote radio observations. The kinetic simulations show that the level of beam-driven Langmuir waves is reduced as the two beams cross. We show that the slow beam induced a strong reduction of the quasilinear relaxation of the fast beam, limiting the amplitude of the generated Langmuir waves. Moreover, in the case of two electron beams, the lack of Langmuir waves coherence reduces the efficiency of the Langmuir parametric decay. We thus conclude that the observed depletion of the type III radio 5 is independent of the radio emission mechanism, as long as it depends on the Langmuir amplitude and coherence.

Beschreibung: The solar minimum period between solar cycles 23 and 24 was the longest since the beginning of space-based measurements, and many manifestations of solar activity were unusually low. Thermospheric neutral density was about 30% lower than during the previous solar minimum, but changes in the ionosphere between the two solar minima are more controversial. Solar radiation, geomagnetic activity, and anthropogenic increases in greenhouse gases, can all play a role in these changes. In this paper, we address the latter of these potential contributions, the degree to which secular change driven by greenhouse gases, primarily CO 2 , could be responsible for the observed changes. New 3D model simulations find a global mean density decrease at 400 km of 5.8% between the two recent solar minima, which is larger than earlier 1D model results, and in better agreement with observations. From these model simulations and from other observational work, we estimate that the contribution of secular change to global mean neutral density decrease between the two recent solar minima is less than ~6%. The contribution of secular change to the global average decrease of F-region ionosphere peak density ( N m F 2 ) and altitude ( h m F 2 ), near mid-day, is estimated to be 1.5% and 1.5 km, respectively. However, secular changes in the ionosphere exhibit large variations with local time, geographic location, and season. The mid-day change of N m F 2 seen in the model simulations ranged between +6% and -9%, and the change of h m F 2 ranged between +11 km and -11 km, depending on geographic location.

Beschreibung: We present a numerical study of nanosecond pulsed dielectric barrier discharge (DBD) actuator operating in quiescent air at atmospheric condition. Our study concentrates on plasma discharge induced fluid dynamics and on exploration of parametric space of interest for voltage pulse in an attempt to shed some light into elucidation of the mechanisms whereby the generated shock wave propagates through and affects the external flow. Specifically, a one-dimensional, self-similar, local ionization kinetic model recently developed to predict key parameters of nanosecond pulsed plasma discharge is coupled with the compressible Navier-Stokes equations possibly for the first time. Within the considered range of parameters of the plasma model which is justified for the modeling of surface nanosecond pulsed discharge at atmospheric pressure, our coupled method is able to provide satisfactory prediction of the shock structure generated by the actuator for comparison with experiment, not only in the qualitative shock wave shape but also in quantitative shock front displacement. We provide a comprehensive analysis of the gas heating, shock wave initiation and evolution processes. For example, the characteristic time of the rapid localized heating responsible for shock wave generation, which is yet to be quantified experimentally, is found to be ∼350 ns. We conduct a parametric investigation by varying the peak voltage from 10 kV to 50 kV and rise time from 5 ns to 150 ns. The pressure wave whose behavior is found to be dominated by input voltage amplitude, introduces highly transient, localized disturbance to the quiescent air. In addition, the vortex induced by the shock passage is relatively weak. The interplay of the induced flows by a few successive plasma discharges operating at continuous mode does not appear to be significant, especially at low voltage amplitude.

Beschreibung: In this paper, we investigate the decay of incompressible homogeneous isotropic turbulence in a variable viscosity fluid. The viscosity coefficient is assumed to depend linearly on a scalar, representing either a temperature or a concentration, and obeying a simple advection-diffusion equation. At high Reynolds numbers, Direct Numerical Simulations (DNS) allow us to confirm the validity of Taylor's postulate that the dissipation is independent from the viscosity and its fluctuations. At low Reynolds numbers, we report the presence of extra energy at small scales due to these variable viscosity effects. This implies that the turbulent kinetic energy decreases less rapidly as a function of time in variable viscosity fluids. In order to explain this phenomenon and quantify its importance on the turbulent flow, we propose a statistical approach based on an eddy-damped quasi-normal Markovian (EDQNM) spectral closure which takes into account the nonlinearity introduced by variable viscosity. It is shown that this latter additional term is of constant sign in the energy spectrum equation and reduces the dissipation of the flow as observed. Also, by assuming the dominance of distant interactions between wave numbers, we can propose a simple formula expressing that variable viscosity effects lead to an effective reduction of the mean viscosity proportional to the variance of viscosity fluctuations.

Beschreibung: Liquid droplets flowing through a rectangular microfluidic channel develop a vortical flow field due to the presence of shear forces from the surrounding fluid. In this paper, we present an experimental and computational study of droplet velocities and internal flow patterns in a rectangular pressure-driven flow for droplet diameters ranging from 0.1 to 2 times the channel height. Our study shows excellent agreement with asymptotic predictions of droplet and interfacial velocities for infinitesimally small droplets. As the droplet diameter nears the size of the channel height, the droplet velocity slows significantly, and the changing external flow field causes a qualitative change in the location of internal vortices. This behavior is relevant for future studies of mass transfer in microfluidic devices.

Beschreibung: Laboratory experiments are performed to examine the formation of a crater in sediment by an impinging vertical turbulent jet. Light attenuation and a “depositometer,” which records conductivity through the bed from an array of electrodes, are used to measure the crater depth as a function of space and time. The onset of crater formation and deepening is best characterized in terms of the Rouse number, Rs (proportional to the particle settling speed divided by the centerline jet speed), rather than Shields number, Sh (proportional to the stress divided by the particle weight per unit area). The critical Rouse number, Rs c , is found to increase with the particle Reynolds number, Re p , as a power law with exponent 0.45 ± 0.03 for Re p ranging between 0.6 and 160. For smaller Rs, the crater is observed to deepen at a near-constant speed, while the crater radius remains constant. Bedload transport, measured in terms of the crater deepening speed, is determined to increase as Re p times the difference between Rs c and Rs.

Beschreibung: In this paper we study the planetary magnetic disturbance during the magnetic storm occurring on 05 April 2010 associated with high speed solar wind stream due to a coronal hole following a CME. We separate the magnetic disturbance associated to the ionospheric disturbance dynamo (Ddyn) from the magnetic disturbance associated to the prompt penetration of magnetospheric electric field (DP2). This event exhibits different responses of ionospheric disturbance dynamo in the different longitude sectors (European-African, Asian and American). The strongest effect is observed in the European-African sector. The Ddyn disturbance reduces the amplitude of the daytime H-component at low latitudes during four consecutive days in agreement with the Blanc and Richmond's model of ionospheric disturbance dynamo. The amplitude of Ddyn decreased with time during the four days. We discuss its diverse worldwide effects. The observed signature of magnetic disturbance process in specific longitude sector is strongly dependent on which Earth's side faces the magnetic storms (i.e., there is a different response depending on which longitude sector is at noon when the SSC hits).Finally we determined an average period of 22 hours for Ddyn using wavelet analysis.

Beschreibung: Complex magnetosphere-ionosphere coupling mechanisms result in high latitude irregularities that are difficult to characterize using only Global Navigation Satellite System (GNSS) scintillation measurements. However, GNSS observations combined with physical parameters derived from modeling can be used to study the physics of these irregularities. We have developed a full three dimensional (3D) electromagnetic (EM) wave propagation model called “Satellite-beacon Ionospheric-scintillation Global Model ofthe upper Atmosphere" (SIGMA), to simulate GNSS scintillations. This model eliminates the most significant approximation made by the previous simulation approaches about the correlation length of the irregularity. Thus, for the first time, using SIGMA wecan accomplish scintillation simulations of significantly high fidelity. While the model is global, it is particularly applicable at high latitudes as it accounts for the complicated geometry of the magnetic field lines in these regions. Using SIGMA wesimulate the spatial and temporal variations in the GNSS signal phase and amplitude on the ground. In this paper, we present the model and results from a study to determine the sensitivity of the SIGMA outputs to different input parameters. We have deduced from our sensitivity study that the peak to peak (P2P) power gets most affected by the spectral index and line of sight (LOS) direction, while the P2P phase and standard deviation of the phase ( σ φ ) are more sensitive to the anisotropy of the irregularity. The sensitivity study of SIGMA narrows the parametric space to investigate when comparing the modeled results to the observations.