A smoothed particle hydrodynamics-based fluid model with a spatially dependent viscosity: application to flow of a suspension with a non-Newtonian fluid matrix

A smoothed particle hydrodynamics approach is utilized to model a non-Newtonian fluid with a spatially varying viscosity. In the limit of constant viscosity, this approach recovers an earlier model for Newtonian fluids of Español and Revenga (Phys Rev E 67:026705, 2003). Results are compared with numerical solutions of the general Navier–Strokes equation using the “regularized” Bingham model of Papanastasiou (J Rheol 31:385–404, 1987) that has a shear-rate-dependent viscosity. As an application of this model, the effect of having a non-Newtonian fluid matrix, with a shear-rate-dependent viscosity in a moderately dense suspension, is examined. Simulation results are then compared with experiments on mono-size silica spheres in a shear-thinning fluid and for sand in a calcium carbonate paste. Excellent agreement is found between simulation and experiment. These results indicate that measurements of the shear viscosity of simple shear-rate-dependent non-Newtonian fluids may be used in simulation to predict the viscosity of concentrated suspensions having the same matrix fluid.     

High shear rate viscometry

We investigate the use of two distinct and complementary approaches in measuring the viscometric properties of low viscosity complex fluids at high shear rates up to 80,000 s^?1. Firstly, we adapt commercial controlled-stress and controlled-rate rheometers to access elevated shear rates by using parallel-plate fixtures with very small gap settings (down to 30 μm). The resulting apparent viscosities are gap dependent and systematically in error, but the data can be corrected—at least for Newtonian fluids—via a simple linear gap correction originally presented by Connelly and Greener, J. Rheol, 29(2):209–226, 1985). Secondly, we use a microfabricated rheometer-on-a-chip to measure the steady flow curve in rectangular microchannels. The Weissenberg–Rabinowitsch–Mooney analysis is used to convert measurements of the pressure-drop/flow-rate relationship into the true wall-shear rate and the corresponding rate-dependent viscosity. Microchannel measurements are presented for a range of Newtonian calibration oils, a weakly shear-thinning dilute solution of poly(ethylene oxide), a strongly shear-thinning concentrated solution of xanthan gum, and a wormlike micelle solution that exhibits shear banding at a critical stress. Excellent agreement between the two approaches is obtained for the Newtonian calibration oils, and the relative benefits of each technique are compared and contrasted by considering the physical processes and instrumental limitations that bound the operating spaces for each device.     

Existence and Stability of Solutions for Steady Flows of Fibre Suspension Flows

We establish existence, uniqueness, convergence and stability of solutions to the equations of steady flows of fibre suspension flows. The existence of a unique steady solution is proven by using an iterative scheme. One of the restrictions imposed on the data confirms a well known fact proven in Galdi and Reddy (J Non-Newtonian Fluid Mech 83:205–230, 1999), Munganga and Reddy (Math Models Methods Appl Sci 12:1177–1203, 2002) and Munganga et al. (J Non-Newtonian fluid Mech 92:135–150, 2000) that the particle number N _p must be less than 35/2. Exact solutions are calculated for Couette and Poiseuille flows. Solutions of Poiseuille flows are shown to be more accurate than those of Couette flow when a time perturbation is considered.     

Experimental Measurement of Local Burning Velocity Within a Rotating Flow

The work presented in this paper details the implementation of a new technique for the measurement of local burning velocity using asynchronous particle image velocimetry. This technique uses the local flow velocity ahead of the flame front to measure the movement of the flame by the surrounding fluid. This information is then used to quantify the local burning velocity by taking into account the translation of the flame via convection. In this paper the developed technique is used to study the interaction between a flame front and a single toroidal vortex for the case of premixed stoichiometric methane and air combustion. This data is then used to assess the impact of vortex structure on flame propagation rates. The burning velocity data demonstrates that there is a significant enhancement to the rate of flame propagation where the flame directly interacts with the rotating vortex. The increases found were significantly higher than expected but are supported by burning velocities (Filatyev et al, Combust Flame 141:1?C21, 2005; Kobayashi et al, Proc Combust Inst 29:1793?C1800, 2002; Shepherd et al. 1998) found in turbulent flames of the same mixture composition. Away from this interaction with the main vortex core, the flame exhibits propagation rates around the value recorded in literature for unperturbed laminar combustion (Tahtouh et al, Combust Flame 159:1735?C1743, 2009; Hassan et al, Combust Flame 115:539?C550, 1998); Halter et al, Proc Combust Inst 30:201?C208, 2005; Coppens et al, Exp Therm Fluid Sci 31:437?C444, 2007).     

