Kernel-independent fast multipole method within the framework of regularized Stokeslets

The method of regularized Stokeslets (MRS) uses a radially symmetric blob function of infinite support to smooth point forces and allows for evaluation of the resulting flow field. This is a common method to study swimmers at zero Reynolds number where the Stokeslet is the fundamental solution corresponding to the kernel of the single layer potential. Simulating the collective motion of N micro-swimmers using the MRS results in at least N² pair-wise interactions. Efficient simulation of a large number of swimmers in free space is observed with the implementation of the kernel-independent fast multipole method (FMM) for radial basis functions. We illustrate the complexity of the algorithm on a simple test case where we study regularized point forces, showing that the method is of order N. Additionally, we explore accuracy in time for the MRS where the swimmers are modeled as Kirchhoff rods and the kernel-independent FMM is compared to the direct calculation using the standard MRS. Optimal hydrodynamic efficiency is also explored for different configurations of swimmers.     

Method of fundamental solutions for partial‐slip fibrous filtration flows

In this study a Stokeslet‐based method of fundamental solutions (MFS) for two‐dimensional low Reynolds number partial‐slip flows has been developed. First, the flow past an infinitely long cylinder is selected as a benchmark. The numerical accuracy is investigated in terms of the location and the number of the Stokeslets. The benchmark study shows that the numerical accuracy increases when the Stokeslets are submerged deeper beneath the cylinder surface, as long as the formed linear system remains numerically solvable. The maximum submergence depth increases with the decrease in the number of Stokeslets. As a result, the numerical accuracy does not deteriorate with the dramatic decrease in the number of Stokeslets. A relatively small number of Stokeslets with a substantial submergence depth is thus chosen for modeling fibrous filtration flows. The developed methodology is further examined by application to Taylor–Couette flows. A good agreement between the numerical and analytical results is observed for no‐slip and partial‐slip boundary conditions. Next, the flow about a representative set of infinitely long cylindrical fibers confined between two planar walls is considered to represent the fibrous filter flow. The obtained flowfield and pressure drop agree very well with the experimental data for this setup of fibers. The developed MFS with submerged Stokeslets is then applied to partial‐slip flows about fibers to investigate the slip effect at fiber–fluid interface on the pressure drop. The numerical results compare qualitatively with the analytical solution available for the limit case of infinite number of fibers. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydroelastic stability of a rectangular plate interacting with a layer of ideal flowing fluid

The three-dimensional formulation of the problem on the natural vibrations and stability of an elastic plate which interacts with a quiescent or flowing fluid is represented and a finite element algorithm of its numerical implementation is proposed. The governing equations, which describe vortex-free ideal fluid dynamics in the case of small perturbations, are written in terms of the perturbation velocity potential and transformed using the Bubnov–Galerkin method. The plate strains are determined on the basis of the Timoshenko theory. The variational principle of virtual displacements which takes into account the work done by inertial forces and the hydrodynamic pressure is used for the mathematical formulation of the dynamic problem of elastic structure. The solution of the problem is reduced to calculations and an analysis of complex eigenvalues of a coupled system of two equations. The effect of the fluid layer height on the eigenfrequencies and the critical velocities of the loss of stability is estimated numerically. It is shown that there exist different types of instability determined by combinations of the kinematic boundary conditions prescribed at the plate edges.     

A unified Chebyshev–Ritz formulation for vibration analysis of composite laminated deep open shells with arbitrary boundary conditions

In this paper, a unified Chebyshev–Ritz formulation is presented to investigate the vibrations of composite laminated deep open shells with various shell curvatures and arbitrary restraints, including cylindrical, conical and spherical ones. The general first-order shear deformation shell theory is employed to include the effects of rotary inertias and shear deformation. Under the current framework, regardless of boundary conditions, each of displacements and rotations of the open shells is invariantly expressed as Chebyshev orthogonal polynomials of first kind in both directions. Then, the accurate solutions are obtained by using the Rayleigh–Ritz procedure based on the energy functional of the open shells. The convergence and accuracy of the present formulation are verified by a considerable number of convergence tests and comparisons. A variety of numerical examples are presented for the vibrations of the composite laminated deep shells with various geometric dimensions and lamination schemes. Different sets of classical constraints, elastic supports as well as their combinations are considered. These results may serve as reference data for future researches. Parametric studies are also undertaken, giving insight into the effects of elastic restraint parameters, fiber orientation, layer number, subtended angle as well as conical angle on the vibration frequencies of the composite open shells.     

