Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Beyond classical dynamic structural plasticity using mesh-free modelling techniques

The problem of modelling the transient response of an elastic-perfectly-plastic cantilever beam, carrying an impulsively loaded tip mass, is,often referred to as the Parkes cantilever problem 25]; The permanent deformation of a cantilever struck transversely at its tip, Proc. R. Soc. A., 288, pp. 462). This paradigm for classical modelling of projectile impact on structures is re-visited and updated using the mesh-free method, smoothed particle hydrodynamics (SPH). The purpose of this study is to investigate further the behaviour of cantilever beams subjected to projectile impact at its tip, by considering especially physically real effects such as plastic shearing close to the projectile, shear deformation, and the variation of the shear strain along the length and across the thickness of the beam. Finally, going beyond macroscopic structural plasticity, a strategy to incorporate physical discontinuity (due to crack formation) in SPH discretization is discussed and explored in the context of tip-severance of the cantilever beam. Consequently, the proposed scheme illustrates the potency for a more refined treatment of penetration mechanics, paramount in the exploration of structural response under ballistic loading. The objective is to contribute to formulating a computational modelling framework within which transient dynamic plasticity and even penetration/failure phenomena for a range of materials, structures and impact conditions can be explored ab initio, this being essential for arriving at suitable tools for the design of armour systems. (C) 2014 Elsevier Ltd. All rights reserved.

Simulations using meshfree methods

In this paper, attempt is made to solve a few problems using the Polynomial Point Collocation Method (PPCM), the Radial Point Collocation Method (RPCM), Smoothed Particle Hydrodynamics (SPH), and the Finite Point Method (FPM). A few observations on the accuracy of these methods are recorded. All the simulations in this paper are three dimensional linear elastostatic simulations, without accounting for body forces.

DOES THE CHEMOTHERMAL INSTABILITY HAVE ANY ROLE IN THE FRAGMENTATION OF PRIMORDIAL GAS

The collapse of the primordial gas in the density regime similar to 10(8)-10(10) cm(-3) is controlled by the three-body H-2 formation process, in which the gas can cool faster than free-fall time-a condition proposed as the chemothermal instability. We investigate how the heating and cooling rates are affected during the rapid transformation of atomic to molecular hydrogen. With a detailed study of the heating and cooling balance in a 3D simulation of Pop III collapse, we follow the chemical and thermal evolution of the primordial gas in two dark matter minihalos. The inclusion of sink particles in modified Gadget-2 smoothed particle hydrodynamics code allows us to investigate the long-term evolution of the disk that fragments into several clumps. We find that the sum of all the cooling rates is less than the total heating rate after including the contribution from the compressional heating (pdV). The increasing cooling rate during the rapid increase of the molecular fraction is offset by the unavoidable heating due to gas contraction. We conclude that fragmentation occurs because H-2 cooling, the heating due to H-2 formation and compressional heating together set a density and temperature structure in the disk that favors fragmentation, not the chemothermal instability.

An improved free surface modeling for incompressible SPH

In Incompressible Smooth Particle Hydrodynamics (ISPH), a pressure Poisson equation (PPE) is solved to obtain a divergence free velocity field. When free surfaces are simulated using this method a Dirichlet boundary condition for pressure at the free surface has to be applied. In existing ISPH methods this is achieved by identifying free surface particles using heuristically chosen threshold of a parameter such as kernel sum, density or divergence of the position, and explicitly setting their pressure values. This often leads to clumping of particles near the free surface and spraying off of surface particles during splashes. Moreover, surface pressure gradients in flows where surface tension is important are not captured well using this approach. We propose a more accurate semi-analytical approach to impose Dirichlet boundary conditions on the free surface. We show the efficacy of the proposed algorithm by using test cases of elongation of a droplet and dam break. We perform two dimensional simulations of water entry and validate the proposed algorithm with experimental results. Further, a three dimensional simulation of droplet splash is shown to compare well with the Volume-of-Fluid simulations. (C) 2014 Elsevier Ltd. All rights reserved.

Formation of Prestellar Cores via Non-Isothermal Gas Fragmentation

