Pure axial flow of viscoelastic fluids in rectangular microchannels under combined effects of electro-osmosis and hydrodynamics

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien–Tanner (PTT) model with the Gordon–Schowalter convected derivative which covers the upper convected Maxwell, Johnson–Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson–Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid system and is verified against 1D analytical solution of the velocity profile with less than 0.06% relative error. Also, a parametric study is carried out to investigate the effect of channel aspect ratio (width to height), wall zeta potential and the Debye–Hückel parameter on 2D velocity profile, volumetric flow rate and the Poiseuille number in the mixed EO/PD flows of viscoelastic fluids with different Weissenberg numbers. Our results show that, for low channel aspect ratios, the previous 1D analytical models underestimate the velocity profile at the channel half-width centerline in the case of favorable pressure gradients and overestimate it in the case of adverse pressure gradients. The results reveal that the inapplicability of the Debye–Hückel approximation at high zeta potentials is more significant for higher Weissenberg number fluids. Also, it is found that, under the specified values of electrokinetic parameters, there is a threshold for velocity scale ratio in which the Poiseuille number is approximately independent of channel aspect ratio.     

Analytical solutions for tapering quadratic and cubic rat-holes in highly frictional granular solids

New exact analytical solutions are presented for both stress and velocity fields for a Coulomb–Mohr granular solid assuming non-dilatant double-shearing theory. The solutions determined apply to highly frictional materials for which the angle of internal friction φ is assumed equal to 90°. This major assumption is made primarily to facilitate exact analytical solutions, and it is discussed at length in the Introduction, both in the context of real materials which exhibit large angles of internal friction, and in the context of using the solutions derived here as the leading term in a regular perturbation solution involving powers of 1−sinφ. The analytical velocity fields so obtained are illustrated graphically by showing the direction of the principal stress as compared to the streamlines. The stress solutions are also exploited to determine the static stress distribution for a granular material contained within vertical boundaries and a horizontal base, which is assumed to have an infinitesimal central outlet through which material flows until a rat-hole of parabolic or cubic profile is obtained, and no further flow takes place. A rat-hole is a stable structure that may form in storage hoppers and stock-piles, preventing any further flow of material. Here we consider the important problems of two-dimensional parabolic rat-holes of profile y=ax², and three-dimensional cubic rat-holes of profile z=ar³, which are both physically realistic in practice. Analytical solutions are presented for both two and three-dimensional rat-holes for the case of a highly frictional granular solid, which is stored at rest between vertical walls and a horizontal rigid plane, and which has an infinitesimal central outlet. These solutions are bona fide exact solutions of the governing equations for a Coulomb–Mohr granular solid, and satisfy exactly the free surface condition along the rat-hole surface, but approximate frictional conditions along the containing boundaries. The analytical solutions presented here constitute the only known solutions for any realistic rat-hole geometry, other than the classical solution which applies to a perfectly vertical cylindrical cavity.     

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.     

A note on start-up and large amplitude oscillatory shear flow of multimode viscoelastic fluids

Start up of plane Couette flow and large amplitude oscillatory shear flow of single and multimode Maxwell fluids as well as Oldroyd-B fluids have been analyzed by analytical or semi-analytical procedures. The result of our analysis indicates that if a single or a multimode Maxwell fluid has a relaxation time comparable or smaller than the rate of change of force imparted on the fluid, then the fluid response is not singular as Elasticity Number (E ). However, if this is not the case, as E , perturbations of single and multimode Maxwell fluids give rise to highly oscillatory velocity and stress fields. Hence, their behavior is singular in this limit. Moreover, we have observed that transients in velocity and stresses that are caused by propagation of shear waves in Maxwell fluids are damped much more quickly in the presence of faster and faster relaxing modes. In addition, we have shown that the Oldroyd-B model gives rise to results quantitatively similar to multimode Maxwell fluids at times larger than the fastest relaxation time of the multimode Maxwell fluid. This suggests that the effect of fast relaxing modes is equivalent to viscous effects at times larger than the fastest relaxation time of the fluid. Moreover, the analysis of shear wave propagation in multimode Maxwell fluids clearly show that the dynamics of wave propagation are governed by an effective relaxation and viscosity spectra. Finally, no quasi-periodic or chaotic flows were observed as a result of interaction of shear waves in large amplitude oscillatory shear flows for any combination of frequency and amplitudes.     

Criteria of negative wake generation behind a cylinder

