A predictor–corrector time integration algorithm for dynamic analysis of nonlinear systems

Nonlinear Dynamics - This paper presents a step-by-step time integration algorithm for efficiently solving second-order nonlinear dynamic problems. The method employs the rewriting of motion as two...     

Friction factor correlations for gas–liquid/liquid–liquid flows with drag-reducing polymers in horizontal pipes

In the present paper, two correlations have been developed to predict the effect of DRP on friction factor of two-phase flow for any pipe diameter. The correlations have been verified using published experimental results of DRP added to air–liquid annular flow and for oil–water flows with any flow pattern at the asymptotic state (maximum drag reduction). Such correlations are not available in literature and considered to be very useful in predicting the drag reduction using DRP and in understanding the most significant parameters that affecting the mechanism of drag reduction by polymers in multiphase flow.     

Some issues related to the modeling of interfacial areas in gas–liquid flows,II. Modeling the source terms for dispersed flows

The first part of this paper dealt with the conceptual issues encountered in the definition of the interfacial area and in the derivation of a transport equation. The second part has the following objectives: 1. to address the closure issues for the source terms appearing in the transport equation when dealing with a bubbly flow in a vertical pipe; 2. to provide a list of open questions to be answered before introducing a transport equation for the interfacial area concentration in a thermal–hydraulic computer code.     

A scaling analysis for point–particle approaches to turbulent multiphase flows

Simple dimensional arguments are used in establishing three different regimes of particle time scale, where explicit expression for particle Reynolds number and Stokes number are obtained as a function of nondimensional particle size (d/η)$(d/\eta )$ and density ratio. From a comparative analysis of the different computational approaches available for turbulent multiphase flows it is argued that the point–particle approach is uniquely suited to address turbulent multiphase flows where the Stokes number, defined as the ratio of particle time scale to Kolmogorov time scale (_τp/_τk)$({\tau }_p/{\tau }_k)$, is greater than 1. The Stokes number estimate has been used to establish parameter range where point–particle approach is ideally suited. The point–particle approach can be extended to handle “finite-sized” particles whose diameter approach that of the smallest resolved eddies. However, new challenges arise in the implementation of Lagrangian–Eulerian coupling between the particles and the carrier phase. An approach where the inter-phase momentum and energy coupling can be separated into a deterministic and a stochastic contribution has been suggested.     

A strongly coupled,embedded-boundary method for fluid–structure interactions of elastically mounted rigid bodies

In the present paper, an embedded-boundary formulation that is applicable to fluid–structure interaction problems is presented. The Navier–Stokes equations for incompressible flow are solved on a Cartesian grid which is not aligned with the boundaries of a body that undergoes large-angle/large-displacement rigid body motions through the fixed grid. A strong-coupling scheme is adopted, where the fluid and the structure are treated as elements of a single dynamical system, and all of the governing equations are integrated simultaneously and interactively in the time domain. A demonstration of the accuracy and efficiency of the method is given for a variety of fluid–structure interaction problems.     

Stroh-like formalism for the coupled stretching–bending analysis of composite laminates

Due to the two-dimensional nature of thin plates, the lamination theory considering the composite laminates with in-plane and plate bending problems coupling each other is treated in this paper by using complex variable formulation. By following the steps of Stroh formalism for two-dimensional linear anisotropic elasticity, a displacement complex variable formalism developed by the other researchers was introduced and re-derived in a different but more Stroh-like way. In addition, a brand-new mixed formalism (mixed use of displacements and stresses as basic functions) is established to compensate the displacement formalism. In order to transfer all the related formulae and mathematical techniques of the Stroh formalism to these two formalisms, the general solutions for the basic equations of lamination theory and their associated eigenrelations have been purposely arranged in the form of Stroh formalism. Moreover, by using the presently developed mixed formalism, the explicit expressions for the fundamental matrix and eigenvectors are obtained first time for the most general composite laminates. Furthermore, letting the coupling stiffness vanish, the formalism has been reduced to the case of symmetric laminates and checked by a recently developed Stroh-like formalism for the plate bending problems. The comparison between Stroh formalism for two-dimensional problem, Stroh-like formalism for plate bending problem, displacement formalism and mixed formalism is then made at the end of this paper, and through their connection some useful relations are obtained.     

Higher order zig-zag theory for fully coupled thermo-electric–mechanical smart composite plates