Interaction between an incompressible flow and elastic cantilevers of circular cross-section

The purpose of this work is to study the deformation of elastic cantilevers due to hydrodynamic forces by coupled fluid–structure interaction simulations. The cantilever is placed in a rectangular duct and the Reynolds number based on bulk velocity and cantilever diameter is 400. Reduced velocities in the range π/4 to 2π are studied, which covers both un-synchronised motion and the initial branch of synchronisation. The cantilever surface is represented by a virtual boundary method which replaces a solid object in flow by additional force distribution to satisfy local boundary condition. The flow field is solved using a Cartesian finite difference code and the deformation of the cylinder a finite element approach using one-dimensional beam elements is used.     

A Variational Approach to Obstacle Problems for Shearable Nonlinearly Elastic Rods

We use variational methods to study obstacle problems for geometrically exact (Cosserat) theories for the planar deformation of nonlinearly elastic rods. These rods can suffer flexure, extension, and shear. There is a marked difference between the behavior of a shearable and an unshearable rod. The set of admissible deformations is not convex, because of the exact geometry used. We first investigate the fundamental question of describing contact forces, which we necessarily treat as vector‐valued Borel measures. Moreover, we introduce techniques for describing point obstacles. Then we prove existence for a very large class of problems. Finally, using nonsmooth analysis for handling the obstacle, we show that the Euler‐Lagrange equations are satisfied almost everywhere. These equations provide very detailed structural information about the contact forces. Accepted June 3, 1996     

Analytical dynamics of continuous medium and its application

综述了国内和国外学者研究连续介质分析动力学问题的进展,阐明了本文主要论述将Lagrange方程应用于连续介质动力学的问题.论文采用Lagrange-Hamilton体系,分别论述了非保守非线性弹性动力学、不可压缩黏性流体动力学、黏弹性动力学、热弹性动力学、刚--弹耦合动力学和刚--液耦合动力学的Lagrange方程及其应用.论述了应用Lagrange方程建立有限元计算模型的问题. 最后,展望了将Lagrange方程应用于连续介质动力学问题的研究前景.     

连续介质分析动力学及其应用

综述了国内和国外学者研究连续介质分析动力学问题的进展,阐明了本文主要论述将Lagrange方程应用于连续介质动力学的问题.论文采用Lagrange-Hamilton体系,分别论述了非保守非线性弹性动力学、不可压缩黏性流体动力学、黏弹性动力学、热弹性动力学、刚--弹耦合动力学和刚--液耦合动力学的Lagrange方程及其应用.论述了应用Lagrange方程建立有限元计算模型的问题. 最后,展望了将Lagrange方程应用于连续介质动力学问题的研究前景.      

Analysis of mode Ⅲ crack perpendicular to the interface between two dissimilar strips

The present work is concerned with the problem of mode Ⅲ crack perpendicular to the interface of a bi-strip composite. One of these strips is made of a functionally graded material and the other of an isotropic material, which contains an edge crack perpendicular to and terminating at the interface. Fourier transforms and asymptotic analysis are employed to reduce the problem to a singular integral equation which is numerically solved using Gauss-Chebyshev quadrature formulae. Furthermore, a parametric study is carried out to investigate the effects of elastic and geometric characteristics of the composite on the values of stress intensity factor.     

Existence of Weak Solutions for the Three-Dimensional Motion of an Elastic Structure in an Incompressible Fluid

We study here the three-dimensional motion of an elastic structure immersed in an incompressible viscous fluid. The structure and the fluid are contained in a fixed bounded connected set Ω. We show the existence of a weak solution for regularized elastic deformations as long as elastic deformations are not too important (in order to avoid interpenetration and preserve orientation on the structure) and no collisions between the structure and the boundary occur. As the structure moves freely in the fluid, it seems natural (and it corresponds to many physical applications) to consider that its rigid motion (translation and rotation) may be large. The existence result presented here has been announced in [4]. Some improvements have been provided on the model: the model considered in [4] is a simplified model where the structure motion is modelled by decoupled and linear equations for the translation, the rotation and the purely elastic displacement. In what follows, we consider on the structure a model which represents the motion of a structure with large rigid displacements and small elastic perturbations. This model, introduced by [15] for a structure alone, leads to coupled and nonlinear equations for the translation, the rotation and the elastic displacement.     

