Simulating the fluid forces and fluid-elastic instabilities of a clamped–clamped cylinder in turbulent axial flow

In this paper, the fluid forces and the dynamics of a flexible clamped–clamped cylinder in turbulent axial flow are computed numerically. In the presented numerical model, there is no need to tune parameters for each specific case or to obtain coefficients from experiments. The results are compared with the dynamics measured in experiments available in the literature. The specific case studied here consists of a silicone cylinder mounted in axial water flow. Computationally it is found that the cylinder loses stability first by buckling. The threshold for buckling is in quantitative agreement with experimental results and weakly nonlinear theory. At higher flow speed a fluttering motion is predicted, in agreement with experimental results. It is also shown that even a small misalignment between the flow and the structure can have a significant impact on the dynamical behavior. To provide insight in the results of these fluid–structure interaction simulations, forces are computed on rigid inclined and curved cylinders, showing the existence of two different flow regimes. Furthermore it is shown that the inlet turbulence state has a non-negligible effect on these forces and thus on the dynamics of the cylinder.     

A theoretical model for gas–liquid slug flow in down inclined ducts

Analytical solutions are obtained for flows in downwardly inclined ducts, partly filled by a liquid and containing finite amplitude moving jumps. A unified theory for both roll waves and periodic slug flows in rounded ducts of arbitrary cross-section is worked out by means of some simplifications. The article is focused on slugs: a set of equations is obtained, which predicts the transition between roll waves and slug regimes and gives access to all flow characteristics without any need of closure laws concerning either the speed of propagation or the slug length. As a result, we gain a new insight on the physical structure of slug flow. The proposed model is valid for sufficient inclination, small pressure gradient along the duct and negligible superficial tension. Owing to assumptions, only main trends and orders of magnitude observed in experiments are to be checked. In this connection the model fits most of the previously published experimental results obtained in ducts of circular cross-section: the domain of occurrence of downwardly propagating slugs is satisfactorily predicted, the limitations in drift velocity and in liquid layer thickness are demonstrated and upwardly propagating slugs are possible.     

An approximate solution for the Couette–Poiseuille flow of the Giesekus model between parallel plates

An approximate analytical solution is derived for the Couette–Poiseuille flow of a nonlinear viscoelastic fluid obeying the Giesekus constitutive equation between parallel plates for the case where the upper plate moves at constant velocity, and the lower one is at rest. Validity of this approximation is examined by comparison to the exact solution during a parametric study. The influence of Deborah number (De) and Giesekus model parameter (α) on the velocity profile, normal stress, and friction factor are investigated. Results show strong effects of viscoelastic parameters on velocity profile and normal stress. In addition, five velocity profile types were obtained for different values of α, De, and the dimensionless pressure gradient (G).     

An improved anisotropy-resolving subgrid-scale model for flows in laminar–turbulent transition region

Some types of mixed subgrid-scale (SGS) models combining an isotropic eddy-viscosity model and a scale-similarity model can be used to effectively improve the accuracy of large eddy simulation (LES) in predicting wall turbulence. Abe (2013) has recently proposed a stabilized mixed model that maintains its computational stability through a unique procedure that prevents the energy transfer between the grid-scale (GS) and SGS components induced by the scale-similarity term. At the same time, since this model can successfully predict the anisotropy of the SGS stress, the predictive performance, particularly at coarse grid resolutions, is remarkably improved in comparison with other mixed models. However, since the stabilized anisotropy-resolving SGS model includes a transport equation of the SGS turbulence energy, k_SGS, containing a production term proportional to the square root of k_SGS, its applicability to flows with both laminar and turbulent regions is not so high. This is because such a production term causes k_SGS to self-reproduce. Consequently, the laminar–turbulent transition region predicted by this model depends on the inflow or initial condition of k_SGS. To resolve these issues, in the present study, the mixed-timescale (MTS) SGS model proposed by Inagaki et al. (2005) is introduced into the stabilized mixed model as the isotropic eddy-viscosity part and the production term in the k_SGS transport equation. In the MTS model, the SGS turbulence energy, k_es, estimated by filtering the instantaneous flow field is used. Since the k_es approaches zero by itself in the laminar flow region, the self-reproduction property brought about by using the conventional k_SGS transport equation model is eliminated in this modified model. Therefore, this modification is expected to enhance the applicability of the model to flows with both laminar and turbulent regions. The model performance is tested in plane channel flows with different Reynolds numbers and in a backward-facing step flow. The results demonstrate that the proposed model successfully predicts a parabolic velocity profile under laminar flow conditions and reduces the dependence on the grid resolution to the same degree as the unmodified model by Abe (2013) for turbulent flow conditions. Moreover, it is shown that the present model is effective at transitional Reynolds numbers. Furthermore, the present model successfully provides accurate results for the backward-facing step flow with various grid resolutions. Thus, the proposed model is considered to be a refined anisotropy-resolving SGS model applicable to laminar, transitional, and turbulent flows.     

