On the spray drying of uniform functional microparticles

Spray drying is a typical method to produce particles in dry powder forms at industrial scale. Most spray-dried products often show a wide range of particle properties even within the same batch. At Monash University, we utilise a microfluidic spray drying approach to generate uniform microparticles with tightly controlled characteristics and sizes in a scalable, almost waste-free process. The technique is useful to correlate the effects of formulation and spray drying conditions on the properties of spray-dried particles, and can be used to test new formulations for targeted applications such as encapsulation and release of active ingredients. The synthesis route can be applied to other self-assembling systems, including mesoporous, crystalline, and hierarchically structured microparticles. As spray drying is commonly used in commercial scales, the understanding of how functional particles are formed in relation to formulations and process conditions could assist in developing a cost effective, energy and material-efficient route to produce powders with better properties and ease of handling for more advanced applications such as selective adsorption and bio-separation.     

Demonstration of a coupled floating offshore wind turbine analysis with high-fidelity methods

This paper presents results of numerical computations for floating off-shore wind turbines using, as an example, a machine of 10-MW rated power. The aerodynamic loads on the rotor are computed using the Helicopter Multi-Block flow solver developed at the University of Liverpool. The method solves the Navier–Stokes equations in integral form using the arbitrary Lagrangian–Eulerian formulation for time-dependent domains with moving boundaries. Hydrodynamic loads on the support platform are computed using the Smoothed Particle Hydrodynamics method, which is mesh-free and represents the water and floating structures by a set of discrete elements, referred to as particles. The motion of the floating offshore wind turbine is computed using a Multi-Body Dynamic Model of rigid bodies and frictionless joints. Mooring cables are modelled as a set of springs and dampers. All solvers were validated separately before coupling, and the results are presented in this paper. The importance of coupling is assessed and the loosely coupled algorithm used is described in detail alongside the obtained results.     

Experimental benchmark of a free plunging wing with imposed flap oscillations

The validation of fluid–structure interaction solvers is difficult since there is a lack of experimental data. Therefore, in this work an aeroelastic experiment is presented. The focus is on the temporal coupling between fluid and structure dynamics. Issues in the spatial coupling are eliminated by using a rigid wing. The wing, with a harmonically actuated 0.2c trailing edge flap, has a degree of freedom in the plunge (vertical) direction. The wing has a chord of 0.5 m and is suspended with springs. The wing motion is constrained by a vertical rail system.For simplicity attached flow is desired and therefore the set angle of attack is α=0°. The Reynolds number is approximately Re=700 000 and the flap deflects over a range of about ±2°. The damped natural frequency of the structure expressed as a reduced frequency is about k=0.194 and measurements are performed for reduced flap frequencies ranging from k=0.1 to k=0.3. Displacements and time dependent aerodynamic forces are measured and for k=0.198 2-D PIV measurements are performed. The planar PIV measurements are used to intrinsically determine the unsteady loads using Noca׳s method.As expected the aeroelastic problem shows similarities with a viscously damped mass–damper–spring, meaning the maximum excursion of the wing is found near the system eigenfrequency. The lift is dominated by the flap motion and the effective angle of attack due to the motion introduces phase shifts of the lift signal with respect to the flap phase angle.The experiment has been set up and executed with the necessary precision, but small ambiguities are found in the lift and drag disqualifying the data for validation. Nevertheless the data set provides a clear insight into typical loads and motions and can be used for comparative studies. It can also be used to (re)design future experiments to improve the quality of the data to the desired level of accuracy for validation.     

Retrieval of aerosol size distribution using improved quantum-behaved particle swarm optimization on spectral extinction measurements

