Boundary-layer non-Newtonian flow over vertical plate in porous medium saturated with nanofluid

The free convective heat transfer to the power-law non-Newtonian flow from a vertical plate in a porous medium saturated with nanofluid under laminar conditions is investigated. It is considered that the non-Newtonian nanofluid obeys the mathematical model of power-law. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis. The partial differential system governing the problem is transformed into an ordinary system via a usual similarity transformation. The numerical solutions of the resulting ordinary system are obtained. These solutions depend on the power-law index n, Lewis number Le, buoyancy-ratio number N _r, Brownian motion number N _b, and thermophoresis number N _t. For various values of n and Le, the effects of the influence parameters on the fluid behavior as well as the reduced Nusselt number are presented and discussed.     

On the thermodynamics of some generalized second-grade fluids

The generalized second-grade fluids, which have been used for modeling the creep of ice and the flow of coal-water and coal-oil slurries, are among the simplest non-Newtonian fluid models that can describe shear-thinning/thickening and exhibit normal stress effects. In this article, we conduct thermodynamic analysis on a class of generalized second-grade fluids, one distinguishing feature of which is the existence of a constitutive function Φ that describes frictional heating. We work within the framework of Serrin’s original formulation of neoclassical thermodynamics, where internal energy and entropy functions, if they exist for a continuous body at all, are to be derived from the classical First Law and (quantitatively reformulated) Second Law of thermodynamics for cycles. For the class of generalized second-grade fluids in question, we show from the First Law that an internal energy density u exists, and we derive the equation of energy balance; from the Second Law, we demonstrate the existence of an entropy density s and derive the Clausius–Duhem inequality that it satisfies. We obtain explicit expressions for u, s and the frictional heating Φ, and derive thermodynamic restrictions on the material functions of temperature μ, α ₁, and α ₂ that appear in the constitutive relation for the Cauchy stress. For the special case of second-grade fluids, our expressions for u and s agree with those which Dunn and Fosdick [6] derived under the theoretical framework of the rational thermodynamics of Coleman and Noll.     

Exact solutions for extended boundary value problems

Summary Extended definition of a stress tensor for a non-Newtonian fluid brings in higher degree derivatives with coefficients as powers of non-Newtonian parameter in the differential equations of motion. Yet, these differential equations need to be solved subject to the same boundary conditions as in the corresponding Newtonian flow problem. A technique is developed to obtain exact solutions for such an extended boundary value problem. Some flow problems forWalters liquidB are considered.     

Computation of confined coflow jets with three turbulence models

A numerical study of confined jets in a cylindrical duct is carried out to examine the performance of two recently proposed turbulence models: an RNG-based K-? model and a realizable Reynolds stress algebraic equation model. The former is of the same form as the standard K-? model but has different model coefficients. The latter uses an explicit quadratic stress-strain relationship to model the turbulent stresses and is capable of ensuring the positivity of each turbulent normal stress. The flow considered involves recirculation with unfixed separation and reatachment points and severe adverse pressure gradients, thereby providing a valuable test of the predictive capability of the models for complex flows. Calculations are performed with a finite volume procedure. Numerical credibility of the solutions is ensured by using second-order-accurate differencing schemes and sufficiently fine grids. Calculations with the standard K-? model are also made for comparison. Detailed comparisons with experiments show that the realizable Reynolds stress algebraic equation model consistently works better than does the standard K-? model in capturing the essential flow features, while the RNG-based K-? model does not seem to give improvements over the standard K-? model under the flow conditions considered.     

A remark on similarity analysis of boundary layer equations of a class of non-Newtonian fluids

A similarity analysis of three-dimensional boundary layer equations of a class of non-Newtonian fluid in which the stress, an arbitrary function of rates of strain, is studied. It is shown that under any group of transformation, for an arbitrary stress function, not all non-Newtonian fluids possess a similarity solution for the flow past a wedge inclined at arbitrary angle except Ostwald-de-Waele power-law fluid. Further it is observed, for non-Newtonian fluids of any model only $90°$ of wedge flow leads to similarity solutions. Our results contain a correction to some flaws in Pakdemirli׳s [14] (1994) paper on similarity analysis of boundary layer equations of a class of non-Newtonian fluids.     

Phenomenological Modeling of Electrorheological Fluids with an Extended Casson-Model

