Compressive tidal heating of a disk galaxy in a rich cluster

The restricted three-body method is used to model the effect of the mean tidal field of a cluster of galaxies on the internal dynamics of a disk galaxy falling into the cluster for the first time. In the model adopted the galaxy experiences a tidal field that is compressive within the core of the cluster. The planar random velocities of all components in the disk increase after the galaxy passes through the core of the cluster. The low-velocity dispersion gas clouds experience a relatively larger increase in random velocity than the hotter stellar components. The increase in planar velocities results in a strong anisotropy between the planar and vertical velocity dispersions. It is argued that this will make the disk unstable to the 'fire-hose instability' which leads to bending modes in the disk and which will thicken the disk slightly. The mean tidal fields in rich clusters were probably stronger during the epoch of cluster formation and relaxation than they are in present-day relaxed clusters.

Implications of rapid footpoint motions of photospheric flux tubes for coronal heating

Some recent observations at Pic-du-Midi (Mulleret al., 1992a) suggest that the photospheric footpoints of coronal magnetic field lines occasionally move rapidly with typical velocities of the order 3 km s–1 for about 3 or 4 min. We argue that such occasional rapid footpoint motions could have a profound impact on the heating of the quiet corona. Qualitative estimates indicate that these occasional rapid motions can account for the entire energy flux needed to heat the quiet corona. We therefore carry out a mathematical analysis to study in detail the response of a vertical thin flux tube to photospheric footpoint motions in terms of a superposition of linear kink modes for an isothermal atmosphere. We find the resulting total energy that is asymptotically injected into an isothermal atmosphere (i.e., an atmosphere without any back reflection). By using typical parameter values for fast and slow footpoint motions, we show that, even if the footpoints spend only 2.5% of the time undergoing rapid motions, still these rapid motions could be more efficient in transporting energy to the corona than the slow motions that take place most of the time.

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]

The scratch hardness and friction of a soft rigid-plastic solid

The paper describes an experimental and analytical study of the normal and scratch hardnesses of a model soft rigid-plastic solid. The material known as ‘Plasticine’, a mixture of dry particles and a mineral oil, has been deformed with a range of rigid conical indentors with included angles of between 30° and 170°. The sliding velocity dependence of the computed scratch hardness and friction has been examined in the velocity range 0.19 mm/s to 7.3 m/s. Data are also described for the time dependence of the normal hardness and also the estimated rate dependence of the intrinsic flow stress. The latter values were estimated from data obtained during the upsetting of right cylinders. Three major conclusions are drawn from these data and the associated analysis. (1) A first-order account of the scratching force may be provided by adopting a model which sums the computed plastic deformation and interfacial sliding contributions to the total sliding work. This is tantamount to the adoption of the two-term non-interacting model of friction. (2) For this system during sliding, at high sliding velocities at least, the interface shear stress which defines the boundary condition is not directly related to the bulk shear stress. The interface rheological characteristics indicate an appreciable dependence on the imposed strain or strain rate. In particular, the relative contributions of the slip and stick boundary conditions appear to be a function of the imposed sliding velocity. (3) The computed normal and scratch hardness values are not simply interrelated primarily because of the evolving boundary conditions which appear to exist in the scratching experiments.

Symmetry Properties of the Large-Deviation Function of the Velocity of a Self-Propelled Polar Particle

A geometrically polar granular rod confined in 2D geometry, subjected to a sinusoidal vertical oscillation, undergoes noisy self-propulsion in a direction determined by its polarity. When surrounded by a medium of crystalline spherical beads, it displays substantial negative fluctuations in its velocity. We find that the large-deviation function (LDF) for the normalized velocity is strongly non-Gaussian with a kink at zero velocity, and that the antisymmetric part of the LDF is linear, resembling the fluctuation relation known for entropy production, even when the velocity distribution is clearly non-Gaussian. We extract an analogue of the phase-space contraction rate and find that it compares well with an independent estimate based on the persistence of forward and reverse velocities.

Unsteady three-dimensional boundary layer flow due to a stretching surface