An improved quantum-behaved particle swarm optimization (IQPSO) algorithm is employed to determine aerosol size distribution (ASD). The direct problem is solved using the anomalous diffraction approximation and Lambert–Beer's Law. Compared with the standard particle swarm optimization algorithm, the stochastic particle size optimization algorithm and the original QPSO, our IQPSO has faster convergence speed and higher accuracy within a smaller number of generations. Optimization parameters for the IQPSO were also evaluated; we recommend using four measurement wavelengths and 50 particles. Size distributions of various aerosol types were estimated using the IQPSO under dependent and independent models. Finally, experimental ASDs at different locations in Harbin were recovered using the IQPSO. All our results confirm that the IQPSO algorithm is an effective and reliable technique for estimating ASD.     

Relating geologic units and mobility system kinematics contributing to Curiosity wheel damage at Gale Crater,Mars

Curiosity landed on plains to the north of Mount Sharp in August 2012. By June 2016 the rover had traversed 12.9 km to the southwest, encountering extensive strata that were deposited in a fluvial-deltaic-lacustrine system. Initial drives across sharp sandstone outcrops initiated an unacceptably high rate of punctures and cracks in the thin aluminum wheel skin structures. Initial damage was found to be related to the drive control mode of the six wheel drive actuators and the kinematics of the rocker-bogie suspension. Wheels leading a suspension pivot were forced onto sharp, immobile surfaces by the other wheels as they maintained their commanded angular velocities. Wheel damage mechanisms such as geometry-induced stress concentration cracking and low-cycle fatigue were then exacerbated. A geomorphic map was generated to assist in planning traverses that would minimize further wheel damage. A steady increase in punctures and cracks between landing and June 2016 was due in part because of drives across the sharp sandstone outcrops that could not be avoided. Wheel lifetime estimates show that with careful path planning the wheels will be operational for an additional ten kilometers or more, allowing the rover to reach key strata exposed on the slopes of Mount Sharp.     

Estimating the tensile strength of super hard brittle materials using truncated spheroidal specimens

New approaches need to be introduced to measure the tensile capacity of super hard materials since the standard methods are not effective. To pursue this objective, a series of laboratory tests were constructed to replicate the fracture mechanism of diamond-based materials. Experiments indicate that under a certain compressive test condition, stresses normal to the axisymmetric line in truncated spheroidal specimens (bullet-shaped specimens) are in tension contributing to the tensile fracture of the material. From experimental and numerical studies, it is concluded that semi-prolate spheroidal specimens can be used to determine precisely the tensile strength of brittle stiff diamond-like composites.     

Stress induced phase transitions in silicon

Silicon has a tremendous importance as an electronic, structural and optical material. Modeling the interaction of a silicon surface with a pointed asperity at room temperature is a major step towards the understanding of various phenomena related to brittle as well as ductile regime machining of this semiconductor. If subjected to pressure or contact loading, silicon undergoes a series of stress-driven phase transitions accompanied by large volume changes. In order to understand the material's response for complex non-hydrostatic loading situations, dedicated constitutive models are required. While a significant body of literature exists for the dislocation dominated high-temperature deformation regime, the constitutive laws used for the technologically relevant rapid low-temperature loading have severe limitations, as they do not account for the relevant phase transitions. We developed a novel finite deformation constitutive model set within the framework of thermodynamics with internal variables that captures the stress induced semiconductor-to-metal ($\mathrm{cd-Si}\to \mathrm{\beta }\mathrm{-Si}$), metal-to-amorphous ($\mathrm{\beta }\mathrm{-Si}\to \mathrm{a-Si}$) as well as amorphous-to-amorphous ($\mathrm{a-Si}\to \mathrm{hda-Si}$, $\mathrm{hda-Si}\to \mathrm{a-Si}$) transitions. The model parameters were identified in part directly from diamond anvil cell data and in part from instrumented indentation by the solution of an inverse problem. The constitutive model was verified by successfully predicting the transformation stress under uniaxial compression and load–displacement curves for different indenters for single loading–unloading cycles as well as repeated indentation. To the authors' knowledge this is the first constitutive model that is able to adequately describe cyclic indentation in silicon.     