Accounting for Polydispersion in the Eulerian Large Eddy Simulation of the Two-Phase Flow in an Aeronautical-type Burner

A major issue for the simulation of two-phase flows in engines concerns the modeling of the liquid disperse phase, either in the Lagrangian or the Eulerian approach. In the perspective of massively parallel computing, the Eulerian approach seems to be more suitable, as it uses the same algorithms as the gaseous phase solver. However taking into account the whole physics of a turbulent spray, especially in terms of polydispersity, requires an additional modeling effort. The Mesoscopic Eulerian Formalism (MEF) [13] accounts for the effect of turbulence on the disperse phase, and was extended to the Large Eddy Simulation framework [41], but is limited to monodisperse flows. In [38], the influence of polydispersity on resolved and unresolved turbulent motions of the disperse phase was highlighted, and a first model was proposed, based on size-conditioned statistics. Starting from this idea, a coupling between the MEF and the Multifluid Approach (MA) [30] is proposed. The MA decomposes the Eulerian phase into several fluid classes called sections, and corresponding to size intervals. Each section uses then size-conditioned closures. The original idea of this work is to use the MEF closures independently in each section, taking into account the mean droplet size of this section. This new approach, called Multifluid Mesoscopic Eulerian Formalism (MMEF), is then able to capture polydispersion with associated size-conditioned turbulent dynamics. First, the importance of polydispersity and the ability of MMEF to capture it are highlighted with a 0D evaporation case and a 2D vortex case, showing its impact on dynamics in both size and physical spaces. Then, the MMEF is applied to the MERCATO configuration of ONERA [18]. Results are compared to monodisperse Eulerian, Lagrangian and experimental results.     

The Two-Dimensional Euler Equations on Singular Domains

We establish the existence of global weak solutions of the two-dimensional incompressible Euler equations for a large class of non-smooth open sets. Loosely, these open sets are the complements (in a simply connected domain) of a finite number of obstacles with positive Sobolev capacity. Existence of weak solutions with L ^p vorticity is deduced from a property of domain continuity for the Euler equations that relates to the so-called γ-convergence of open sets. Our results complete those obtained for convex domains in Taylor (Progress in Nonlinear Differential Equations and their Applications, Vol. 42, 2000), or for domains with asymptotically small holes (Iftimie et al. in Commun Partial Differ Equ 28(1–2), 349–379, 2003; Lopes Filho in SIAM J Math Anal 39(2), 422–436, 2007).     

Global Solvability of a Free Boundary Three-Dimensional Incompressible Viscoelastic Fluid System with Surface Tension

Motivated by Beale (Commun Pure Appl Math 34:359–392, 1981; Arch Ration Mech Anal 84:307–352, 1983/1984), we investigate the global well-posedness of a free boundary problem of a three-dimensional incompressible viscoelastic fluid system in an infinite strip and with surface tension on the upper free boundary, provided that the initial data is sufficiently close to the equilibrium state.     

On the Formulation of Continuum Thermodynamic Models for Solids as General Equations for Non-equilibrium Reversible-Irreversible Coupling