Sheet-like clouds are common in turbulent gas and perhaps form via collisions between turbulent gas flows. Having examined the evolution of an isothermal shocked slab in an earlier contribution, in this work we follow the evolution of a sheet-like cloud confined by (thermal) pressure and gas in it is allowed to cool. The extant purpose of this endeavour is to study the early phases of core-formation. The observed evolution of this cloud supports the conjecture that molecular clouds themselves are three-phase media (comprising viz. a stable cold and warm medium, and a third thermally unstable medium), though it appears, clouds may evolve in this manner irrespective of whether they are gravitationally bound. We report, this sheet fragments initially due to the growth of the thermal instability (TI) and some fragments are elongated, filament-like. Subsequently, relatively large fragments become gravitationally unstable and sub-fragment into smaller cores. The formation of cores appears to be a three stage process: first, growth of the TI leads to rapid fragmentation of the slab; second, relatively small fragments acquire mass via gas-accretion and/or merger and third, sufficiently massive fragments become susceptible to the gravitational instability and sub-fragment to form smaller cores. We investigate typical properties of clumps (and smaller cores) resulting from this fragmentation process. Findings of this work support the suggestion that the weak velocity field usually observed in dense clumps and smaller cores is likely seeded by the growth of dynamic instabilities. Simulations were performed using the smooth particle hydrodynamics algorithm.

Rheology of particle-loaded semi-dilute polymer solutions

We present measurements of the rheology of suspensions of rigid spheres in a semi-dilute polymer solution from experiments of steady and oscillatory shear. For a given value of the shear rate gamma, addition of particles enhances the viscosity and the first normal stress difference but decreases the magnitude of the second normal stress difference. The viscosity eta exhibits a power law variation in gamma for a range of gamma that grows with phi. The first normal stress N-1 is positive and its value grows with phi; it exhibits a clear power law variation for the entire range of gamma that was studied. The second normal stress difference N-2 is negative for the pure polymer solution and much smaller in magnitude than N-1; on addition of particles, its magnitude further decreases, and it appears to change sign at large phi. The behavior of N-1 and N-2 is at odds with the findings of recent studies on particle-loaded dilute polymer solutions and polymer melts. The small-amplitude oscillatory shear experiments show the linear viscoelastic properties, G(') and G('), increasing with phi at a given value of the angular frequency omega. The dynamic viscosity of the suspension differs substantially from its steady shear viscosity, and the difference increases as gamma, omega -> 0.

Molecular theory for the effects of specific solute-solvent interaction on the diffusion of a solute particle in a molecular liquid

An understanding of the effect of specific solute-solvent interactions on the diffusion of a solute probe is a long standing problem of physical chemistry. In this paper a microscopic treatment of this effect is presented. The theory takes into account the modification of the solvent structure around the solute due to this specific interaction between them. It is found that for strong, attractive interaction, there is an enhanced coupling between the solute and the solvent dynamic modes (in particular, the density mode), which leads to a significant increase in the friction on the solute. The diffusion coefficient of the solute is found to depend strongly and nonlinearly on the magnitude of the attractive interaction. An interesting observation is that specific solute-solvent interaction can induce a crossover from a sliplike to a sticklike diffusion. In the limit of strong attractive interaction, we recover a dynamic version of the solvent-berg picture. On the other hand, for repulsive interaction, the diffusion coefficient of the solute increases. These results are in qualitative agreement with recent experimental observations.

Particle field imaging using digital in-line holography

Digital holography is the direct recording of holograms using a CCD camera and is an alternative to the use of a film or a plate. In this communication in-line digital holographic microscopy has been explored for its application in particle imaging in 3D. Holograms of particles of about 10 mu m size have been digitally reconstructed. Digital focusing was done to image the particles in different planes along the depth of focus. Digital holographic particle imaging results were compared with conventional optical microscope imaging. A methodology for dynamic analysis of microparticles in 3D using in-line digital holography has been proposed.

Soft-particle model analysis of effect of LPS on electrophoretic softness of Acidithiobacillus ferrooxidans grown in presence of different metal ions