In this paper we propose a phenomenological theory for electrorheological fluids. In general these are suspensions which undergo dramatic changes in their material properties if they are exposed to an electric field. In the context of continuum mechanics these fluids can be modeled as non-Newtonian fluids. Recalling the governing equations of rational thermodynamics and electrodynamics of moving media (Maxwell-Minkowski-equations), we derive suitable governing equations of electrorheology using essentially two assumptions concerning magnetic quantities. Furthermore we introduce a 3-dimensional nonlinear constitutive equation for the Cauchy stress tensor which is an extension of the model proposed by Ružička (see [14]). Assuming a viscometric flow, we compare the shear stress of our model with other well known models and fit the parameters by using measurements that were obtained in a rotational viscometer. Excellent agreement between model and measurements is achieved. On the basis of these results we propose a 3-dimensional model, the so-called extended Casson -model. This model is investigated further for a channel flow configuration with a homogeneous electric field. We determine analytical solutions for the electric field, the velocity and the volumetric flow rate and illustrate the velocity profiles and the predicted pressure drop. The velocity profiles are flattened compared to parabolic profiles and become more flat if the electric field increases. Received March 21, 2000     

An extended thermodynamic approach to non-Newtonian fluids and related results in Marangoni instability problem

It is shown that extended irreversible thermodynamics provides a simple and coherent modelling of non-Newtonian fluids. The basic hypothesis underlyig the present formalism is to raise the stress tensor to the status of independent variable. Restrictions are placed on the constitutive equations and the material coefficients by the second law of thermodynamics, stability of equilibrium and the objectivity principle. The general procedure is applied to derive the Reiner-Rivlin and Rivlin-Ericksen second-order fluids. As an illustration, Marangoni convection in a thin horizontal layer of a non-Newtonian fluid submitted to a temperature gardient and placed in a microgravity environment is treated.     

Weibull modulus for diverse strength due to sample-specificity

By sample specificity it is meant that specimens with the same nominal material parameters and tested under the same environmental conditions may exhibit different behavior with diversified strength. Such an effect has been widely observed in the testing of material failure and is usually attributed to the heterogeneity of material at the mesoscopic level. The degree with which mesoscopic heterogeneity affects macroscopic failure is still not clear. Recently, the problem has been examined by making use of statistical ensemble evolution of dynamical system and the mesoscopic stress re-distribution model (SRD). Sample specificity was observed for non-global mean stress field models, such as the cluster mean field model, stress concentration at tip of microdamage, etc. Certain heterogeneity of microdamage could be sensitive to particular SRD leading to domino type of coalescence. Such an effect could start from the microdamage heterogeneity and then be magnified to other scale levels. This trans-scale sensitivity is the origin of sample specificity. The sample specificity leads to a failure probability Φ_N with a transitional region 0<Φ_N<1, so that globally stable and catastrophic modes could co-exist. It is found that the scatter in strength can fit the Weibull distribution very well. Hence, the Weibull modulus is indicative of sample specificity. Numerical results obtained from the SRD for different non-global mean stress fields show that Weibull modulus increases with increasing sample span and influence region of microdamage.     

A Theoretical Model of Hysteresis and Dynamic Effects in the Capillary Relation for Two-phase Flow in Porous Media

It is well known that the relationship between capillary pressure and saturation, in two-phase flow problems demonstrates memory effects and, in particular, hysteresis. Explicit representation of full hysteresis with a myriad of scanning curves in models of multiphase flow has been a difficult problem. A second complication relates to the fact that P ^c–S relationships, determined under static conditions, are not necessarily valid in dynamics. There exist P ^c–S relationships which take into account dynamic effects. But the combination of hysteretic and dynamic effects in the capillary relationship has not been considered yet. In this paper, we have developed new models of capillary hysteresis which also include dynamic effects. In doing so, thermodynamic considerations are employed to ensure the admissibility of the new relationships. The simplest model is constructed around main imbibition and drainage curves and assumes that all scanning curves are vertical lines. The dynamic effect is taken into account by introducing a damping coefficient in P ^c–S equation. A second-order model of hysteresis with inclined scanning curves is also developed. The simplest version of proposed models is applied to two-phase incompressible flow and an example problem is solved.     

Evaluation of advanced anisotropic models with mixed hardening for general associated and non-associated flow metal plasticity

