Simulation of particle dynamics in a rapidly varying viscous flow

A method for the numerical simulation of the dynamics of particles in a rapidly varying viscous flow has been developed and implemented as a software package. The frequency of variations in the fluid velocity is assumed to be such that the nonlinear terms in the equations of motion can be neglected in comparison with the nonstationary terms. The hydrodynamic interaction of the particles is taken into account. The velocities of the particles and their trajectories are computed. It is found that the trajectories of the particles depend substantially on the ratio of their radii. In the case of dipole particles in a rapidly varying external magnetic field, the hydrodynamic interaction is shown to prevent the particles from approaching each other under the influence of dipole-dipole interaction forces.     

Dynamical models of molecular chains and efficient integration algorithms

Chains of point masses and chains of rigid bodies are used to model biological polymers. To investigate their dynamics we propose a method which allows an efficient realization of the constraints jointly with a simple and accurate integration of the free rigid body motion. The method is quite effective to evolute the geodesic flow of a rigid body chain and the global performance depends on the computational complexity of the algorithms used to compute the interaction forces. Our approach is suitable to describe a chain of rigid bodies immersed in a thermal bath. In the method we propose, the constraints are realized by hard springs whose elastic constant is set to maximize the energy dissipation rate of a Runge–Kutta integrator scheme. Moreover the use of local Lagrangian coordinates is introduced using the possibility of a continuous change of chart, such that the distance from the coordinate singularities is the highest possible. For a chain of point masses the numerical results are checked with another method where the constraints are exactly realized by means of Lagrangian coordinates. When the chain is subject to regular interactions potentials plus a thermal bath the exact and approximate constraints realization provide comparable results.     

Simulation of particle dynamics in a viscous fluid near a plane wall

A method was developed (and implemented as a software package) for the numerical simulation of the dynamics of a large number of particles in a viscous flow near a plane wall. The fluid motion is considered in the Stokes approximation with allowance for the hydrodynamic interaction of all the particles. The velocities and trajectories of the particles are computed, and fluid streamlines are constructed. In the case of even two particles falling onto the wall, it is found that a complex vortex flow is formed between them. This gives rise to the particles having a velocity component along the wall. The dynamics of a group of several dozen particles moving parallel and perpendicular to the wall are simulated.     

Regular and chaotic recurrence in FPU cell-chains

In contrast to the classical Fermi–Pasta–Ulam (FPU) chain, the inhomogeneous FPU chain shows nearly all the principal resonances. Using this fact, we can construct a periodic FPU chain of low dimension, for instance a FPU cell of four degrees-of-freedom, that can be used as a building block for a FPU cell-chain. This is a new type of chain. Differences between chains in nearest-neighbor interaction and those in overall interaction are caused by symmetry. We will show some striking results on the dynamics of FPU cell-chains where near-integrable behavior or chaos plays a part near stable equilibrium.     

供应链交互作用对牛鞭效应的影响分析

考虑两平行供应链系统,建立了需求依赖于两种产品价格的需求函数模型,分析了平行供应链交互作用对牛鞭效应的影响。研究表明:(1)供应链交互作用可能增加或减弱牛鞭效应。(2)对于具有产品可替代性的两竞争型供应链系统,若产品价格交互敏感性不强,则较大协方差的引入可以抑制牛鞭效应。(3)对于具有产品互补性的两合作型供应链系统,若产品价格交互敏感性较强,则较小协方差的引入可以抑制牛鞭效应。     

泥石流固液分相流速计算方法研究

泥石流固液分相流速是泥石流对岸坡、防治结构冲击、磨损机理的核心问题．将泥石流体简化为具有相同粒径的固相和具有相同力学性质的液相,基于泥石流体为沿流动方向的一维两相流体,运用两相流理论建立了泥石流固液分相流速控制方程．构建了泥石流平均压力、彻体力及平均表面力的计算方法,尤其通过浆体的Binhanm体流变方程、Bagnold颗粒相互作用试验成果建立了控制体平均表面力计算方法;建立了固液两相流速比例系数,以及理论固相流速与实际流速的比例系数．据此求解控制方程得到了固液分相流速计算方法,该方法既可同时适用于粘性泥石流和稀性泥石流,也可在泥石流爆发以后通过现场采集沉积物分析反求泥石流爆发期间的分相流速．工程实例分析显示,该方法计算结果与实测结果吻合较好．     

