Jamming for a 2D granular material

This paper focuses on the nature of jamming, as seen in two-dimensional frictional granular systems consisting of photoelastic particles. The photoelastic technique is unique at this time, in its capability to provide detailed particle-scale information on forces and kinematic quantities such as particle displacements and rotations. These experiments first explore isotropic stress states near point J through measurements of the mean contact number per particle, Z, and the pressure, P as functions of the packing fraction, . In this case, the experiments show some but not all aspects of jamming, as expected on the basis of simulations and models that typically assume conservative, hence frictionless, forces between particles. Specifically, there is a rapid growth in Z, at a reasonable which we identify with as c. It is possible to fit Z and P, to power law expressions in - c above c, and to obtain exponents that are in agreement with simulations and models. However, the experiments differ from theory on several points, as typified by the rounding that is observed in Z and P near c. The application of shear to these same 2D granular systems leads to phenomena that are qualitatively different from the standard picture of jamming. In particular, there is a range of packing fractions below c, where the application of shear strain at constant leads to jammed stress-anisotropic states, i.e. they have a non-zero shear stress, τ. The application of shear strain to an initially isotropically compressed (hence jammed) state, does not lead to an unjammed state per se. Rather, shear strain at constant first leads to an increase of both τ and P. Additional strain leads to a succession of jammed states interspersed with relatively localized failures of the force network leading to other stress-anisotropic states that are jammed at typically somewhat lower stress. The locus of jammed states requires a state space that involves not only and τ, but also P. P, τ, and Z are all hysteretic functions of shear strain for fixed . However, we find that both P and τ are roughly linear functions of Z for strains large enough to jam the system. This implies that these shear-jammed states satisfy a Coulomb like-relation, τ = μP. © 2010 The Royal Society of Chemistry.

Resonant and Secular Orbital Interations

In stable solar systems, planets remain in nearly elliptical orbits around their stars. Over longer timescales, however, their orbital shapes and sizes change due to mutual gravitational perturbations. Orbits of satellites around a planet vary for the same reason. Because of their interactions, the orbits of planets and satellites today are different from what they were earlier. In order to determine their original orbits, which are critical constraints on formation theories, it is crucial to understand how orbits evolve over the age of the Solar System. Depending on their timescale, we classify orbital interactions as either short-term (orbital resonances) or long-term (secular evolution). My work involves examples of both interaction types. Resonant history of the small Neptunian satellites In satellite systems, tidal migration brings satellite orbits in and out of resonances. During a resonance passage, satellite orbits change dramatically in a very short period of time. We investigate the resonant history of the six small Neptunian moons. In this unique system, the exotic orbit of the large captured Triton (with a circular, retrograde, and highly tilted orbit) influences the resonances among the small satellites very strongly. We derive an analytical framework which can be applied to Neptune's satellites and to similar systems. Our numerical simulations explain the current orbital tilts of the small satellites as well as constrain key physical parameters of both Neptune and its moons. Secular orbital interactions during eccentricity damping Long-term periodic changes of orbital shape and orientation occur when two or more planets orbit the same star. The variations of orbital elements are superpositions of the same number of fundamental modes as the number of planets in the system. We investigate how this effect interacts with other perturbations imposed by external disturbances, such as the tides and relativistic effects. Through analytical studies of a system consisting of two planets, we find that an external perturbation exerted on one planet affects the other indirectly. We formulate a general theory for how both orbits evolve in response to an arbitrary externally-imposed slow change in eccentricity.

FSBS resonances observed in a standard highly nonlinear fiber.

Forward stimulated Brillouin scattering (FSBS) is observed in a standard 2-km-long highly nonlinear fiber. The frequency of FSBS arising from multiple radially guided acoustic resonances is observed up to gigahertz frequencies. The tight confinement of the light and acoustic field enhances the interaction and results in a large gain coefficient of 34.7 W(-1) at a frequency of 933.8 MHz. We also find that the profile on the anti-Stokes side of the pump beam have lineshapes that are asymmetric, which we show is due to the interference between FSBS and the optical Kerr effect. The measured FSBS resonance linewidths are found to increase linearly with the acoustic frequency. Based on this scaling, we conclude that dominant contribution to the linewidth is from surface damping due to the fiber jacket and structural nonuniformities along the fiber.

Universal quantum viscosity in a unitary Fermi gas.

