Wave body interaction in 2D using smoothed particle hydrodynamics (SPH) with variable particle mass

Wave interaction with bodies is an important practical application for smoothed particle hydrodynamics (SPH) which in principle applies to steep and breaking waves without special treatment. However, few detailed tests have been undertaken even with small amplitude waves. In order to reduce computer time a variable particle mass distribution is tested here with fine resolution near the body and coarse resolution further away, while maintaining a uniform kernel size. We consider two well‐defined test cases, in two dimensions, of waves generated by a heaving semi‐immersed cylinder and progressive waves interacting with a fixed cylinder. But first, still water with hydrostatic pressure is tested. The open‐source code SPHysics ( http://www.sphysics.org )^{§Update made here after initial online publication.} is used with a Riemann solver in an Arbitrary Lagrangian–Eulerian formulation. For the heaving cylinder, SPH results for far field wave amplitude and cylinder force show good agreement with the data of Yu and Ursell (J. Fluid Mech. 1961; 11 :529–551). For wave loading on a half‐submerged cylinder the agreement with the experimental data of Dixon et al. (J. Waterway Port Coastal Ocean Div. 1979; 105 :421–438) for the root mean square force is within 2%. For more submerged cases, the results show some discrepancy, but this was also found with other modelling approaches. The sensitivity of results to the value of the slope limiter used in the MUSCL‐based Riemann solver is demonstrated. The variable mass distribution leads to a computer run speedup of nearly 200% in these cases. Copyright © 2011 John Wiley & Sons, Ltd.     

Free convection of a dusty‐gas flow along a semi‐infinite vertical cylinder

An analysis is performed to study the free convection of a dusty‐gas flow along a semi‐infinite isothermal vertical cylinder. The governing equations of the flow problem are transformed into non‐dimensional form and the resulting nonlinear, coupled parabolic partial differential equations have been solved numerically using an implicit finite difference scheme of Crank–Nicholson type. The flow variables such as gas–velocity, dust‐particle velocity and temperature, shearing stress and heat transfer coefficients are calculated numerically for various parameters occurring in the problem. It is observed that due to the presence of dust particles, the gas velocity is found to decrease. Copyright © 2009 John Wiley & Sons, Ltd.     

Start‐up of electrophoresis of an arbitrarily oriented dielectric cylinder

An analytical study is presented for the transient electrophoretic response of a circular cylindrical particle to the step application of an electric field. The electric double layer adjacent to the particle surface is thin but finite compared with the radius of the particle. The time‐evolving electroosmotic velocity at the outer boundary of the double layer is utilized as a slip condition so that the transient momentum conservation equation for the bulk fluid flow is solved. Explicit formulas for the unsteady electrophoretic velocity of the particle are obtained for both axially and transversely applied electric fields, and can be linearly superimposed for an arbitrarily‐oriented applied field. If the cylindrical particle is neutrally buoyant in the suspending fluid, the transient electrophoretic velocity is independent of the orientation of the particle relative to the applied electric field and will be in the direction of the applied field. If the particle is different in density from the fluid, then the direction of electrophoresis will not coincide with that of the applied field until the steady state is attained. The growth of the electrophoretic mobility with the elapsed time for a cylindrical particle is substantially slower than for a spherical particle.     

A ghost fluid Lattice Boltzmann method for complex geometries

A ghost fluid Lattice Boltzmann method (GF‐LBM) is developed in this study to represent complex boundaries in Lattice Boltzmann simulations of fluid flows. Velocity and density values at the ghost points are extrapolated from the fluid interior and domain boundary via obtaining image points along the boundary normal inside the fluid domain. A general bilinear interpolation algorithm is used to obtain values at image points which are then extrapolated to ghost nodes thus satisfying hydrodynamic boundary conditions. The method ensures no‐penetration and no‐slip conditions at the boundaries. Equilibrium distribution functions at the ghost points are computed using the extrapolated values of the hydrodynamic variables, while non‐equilibrium distribution functions are extrapolated from the interior nodes. The method developed is general, and is capable of prescribing Dirichlet as well as Neumann boundary conditions for pressure and velocity. Consistency and second‐order accuracy of the method are established by running three test problems including cylindrical Couette flow, flow between eccentric rotating cylinders and flow over a cylinder in a confined channel. Copyright © 2011 John Wiley & Sons, Ltd.     

Comparison of adjustable permanent magnetic field sources

A permanent magnet assembly in which the flux density can be altered by a mechanical operation is often significantly smaller than comparable electromagnets and also requires no electrical power to operate. In this paper five permanent magnet designs in which the magnetic flux density can be altered are analyzed using numerical simulations, and compared based on the generated magnetic flux density in a sample volume and the amount of magnet material used. The designs are the concentric Halbach cylinder, the two half Halbach cylinders, the two linear Halbach arrays and the four and six rod mangle. The concentric Halbach cylinder design is found to be the best performing design, i.e. the design that provides the most magnetic flux density using the least amount of magnet material. A concentric Halbach cylinder has been constructed and the magnetic flux density, the homogeneity and the direction of the magnetic field are measured and compared with numerical simulation and a good agrement is found.     

Optical system and lens design for Blu-ray disc optical pick-up

During the design process of Blu-ray disc optical pick-up (BD OPU), the optical system and the collimator lens design is especially important. This paper designs an optical system and some lenses for the BD OPU, including collimator lens, beam shaping lens for laser beam shaping, and cylinder lens for signal detecting. In this OPU, we use a triplet lens to collimate the laser beam. At the same time, we build a series of assembly jigs to make sure that each lens can be put into the OPU basement properly. At last, we get the reading spot image and S-curve photo of OPU, which can be used to read Blu-ray disc (BD).     

Homotopy Structures for Algebras over a Monad

Let T be a monad over a category A. Then a homotopy structure for A, defined by a cocylinder P : A A, or path-endofunctor, can be lifted to the category A ^T of Eilenberg–Moore algebras over T, provided that P is consistent with T in a natural sense, i.e. equipped with a natural transformation : T P P T satisfying some obvious axioms. In this way, homotopy can be lifted from well-known, basic situations to various categories of algebras for instance, from topological spaces to topological semigroups, or spaces over a fixed space (fibrewise homotopy), or actions of a fixed topological group (equivariant homotopy); from categories to strict monoidal categories; from chain complexes to associative chain algebras. The interest is given by the possibility of lifting the homotopy operations (as faces, degeneracy, connections, reversion, interchange, vertical composition, etc.) and their axioms from A to A ^T , just by verifying the consistency between these operations and : T P P T. When this holds, the structure we obtain on our category of algebras is sufficiently powerful to ensure the main general properties of homotopy.     

Slow magnetohydrodynamic flow past a circular cylinder

Oseen’ approximations are used to study the slow motion of a viscous, incompressible, electrically conducting fluid past a circular cylinder in the presence of a uniform aligned magnetic field. Using series truncation method, the analytical solutions for the first three terms in the Fourier sine series expansion of the stream function are obtained. Numerical values of the tangential drag for different values of magnetic interaction parameter and viscous Reynolds number are calculated.     

测定圆柱体杨氏模量的原理分析

分析了测量圆柱体弯曲形变杨氏模量的实验原理,得出了理论公式,将测量拉压形变、弯曲形变的杨氏模量以及扭转形变的剪切模量3个实验中样品形状统一成圆柱体.     

垂直入射时无限长分层柱电磁散射的改进算法及应用

结合非均匀球粒子对平面波散射的散射场计算的改进算法,提出了平面波垂直入射无限长分层圆柱散射场快速稳定而有效的改进电磁散射算法。与已有算法相比,改进算法所能计算的无限长非均匀介质圆柱的尺寸参量突破10000,计算层数达到106,并且计算时间很短,最多仅为几秒。该算法可以用于不同的波段以及不同领域的任意无耗或吸收无限长圆柱体散射场的计算。最后将该算法应用于梯度折射率聚合物光纤(GI-POF)散射特性的研究,为非接触、在线测量聚合物光纤折射率分布提供了理论依据。