The purpose of this work is the comparison of some aspects of the formulation of material models in the context of continuum thermodynamics (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997) with their formulation in the form of a General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC: e.g., Grmela and Öttinger in Phys. Rev. E 56: 6620–6632, 1997; Öttinger and Grmela in Phys. Rev. E 56: 6633–6655, 1997; Öttinger in Beyond equilibrium thermodynamics, Wiley, New York, 2005; Grmela in J. Non-Newton. Fluid Mech. 165: 980–998, 2010). A GENERIC represents a generalization of the Ginzburg-Landau model for the approach of non-equilibrium systems to thermodynamic equilibrium. Originally developed to formulate non-equilibrium thermodynamic models for complex fluids, it has recently been applied to anisotropic inelastic solids in a Eulerian setting (Hütter and Tervoort in J. Non-Newton. Fluid Mech. 152: 45–52, 2008; 53–65, 2008; Adv. Appl. Mech. 42: 254–317, 2009) as well as to damage mechanics (Hütter and Tervoort in Acta Mech. 201: 297–312, 2008). In the current work, attention is focused for simplicity on the case of thermoelastic solids with heat conduction and viscosity in a Lagrangian setting (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997, Chaps. 9–12). In the process, the relation of the two formulations to each other is investigated in detail. A particular point in this regard is the concept of dissipation and its model representation in both contexts.     

Numerical investigation of pyrolysis of a Loy Yang coal in a lab-scale furnace at elevated pressures

A computational fluid dynamics (CFD) model of the pyrolysis of a Loy Yang low-rank coal in a pressurised drop tube furnace (pdtf) was undertaken evaluating Arrhenius reaction rate constants. The paper also presents predictions of an isothermal flow through the drop tube furnace. In this study, a pdtf reactor operated at pressures up to 15 bar and at a temperature of 1,173 K with particle heating rates of approximately 10⁵ K s^?1 was used. The CFD model consists of two geometrical sections; flow straightner and injector. The single reaction and two competing reaction models were employed for this numerical investigation of the pyrolysis process. The results are validated against the available experimental data in terms of velocity profiles for the drop tube furnace and the particle mass loss versus particle residence times. The isothermal flow results showed reasonable agreement with the available experimental data at different locations from the injector tip. The predicted results of both the single reaction and competing reaction modes showed slightly different results. In addition, several reaction rate constants were tested and validated against the available experimental data. The most accurate results were being Badzioch and Hawksley (Ind Eng Chem Process Des Dev 9:521–530, 1970) with a single reaction model and Ubhayakar et al. (Symp (Int) Combust 16:427–436, 1977) for two competing reactions. These numerical results can provide useful information towards future modelling of the behaviour of Loy Yang coal in a full scale tangentially-fired furnace.     

Flow temporal reconstruction from non-time-resolved data part I: mathematic fundamentals

At least two circumstances point to the need of postprocessing techniques to recover lost time information from non-time-resolved data: the increasing interest in identifying and tracking coherent structures in flows of industrial interest and the high data throughput of global measuring techniques, such as PIV, for the validation of computational fluid dynamics (CFD) codes. This paper offers the mathematic fundamentals of a space--time reconstruction technique from non-time-resolved, statistically independent data. An algorithm has been developed to identify and track traveling coherent structures in periodic flows. Phase-averaged flow fields are reconstructed with a correlation-based method, which uses information from the Proper Orthogonal Decomposition (POD). The theoretical background shows that the snapshot POD coefficients can be used to recover flow phase information. Once this information is recovered, the real snapshots are used to reconstruct the flow history and characteristics, avoiding neither the use of POD modes nor any associated artifact. The proposed time reconstruction algorithm is in agreement with the experimental evidence given by the practical implementation proposed in the second part of this work (Legrand et al. in Exp Fluids, 2011), using the coefficients corresponding to the first three POD modes. It also agrees with the results on similar issues by other authors (Ben Chiekh et al. in 9 Congrès Francophone de Vélocimétrie Laser, Bruxelles, Belgium, 2004; Van Oudheusden et al. in Exp Fluids 39-1:86?C98, 2005; Meyer et al. in 7th International Symposium on Particle Image Velocimetry, Rome, Italy, 2007a; in J Fluid Mech 583:199?C227, 2007b; Perrin et al. in Exp Fluids 43-2:341?C355, 2007). Computer time to perform the reconstruction is relatively short, of the order of minutes with current PC technology.     

Analytical Validation of a Continuum Model for Epitaxial Growth with Elasticity on Vicinal Surfaces