The effect of Surface lipopolysaccharides (LPS) on the electrophoretic softness and fixed charge density in the ion-penetrable layer of Acidithiobacillus ferrooxidans cells grown in presence of copper or arsenic ions have been discussed, The electrophoretic mobility data were analyzed using the soft-particle electrophoresis theory. Cell surface potentials of all the strains based on soft-particle theory were lower than those estimated using the conventional Smoluchowski theory, Exposure to metal ions increased the Surface electrophoretic softness with decrease in the fixed charge density. Effect of cell surface lipopolysaccharides on the model parameters are investigated and discussed.

Relativistic particle models with separable interactions

We present relativistic, classical particle models that possess Poincaré invariance, invariant world lines, particle interaction, and separability.

Separability in relativistic Hamiltonian particle dynamics

The problem of separability in recent models of classical relativistic interacting particles is examined. This physical requirement is shown to be more subtle than naive separability of all the constraints defining the system: it is adequate to be able to canonically transform the time-fixing constraints from an unseparated to a separated form when clusters emerge. Viewing separability in this way, and within a specific framework, we are led to a new no-interaction theorem which states the incompatibility of nontrivial interaction with relativistic invariance, separability, and invariant world lines for more than two particles.

Fluidity of mica particle dispersed aluminium alloy

Cast aluminium alloy-mica particle composites were made by dispersing mica particles in a vortex produced by stirring the liquid Al-4 wt% Cu-1.5 wt% Mg alloy and then casting the melt containing the suspended particles into permanent moulds. Spiral fluidity and casting fluidity of the alloy containing mica particles in suspension were determined. Both the spiral fluidity and the casting fluidity of the base alloy were found to decrease with an increase in volume or weight percent of mica particles (of a given size), and with a decrease in particle size (for a given amount of particles). The fluidities of Al-4 wt% Cu-1.5 wt% Mg alloys containing suspended mica particles were found to correlate very well with the surface area of suspended mica particles. The regression equation for spiral fluidity Y (cm) as a function of surface area of mica particles per gram of spiral X (cm2 g–1) at 700° C was found to be Y=42.62–0.42X with a correlation coefficient of 0.9634. The regression equations for casting fluidity Yprime (cm) as a functiono of surface area of mica particles per gram of fluidity test piece Xprime (cm2 g–1) at 710 and 670° C were found to be Yprime=19.71–0.17Xprime and Yprime=13.52–0.105Xprime with correlation coefficients of 0.9194 and 0.9612 respectively. The percentage decrease in casting fluidity of composite melts containing up to 2.5 wt% mica with a drop in temperature is quite similar to the corresponding decrease in the casting fluidity of base alloy melts (without mica). The change in fluidity due to mica dispersions has been discussed in terms of changes in viscosity of the composite melts. However, the fluidities of these composite alloys containing up to 2.5 wt% mica are adequate for making a variety of simple castings including bearings for which these alloys have been developed.

"Small-particle limit" the second harmonic generation from noble metal nanoparticles

In this paper, we have probed the origin of SHG in copper nanoparticles by polarization-resolved hyper-Rayleigh scattering (HRS). Results obtained with various sizes of copper nanoparticles at four different wavelengths covering the wavelength range 738-1907 nm reveal that the origin of second harmonic generation (SHG) in these particles is purely dipolar in nature as long as the size (d) of the particles remains smaller compared to the wavelength (;.) of light ("small-particle limit"). However, contribution of the higher order multipoles coupled with retardation effect becomes apparent with an increase in the d/lambda ratio. We have identified the "small-particle limit" in the second harmonic generation from noble metal nanoparticles by evaluating the critical d/lambda ratio at which the retardation effect sets in the noble metal nanoparticles. We have found that the second-order nonlinear optical property of copper nanoparticles closely resembles that of gold, but not that of silver. (C) 2009 Elsevier B.V. All rights reserved.

Numerical simulation of conductive particle behaviour at the surface of a plate electrode affected by a DC corona field

The electric field in certain electrostatic devices can be modeled by a grounded plate electrode affected by a corona discharge generated by a series of parallel wires connected to a DC high-voltage supply. The system of differential equations that describe the behaviour (i.e., charging and motion) of the conductive particle in such an electric field has been numerically solved, using several simplifying assumptions. Thus, it was possible to investigate the effect of various electrical and mechanical factors on the trajectories of conductive particles. This model has been employed to study the behaviour of coalparticles in fly-ash corona separators.