Effect of spherically symmetric mass flow from the surface of a particle on the force of interaction with a plane surface

A stationary velocity field of the flow of a gaseous medium generated by uniform radial injection from the surface of a spherical particle near a wall is considered in the Stokes' approximation. Bispherical coordinates are used to write the expression for the stream function. A formula is obtained for the force acting on the spherical particle when there is an arbitrary mass flow from its surface, generalizing earlier results /1, 2/. An expression for the force acting on the particle is obtained for the case of spherically symmetric injection from the surface of the particle, and asymptotic formulas at short and long distances from the wall are studied.

An analogous problem concerning the forces of interaction between two spherical particles of the same radius, when uniform injection of equal intensity takes place from their surfaces, is discussed. This is equivalent to the problem of the interaction of a spherical particle with a free surface. A general expression for the force of interaction, and its asymptotic forms for short and long distances, are obtained.     

Fluid-Structure Port-Hamiltonian Model for Incompressible Flows in Tubes with Time Varying Geometries

ABSTRACT

A simple and scalable finite-dimensional model based on the port-Hamiltonian framework is proposed to describe the fluid–structure interaction in tubes with time-varying geometries. For this purpose, the moving tube wall is described by a set of mass-spring-damper systems while the fluid is considered as a one-dimensional incompressible flow described by the average momentum dynamics in a set of incompressible flow sections. To couple these flow sections small compressible volumes are defined to describe the pressure between two adjacent fluid sections. The fluid-structure coupling is done through a power-preserving interconnection between velocities and forces. The resultant model includes external inputs for the fluid and inputs for external forces over the mechanical part that can be used for control or interconnection purposes. Numerical examples show the accordance of this simplified model with finite-element models reported in the literature.     

Turbulent particle dispersion in an electrostatic precipitator

The behaviour of charged particles in turbulent gas flow in electrostatic precipitators (ESPs) is crucial information to optimise precipitator efficiency. This paper describes a strongly coupled calculation procedure for the rigorous computation of particle dynamics during ESP taking into account the statistical particle size distribution. The turbulent gas flow and the particle motion under electrostatic forces are calculated by using the commercial computational fluid dynamics (CFD) package FLUENT linked to a finite volume solver for the electric field and ion charge. Particle charge is determined from both local electrical conditions and the cell residence time which the particle has experienced through its path. Particle charge density and the particle velocity are averaged in a control volume to use Lagrangian information of the particle motion in calculating the gas and electric fields. The turbulent particulate transport and the effects of particulate space charge on the electrical current flow are investigated. The calculated results for poly-dispersed particles are compared with those for mono-dispersed particles, and significant differences are demonstrated.     

The effect of environmental parameters on product recovery

Economical and environmental issues are the main driving forces for the development of closed-loop supply chains. This paper examines the impact of environmental issues on long-term behaviour of a single product supply chain with product recovery. The environmental issues examined are the firm's `green image' effect on customer demand, the take back obligation imposed by legislation, and the state campaigns for proper disposal of used products. The behaviour of the system is analyzed through a dynamic simulation model based on the principles of the system dynamics (SD) methodology. This model includes all major inventories of new, used and recovered products and the flows among them. Inventory levels and flow rates are linked through differential equations. The dynamic model provides an experimental simulation tool, which can be used to evaluate the effect of environmental issues on long-term decision making in collection and remanufacturing activities and on product demand. Numerical analysis illustrates the potential uses of the methodology.     

The Invariant Measures of Markov Chains on Product Spaces and a Measurement of Dependence

For probability measures on product spaces, we define a notion of dependence among each coordinate. We study Markov chains on product spaces in which the time developement of each coordinate corresponds to the movement of a particle in an interacting particle system described by the Markov chain. If there is no interaction among these particles, the dependence of them decreases monotonically. We establish an inequality which states that if the interaction among each particle is small, then, in stationarity, the dependence among each particle is small.     