In this numerical study, the unsteady laminar incompressible boundary-layer flow over a continuously stretching surface has been investigated when the velocity of the stretching surface varies arbitrarily with time. Both the nodal and the saddle point regions of flow have been considered for the analysis. Also, constant wall temperature/concentration and constant heat/mass flux at the stretching surface have been taken into account. The quasilinearisation method with an implicit finite-difference scheme is used in the nodal point region (0 less-than-or-equal-to c less-than-or-equal-to 1) where c denotes the stretching ratio. This method fails in the saddle point region (-1 less-than-or-equal-to c less-than-or-equal-to 0) due to the occurrence of reverse flow in the y-component of velocity. In order to overcome this difficulty, the method of parametric differentiation with an implicit finite-difference scheme is used, where the values at c = 0 are taken as starting values. Results have been obtained for the stretching velocities which are accelerating and decelerating with time. Results show that the skin friction, the heat transfer and the mass transfer parameters respond significantly to the time dependent stretching velocities. Suction (A > 0) is found to be an important parameter in obtaining convergent solution in the case of the saddle point region of flow. The Prandtl number and the Schmidt number strongly affect the heat and mass transfer of the diffusing species, respectively.

Measurement of complex optical flow with use of an augmented generalized gradient scheme

We present a method for measuring the local velocities and first-order variations in velocities in a timevarying image. The scheme is an extension of the generalized gradient model that encompasses the local variation of velocity within a local patch of the image. Motion within a patch is analyzed in parallel by 42 different spatiotemporal filters derived from 6 linearly independent spatiotemporal kernels. No constraints are imposed on the image structure, and there is no need for smoothness constraints on the velocity field. The aperture problem does not arise so long as there is some two-dimensional structure in the patch being analyzed. Among the advantages of the scheme is that there is no requirement to calculate second or higher derivatives of the image function. This makes the scheme robust in the presence of noise. The spatiotemporal kernels are of simple form, involving Gaussian functions, and are biologically plausible receptive fields. The validity of the scheme is demonstrated by application to both synthetic and real video images sequences and by direct comparison with another recently published scheme Biol. Cybern. 63, 185 (1990)] for the measurement of complex optical flow.

Measurement of complex optical flow with use of an augmented generalized gradient scheme

We present a method for measuring the local velocities and first-order variations in velocities in a time-varying image. The scheme is an extension of the generalized gradient model that encompasses the local variation of velocity within a local patch of the image. Motion within a patch is analyzed in parallel by 42 different spatiotemporal filters derived from 6 linearly independent spatiotemporal kernels. No constraints are imposed on the image structure, and there is no need for smoothness constraints on the velocity field. The aperture problem does not arise so long as there is some two-dimensional structure in the patch being analyzed. Among the advantages of the scheme is that there is no requirement to calculate second or higher derivatives of the image function. This makes the scheme robust in the presence of noise. The spatiotemporal kernels are of simple form, involving Gaussian functions, and are biologically plausible receptive fields. The validity of the scheme is demonstrated by application to both synthetic and real video images sequences and by direct comparison with another recently published scheme [Biol. Cybern. 63, 185 (1990)] for the measurement of complex optical flow.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Kinematics of the Cometary Globules in the Gum Nebula

We report a study of the kinematics of the cometary globules in the Gum Nebula using the J = 1 yields 0 transition line of (CO-12)O. A morphological center for the system with which 60 percent of the globules are associated is identified. It is shown that the observed radial velocities of the heads of the globules are consistent with an expansion of the system. Systematic velocity gradients are present along some of the tails. The estimated expansion age and the tail stretching age are both about a few million years, suggesting a common origin for the expansion and the formation of the tails. The presence of young stars of similar ages in some of the globules points to star formation triggered by the same cause. Possible scenarios are briefly discussed.

Translational velocity oscillations of piston generated vortex rings

Experimental results are presented that show that the translational velocities of piston generated vortex rings often undergo oscillations, similar to those recently discovered for drop generated rings. An attempt has been made to minimize uncertainties by utilizing both dye and hydrogen bubbles for visualization and carefully repeating measurements on the same ring and on different realizations under the same nominal piston conditions. The results unambiguously show that under most conditions, both for laminar and turbulent rings and for rings generated from pipes and orifices, the oscillations are present. The present results, together with the earlier results on drop generated rings, give support to the view that translational velocity oscillations are probably an inherent feature of translating vortex ring fields. (C) 1995 American Institute of Physics.

Vibration and Stability of Fluid Conveying Pipes with Stochastic Parameters