Within the context of heteroepitaxial growth of a film onto a substrate, terraces and steps self-organize according to misfit elasticity forces. Discrete models of this behavior were developed by Duport et al. (J Phys I 5:1317–1350, 1995) and Tersoff et al. (Phys Rev Lett 75:2730–2733, 1995). A continuum limit of these was in turn derived by Xiang (SIAM J Appl Math 63:241–258, 2002) (see also the work of Xiang and Weinan Phys Rev B 69:035409-1–035409-16, 2004; Xu and Xiang SIAM J Appl Math 69:1393–1414, 2009). In this paper we formulate a notion of weak solution to Xiang’s continuum model in terms of a variational inequality that is satisfied by strong solutions. Then we prove the existence of a weak solution.     

Particle-based numerical modeling of AE statistics in disordered materials

We present some numerical simulations of AE due to damage propagation in disordered materials under compression and bending. To this purpose, the AE cumulative number, the time frequency analysis and the statistical properties of AE time series will be numerically simulated adopting the so-called “particle method strategy” (Cundall in Proceedings ISRM Symp., Nancy, France, vol. 2, pp. 129–136, 1971). The method provides the velocity of particles in a set simulating the behavior of a granular system and, therefore, is suitable to model the compressive wave propagation and acoustic emission (corresponding to cracking) in a solid body. The numerical simulations (Abe et al. in Pure Appl. Geophys. 161:2265–2277, 2004) correctly describe the compression test in terms of mean stress-strain response and crack pattern (Invernizzi et al. in Proceedings of the SEM annual conference, Society for Experimental Mechanics Inc., USA, SEM Annual Conference, Indianapolis, Indiana, USA, June 7–10, 2010). The size effects on the peak compressive strength and on the AE count are correctly reproduced. In addition, the amplitude distribution (b-value) and temporal evolution of AE events due to cracking, crucial for the evaluation of damage and remaining lifetime, were simulated and result in agreement with the experimental evidences.     

Scale effect on unsteady cloud cavitation

No experiment was conducted, yet, to investigate the scale effects on the dynamics of developed cavitating flow with periodical cloud shedding. The present study was motivated by the unclear results obtained from the experiments in a Venturi-type section that was scaled down 10 times for the purpose of measurements by ultra-fast X-ray imaging (Coutier-Delgosha et al. 2009). Cavitation in the original size scale section (Stutz and Reboud in Exp Fluids 23:191–198, 1997, Exp Fluids 29:545–552 2000) always displays unsteady cloud separation. However, when the geometry was scaled down, the cavitation became quasi steady although some oscillations still existed. To investigate this phenomenon more in detail, experiments were conducted in six geometrically similar Venturi test sections where either width or height or both were scaled. Various types of instabilities are obtained, from simple oscillations of the sheet cavity length to large vapor cloud shedding when the size of the test section is increased. It confirms that small scale has a significant influence on cavitation. Especially the height of the test section plays a major role in the dynamics of the re-entrant jet that drives the periodical shedding observed at large scale. Results suggest that the sheet cavity becomes stabile when the section is scaled down to a certain point because re-entrant jet cannot fully develop.     

On Co-Circular Central Configurations in the Four and Five Body-Problems for Homogeneous Force Law

We study the central configurations (cc for short) for four masses arranged on a common circle (called co-circular cc) in two different situations, namely with no mass inside and later adding a fifth mass at the center of the circle. In the former, we focus the kite shape configurations by proving the existence of a one-parameter family of cc which goes from the kite containing an equilateral triangle up to the square shape. After, by putting a fifth mass at the center, we feature the planar cc of five bodies as a tensor of corange two see, “Albouy and Chenciner (Invent Math 131:151–184, 1998)” and we prove that cc is stacked see, “Hampton (Nonlinearity 18:2299–2304, 2005b)” in a such way that the center of mass of the four bodies should be the center of the circle. We emphasize that our approach includes not only the Newtonian force law, but the homogeneous ones with exponent $a\le -1$ a ≤ ? 1 .     