Features of the Formation of Composite Materials in Anisotropic Liquid Systems with Elongated Suspended Particles

The equations of rotational motion of nondeformable spherical and axisymmetric elongated particles and a rheologic equation for stresses in arbitrary gradient flows of dilute suspensions of such particles in an anisotropic carrying fluid are obtained within the framework of a structural-phenomenological approach. As a rheologic model of the suspension-carrying fluid and a hydrodynamic model of the suspended particles, we use the Ericksen simple anisotropic fluid and a symmetric triaxial dumbbell, respectively. The constitutive equations obtained are used to study the effect of anisotropy of the carrying fluid on the dynamics of suspended particles and on the rheologic properties of suspensions in a simple shear flow. A stationary orientation of the elongated suspended particles under the action of hydrodynamic forces is discovered. The possibility of applying this phenomenon to the formation of composite materials is discussed.     

What are causes of cash flow bullwhip effect in centralized and decentralized supply chains?

Bullwhip effect in supply chain is a phenomenon which can emerge in both inventory levels and replenishment orders. Bullwhip effect causes variations in cash conversion cycle (CCC) across cash flow of supply chain. As a result, it can lead to inefficiencies such as cash flow bullwhip (CFB). Due to negative impact of CFB on cash flow of supply chain, it can lead to a decrease in efficiency of supply chain management (SCM). That is why supply chain modeling is a proper start point for effective management and control of the CFB. This paper aims to analyze concurrent impact of causes of inventory bullwhip effect and effect of their interactions on CFB based on generalized OUT policy from aspect of CCC variance. To this end, first we develop system dynamics structure of beer distribution game as simulation model which includes multi-stage supply chain under both centralized and decentralized supply chains. Then, in order to develop CFB function, we design experiments in developed simulation model using response surface methodology (RSM). Results demonstrate that if each chain member uses generalized OUT policy as replenishment model, there still exists CFB in both chains and CFB largely stems from rationing and shortage gaming in both centralized and decentralized supply chain. In addition, when information on ordering parameters are not shared among members, parameters of downstream stage (i.e. retailer) are more important than parameters of upstream stage (i.e. manufacturer) in reducing CFB function.     

A multi-structural framework for adaptive supply chain planning and operations control with structure dynamics considerations

A trend in up-to-date developments in supply chain management (SCM) is to make supply chains more agile, flexible, and responsive. In supply chains, different structures (functional, organizational, informational, financial, etc.) are (re)formed. These structures interrelate with each other and change in dynamics. The paper introduces a new conceptual framework for multi-structural planning and operations of adaptive supply chains with structure dynamics considerations. We elaborate a vision of adaptive supply chain management (A-SCM), a new dynamic model and tools for the planning and control of adaptive supply chains. SCM is addressed from perspectives of execution dynamics under uncertainty. Supply chains are modelled in terms of dynamic multi-structural macro-states, based on simultaneous consideration of the management as a function of both states and structures. The research approach is theoretically based on the combined application of control theory, operations research, and agent-based modelling. The findings suggest constructive ways to implement multi-structural supply chain management and to transit from a “one-way” partial optimization to the feedback-based, closed-loop adaptive supply chain optimization and execution management for value chain adaptability, stability and crisis-resistance. The proposed methodology enhances managerial insight into advanced supply chain management.     

Numerical investigation of blood flow. Part I: In microvessel bifurcations