Assessment of the kinetic–frictional model for dense granular flow

This paper aims to quantitatively assess the application of kinetic–frictional model to simulate the motion of dry granular materials in dense condition, in particular, the annular shearing in Couette configuration. The weight of frictional stress was varied to study the contribution of the frictional stress in dense granular flows. The results show that the pure kinetic-theory-based computational fluid dynamics (CFD) model (without frictional stress) over-predicts the dominant solids motion of dense granular flow while adding frictional stress [Schaeffer, D. G. (1987). Instability in the evolution equations describing incompressible granular flow. Journal of Differential Equations, 66(1), 19–50] with the solids pressure of [Lun, C. K. K., Savage, S. B., Jeffrey, D. J., & Chepurniy, N. (1984). Kinetic theories for granular flow: Inelastic particles in Couette flow and slightly inelastic particles in a general flow field. Journal of Fluid Mechanics, 140, 223–256] in the CFD model improves the simulation to better conform available experimental results. The results also suggest that frictional stress transmission plays an important role in dense granular flow and should not be neglected in granular flow simulations. Compatible simulation results to the experimental data are seen by increasing the weight of frictional stress to a factor of 1.25–1.5. These improved simulation results suggest the current constitutive relations (kinetic–frictional model) need to be improved in order to better reflect the real dense granular flow.     

Criticality of bifurcation in the tuned axial–torsional rotary drilling model

We study the bifurcation characteristics of a lumped-parameter model of rotary drilling with 1:1 internal resonance between the axial and the torsional modes which leads to the largest stability thresholds. For this special case, the two-degree-of-freedom model for the drill-string reduces to an effectively single-degree-of-freedom system facilitating further analysis. The regenerative effect of the cutting action due to the axial vibrations is incorporated through a delayed term in the cutting force with the delay depending on the torsional oscillations. This state dependency of the delay introduces nonlinearity in the current model. Steady drilling loses stability via a Hopf bifurcation, and the nature of the bifurcation is determined by an analytical study using the method of multiple scales. We find that both subcritical and supercritical Hopf bifurcations are present in this system depending on the choice of operating parameters. Hence, the nonlinearity due to the state-dependent delay term could both be stabilizing or destabilizing in nature, and the self-interruption nonlinearity is essential to capture the global behavior. Numerical bifurcation analysis of a global axial–torsional model of rotary drilling further confirms the analytical results from the method of multiple scales. Further exploration of the rotary drilling dynamics unravels more complex phenomena including grazing bifurcations and possibly chaotic solutions.     

The flow of magnetohydrodynamic Maxwell nanofluid over a cylinder with Cattaneo–Christov heat flux model

A theoretical analysis is performed for studying the flow and heat and mass transfer characteristics of Maxwell fluid over a cylinder with Cattaneo–Christov and non-uniform heat source/sink. The Brownian motion and thermophoresis parameters also considered into account. Numerical solutions are carried out by using Runge–Kutta-based shooting technique. The effects of various governing parameters on the flow and temperature profiles are demonstrated graphically. We also computed the friction factor coefficient, local Nusselt and Sherwood numbers for the permeable and impermeable flow over a cylinder cases. It is found that the rising values of Biot number, non-uniform heat source/sink and thermophoresis parameters reduce the rate of heat transfer. It is also found that the friction factor coefficient is high in impermeable flow over a cylinder case when compared with the permeable flow over a cylinder case.     