A Fermi gas of atoms with resonant interactions is predicted to obey universal hydrodynamics, in which the shear viscosity and other transport coefficients are universal functions of the density and temperature. At low temperatures, the viscosity has a universal quantum scale ħ n, where n is the density and ħ is Planck's constant h divided by 2π, whereas at high temperatures the natural scale is p(T)(3)/ħ(2), where p(T) is the thermal momentum. We used breathing mode damping to measure the shear viscosity at low temperature. At high temperature T, we used anisotropic expansion of the cloud to find the viscosity, which exhibits precise T(3/2) scaling. In both experiments, universal hydrodynamic equations including friction and heating were used to extract the viscosity. We estimate the ratio of the shear viscosity to the entropy density and compare it with that of a perfect fluid.

Evaluación del desempeño de tres cuchillas de corte y remoción de una sembradora para siembra directa sobre un suelo vertisol

El presente trabajo tuvo como objetivo evaluar el desempeno de tres cuchillas de corte y remoción en un suelo vertisol bajo siembra directa mediante la eficiencia de implantación y la velocidad de emergencia de un cultivo de soja. Las cuchillas utilizadas fueron: lisa, rippled y dura flute. La hipotesis de trabajo fue: la eficiencia de implantación y el dia medio a emergencia del cultivo de soja bajo siembra directa y en un suelo Vertisol son dependientes del tipo de cuchilla labrasurco que compone el tren de siembra de la sembradora. El ensayo se llevó a cabo en un establecimiento ubicado en el departamento de Gualeguaychu, sur de la Provincia de Entre Ríos, siendo sus coordenadas 32 grados 50 minutos 8.18 segundos latitud sur y 58 grados 40 minutos 05.69 segundos longitud oeste. Los principales resultados obtenidos señalan que el máximo logro, en cuanto al porcentaje de emergencia, se observa a los 13 días después de la siembra y luego hay una declinación en el stand de plantas hasta el dia 25 posterior a la siembra, siendo la cuchilla lisa la que mas se vio afectada. Las condiciones de alta humedad debido a las continuas precipitaciones sucedidas en el transcurso de la emergencia del cultivo afectaron el establecimiento de las plantulas de soja. Las conclusiones fueron: a) la cuchilla rippled produjo la mayor eficiencia de implantacion, b) la velocidad de emergencia de las plantas de soja fue superior en aquellas cuchillas que, ademas de cumplir con la función de cortar el rastrojo, generan una remoción de la línea de siembra, c) el día medio a emergencia resultó ser menor con la cuchilla dura flute, resultando ser mayor con la cuchilla lisa, d) donde se utilizó la cuchilla labrasurco lisa hubo una mayor muerte de plántulas por damping off.

An experimental and numerical CFD study of turbulence in a tundish container

This paper describes work performed at IRSID/USINOR in France and the University of Greenwich, UK, to investigate flow structures and turbulence in a water-model container, simulating aspects typical of metal tundish operation. Extensive mean and fluctuating velocity measurements were performed at IRSID using LDA to determine the flow field and these form the basis for a numerical model validation. This apparently simple problem poses several difficulties for the CFD modelling. The flow is driven by the strong impinging jet at the inlet. Accurate description of the jet is most important and requires a localized fine grid, but also a turbulence model that predicts the correct spreading rates of jet and impinging wall boundary layers. The velocities in the bulk of the tundish tend to be (indeed need to be) much smaller than those of the jet, leading to damping of turbulence, or even laminar flow. The authors have developed several low-Reynolds number (low-Re) k–var epsilon model variants to compute this flow and compare against measurements. Best agreement is obtained when turbulence damping is introduced to account not only for walls, but also for low-Re regions in the bulk – the k–var epsilon model otherwise allows turbulence to accumulate in the container due to the restricted outlet. Several damping functions are tested and the results reported here. The k–ω model, which is more suited to transitional flow, also seems to perform well in this problem.

Dynamic solid mechanics using finite volume methods

A procedure for evaluating the dynamic structural response of elastic solid domains is presented. A prerequisite for the analysis of dynamic fluid–structure interaction is the use of a consistent set of finite volume (FV) methods on a single unstructured mesh. This paper describes a three-dimensional (3D) FV, vertex-based method for dynamic solid mechanics. A novel Newmark predictor–corrector implicit scheme was developed to provide time accurate solutions and the scheme was evaluated on a 3D cantilever problem. By employing a small amount of viscous damping, very accurate predictions of the fundamental natural frequency were obtained with respect to both the amplitude and period of oscillation. This scheme has been implemented into the multi-physics modelling software framework, PHYSICA, for later application to full dynamic fluid structure interaction.