Numerical simulation of viscous flow interaction with an elastic membrane

A numerical fluid–structure interaction model is developed for the analysis of viscous flow over elastic membrane structures. The Navier–Stokes equations are discretized on a moving body‐fitted unstructured triangular grid using the finite volume method, taking into account grid non‐orthogonality, and implementing the SIMPLE algorithm for pressure solution, power law implicit differencing and Rhie–Chow explicit mass flux interpolations. The membrane is discretized as a set of links that coincide with a subset of the fluid mesh edges. A new model is introduced to distribute local and global elastic effects to aid stability of the structure model and damping effects are also included. A pseudo‐structural approach using a balance of mesh edge spring tensions and cell internal pressures controls the motion of fluid mesh nodes based on the displacements of the membrane. Following initial validation, the model is applied to the case of a two‐dimensional membrane pinned at both ends at an angle of attack of 4° to the oncoming flow, at a Reynolds number based on the chord length of 4 × 10³. A series of tests on membranes of different elastic stiffness investigates their unsteady movements over time. The membranes of higher elastic stiffness adopt a stable equilibrium shape, while the membrane of lowest elastic stiffness demonstrates unstable interactions between its inflated shape and the resulting unsteady wake. These unstable effects are shown to be significantly magnified by the flexible nature of the membrane compared with a rigid surface of the same average shape. Copyright © 2007 John Wiley & Sons, Ltd.     

Generalized Euler–Lagrange equation for nonsmooth calculus of variations

In this paper, we first introduce a novel generalized derivative and obtain the generalized first-order Taylor expansion of the nonsmooth functions. Then we derive the generalized Euler–Lagrange equation for the nonsmooth calculus of variations and solve this equation by using Chebyshev pseudospectral method, approximately. Finally, the optimal solutions of some problems in the nonsmooth calculus of variations are approximated.     

三维裂纹问题的高精度数值解法

提出了一种求解三维均质弹性体中任意形状平片裂纹问题超奇异积分方程组的Chebyshev多项式数值解法。数值计算结果表明：文中方法不仅收敛快，而且精度高。     

A binary-tree based model for rate-independent polycrystals

A model of a rigid-plastic rate-independent polycrystalline aggregate wherein sub-aggregates are represented as the nodes of a binary tree is proposed. The lowest nodes of the binary tree represent grains. Higher binary tree nodes represent increasingly larger sub-aggregates of grains, culminating with the root of the tree, which represents the entire polycrystalline aggregate. Planar interfaces are assumed to separate the sub-aggregates represented by nodes in the binary tree. Equivalence between the governing equations of the model and a standard linear program is established. The objective function of the linear program is given by the plastic power associated with polycrystal deformation and the linear constraints are given by compatibility requirements between the sub-aggregates represented by sibling nodes in the binary tree. The deviatoric part of the Cauchy stress in each sub-aggregate is deduced as linear combinations of the Lagrange multipliers associated with the constraints. It is shown that the present model allows a generalization of Taylor’s principle to polycrystals. The proposed model is applied to simulate tensile, compressive, torsional, and plane-strain deformation of copper polycrystals. The predicted macroscopic response is in good agreement with published experimental data. The effect of the initial distribution of the planar interfaces separating the sub-aggregates represented by the binary tree on the predicted mechanical response in tension, compression and torsion is studied. Also, the role of constraints relaxation in simulations of plane strain compression is investigated in detail.     

Finite element computations for unsteady fluid and elastic membrane interaction problems

The interaction between the hydrodynamic forces of a flow field and the elastic forces of adjacent deformable boundaries is described by elastohydrodynamics, a coupled fluid–elastic membrane problem. Direct numerical solution of the unsteady, highly non-linear equations requires that the dynamic evolution of both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. This paper describes a numerical algorithm based on the deformable spatial domain space–time (DSD/ST) finite element method for the unsteady motion of an incompressible, viscous fluid with elastic membrane interaction. The unsteady Navier–Stoke and elastic membrane equations are solved separately using an iterative procedure by the GMRES technique with an incomplete lower-upper (ILU) decomposition at every time instant. One-dimensional, two-dimensional and deformable domain model problems are used to demonstrate the capabilities and accuracy of the present algorithm. Both steady state and transient problems are studied. © 1997 John Wiley & Sons, Ltd.     

Viscoelastic flow analysis using the software OpenFOAM and differential constitutive equations