Flexible cantilever pipes conveying fluids with high velocity are analysed for their dynamic response and stability behaviour. The Young's modulus and mass per unit length of the pipe material have a stochastic distribution. The stochastic fields, that model the fluctuations of Young's modulus and mass density are characterized through their respective means, variances and autocorrelation functions or their equivalent power spectral density functions. The stochastic non self-adjoint partial differential equation is solved for the moments of characteristic values, by treating the point fluctuations to be stochastic perturbations. The second-order statistics of vibration frequencies and mode shapes are obtained. The critical flow velocity is-first evaluated using the averaged eigenvalue equation. Through the eigenvalue equation, the statistics of vibration frequencies are transformed to yield critical flow velocity statistics. Expressions for the bounds of eigenvalues are obtained, which in turn yield the corresponding bounds for critical flow velocities.

Lopsidedness in WHISP galaxies I. Rotation curves and kinematic lopsidedness

The frequently observed lopsidedness of the distribution of stars and gas in disc galaxies is still considered as a major problem in galaxy dynamics. It is even discussed as an imprint of the formation history of discs and the evolution of baryons in dark matter haloes. Here, we analyse a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. The Hi data allow us to follow the morphology and the kinematics out to very large radii. In the present paper, we present the rotation curves and study the kinematic asymmetry. We extract the rotation curves of the receding and approaching sides separately and show that the kinematic behaviour of disc galaxies can be classified into five different types: symmetric velocity fields where the rotation curves of the receding and approaching sides are almost identical; global distortions where the rotation velocities of the receding and approaching sides have an offset that is constant with radius; local distortions leading to large deviations in the inner and negligible deviations in the outer parts (and vice versa); and distortions that divide the galaxies into two kinematic systems that are visible in terms of the different behaviour of the rotation curves of the receding and approaching sides, which leads to a crossing and a change in side. The kinematic lopsidedness is measured from the maximum rotation velocities, averaged over the plateau of the rotation curves. This gives a good estimate of the global lopsidedness in the outer parts of the sample galaxies. We find that the mean value of the perturbation parameter denoting the lopsided potential as obtained from the kinematic data is 0.056. Altogether, 36% of the sample galaxies are globally lopsided, which can be interpreted as the disc responding to a halo that was distorted by a tidal encounter. In Paper II, we study the morphological lopsidedness of the same sample of galaxies.

Velocity distribution function for a dilute granular material in shear flow

The velocity distribution function for the steady shear flow of disks (in two dimensions) and spheres (in three dimensions) in a channel is determined in the limit where the frequency of particle-wall collisions is large compared to particle-particle collisions. An asymptotic analysis is used in the small parameter epsilon, which is naL in two dimensions and na(2)L in three dimensions, where; n is the number density of particles (per unit area in two dimensions and per unit volume in three dimensions), L is the separation of the walls of the channel and a is the particle diameter. The particle-wall collisions are inelastic, and are described by simple relations which involve coefficients of restitution e(t) and e(n) in the tangential and normal directions, and both elastic and inelastic binary collisions between particles are considered. In the absence of binary collisions between particles, it is found that the particle velocities converge to two constant values (u(x), u(y)) = (+/-V, O) after repeated collisions with the wall, where u(x) and u(y) are the velocities tangential and normal to the wall, V = (1 - e(t))V-w/(1 + e(t)), and V-w and -V-w, are the tangential velocities of the walls of the channel. The effect of binary collisions is included using a self-consistent calculation, and the distribution function is determined using the condition that the net collisional flux of particles at any point in velocity space is zero at steady state. Certain approximations are made regarding the velocities of particles undergoing binary collisions :in order to obtain analytical results for the distribution function, and these approximations are justified analytically by showing that the error incurred decreases proportional to epsilon(1/2) in the limit epsilon --> 0. A numerical calculation of the mean square of the difference between the exact flux and the approximate flux confirms that the error decreases proportional to epsilon(1/2) in the limit epsilon --> 0. The moments of the velocity distribution function are evaluated, and it is found that [u(x)(2)] --> V-2, [u(y)(2)] similar to V-2 epsilon and -[u(x)u(y)] similar to V-2 epsilon log(epsilon(-1)) in the limit epsilon --> 0. It is found that the distribution function and the scaling laws for the velocity moments are similar for both two- and three-dimensional systems.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.