Modelling electromagnetically levitated liquid droplet oscillations

This work comprises accurate computational analysis of levitated liquid droplet oscillations in AC and DC magnetic fields. The AC magnetic field interacting with the induced electric current within the liquid metal droplet generates intense fluid flow and the coupled free surface oscillations. The pseudo-spectral technique is used to solve the turbulent fluid flow equations for the continuously dynamically transformed axisymmetric fluid volume. The volume electromagnetic force distribution is updated with the shape and position change. We start with the ideal fluid test case for undamped Rayleigh frequency oscillations in the absence of gravity, and then add the viscous and the DC magnetic field damping. The oscillation frequency spectra are further analysed for droplets levitated against gravity in AC and DC magnetic fields at various combinations. In the extreme case electrically poorly conducting, diamagnetic droplet (water) levitation dynamics are simulated. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.

Maximising heat transfer efficiency in the cold crucible induction melting process

Induction heating is an efficient method used to melt electrically conductive materials, particularly if melting takes place in a ceramic crucible. This form of melting is particularly good for alloys, as electromagnetic forces set up by the induction coil lead to vigorous stirring of the melt ensuring homogeneity and uniformity in temperature. However, for certain reactive alloys, or where high purity is required, ceramic crucibles cannot be used, but a water-cooled segmented copper crucible is employed instead. Water cooling prevents meltdown or distortion of the metal wall, but much of the energy goes into the coolant. To reduce this loss, the electromagnetic force generated by the coil is used to push the melt away from the walls and so minimise contact with water-cooled surfaces. Even then, heat is lost through the crucible base where contact is inevitable. In a collaborative programme between Greenwich and Birmingham Universities, computer modelling has been used in conjunction with experiments to improve the superheat attainable in the melt for a,number of alloys, especially for y-TiAl intermetallics to cast aeroengine turbine blades. The model solves the discretised form of the turbulent Navier-Stokes, thermal energy conservation and Maxwell equations using a Spectral Collocation technique. The time-varying melt envelope is followed explicitly during the computation using an adaptive mesh. This paper briefly describes the mathematical model used to represent the interaction between the magnetic field, fluid flow, heat transfer and change of phase in the crucible and identifies the proportions of energy used in the melt, lost in the crucible base and in the crucible walls. The role of turbulence is highlighted as important in controlling heat losses and turbulence damping is introduced as a means of improving superheat. Model validation is against experimental results and shows good agreement with measured temperatures and energy losses in the cooling fluid throughout the melting cycle.

Identifying the damage mechanisms in paper using the acoustic emission monitoring technique

This paper investigates the use of the acoustic emission (AE) monitoring technique for use in identifying the damage mechanisms present in paper associated with its production process. The microscopic structure of paper consists of a random mesh of paper fibres connected by hydrogen bonds. This implies the existence of two damage mechanisms, the failure of a fibre-fibre bond and the failure of a fibre. This paper describes a hybrid mathematical model which couples the mechanics of the mass-spring model to the acoustic wave propagation model for use in generating the acoustic signal emitted by complex structures of paper fibres under strain. The derivation of the mass-spring model can be found in [1,2], with details of the acoustic wave equation found in [3,4]. The numerical implementation of the vibro-acoustic model is discussed in detail with particular emphasis on the damping present in the numerical model. The hybrid model uses an implicit solver which intrinsically introduces artificial damping to the solution. The artificial damping is shown to affect the frequency response of the mass-spring model, therefore certain restrictions on the simulation time step must be enforced so that the model produces physically accurate results. The hybrid mathematical model is used to simulate small fibre networks to provide information on the acoustic response of each damage mechanism. The simulated AEs are then analysed using a continuous wavelet transform (CWT), described in [5], which provides a two dimensional time-frequency representation of the signal. The AEs from the two damage mechanisms show different characteristics in the CWT so that it is possible to define a fibre-fibre bond failure by the criteria listed below. The dominant frequency components of the AE must be at approximately 250 kHz or 750 kHz. The strongest frequency component may be at either approximately 250 kHz or 750 kHz. The duration of the frequency component at approximately 250 kHz is longer than that of the frequency component at approximately 750 kHz. Similarly, the criteria for identifying a fibre failure are given below. The dominant frequency component of the AE must be greater than 800 kHz. The duration of the dominant frequency component must be less than 5.00E-06 seconds. The dominant frequency component must be present at the front of the AE. Essentially, the failure of a fibre-fibre bond produces a low frequency wave and the failure of a fibre produces a high frequency pulse. Using this theoretical criteria, it is now possible to train an intelligent classifier such as the Self-Organising Map (SOM) [6] using the experimental data. First certain features must be extracted from the CWTs of the AEs for use in training the SOM. For this work, each CWT is divided into 200 windows of 5E-06s in duration covering a 100 kHz frequency range. The power ratio for each windows is then calculated and used as a feature. Having extracted the features from the AEs, the SOM can now be trained, but care is required so that the both damage mechanisms are adequately represented in the training set. This is an issue with paper as the failure of the fibre-fibre bonds is the prevalent damage mechanism. Once a suitable training set is found, the SOM can be trained and its performance analysed. For the SOM described in this work, there is a good chance that it will correctly classify the experimental AEs.