A higher order zig-zag plate theory is developed to refine the prediction of the mechanical, thermal, and electric behaviors fully coupled. Both in-plane displacement and temperature fields through the thickness are constructed by superimposing linear zig-zag field to the smooth globally cubic varying field. Smooth parabolic distribution through the thickness is assumed in the out-of-plane displacement field in order to consider transverse normal deformation. Linear zig-zag form is adopted in the electric potential. The layer-dependent degrees of freedom of displacement and temperature fields are expressed in terms of reference primary degrees of freedom by applying interface continuity conditions as well as bounding surface conditions of transverse shear stresses and transverse heat fluxes. Thus the proposed theory is not only accurate but also efficient. Through the numerical examples of coupled and uncoupled analysis, the accuracy and efficiency of the present theory are demonstrated. The present theory is suitable in the predictions of fully coupled behaviors of thick smart composite plate under mechanical, thermal, and electric loads combined.     

A priori study for the modelling of velocity–interface correlations in the stratified air–water flows

This work concerns the modelling of stratified two-phase turbulent flows with interfaces. We consider an equation for an intermittency function $\alpha (\mathbf{x}\text{,}t)$ which denotes the probability of finding an interface at a given time t and a given point $\mathbf{x}$. In Wacławczyk and Oberlack (2011) a model for the unclosed terms in this equation was proposed. Here, we investigate the performance of this model by a priori tests, and finally, based on the a priori data discuss its possible modification and improvements.     

Differentially weighted direct simulation Monte Carlo method for particle collision in gas–solid flows

In gas–solid flows, particle–particle interaction (typical, particle collision) is highly significant, despite the small particles fractional volume. Widely distributed polydisperse particle population is a typical characteristic during dynamic evolution of particles (e.g., agglomeration and fragmentation) in spite of their initial monodisperse particle distribution. The conventional direct simulation Monte Carlo (DSMC) method for particle collision tracks equally weighted simulation particles, which results in high statistical noise for particle fields if there are insufficient simulation particles in less-populated regions. In this study, a new differentially weighted DSMC (DW-DSMC) method for collisions of particles with different number weight is proposed within the framework of the general Eulerian–Lagrangian models for hydrodynamics. Three schemes (mass, momentum and energy conservation) were developed to restore the numbers of simulation particle while keeping total mass, momentum or energy of the whole system unchanged respectively. A limiting case of high-inertia particle flow was numerically simulated to validate the DW-DSMC method in terms of computational precision and efficiency. The momentum conservation scheme which leads to little fluctuation around the mass and energy of the whole system performed best. Improved resolution in particle fields and dynamic behavior could be attained simultaneously using DW-DSMC, compared with the equally weighted DSMC. Meanwhile, computational cost can be largely reduced in contrast with direct numerical simulation.     

Lie symmetry analysis,soliton and numerical solutions of boundary value problem for variable coefficients coupled KdV–Burgers equation

Nonlinear Dynamics - In this article, an initial and boundary value problem for variable coefficients coupled KdV–Burgers equation is considered. With the help of Lie group approach, initial...     

Numerical evaluation of turbulence models for dense to dilute gas–solid flows in vertical conveyor

A two-fluid model (TFM) of multiphase flows based on the kinetic theory and small frictional limit boundary condition of granular flow was used to study the behavior of dense to dilute gas–solid flows in vertical pneumatic conveyor. An axisymmetric 2-dimensional, vertical pipe with 5.6 m length and 0.01 m internal diameter was chosen as the computation domain, same to that used for experimentation in the literature. The chosen particles are spherical, of diameter 1.91 mm and density 2500 kg/m³. Turbulence interaction between the gas and particle phases was investigated by Simonin's and Ahmadi's models and their numerical results were validated for dilute to dense conveying of particles. Flow regimes transition and pressure drop were predicted. Voidage and velocity profiles of each phase were calculated in radial direction at different lengths of the conveying pipe. It was found that the voidage has a minimum, and gas and solid velocities have maximum values along the center line of the conveying pipe and pressure drop has a minimum value in transition from dense slugging to dilute stable flow regime. Slug length and pressure fluctuation reduction were predicted with increasing gas velocity, too. It is shown that solid phase turbulence plays a significant role in numerical prediction of hydrodynamics of conveyor and the capability of particles turbulence models depends on tuning parameters of slip-wall boundary condition.     

Vector potential–vorticity relationship for the Stokes flows: application to the Stokes eigenmodes in 2D/3D closed domain

The unsteady dynamics of the Stokes flows, where , is shown to verify the vector potential–vorticity ( ) correlation , where the field is the pressure-gradient vector potential defined by . This correlation is analyzed for the Stokes eigenmodes, , subjected to no-slip boundary conditions on any two-dimensional (2D) closed contour or three-dimensional (3D) surface. It is established that an asymptotic linear relationship appears, verified in the core part of the domain, between the vector potential and vorticity, , where is a constant offset field, possibly zero.     

Variational methods for the steady state response of elastic–plastic solids subjected to cyclic loads