Effect of adhesive thickness,adhesive type and scarf angle on the mechanical properties of scarf adhesive joints

The effects of adhesive thickness, adhesive type and scarf angle, which are determined as the main control parameters by the dimensional analysis, on the mechanical properties of a scarf adhesive joint (SJ) subjected to uniaxial tensile loading are examined using a mixed-mode cohesive zone model (CZM) with a bilinear shape to govern the interface separation. Particularly, the adhesive-dependence of the vital cohesive parameters of CZM, which mainly include initial stiffness, total fracture energy and separation strength, is introduced emphatically. The numerical results demonstrate that the ultimate tensile loading increases as the adhesive thickness decreases. Cross the ultimate tension, the joint loses the load-bearing capacity when adopting the brittle adhesive but sustains partial load-bearing capacity while selecting the ductile adhesive. In addition, for the joint with the ductile adhesive, the maximum applied displacement until the complete failure of it is directly proportional to the adhesive thickness, which is different from the case using the brittle adhesive. Taking the combination of the ultimate loading and applied displacement into account, failure energy is employed to evaluate the joint performances. The results show that the failure energy of the joint with the brittle adhesive increases as the adhesive thickness decreases. Conversely, the situation of the joint using the ductile adhesive is vice versa. Moreover, the effect of the adhesive thickness becomes more noticeable with decreasing the scarf angle owing to the variation of the proportion of each component of the mixed-mode. Furthermore, all the characteristic parameters (the ultimate tensile loading, the maximum applied displacement and the failure energy) that adopted to describe the performances of SJ increase as the scarf angle decreases. Finally, the numerical method employed in this study is validated by comparing with existing experimental results.     

Modeling dynamic brittle behavior of materials with circular flaws or pores

Compressive failure of brittle materials is driven primarily by crack growth from pre-existing flaws in the material. These flaws, such as grain boundaries, pores, preexisting cracks, inclusions and missing grains, are randomly spaced and have a range of possible shapes and sizes. The current work proposes a micromechanics-based model for compressive dynamic failure of brittle materials with circular pore flaws, which incorporates both the number density and the size distribution of flaws. Results show that the distribution of flaw sizes is very important, particularly at moderate strain rate, since analyses based solely on the mean flaw size overpredict strength. Therefore, in order to increase dynamic strength at low to moderate strain rates, it is most effective to control the presence of large flaws. At very high strain rates, however, crack growth is activated even in small flaws and therefore controlling the total number density rather than the size of the flaws is effective for increasing dynamic strength. Finally, the model shows that neglecting very small flaws in the pore population may not have significant effects on the results in many cases, suggesting that the model is a useful tool for identifying a minimum resolution required for experimental characterization of microstructure.     

NUMERICAL SIMULATION OF WIND-BLOWN SAND MOVEMENT BASED ON SPH

运用光滑粒子流体动力学(smoothed particles hydrodynamics, SPH)方法对沙粒和气流的相互耦合运动特性进行了分析,研究提出了风沙流的SPH数值方法并进行了数值模拟. 首先提出了风沙流的SPH建模方法和基本理论,建立了风沙流动的SPH数值模拟平台. 其次通过建立风沙流的SPH 模型并施加边界条件,对自然风作用下沙粒的运动情况进行了数值模拟,详细分析了沙粒运动轨迹及特性,最后通过与相关研究成果对比分析,验证了完善后的SPH方法有效性. 通过考虑气流场的可变性,在风沙流SPH计算模型中引入了加载(起风)和卸载(停风)方式,观察并对比分析了沙粒的运动轨迹和特性. 为进一步研究风沙流的实时动态非线性行为提供了SPH理论基础和数值分析方法.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

固体介质中SPH方法的拉伸不稳定性问题研究进展