It was demonstrated by simulation in our previous study that both the normal stress and its gradient are responsible for the negative wake generation (overshoot in the axial velocity) and streamline shifting. Extensional properties of the fluids dominate the generation of the negative wake, while other factors strengthen or weaken the formation of velocity overshoot. In this study, the criteria for the negative wake generation are discussed in detail for various fluid models, including the PTT, the FENE-CR, the FENE-P, and the Giesekus models. With the FENE-CR fluid, it is easier to generate negative wake than with the FENE-P fluid. This confirms that the constant shear viscosity FENE-CR fluid enhances the velocity overshoot, and that the shear-thinning viscosity FENE-P fluid delays the negative wake generation. The Giesekus fluid has a similar behaviour to the PTT fluid with regarding to the critical conditions of negative wake generation when appropriate fluid parameters are selected. The mechanism of wall proximity in enhancing the negative wake generation is also demonstrated with the analysis for the first time.     

The exact solution of stress distribution in fillet welds

I.Intr0ducti0nInanearlierpaper['],theexactsolutionofstressdistributioninweldsundertheactionofconcentratedforcePwasobtained.Thedesignmethoduptonowaretheutilizationofaveragestressactingontherelatingthroatsectionofweldsandthestressdistribution0btainedbythisa…     

A semi-analytical solution for linearized multicomponent cation exchange reactive transport in groundwater

Cation exchange in groundwater is one of the dominant surface reactions. Mass transfer of cation exchanging pollutants in groundwater is highly nonlinear due to the complex nonlinearities of exchange isotherms. This makes difficult to derive analytical solutions for transport equations. Available analytical solutions are valid only for binary cation exchange transport in 1-D and often disregard dispersion. Here we present a semi-analytical solution for linearized multication exchange reactive transport in steady 1-, 2- or 3-D groundwater flow. Nonlinear cation exchange mass–action–law equations are first linearized by means of a first-order Taylor expansion of log-concentrations around some selected reference concentrations and then substituted into transport equations. The resulting set of coupled partial differential equations (PDEs) are decoupled by means of a matrix similarity transformation which is applied also to boundary and initial concentrations. Uncoupled PDE’s are solved by standard analytical solutions. Concentrations of the original problem are obtained by back-transforming the solution of uncoupled PDEs. The semi-analytical solution compares well with nonlinear numerical solutions computed with a reactive transport code (CORE^2D) for several 1-D test cases involving two and three cations having moderate retardation factors. Deviations of the semi-analytical solution from numerical solutions increase with increasing cation exchange capacity (CEC), but do not depend on Peclet number. The semi-analytical solution captures the fronts of ternary systems in an approximate manner and tends to oversmooth sharp fronts for large retardation factors. The semi-analytical solution performs better with reference concentrations equal to the arithmetic average of boundary and initial concentrations than it does with reference concentrations derived from the arithmetic average of log-concentrations of boundary and initial waters.     

Exact integration of the von Mises elastoplasticity model with combined linear isotropic-kinematic hardening

The main purpose of this work is to present two semi-analytical solutions for the von Mises elastoplasticity model governed by combined linear isotropic-kinematic hardening. The first solution (SOLε) corresponds to strain-driven problems with constant strain rate assumption, whereas the second one (SOLσ) is proposed for stress-driven problems using constant stress rate assumption. The formulas are derived within the small strain theory Besides the new analytical solutions, a new discretized integration scheme (AMε) based on the time-continuous SOLε is also presented and the corresponding algorithmically consistent tangent tensor is provided. A main advantage of the discretized stress updating algorithm is its accuracy; it renders the exact solution if constant strain rate is assumed during the strain increment, which is a commonly adopted assumption in the standard finite element calculations. The improved accuracy of the new method (AMε) compared with the well-known radial return method (RRM) is demonstrated by evaluating two simple examples characterized by generic nonlinear strain paths.     

Functional and generalized separable solutions to unsteady Navier–Stokes equations

The paper studies unsteady Navier–Stokes equations with two space variables. It shows that the non-linear fourth-order equation for the stream function with three independent variables admits functional separable solutions described by a system of three partial differential equations with two independent variables. The system is found to have a number of exact solutions, which generate new classes of exact solutions to the Navier–Stokes equations. All these solutions involve two or more arbitrary functions of a single argument as well as a few free parameters. Many of the solutions are expressed in terms of elementary functions, provided that the arbitrary functions are also elementary; such solutions, having relatively simple form and presenting significant arbitrariness, can be especially useful for solving certain model problems and testing numerical and approximate analytical hydrodynamic methods. The paper uses the obtained results to describe some model unsteady flows of viscous incompressible fluids, including flows through a strip with permeable walls, flows through a strip with extrusion at the boundaries, flows onto a shrinking plane, and others. Some blow-up modes, which correspond to singular solutions, are discussed.     

Analysis of functionally graded rotating disks with variable thickness

Elastic solutions for axisymmetric rotating disks made of functionally graded material with variable thickness are presented. The material properties and disk thickness profile are assumed to be represented by two power-law distributions. In the case of hollow disk, based on the form of the power-law distribution for the mechanical properties of the constituent components and the thickness profile function, both analytical and semi-analytical solutions are given under free–free and fixed-free boundary conditions. For the solid disk, only semi-analytical solution is presented. The effects of the material grading index and the geometry of the disk on the stresses and displacements are investigated. It is found that a functionally graded rotating disk with parabolic or hyperbolic convergent thickness profile has smaller stresses and displacements compared with that of uniform thickness. It is seen that the maximum radial stress for the solid functionally graded disk with parabolic thickness profile is not at the centre like uniform thickness disk. Results of this paper suggest that a rotating functionally graded disk with parabolic concave or hyperbolic convergent thickness profile can be more efficient than the one with uniform thickness.     

Transient wave propagation of isotropic plates using a higher-order plate theory

Transient wave propagation of isotropic thin plates using a higher-order plate theory is presented in this paper. The aim of this investigation is to assess the applicability of the higher-order plate theory in describing wave behavior of isotropic plates at higher frequencies. Both extensional and flexural waves are considered. A complete discussion of dispersion of isotropic plates is first investigated. All the wave modes and wave behavior for each mode in the low and high-frequency ranges are provided in detail. Using the dispersion relation and integral transforms, exact integral solutions for an isotropic plate subjected to pure impulse load and a number of wave excitations based on the higher-order theory are obtained and asymptotic solutions which highlight the physics of waves are also presented. The axisymmetric three-dimensional analytical solutions of linear wave equations are also presented for comparison. Results show that the higher-order theory can predict the wave behavior closely with exact linear wave solutions at higher frequencies.     

An exact analytical solution for creeping Dean flow of Bingham plastics through curved rectangular ducts

In this paper, an exact analytical solution for creeping flow of Bingham plastic fluid passing through curved rectangular ducts is presented for the first time. The closed form of axial velocity distribution, flow resistance ratio, and wall shear stress are derived using bounded Fourier transformation. An extensive investigation on mutual effects of Hedstrom number, curvature ratio, and aspect ratio is conducted. The results indicate that a drag reduction is caused in the flow field by increasing the Hedstrom number. It is shown that unlike the Newtonian creeping Dean flow, the critical aspect ratio (an aspect ratio in which the flow resistance ratio is independent from curvature ratio) does not exist at large enough Hedstrom numbers. Analytical solution also indicated that as Hedstrom number is increased, the value of Poiseuille number is enhanced, and unlike the Newtonian flows, the value of Poiseuille number is not zero at edges of cross section.     

An analytical solution for the laminar flow over a flat plate with similarity preserving suction

An exact analytical solution is presented for the laminar boundary-layer flow over a semi-infinite flat plate subjected to a type of similarity preserving suction. The solution is developed for the case of a plate immersed in either a uniform compressible stream with viscosity proportional to temperature or a uniform incompressible stream with constant viscosity. The problem is formulated in Crocco's variables. It is described by a second-order, non-linear, ordinary differential (and singular) boundary-value problem for the shear stress as a function of the velocity in the boundary layer. A unique solution is shown to exist and to possess a power series representation for all magnitudes of suction. The series is constructed explicitly and provides a transcendental equation for the shear stress at the plate (the important skin friction) which can be solved to any desired accuracy. Examples of upper and lower bounds for the wall shear are presented for several magnitudes of suction and confirm the reasonable accuracy of results obtained heretofore only by numerical solutions of the problem. In addition to the intrinsic value of the technique developed, it can be the basis of accurate checks for the numerical solution of more complex problems.     

Practical comparison of differential viscoelastic constitutive equations in finite element analysis of planar 4:1 contraction flow

Finite element modeling of planar 4:1 contraction flow (isothermal incompressible and creeping) around a sharp entrance corner is performed for favored differential constitutive equations such as the Maxwell, Leonov, Giesekus, FENE-P, Larson, White-Metzner models and the Phan Thien-Tanner model of exponential and linear types. We have implemented the discrete elastic viscous stress splitting and streamline upwinding algorithms in the basic computational scheme in order to augment stability at high flow rate. For each constitutive model, we have obtained the upper limit of the Deborah number under which numerical convergence is guaranteed. All the computational results are analyzed according to consequences of mathematical analyses for constitutive equations from the viewpoint of stability. It is verified that in general the constitutive equations proven globally stable yield convergent numerical solutions for higher Deborah number flows. Therefore one can get solutions for relatively high Deborah number flows when the Leonov, the Phan Thien-Tanner, or the Giesekus constitutive equation is employed as the viscoelastic field equation. The close relationship of numerical convergence with mathematical stability of the model equations is also clearly demonstrated.     

Computation of flow of viscoelastic fluids by parameter differentiation

A technique combining the features of parameter differentiation and finite differences is presented to compute the flow of viscoelastic fluids. Two flow problems are considered: (i) three-dimensional flow near a stagnation point and (ii) axisymmetric flow due to stretching of a sheet. Both flows are characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. The exact numerical solutions are obtained using the technique described in the paper. Also, the first-order perturbation solutions (in terms of the viscoelastic fluid parameter) are derived. A comparison of the results shows that the perturbation method is inadequate in predicting some of the vital characteristic features of the flows, which can possibly be revealed only by the exact numerical solution.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

Semi-analytical solution to one-dimensional consolidation in unsaturated soils

In this paper, a series of semi-analytical solutions to one-dimensional consolidation in unsaturated soils are obtained. The air governing equation by Fredlund for unsaturated soils consolidation is simplified. By applying the Laplace transform and the Cayley-Hamilton theorem to the simplified governing equations of water and air, Darcy＇s law, and Fick＇s law, the transfer function between the state vectors at top and at any depth is then constructed. Finally, by the boundary conditions, the excess pore-water pressure, the excess pore-air pressure, and the soil settlement are obtained under several kinds of boundary conditions with the large-area uniform instantaneous loading. By the Crump method, the inverse Laplace transform is performed, and the semi-analytical solutions to the excess pore-water pressure, the excess pore-air pressure, and the soils settlement are obtained in the time domain. In the case of one surface which is permeable to air and water, comparisons between the semi-analytical solutions and the analytical solutions indicate that the semi-analytical solutions are correct. In the case of one surface which is permeable to air but impermeable to water, comparisons between the semi-analytical solutions and the results of the finite difference method are made, indicating that the semi-analytical solution is also correct.     

Start-up flows of second grade fluids in domains with one finite dimension

A number of unidirectional transient flows of a second grade fluid in a domain with one finite dimension are studied. The method of integral transforms (Fourier, Hankel or Laplace) is applied to obtain exact solutions. A general theorem on start-up flows for second grade fluids is presented that allows us to determine unidirectional flows of second grade fluids once the corresponding solution is known within the context of the Navier-Stokes theory. In the process of obtaining solutions for the fluid of second grade, we find several new exact solutions within the context of the classical Navier-Stokes theory.     

Convergence and exact solutions of spline finite strip method using unitary transformation approach

The spline finite strip method(PSM) is one of the most popular numerical methods for analyzing prismatic structures.Efficacy and convergence of the method have been demonstrated in previous studies by comparing only numerical results with analytical results of some benchmark problems.To date,no exact solutions of the method or its explicit forms of error terms have been derived to show its convergence analytically. As such,in this paper,the mathematical exact solutions of spline finite strips in the plat...     

A new finite element scheme using the Lagrangian framework for simulation of viscoelastic fluid flows

A new numerical scheme for simulation of viscoelastic fluid flows was designed, making use of finite element algorithms generally regarded as advantageous for tackling the problem. This includes the Lagrangian approach for the solution of viscoelastic constitutive equation using the co-deformational frame of reference with a possibility of analytically solving the equation along the particles trajectories, which in turn allowed eluding the solution of any system of linear equations for the stress. Then, the full ellipticity of the momentum conservation equation was utilised thanks to a possibility of accurate determination of the stress tensor independently of the velocity field at the current stage of computation. The needed independent stress was calculated at each time step on the basis of the past deformation history, which in turn was determined on the basis of the past velocity fields, all incorporated into a modified Euler time stepping algorithm. Owing to explicit inclusion of the full viscous term from the viscoelastic model into the momentum conservation equation, no stress splitting was necessary. The trajectory feet tracking was done accurately using a semi-analytic solution of the displacement gradient evolution equation and a weak formulation of the kinematics equation, the latter at the expense of solving an extra symmetric system of linear equations.The error expressed in the form of the Sobolev norms was determined using a comparison with available analytical solution for UCM fluid in the transient regime or numerically obtained steady-state stress values for the PTT fluid in Couette flow. The implementation of the PTT fluid model was done by modifying the relative displacement gradient tensor so that a new convective frame was defined.The stability of the algorithm was assessed using the well-known benchmark problem of a sphere sedimenting in a tube with viscoelastic fluid. The stable numerical results were obtained at high Weissenberg numbers, with the limit of convergence Wi=6.6, exceeding any previously reported values. The robustness of the code was proven by simulation of the Weissenberg effect (the rod-climbing phenomenon) with the use of PTT fluid.