Solids (or structures) of elastic–plastic internal variable material models and subjected to cyclic loads are considered. A minimum net resistant power theorem, direct consequence of the classical maximum intrinsic dissipation theorem of plasticity theory, is envisioned which describes the material behavior by determining the plastic flow mechanism (if any) corresponding to a given stress/hardening state. A maximum principle is provided which characterizes the optimal initial stress/hardening state of a cyclically loaded structure as the one such that the plastic strain and kinematic internal variable increments produced over a cycle are kinematically admissible. A steady cycle minimum principle, integrated form of the aforementioned minimum net resistant power theorem, is provided, which characterizes the structure’s steady state response (steady cycle) and proves to be an extension to the present context of known principles of perfect plasticity. The optimality equations of this minimum principle are studied and two particular cases are considered: (i) loads not exceeding the shakedown limit (so recovering known results of shakedown theory) and (ii) specimen under uniform cyclic stress (or strain). Criteria to assess the structure’s ratchet limit loads are given. These, together with some insensitivity features of the structure’s alternating plasticity state, provide the basis to the ratchet limit load analysis problem, for which solution procedures are discussed.     

Finite amplitude oscillations of flanged laminas in viscous flows: Vortex–structure interactions for hydrodynamic damping control

In this paper, we study the problem of harmonic oscillations of a flanged lamina in a quiescent Newtonian incompressible viscous fluid. We conduct a comprehensive fluid–structure interaction investigation with the goal of assessing the effect of the presence of the flanges on the added mass and hydrodynamic damping experienced by the oscillating solid. We determine the complex nonlinear hydrodynamic function incorporating these effects via its real and imaginary parts, respectively, and its dependence on three nondimensional parameters that govern the flow evolution. We further investigate in detail the flow physics and the effects of nonlinearities on vortex shedding, convection, and diffusion in the vicinity of the oscillating structure. We find that the added mass effect is relatively independent of the oscillation amplitude and increases with the flange size. On the other hand, the hydrodynamic damping effect is remarkably affected by the interplay of geometry and dynamic parameters resulting into a peculiar non-monotonic behavior. We show the existence of a minimum in the hydrodynamic damping which can be attained via specific control of vortex–structure interaction dynamics and discuss its properties and significance from a physical perspective through analysis of the relevant flow fields. This novel finding has potential application for damping reduction in elastic systems where reduction of energy losses and increase of oscillation quality factor are desired.     

Performance of a combined three-hole conductivity probe for void fraction and velocity measurement in air–water flows

The development of a three-hole pressure probe with back-flushing combined with a conductivity probe, used for measuring simultaneously the magnitude and direction of the velocity vector in complex air–water flows, is described in this paper. The air–water flows envisaged in the current work are typically those occurring around the rotors of impulse hydraulic turbines (like the Pelton and Cross-Flow turbines), where the flow direction is not known prior to the data acquisition. The calibration of both the conductivity and three-hole pressure components of the combined probe in a rig built for the purpose, where the probe was placed in a position similar to that adopted for the flow measurements, will be reported. After concluding the calibration procedure, the probe was utilized in the outside region of a Cross-Flow turbine rotor. The experimental results obtained in the present study illustrate the satisfactory performance of the combined probe, and are encouraging toward its use for characterizing the velocity field of other complex air–water flows.     

Effect of boundary representation on viscous,separated flows in a discontinuous-Galerkin Navier–Stokes solver

The effect of curved-boundary representation on the physics of the separated flow over a NACA 65(1)-412 airfoil is thoroughly investigated. A method is presented to approximate curved boundaries with a high-order discontinuous-Galerkin spectral element method for the solution of the Navier–Stokes equations. Multiblock quadrilateral element meshes are constructed with the grid generation software GridPro. The boundary of a NACA 65(1)-412 airfoil, defined by a cubic natural spline, is piecewise-approximated by isoparametric polynomial interpolants that represent the edges of boundary-fitted elements. Direct numerical simulation of the airfoil is performed on a coarse mesh and fine mesh with polynomial orders ranging from four to twelve. The accuracy of the curve fitting is investigated by comparing the flows computed on curved-sided meshes with those given by straight-sided meshes. Straight-sided meshes yield irregular wakes, whereas curved-sided meshes produce a regular Karman street wake. Straight-sided meshes also produce lower lift and higher viscous drag as compared with curved-sided meshes. When the mesh is refined by reducing the sizes of the elements, the lift decrease and viscous drag increase are less pronounced. The differences in the aerodynamic performance between the straight-sided meshes and the curved-sided meshes are concluded to be the result of artificial surface roughness introduced by the piecewise-linear boundary approximation provided by the straight-sided meshes.