光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种基于核估计的无网格Lagrange数值方法.它用粒子方程离散流体动力学的连续方程, 既可以处理有限元难于处理的大变形和严重扭曲问题, 又可以处理有限差分法不易处理的自由边界和材料界面的问题, 在固体力学中的冲击、爆炸和裂纹模拟中具有广阔的发展前景.但是, 该算法的拉伸不稳定性(tensile instability)问题是它在固体力学领域中应用的最大障碍.对SPH稳定性分析表明, 算法不稳定性的条件仅与应力状态和核函数的2阶导数有关.目前, 应力点法(stress points)、Lagrange核函数法、人工应力法(artificialstress)、修正光滑粒子法(corrective smoothed particle method, CSPM)和守恒光滑法(conservativesmoothing)以及其他一些方法成功地改善了SPH的拉伸不稳定性, 但是每一种方法都不能彻底解决SPH的拉伸不稳定性问题.本文介绍了SPH法的方程和Von Neumann稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.      

含液颗粒材料液固耦合分析的离散颗粒模型及特征线SPH法

对含液颗粒材料流固耦合分析建议了一个基于离散颗粒模型与特征线SPH法的显式拉格朗日-欧拉无网格方案。在已有的用以模拟固体颗粒集合体的离散颗粒模型[1]基础上,将颗粒间间隙内的流体模型化为连续介质,对其提出并推导了基于特征线的SPH法。数值例题显示了所建议方案在模拟颗粒材料与间隙流相互作用的能力和性能以及间隙流体对颗粒结构承载能力及变形的影响。     

SPH-based simulation of granular collapse on an inclined bed

The paper presents a numerical model for simulating a granular flow and its deposition on an inclined bed. A granular material is described as an elastic–plastic continuum and its constitutive law, namely Hooke's law, is discretized on the basis of the Smoothed Particle Hydrodynamics (SPH) method. In the equation of motion, however, the artificial viscosity, which is widely used in SPH, is not applied. The diffusive term derived from Hooke's law is introduced with a diffusion coefficient that varies depending on the stress and strain rate based on the Drucker–Prager yield function. The model is verified and validated through two numerical tests. It is shown that the basic elastic–perfectly plastic characteristics are reproduced with a simple shearing test. The effects of the diffusion coefficient and spatial resolution are investigated to show the validity of the model. In the simulation of the gravitational collapse of a granular column on an inclined bed, the performance of the model from the final deposition profile, the time history of the front position of the granular flow, the maximum runout distance, and the velocity profile are investigated for several cases of basal inclinations. The calculated results show good agreement with the experimental results.     

Simulating bubble dynamics in a buoyant system

Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two-dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three-dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.     

An improved KGF‐SPH with a novel discrete scheme of Laplacian operator for viscous incompressible fluid flows

The kernel gradient free (KGF) smoothed particle hydrodynamics (SPH) method is a modified finite particle method (FPM) which has higher order accuracy than the conventional SPH method. In KGF‐SPH, no kernel gradient is required in the whole computation, and this leads to good flexibility in the selection of smoothing functions and it is also associated with a symmetric corrective matrix. When modeling viscous incompressible flows with SPH, FPM or KGF‐SPH, it is usual to approximate the Laplacian term with nested approximation on velocity, and this may introduce numerical errors from the nested approximation, and also cause difficulties in dealing with boundary conditions. In this paper, an improved KGF‐SPH method is presented for modeling viscous, incompressible fluid flows with a novel discrete scheme of Laplacian operator. The improved KGF‐SPH method avoids nested approximation of first order derivatives, and keeps the good feature of ‘kernel gradient free’. The two‐dimensional incompressible fluid flow of shear cavity, both in Euler frame and Lagrangian frame, are simulated by SPH, FPM, the original KGF‐SPH and improved KGF‐SPH. The numerical results show that the improved KGF‐SPH with the novel discrete scheme of Laplacian operator are more accurate than SPH, and more stable than FPM and the original KGF‐SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of soil–water interaction using smoothed particle hydrodynamics (SPH) method