The main objective of this paper is to develop a generalized finite element formulation of stress integration method for non-quadratic yield functions and potentials with mixed nonlinear hardening under non-associated flow rule. Different approaches to analyze the anisotropic behavior of sheet materials were compared in this paper. The first model was based on a non-associated formulation with both quadratic yield and potential functions in the form of Hill’s (1948). The anisotropy coefficients in the yield and potential functions were determined from the yield stresses and r-values in different orientations, respectively. The second model was an associated non-quadratic model (Yld2000-2d) proposed by Barlat et al. (2003). The anisotropy in this model was introduced by using two linear transformations on the stress tensor. The third model was a non-quadratic non-associated model in which the yield function was defined based on Yld91 proposed by Barlat et al. (1991) and the potential function was defined based on Yld89 proposed by Barlat and Lian (1989). Anisotropy coefficients of Yld91 and Yld89 functions were determined by yield stresses and r-values, respectively. The formulations for the three models were derived for the mixed isotropic-nonlinear kinematic hardening framework that is more suitable for cyclic loadings (though it can easily be derived for pure isotropic hardening). After developing a general non-associated mixed hardening numerical stress integration algorithm based on backward-Euler method, all models were implemented in the commercial finite element code ABAQUS as user-defined material subroutines. Different sheet metal forming simulations were performed with these anisotropic models: cup drawing processes and springback of channel draw processes with different drawbead penetrations. The earing profiles and the springback results obtained from simulations with the three different models were compared with experimental results, while the computational costs were compared. Also, in-plane cyclic tension–compression tests for the extraction of the mixed hardening parameters used in the springback simulations were performed for two sheet materials.     

General relations of strain-gradient stress analysis

Some particular fields of stress gradients in plates are investigated analytically and experimentally. Carriers of empirical information are light beam deflections caused by stress or strain gradients in homogeneous beams subjected to a particular case of flexure with shear. This study is based on theories and techniques of the strain-gradient method that was recently introduced by the authors. This is a generalization of prior analytical-experimental examination of strain-gradient light deflections produced in stressed plates, which had concentrated on the simplest case where information of interest is collected along a line of symmetry of the stress field and where both the information carrying light beams are deflected in the plane of symmetry only. The developed relations were tested experimentally, by using an S-beam as representative example of general plane stress field. The main purpose of the present investigation is to document the efficacy of the strain-gradient method in analysis of the general case of stress state and to test the ranges of applicability of the accepted mathematical models and of the subsequently derived relations. In this respect, the most interesting stress state is that in a beam subjected to the Saint-Venant bending, where the transversal and longitudinal axes of the beam are in pure shear. The obtained results are compared, at each step, with the predictions of the developed analytical models and with the predictions of Filon's stress function. The results of comparison are satisfactory. The procedures of evaluating the photoelastic and material coefficients using strain-gradient techniques were tested positively. The developed method can yield valuable information on the actual features of stress states in fracture mechanics. Part of the problems and first results of this investigation will be presented in the authors' paper “Strain-gradient stress analysis in Saint-Venant bending”.     

Theoretical prediction of the mass transfer coefficients in bubble columns operating in churn-turbulent flow regime. Study in Newtonian and non-Newtonian fluids under different operation conditions

 A comprehensive experimental study of the volumetric transfer coefficient k _L a with Newtonian and non-Newtonian fluids in bubble columns using CO₂ as gas phase is the objective of this work. The evaluation of the hydrodynamic characteristics of the bubble columns and delineated the different hydrodynamic regimes considering column geometry, gas flow, liquid height and type of fluid (Newtonian and non-Newtonian) suggest a general applicability of the proposed model. An explanation about of the k _L a values in non-Newtonian fluid is offered take into account shear rate, column geometry, viscosity and results reported in the literature previously. Received on 31 July 1999     

PROBABILISTIC MODELS FOR LONG FATIGUE CRACK GROWTH RATES OF LZ50 AXLE STEEL

IntroductionStudy on the model of fatigue crack growth rate (da/dN) and measurement of theprobabilistic curves of the model are a basis of the fatigue damage tolerate design and theprobabilistic safety assessment. Some benefit efforts have been made in th…     

A turbulence model with variable coefficients for calculating mixing layers and jets

Currently available turbulence models for calculating mixing layers and jets are considered. The reasons for inadequate modeling of initial jet regions are analyzed. The problem of the sensitivity of certain turbulence models to the freestream value of ω is studied. The necessity of introducing variable turbulent-transport coefficients and a correction for axial symmetry is shown. A differential model with variable coefficients for the Reynolds stresses is proposed which makes it possible to obtain adequate solutions in the self-similar regions of mixing layers and jets.     

Dynamic coefficient evaluation for an algebraic subgrid stress model