Spatiotemporal dynamics of a predator–prey model

Spatial component of ecological interactions has been identified as an important factor in how ecological communities are shaped. In this paper, we consider a Holling?CTanner model with spatial diffusion. Choosing appropriate parameter values in parameter spaces, we obtain rich patterns, including spotted, black-eye, and labyrinthine patterns. The numerical results show that predator?Cprey system can exhibit complicated behavior.     

An improved three-dimensional coupled fluid–structure model for Coriolis flowmeters

The paper presents a coupled numerical model built to simulate the operation of Coriolis flowmeters, which exploit the alteration of the vibration mode shape of the measuring tube for the mass flow rate measurement. The explained measuring effect is a consequence of the interaction between the motion of the tube, vibrating at its natural frequency, and the fluid flow in it. The numerical model is realized by coupling of a finite volume (FV) code for fluid flow analysis with a finite element (FE) code for structural analysis using the conventional staggered solution procedure, with added inner iterations to achieve strong coupling. The simulation algorithm is divided into two steps. A free vibration of the measuring tube considered in the first step is complemented in the second step, after the numerical free vibration response is properly stabilized, with the harmonic excitation force actuating the measuring tube at its resonant frequency of several hundreds of Hertz to resemble the operation of actual Coriolis flowmeters. Different scenarios using zero-order or three-point fluid load predictor and soft application of the fluid load in the initial stages of the simulation are compared to yield a simulation strategy, which will minimize the time needed to obtain the stabilized steady-state response of the vibrating measuring tube. The proposed simulation procedure was applied on a straight-tube Coriolis flowmeter and used for the estimation of the velocity profile effect. The results exhibit sufficient stability (low scatter) to be used for the estimation of sensitivity variations of order of magnitude around tenths of a percent.     

Validation of a flow–structure-interaction computation model of phonation

Computational models of vocal fold (VF) vibration are becoming increasingly sophisticated, their utility currently transiting from exploratory research to predictive research. However, validation of such models has remained largely qualitative, raising questions over their applicability to interpret clinical situations. In this paper, a computational model with a segregated implementation is detailed. The model is used to predict the fluid–structure interaction (FSI) observed in a physical replica of the VFs when it is excited by airflow. Detailed quantitative comparisons are provided between the computational model and the corresponding experiment. First, the flow model is separately validated in the absence of VF motion. Then, in the presence of flow-induced VF motion, comparisons are made of the flow pressure on the VF walls and of the resulting VF displacements. Self-similarity of spatial distributions of flow pressure and VF displacements is highlighted. The self-similarity leads to normalized pressure and displacement profiles. It is shown that by using linear superposition of average and fluctuation components of normalized computed displacements, it is possible to determine displacements in the physical VF replica over a range of VF vibration conditions. Mechanical stresses in the VF interior are related to the VF displacements, thereby the computational model can also determine VF stresses over a range of phonation conditions.     

Kinematics of Bulkley–Herschel fluid flow with a free surface during the filling of a channel

The non-Newtonian fluid flow with a free surface occurring during the filling of a plane channel in the gravity field is modeled. The mathematical formulation of the problem using the rheological Bulkley–Herschel model is presented. A numerical finite-difference algorithm for solving this problem is developed. A parametric investigation of the main characteristics of the process as functions of the control parameters is performed. The effect of the rheological parameters of the fluid on the distribution of the quasisolid motion zones is demonstrated.     

Rich spatial–temporal dynamics in a diffusive population model for pioneer–climax species

Nonlinear Dynamics - A general diffusive population model for interactions of pioneer and climax species subject to the no-flux boundary condition is considered. Local and global steady-state...     

A coupled model for simulation of the gas–liquid two-phase flow with complex flow patterns