In some diseases there is a focal pattern of velocity in regions of bifurcation, and thus the dynamics of bifurcation has been investigated in this work. A computational model of blood flow through branching geometries has been used to investigate the influence of bifurcation on blood flow distribution. The flow analysis applies the time-dependent, three-dimensional, incompressible Navier–Stokes equations for Newtonian fluids. The governing equations of mass and momentum conservation were solved to calculate the pressure and velocity fields. Movement of blood flow from an arteriole to a venule via a capillary has been simulated using the volume of fluid (VOF) method. The proposed simulation method would be a useful tool in understanding the hydrodynamics of blood flow where the interaction between the RBC deformation and blood flow movement is important. Discrete particle simulation has been used to simulate the blood flow in a bifurcation with solid and fluid particles. The fluid particle method allows for modeling the plasma as a particle ensemble, where each particle represents a collective unit of fluid, which is defined by its mass, moment of inertia, and translational and angular momenta. These kinds of simulations open a new way for modeling the dynamics of complex, viscoelastic fluids at the micro-scale, where both liquid and solid phases are treated with discrete particles.     

A chain of interacting particles under strain

We investigate the behaviour of a chain of interacting Brownian particles with one end fixed and the other end moving away at slow speed ε>0, in the limit of small noise. The interaction between particles is through a pairwise potential U with finite range b>0. We consider both overdamped and underdamped dynamics.     

Non-linearity of multibody dynamic equations with respect to Lagrange multipliers: application to railway dynamics

In this paper, the dynamics of multibody systems with closed kinematical chains of bodies is considered. The main focus is set on non-linearity of the multibody equations with respect to the Lagrange multipliers. When closed chains are considered, loop cutting procedure is a solution to express the constraint equations associated with the loops. Dynamic equations of the multibody tree-like structure are thus completed with the constraint forces via the Lagrange multipliers. In the considered case of railway vehicles, constraints arise from the contact between the rigid wheels and the rails. Corresponding contact forces applied to the wheels appears via the Lagrange multipliers λ and the tangent creep forces as well. Resulting differential-algebraic equations can be transformed into an ODE system and then time-integrated using the coordinate partitioning method [3], when the system is linear with respect to λ. This paper presents an algorithm allowing us to solve this system in case of nonlinearities with respect to λ, which is typical of wheel/rail contact force models. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Challenges of living in the harsh environments: A mathematical modeling study

We propose and analyze a mathematical model, which mimics community dynamics of plants and animals in harsh environments. The mathematical model exploits type IV functional responses whose idiosyncrasies have been recognized only in recent years. The interaction of the middle predator with the top predator is cast into Leslie-Gower scheme. Linear and non-linear stability analyses are performed to get an idea of the stability behavior of the model food chain. It turns out that carrying capacity of the prey and the immunity parameter of the middle predator are two crucial parameters governing the model. Availability of alternative food options to the generalist predator also plays a key role in deciding the model dynamics.Simulation runs performed on this model provide insight into population dynamics of monkeys of macaque family found in northern Japan. These monkeys are social animals which reproduce sexually. The characteristic feature of the model dynamics is that the generalist predator (macaque monkeys) is able to avoid impending extinction frequently and recovers at a rate which falsify threats from exogenous external forces; extreme weather conditions, etc.     

Modeling and simulation of the motion of nanoparticles in cylindrical capillaries allowing particle‐to‐wall interactions

The process of transporting nanoparticles at the blood vessels level stumbles upon various physical and physiological obstacles; therefore, a Mathematical modeling will provide a valuable means through which to understand better this complexity. In this paper, we consider the motion of nanoparticles in capillaries having cylindrical shapes (i.e., tubes of finite size). Under the assumption that these particles have spherical shapes, the motion of these particles reduces to the motion of their centers. Under these conditions, we derive the mathematical model, to describe the motion of these centers, from the equilibrium of the gravitational force, the hemodynamic force and the van der Waals interaction forces. We distinguish between the interaction between the particles and the interaction between each particle and the walls of the tube. Assuming that the minimum distance between the particles is large compared with the maximum radius R of the particles and hence neglecting the interactions between the particles, we derive simpler models for each particle taking into account the particles‐to‐wall interactions. At an error of order O(R) or O(R³)(depending if the particles are 'near' or 'very near' to the walls), we show that the horizontal component of each particle's displacement is solution of a nonlinear integral equation that we can solve via the fixed point theory. The vertical components of the displacement are computable in a straightforward manner as soon as the horizontal components are estimated. Finally, we support this theory with several numerical tests. Copyright © 2016 John Wiley & Sons, Ltd.