Static model coefficients for an algebraic subgrid stress (SGS) model are determined using a dynamic approach, based on results from simulations of isotropic decaying turbulence. The study was motivated by the discrepancies in energy transfer predictions using the previously documented coefficients (Bhushan and Warsi, Int. J. Numer. Meth. Fluids 2005; 49 : 489–519). The discrepancies are identified to be due to inconsistent filter functions used in the analytic estimates and the simulations. The study emphasizes that SGS model development should use filter functions compatible with those inherent in CFD application solvers. The dynamic approach predicts consistent model and transfer coefficients for different grid resolutions and is judged to be a reliable basis for model coefficient adjustments. The predicted Leonard's stress coefficient and associated energy transfer coefficients agree very well with the analytic estimates using a Gaussian/cutoff combination filter. This suggests that the modeling of Leonard's stress term using a truncated Taylor series expansion is robust and may not benefit significantly from dynamic modeling. Validation simulations were performed for turbulent channel flow at Re_τ = 180 and 590. The dynamic approach was found to be reliable only for the lower log‐layer of the Re_τ = 590 case, where the scale invariance condition was satisfied. Nonetheless, in this narrow range, the model and transfer coefficients compare well with the isotropic case. The static coefficient algebraic model with new adjusted coefficients shows improved predictions compared with the previous coefficients, for both isotropic decaying turbulence and channel flow cases. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparison of optical and mechanical measurements of second normal stress difference relaxation following step strain

Relaxation of the second normal stress difference (N ₂) following step strain of a concentrated monodisperse polystyrene solution has been studied using mechanical and optical rheometry. Measurements of normal thrust in a parallel plate geometry are corrected for strain inhomogeneity and combined with independent measurements of the first normal stress difference (N ₁) to determine N ₂. Optical experiments were performed using a novel configuration where flow birefringence data collected using multiple light paths within the shear plane are combined with the stress-optical law to determine all three independent stress components for shearing deformations. This technique eliminates end effects, and provides an opportunity to oversample the stress tensor and develop consistency checks of experimental data. N ₂ is found to be nonzero at all accessible times, and relaxes in roughly constant proportion to N ₁. This reflects nonaffine distribution of chain segments, even well within the regime of chain retraction at short times. Data collected with the two techniques are reasonably consistent with each other, and with results of previous studies, generally lying between the predictions of the Doi-Edwards model with and without the independent alignment approximation. The normal stress ratio –N ₂/N ₁ shows pronounced strain thinning in the nonlinear regime.     

Limiting Chain Extensibility Constitutive Models of Valanis–Landel Type

A new general constitutive model in terms of the principal stretches is proposed to reflect limiting chain extensibility resulting in severe strain-stiffening for incompressible, isotropic, homogeneous elastic materials. The strain-energy density involves the logarithm function and has the general Valanis–Landel form. For specific functions in the Valanis–Landel representation, we obtain particular strain-energies, some of which have been proposed in the recent literature. The stress–stretch response in some basic homogeneous deformations is described for these particular strain-energy densities. It is shown that the stress response in these deformations is similar to that predicted by the Gent model involving the first invariant of the Cauchy–Green tensor. The models discussed here depend on both the first and second invariants.      

On non-Newtonian Fluids with Energy Transfer

The model combining incompressible Navier–Stokes’ equation in a non-Newtonian p-power-law modification and the nonlinear heat equation is considered. Existence of its (very) weak solutions is proved for p > 11/5 under mild assumptions of the temperature-dependent stress tensor by careful successive limit passage in a Galerkin approximation.      

A simple sliding‐mesh interface procedure and its application to the CFD simulation of a tidal‐stream turbine

An effective way of using computational fluid dynamics (CFD) to simulate flow about a rotating device—for example, a wind or marine turbine—is to embed a rotating region of cells inside a larger, stationary domain, with a sliding interface between. This paper describes a simple but effective method for implementing this as an internal Dirichlet boundary condition, with interfacial values obtained by interpolation from halo nodes. The method is tested in two finite‐volume codes: one using block‐structured meshes and the other unstructured meshes. Validation is performed for flow around simple, isolated, rotating shapes (cylinder, sphere and cube), comparing, where possible, with experiment and the alternative CFD approach of fixed grid with moving walls. Flow variables are shown to vary smoothly across the sliding interface. Simulations of a tidal‐stream turbine, including both rotor and support, are then performed and compared with towing‐tank experiments. Comparison between CFD and experiment is made for thrust and power coefficients as a function of tip‐speed ratio (TSR) using Reynolds‐averaged Navier–Stokes turbulence models and large‐eddy simulation (LES). Performance of most models is good near the optimal TSR, but simulations underestimate mean thrust and power coefficients in off‐design conditions, with the standard k– ? turbulence model performing noticeably worse than shear stress transport k– ω and Reynolds‐stress‐transport closures. LES gave good predictions of mean load coefficients and vital information about wake structures but at substantial computational cost. Grid‐sensitivity studies suggest that Reynolds‐averaged Navier–Stokes models give acceptable predictions of mean power and thrust coefficients on a single device using a mesh of about 4 million cells. Copyright © 2013 John Wiley & Sons, Ltd.