An application of smoothed particle hydrodynamics (SPH) to simulation of soil–water interaction is presented. In this calculation, water is modeled as a viscous fluid with week compressibility and soil is modeled as an elastic–perfectly plastic material. The Mohr–Coulomb failure criterion is applied to describe the stress states of soil in the plastic flow regime. Dry soil is modeled by one-phase flow while saturated soil is modeled by separate water and soil phases. Interaction between soil and water is taken into account by means of pore water pressure and seepage force. Simulation tests of soil excavation by a water jet are calculated as a challenging example to verify the broad applicability of the SPH method. The excavations are carried out in two different soil models, one is dry soil and the other is fully saturated soil. Numerical results obtained in this paper have shown that the gross discontinuities of soil failure can be simulated without any difficulties. This supports the feasibility and attractiveness of this a new approach in geomechanics applications. Advantages of the method are robustness, conceptual simplicity and relative ease of incorporating new physics.     

A coupled WC-TL SPH method for simulation of hydroelastic problems

A coupled weakly compressible (WC) and total Lagrangian (TL) smoothed particle hydrodynamics (SPH) method is developed for simulating hydroelastic problems. The fluid phase is simulated using WCSPH method, while the structural dynamics are solved using TLSPH method. Fluid and solid components of the method are validated separately. A sloshing water tank problem is solved to test the WCSPH method while oscillation of a thin plate and large deformation of a cantilever beam are simulated to test the TLSPH method. After validating each component, the coupled WC-TL SPH scheme is used to simulate two benchmark hydroelastic problems. The first test case shows the evolution of water column with an elastic boundary gate, and the second one investigates the breaking water column impact on elastic structures. The agreement between WC-TL SPH results and literature data shows the ability of the proposed method in simulating hydroelastic phenomena.     

基于拟流体模型的SPH新方法及其在弹丸超高速碰撞薄板中的应用

引入颗粒动力学理论(拟流体模型)建立了适用于超高速碰撞的SPH新方法。将超高速碰撞中处于损伤状态的碎片等效为拟流体,在描述其运动过程中引入了碎片间相互作用和气体相对碎片的作用。采用该方法对球形弹丸超高速碰撞薄板形成碎片云的过程进行了数值模拟,得到了弹坑直径、外泡碎片云和内核碎片云的形状、分布,并与使用传统SPH方法、自适应光滑粒子流体动力学(ASPH)方法的模拟结果进行对比,结果显示:新方法在内核碎片云形状和分布上计算结果更加准确。同时对Whipple屏超高速碰撞问题进行了研究,分析了不同撞击速度下防护屏弹坑尺寸及舱壁损伤特性等特性,计算结果与实验吻合较好且符合Whipple防护结构的典型撞击极限曲线。     

A corrected symmetric SPH method to simulate viscoelastic free surface flows based on the PTT model

In this work, a corrected symmetric and periodic density reinitialized SPH (CSPDR‐SPH) method is proposed and extended to simulate the viscoelastic free surface flows based on the Phan–Thien–Tanner model. The improvements mainly lie in deriving a corrected symmetric kernel gradient, and combining it with a periodic density reinitialization procedure. In addition, a simple artificial viscosity and a simple artificial stress form are adopted. Thus, the CSPDR‐SPH method has higher accuracy and better stability than the SPH method, and conserves both linear and angular momentums. The consistency and convergence of the CSPDR‐SPH method are justified by approximating a function in one and two dimensions. The merits of CSPDR‐SPH method are demonstrated by several benchmarks. The simple flow in a two‐dimensional channel is investigated to show the capability of the CSPDR‐SPH method to simulate the viscoelastic free surface flow. Then the CSPDR‐SPH method is extended to simulate the impacting drop problem. Numerical results show that the CSPDR‐SPH method can precisely capture the viscoelastic free surface. The Reynolds number, Weissenberg number and elongation parameter have remarkable effect on the flows. Copyright © 2012 John Wiley & Sons, Ltd.