A new model coupling two basic models, the model based on interface tracking method and the two-fluid model, for simulating gas–liquid two-phase flow is presented. The new model can be used to simulate complex multiphase flow in which both large-length-scale interface and small-length-scale gas–liquid interface coexist. By the physical state and the length scale of interface, three phases are divided, including the liquid phase, the large-length-scale-interface phase (LSI phase) and the small-length-scale-interface phase (SSI phase). A unified solution framework shared by the two basic models is built, which makes it convenient to perform the solution process. Based on the unified solution framework, the modified MCBA–SIMPLE algorithm is employed to solve the Navier–Stokes equations for the proposed model. A special treatment called “volume fraction redistribution” is adopted for the special grids containing all three phases. Another treatment is proposed for the advection of large-length-scale interface when some portion of SSI phase coalesces into LSI phase. The movement of the large-length-scale interface is evaluated using VOF/PLIC method. The proposed model is equivalent to the two-fluid model in the zone where only the liquid phase and the SSI phase are present and to the model based on interface tracking method in the zone where only the liquid phase and the LSI phase are present. The characteristics of the proposed model are shown by four problems.     

Spatio-temporal dynamics of a reaction-diffusion system for a predator–prey model with hyperbolic mortality

We investigate the effects of diffusion on the spatial dynamics of a predator–prey model with hyperbolic mortality in predator population. More precisely, we aim to study the formation of some elementary two-dimensional patterns such as hexagonal spots and stripe patterns. Based on the linear stability analysis, we first identify the region of parameters in which Turing instability occurs. When control parameter is in the Turing space, we analyse the existence of stable patterns for the excited model by the amplitude equations. Then, for control parameter away from the Turing space, we numerically investigate the initial value-controlled patterns. Our results will enrich the pattern dynamics in predator–prey models and provide a deep insight into the dynamics of predator–prey interactions.     

Deformations of a clamped–clamped elastica inside a circular channel with clearance

In this paper, we study the planar deformations of an elastica inside a circular channel with clearance. One end of the elastica is fully clamped, while the other end is partially clamped in the lateral direction and is subject to a pushing force longitudinally. In the experiment we first observe various deformation patterns after pushing the elastica through the partial clamp. Both symmetric and asymmetric deformations are recorded. Special attention is focused on the contact conditions between the elastica and the circular channel. In order to analyze the elastica deformation theoretically, we first divide the elastica into several elementary sub-domains depending on the contact condition between the elastica and the circular channel. In each sub-domain the elastica is either loaded only at the ends or in full contact with the outer wall. Armed with these basic equilibrium analyses, we proceed to calculate and classify the loaded elastica into several deformation patterns. Finally, we present the load-deflection curves, both theoretically and experimentally, which relate the longitudinal forces at both ends to the elastica length increase inside the channel. The branching phenomena predicted theoretically agree fairly well quantitatively with the experimental measurements.     

Hidden dynamics in a fractional-order memristive Hindmarsh–Rose model

Nonlinear Dynamics - By showing that there is no any memristor in the integer-order Hindmarsh–Rose (H–R) model, the slow ion channel is remodeled by a fractional-order memristor, where...     

Stress boundary layers in the viscoelastic flow past a cylinder in a channel： limiting solutions

The flow past a cylinder in a channel with the aspect ratio of 2:1 for the upper convected Maxwell (UCM) fluid and the Oldroyd-B fluid with the viscosity ratio of 0.59 is studied by using the Galerkin/Least-square finite element method and a p-adaptive refinement algorithm. A posteriori error estimation indicates that the stress-gradient error dominates the total error. As the Deborah number, D_e, approaches 0.8 for the UCM fluid and 0.9 for the Oldroyd-B fluid, strong stress boundary layers near the rear stagnation point are forming, which are characterized by jumps of the stress-profiles on the cylinder wall and plane of symmetry, huge stress gradients and rapid decay of the gradients across narrow thicknesses. The origin of the huge stress-gradients can be traced to the purely elongational flow behind the rear stagnation point, where the position at which the elongation rate is of 1/2D_e approaches the rear stagnation point as the Deborah number approaches the critical values. These observations imply that the cylinder problem for the UCM and Oldroyd-B fluids may have physical limiting Deborah numbers of 0.8 and 0.9, respectively.The project supported by the National Natural Science Foundation of China (50335010 and 20274041) and the MOLDFLOW Comp. Australia.