Viscoelastic fluids are of great importance in many industrial sectors, such as in food and synthetic polymers industries. The rheological response of viscoelastic fluids is quite complex, including combination of viscous and elastic effects and non-linear phenomena. This work presents a numerical methodology based on the split-stress tensor approach and the concept of equilibrium stress tensor to treat high Weissenberg number problems using any differential constitutive equations. The proposed methodology was implemented in a new computational fluid dynamics (CFD) tool and consists of a viscoelastic fluid module included in the OpenFOAM, a flexible open source CFD package. Oldroyd-B/UCM, Giesekus, Phan-Thien–Tanner (PTT), Finitely Extensible Nonlinear Elastic (FENE-P and FENE-CR), and Pom–Pom based constitutive equations were implemented, in single and multimode forms. The proposed methodology was evaluated by comparing its predictions with experimental and numerical data from the literature for the analysis of a planar 4:1 contraction flow, showing to be stable and efficient.     

Nonlinear normal modes in multi-mode models of an inertially coupled elastic structure

Nonlinear normal modes for elastic structures have been studied extensively in the literature. Most studies have been limited to small nonlinear motions and to structures with geometric nonlinearities. This work investigates the nonlinear normal modes in elastic structures that contain essential inertial nonlinearities. For such structures, based on the works of Crespo da Silva and Meirovitch, a general methodology is developed for obtaining multi-degree-of-freedom discretized models for structures in planar motion. The motion of each substructure is represented by a finite number of substructure admissible functions in a way that the geometric compatibility conditions are automatically assured. The multi degree-of-freedom reduced-order models capture the essential dynamics of the system and also retain explicit dependence on important physical parameters such that parametric studies can be conducted. The specific structure considered is a 3-beam elastic structure with a tip mass. Internal resonance conditions between different linear modes of the structure are identified. For the case of 1:2 internal resonance between two global modes of the structure, a two-mode nonlinear model is then developed and nonlinear normal modes for the structure are studied by the method of multiple time scales as well as by a numerical shooting technique. Bifurcations in the nonlinear normal modes are shown to arise as a function of the internal mistuning that represents variations in the tip mass in the structure. The results of the two techniques are also compared.     

Dynamic response of tower crane induced by the pendulum motion of the payload

Dynamic response of tower cranes coupled with the pendulum motions of the payload is studied in this paper. A simple perturbation scheme and the assumption of small pendulum angle are applied to simplify the governing equation. The tower crane is modeled by the finite element method, while the pendulum motion is represented as rigid-body kinetics. Integrated governing equations for the coupled dynamics problem are derived based on Lagrange’s equations including the dissipation function. Dynamics of a real luffing crane model with the spherical and planar pendulum motions is analyzed using the proposed formulations and computational method. It is found that the dynamic responses of the tower crane are dominated by both the first few natural frequencies of crane structure and the pendulum motion of the payload. The dynamic amplification factors generally increase with the increase of the initial pendulum angle and the changes are just slightly nonlinear for the planar pendulum motion.     

Existence of a Weak Solution to a Nonlinear Fluid–Structure Interaction Problem Modeling the Flow of an Incompressible, Viscous Fluid in a Cylinder with Deformable Walls

We study a nonlinear, unsteady, moving boundary, fluid–structure interaction (FSI) problem arising in modeling blood flow through elastic and viscoelastic arteries. The fluid flow, which is driven by the time-dependent pressure data, is governed by two-dimensional incompressible Navier–Stokes equations, while the elastodynamics of the cylindrical wall is modeled by the one-dimensional cylindrical Koiter shell model. Two cases are considered: the linearly viscoelastic and the linearly elastic Koiter shell. The fluid and structure are fully coupled (two-way coupling) via the kinematic and dynamic lateral boundary conditions describing continuity of velocity (the no-slip condition), and the balance of contact forces at the fluid–structure interface. We prove the existence of weak solutions to the two FSI problems (the viscoelastic and the elastic case) as long as the cylinder radius is greater than zero. The proof is based on a novel semi-discrete, operator splitting numerical scheme, known as the kinematically coupled scheme, introduced in Guidoboni et al. (J Comput Phys 228(18):6916–6937, 2009) to numerically solve the underlying FSI problems. The backbone of the kinematically coupled scheme is the well-known Marchuk–Yanenko scheme, also known as the Lie splitting scheme. We effectively prove convergence of that numerical scheme to a solution of the corresponding FSI problem.