A hybrid stabilization technique for simulating water wave - Structure interaction by incompressible Smoothed Particle Hydrodynamics (ISPH) method

The Smoothed Particle Hydrodynamics (SPH) method is emerging as a potential tool for studying water wave related problems, especially for violent free surface flow and large deformation problems. The incompressible SPH (ISPH) computations have been found not to be able to maintain the stability in certain situations and there exist some spurious oscillations in the pressure time history, which is similar to the weakly compressible SPH (WCSPH). One main cause of this problem is related to the non-uniform and clustered distribution of the moving particles. In order to improve the model performance, the paper proposed an efficient hybrid numerical technique aiming to correct the ill particle distributions. The correction approach is realized through the combination of particle shifting and pressure gradient improvement. The advantages of the proposed hybrid technique in improving ISPH calculations are demonstrated through several applications that include solitary wave impact on a slope or overtopping a seawall, and regular wave slamming on the subface of open-piled structure

Mesoscale Calculations of the Dynamic Behavior of a Granular Ceramic

Mesoscale calculations have been conducted in order to gain further insight into the dynamic compaction characteristics of granular ceramics. The primary goals of this work are to numerically determine the shock response of granular tungsten carbide and to assess the feasibility of using these results to construct the bulk material Hugoniot. Secondary goals include describing the averaged compaction wave behavior as well as characterizing wave front behavior such as the strain rate versus stress relationship and statistically describing the laterally induced velocity distribution. The mesoscale calculations were able to accurately reproduce the experimentally determined Hugoniot slope but under predicted the zero pressure shock speed by 12%. The averaged compaction wave demonstrated an initial transient stress followed by asymptotic behavior as a function of grain bed distance. The wave front dynamics demonstrate non-Gaussian compaction dynamics in the lateral velocity distribution and a power-law strain rate–stress relationship

An integrated study of parallel valveless micropumps

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.We describe an analytical, computational and experimental study of parallel valveless micropumps. A one dimensional model of a parallel micropump is presented and compared with available experimental data. The model confirms the linear decrease of the volume flux with pressure rise which is consistent with the experiments. The computational study showed a similar linear decrease but highlighted the effect of turbulence closures on the rectified mean flow, with the experimental data sitting between the turbulent and laminar closure regimes. The experimental study confirmed the importance of the displacement distance of fluid through the nozzle compared to nozzle length in the setting whether the flow regime is streaming or rectified. General conclusions are made about how to improve the pumping efficiency of micropumps.This study is supported by the Dorothy Hodgkin Postgraduate Award (DHPA) of the Engineering and Physical Sciences Research Council (EPSRC) of United Kingdom and Ebara Research Co. Ltd of Japan

Abundance of unknots in various models of polymer loops

A veritable zoo of different knots is seen in the ensemble of looped polymer chains, whether created computationally or observed in vitro. At short loop lengths, the spectrum of knots is dominated by the trivial knot (unknot). The fractional abundance of this topological state in the ensemble of all conformations of the loop of $N$ segments follows a decaying exponential form, $\sim \exp (-N/N_0)$, where $N_0$ marks the crossover from a mostly unknotted (ie topologically simple) to a mostly knotted (ie topologically complex) ensemble. In the present work we use computational simulation to look closer into the variation of $N_0$ for a variety of polymer models. Among models examined, $N_0$ is smallest (about 240) for the model with all segments of the same length, it is somewhat larger (305) for Gaussian distributed segments, and can be very large (up to many thousands) when the segment length distribution has a fat power law tail.Comment: 13 pages, 6 color figure

Virtual Delamination Testing through Non-Linear Multi-Scale Computational Methods: Some Recent Progress

This paper deals with the parallel simulation of delamination problems at the meso-scale by means of multi-scale methods, the aim being the Virtual Delamination Testing of Composite parts. In the non-linear context, Domain Decomposition Methods are mainly used as a solver for the tangent problem to be solved at each iteration of a Newton-Raphson algorithm. In case of strongly nonlinear and heterogeneous problems, this procedure may lead to severe difficulties. The paper focuses on methods to circumvent these problems, which can now be expressed using a relatively general framework, even though the different ingredients of the strategy have emerged separately. We rely here on the micro-macro framework proposed in (Ladev\`eze, Loiseau, and Dureisseix, 2001). The method proposed in this paper introduces three additional features: (i) the adaptation of the macro-basis to situations where classical homogenization does not provide a good preconditioner, (ii) the use of non-linear relocalization to decrease the number of global problems to be solved in the case of unevenly distributed non-linearities, (iii) the adaptation of the approximation of the local Schur complement which governs the convergence of the proposed iterative technique. Computations of delamination and delamination-buckling interaction with contact on potentially large delaminated areas are used to illustrate those aspects

Arbitrary-Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes

In this article we present a new class of high order accurate Arbitrary-Eulerian-Lagrangian (ALE) one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes. A WENO reconstruction algorithm is used to achieve high order accuracy in space and a high order one-step time discretization is achieved by using the local space-time Galerkin predictor. For that purpose, a new element--local weak formulation of the governing PDE is adopted on moving space--time elements. The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes. Moreover, a polynomial mapping defined by the same local space-time basis functions as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element. To maintain algorithmic simplicity, the final ALE one-step finite volume scheme uses moving triangular meshes with straight edges. This is possible in the ALE framework, which allows a local mesh velocity that is different from the local fluid velocity. We present numerical convergence rates for the schemes presented in this paper up to sixth order of accuracy in space and time and show some classical numerical test problems for the two-dimensional Euler equations of compressible gas dynamics.Comment: Accepted by "Communications in Computational Physics