Laminar,transitional and turbulent annular flow of drag-reducing polymer solutions

Mean and rms axial velocity-profile data obtained using laser Doppler anemometry are presented together with pressure-drop data for the flow through a concentric annulus (radius ratio κ = 0.506) of a Newtonian (a glycerine–water mixture) and non-Newtonian fluids—a semi-rigid shear-thinning polymer (a xanthan gum) and a polymer known to exhibit a yield stress (carbopol). A wider range of Reynolds numbers for the transitional flow regime is observed for the more shear-thinning fluids. In marked contrast to the Newtonian fluid, the higher shear stress on the inner wall compared to the outer wall does not lead to earlier transition for the non-Newtonian fluids where more turbulent activity is observed in the outer wall region. The mean axial velocity profiles show a slight shift (～5%) of the location of the maximum velocity towards the outer pipe wall within the transitional regime only for the Newtonian fluid.     

Quasi-linear versus potential-based formulations of force–flux relations and the GENERIC for irreversible processes: comparisons and examples

An essential part in modeling out-of-equilibrium dynamics is the formulation of irreversible dynamics. In the latter, the major task consists in specifying the relations between thermodynamic forces and fluxes. In the literature, mainly two distinct approaches are used for the specification of force–flux relations. On the one hand, quasi-linear relations are employed, which are based on the physics of transport processes and fluctuation–dissipation theorems (de Groot and Mazur in Non-equilibrium thermodynamics, North Holland, Amsterdam, 1962, Lifshitz and Pitaevskii in Physical kinetics. Volume 10, Landau and Lifshitz series on theoretical physics, Pergamon Press, Oxford, 1981). On the other hand, force–flux relations are also often represented in potential form with the help of a dissipation potential (?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997). We address the question of how these two approaches are related. The main result of this presentation states that the class of models formulated by quasi-linear relations is larger than what can be described in a potential-based formulation. While the relation between the two methods is shown in general terms, it is demonstrated also with the help of three examples. The finding that quasi-linear force–flux relations are more general than dissipation-based ones also has ramifications for the general equation for non-equilibrium reversible–irreversible coupling (GENERIC: e.g., Grmela and Öttinger in Phys Rev E 56:6620–6632, 6633–6655, 1997, Öttinger in Beyond equilibrium thermodynamics, Wiley Interscience Publishers, Hoboken, 2005). This framework has been formulated and used in two different forms, namely a quasi-linear (Öttinger and Grmela in Phys Rev E 56:6633–6655, 1997, Öttinger in Beyond equilibrium thermodynamics, Wiley Interscience Publishers, Hoboken, 2005) and a dissipation potential–based (Grmela in Adv Chem Eng 39:75–129, 2010, Grmela in J Non-Newton Fluid Mech 165:980–986, 2010, Mielke in Continuum Mech Therm 23:233–256, 2011) form, respectively, relating the irreversible evolution to the entropy gradient. It is found that also in the case of GENERIC, the quasi-linear representation encompasses a wider class of phenomena as compared to the dissipation-based formulation. Furthermore, it is found that a potential exists for the irreversible part of the GENERIC if and only if one does for the underlying force–flux relations.     

A phenomenological model for microstructure-dependent inelasticity in shape-memory alloys

Degradation in shape-memory alloy response is a crucial concern for a variety of innovative applications. Under cyclic loadings, these materials generally experience permanent inelastic deformations. The onset of plasticization is known to be very sensitive to the microstructure of the polycrystalline specimen. Moving from recent experimental findings (Malard et al. in Funct Mater Lett 2:45–54, 2009; Acta Mater 59:1542–1556, 2011), we present a phenomenological model for permanent inelastic effects in shape-memory alloys taking into account the polycrystalline microstructure. In particular, the mechanical response under cyclic loadings is investigated in connection with the mean crystal grain size. Formulated within the variational frame of generalized standard materials, the model consists in an extension of the model in Auricchio et al. (Int J Plast 23:207–226, 2007) to the case of microstructure-dependent parameters. The mathematical setting is discussed and numerical simulations showing the capability of the model to reproduce experiments are presented.