Droplet oscillations in high gradient magnetic field

In high intensity and high gradient magnetic fields the volumetric force on diamagnetic material, such as water, leads to conditions very similar to microgravity in a terrestrial laboratory. In principle, this opens the possibility to determine material properties of liquid samples without wall contact, even for electrically non-conducting materials. In contrast, AC field levitation is used for conductors, but then terrestrial conditions lead to turbulent flow driven by Lorentz forces. DC field damping of the flow is feasible and indeed practiced to allow property measurements. However, the AC/DC field combination acts preferentially on certain oscillation modes and leads to a shift in the droplet oscillation spectrum.What is the cause? A nonlinear spectral numerical model is presented, to address these problems

Turbulence model performance in continuous casting simulations

This paper concerns the development and validation (using an oil/water system) of a finite volume computer model of the continuous casting process for steel flat products. The emphasis is on hydrodynamic aspects and in particular the dynamic behaviour of the metal/slag interface. Instability and wave action encourage the entrainment of inclusions into the melt affecting product quality. To track the interface between oil and water a new implicit algorithm was developed, called the Counter Diffusion Method. To prevent excessive damping, a time-filtered version of the k-e model, was found necessary, with appropriate density stratification terms representing interface turbulence damping.

Droplet oscillations in high gradient static magnetic field

In high intensity and high gradient magnetic fields the volumetric force on diamagnetic material, such as water, leads to conditions very similar to microgravity in a terrestrial laboratory. In principle, this opens the possibility to determine material properties of liquid samples without wall contact, even for electrically non-conducting materials. In contrast, AC field levitation is used for conductors, but then terrestrial conditions lead to turbulent flow driven by Lorentz forces. DC field damping of the flow is feasible and indeed practiced to allow property measurements. However, the AC/DC field combination acts preferentially on certain oscillation modes and leads to a shift in the droplet oscillation spectrum.What is the cause? A nonlinear spectral numerical model is presented, to address these problems.

A Note On Tidal Asymmetry In The Severn Estuary

A one-dimensional, non-linear numerical model is used to investigate the tidally averaged frictional stress and set-up of water level due to tidal asymmetry in the Severn Estuary; these quantities depend on the overtides in the region. A linearized model of the overtides is applied to calculations of the M4 currents in order to delineate the mechanisms responsible for their generation. The relative importance of individual non-linear mechanisms to the tidally averaged stress and set-up is determined; these mechanisms are interactions between tidal flow and changes in depth or breadth over a cross-section, frictional interaction between the tidal flow and Stokes drift, interaction between the tidal fluctuations in water depth and frictional retardation and non-linear advection.

Self-lubricating silicone elastomer biomaterials

Silicone has a relatively high coefficient of friction and silicone medical devices therefore lack inherent lubricity, leading to pain on device insertion and potential tissue trauma. In this study, higher molecular weight tetra(alkoxy) silanes, particularly tetra(oleyloxy) silane, have been used as crosslinkers in the condensation cure of a hydroxy end-functionalised linear poly(dimethylsiloxane). The resulting elastomers displayed a persistent lubricous surface of oleyl alcohol, and coefficients of friction (static and dynamic) approaching zero. Chemical structures of the synthesised silanes and surface alcohol exudate were confirmed by nuclear magnetic resonance spectroscopy. Mechanical properties of the elastomers, which were chemically identical to conventionally cured systems, suggested that an 80/20 mixture of tetra(oleyloxy) silane and tetra(propoxysilane) gave the best compromise between desirable mechanical and frictional properties.