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.     

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.     

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).     

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.     

A New Divergence Free Synthetic Eddy Method for the Reproduction of Inlet Flow Conditions for LES

This paper describes a recent development of the Synthetic Eddy Method (SEM) proposed by Jarrin et al. (Int J Heat Fluid Flow 30(3):435–442, 2009) for generation of synthetic turbulence. The present scheme is designed to produce a divergence-free turbulence field that can reproduce almost all possible states of Reynolds stress anisotropy. This improved representation, when used to provide inlet conditions for an LES, leads to reduced near-inlet pressure fluctuations in the LES and to a reduced development length, both of which lead to lower computer resource requirements. An advantage of this method with respect to forcing approaches (which require an iterative approach) is the suitability for direct usage with embedded LES. Results for a turbulent channel flow are reported here and compared to those from the original SEM, and other direct approaches such as the VORTEX method of Sergent (2002) and the Synthesized Turbulence approach of Davidson and Billson (Int J Heat Fluid Flow 27(6):1028–1042, 2006), showing overall improved performance and a more accurate representation of turbulence structures immediately downstream of the inlet.     

Measurements of an unsteady liquid metal flow during spin-up driven by a rotating magnetic field

A cylindrical cavity with an aspect ratio of unity is filled with liquid metal and suddenly exposed to an azimuthal body force generated by a rotating magnetic field (RMF). This experimental study is concerned with the secondary meridional flow during the time, if the fluid spins up from rest. Vertical profiles of the axial velocity have been measured by means of the ultrasound Doppler velocimetry. The flow measurements confirm the spin-up concept by Ungarish (J Fluid Mech 347:105–118, 1997) and the continuative study by Nikrityuk et al. (Phys Fluids 17:067101, 2005) who suggested the existence of two stages during the RMF-driven spin-up, in particular the so-called initial adjustment phase followed by an inertial phase which is dominated by inertial oscillations of the secondary flow. Evolving instabilities of the double-vortex structure of the secondary flow have been detected at a Taylor number of 1.24 × 10⁵ verifying the predictions of Grants and Gerbeth (J Fluid Mech 463:229–240, 2002). Perturbations in form of Taylor–Görtler vortices have been observed just above the instability threshold.     

Scaling the normal stresses in concentrated non-colloidal suspensions of spheres

We have collected six sets of results for the ratio of the second normal stress difference to the shear stress $(N_{2}/\tau )$ in non-colloidal suspensions of spheres in Newtonian matrices. They all show a near-cubic dependence on the volume fraction $\varphi $ in the range $0.1 < \varphi < 0.5$ , in contrast to the square law predictions of Brady and Morris (J Fluid Mech 348:103–139, 1997) for dilute suspensions. We suggest that the difference can be resolved by using a dependence on the square of the effective volume fraction $\varphi _{\textrm e}$ , and good agreement is then found.     

Capillary Floating and the Billiard Ball Problem

In a study of capillary floating, Finn (J Math Fluid Mech 11:443?C458, 2009) described a procedure for determining cross-sections of non-circular, infinite convex cylinders that float horizontally on a liquid surface in every orientation with contact angle ??/2. Finn??s procedure yielded incomplete results for other contact angles; he raised the question as to whether an analogous construction would be feasible in that case. In the note, Finn (J Math Fluid Mech 11:464?C465, 2009) pointed out a connection with an independent problem on billiard caustics citing the unpublished work (Gutkin in Proceedings of the Workshop on Dynamics and Related Questions, PennState University, 1993) of the present author. Here we present a solution of the billiard problem in full detail, thus settling Finn??s question in a surprising way. In particular, we show that such floating cylinders exist if and only if the contact angle lies in a certain, explicitly described countably dense set. Moreover, for each element ?? in this set we exhibit a family of convex, non-circular cylinders that float in every orientation with contact angle ??. Our discussion contains other material of independent interest for the billiard ball problem.     

On Near-Wall Treatment in (U)RANS-Based Closure Models

The present work is concerned with computational evaluation of a recently formulated near-wall relationship providing the value of the dissipation rate ε of the kinetic energy of turbulence k through its exact dependence on the Taylor microscale λ: ε = 10νk/λ ², (Jakirli? and Jovanovi?, J. Fluid Mech. 656:530–539, 2010). Dissipation rate determination benefits from the asymptotic behavior of the Taylor microscale resulting in its linear variation in terms of the wall distance (λ?∝?y) being valid throughout entire viscous sublayer. Accordingly, it can be applied as a unified near-wall treatment in all computational frameworks relying on a RANS-based model of turbulence (including also hybrid LES/RANS schemes) independent of modeling level—both main modeling concepts eddy-viscosity and Reynolds stress models can be employed. Presently, the feasibility of the proposed formulation was demonstrated by applying a conventional near-wall second-moment closure model based on the homogeneous dissipation rate ε _h ( ${\varepsilon_h =\varepsilon -0.5\partial \left( {{\nu \partial k}/ {\partial x_j }} \right)} / {\partial x_j }$ ; Jakirli? and Hanjali?, J. Fluid Mech. 539:139–166, 2002) and its instability-sensitive version, modeled in terms of the inverse turbulent time scale ω _h (ω _h ?=?ε _h /k; Maduta and Jakirli?, 2011), to a fully-developed channel flow with both flat walls and periodic hill-shaped constrictions mounted on the bottom wall in a Reynolds number range. The latter configuration is subjected to boundary layer separation from a continuous curved wall. The influence of the near-wall resolution lowering with respect to the location of the wall-closest computational node, coarsened even up to the viscous sublayer edge situated at $y_P^+ \approx 5$ in equilibrium flows, is analyzed. The results obtained follow closely those pertinent to the conventional near-wall integration with the wall-next node positioned at $y_P^+ \le 0.5$ .     

Local Controllability to Trajectories of the Magnetohydrodynamic Equations

In this paper we prove the local controllability to trajectories of the three dimensional magnetohydrodynamic equations by means of two internal controls, one in the velocity equations and the other in the magnetic field equations and both localized in an arbitrary small subset with not empty interior of the domain. This paper improves the previous results (Barbu et al. in Comm Pure Appl Math 56:732–783, 2003; Barbu et al. in Adv Differ Equ 10:481–504, 2005; Havârneanu et al. in Adv Differ Equ 11:893–929, 2006; Havârneanu, in SIAM J Control Optim 46:1802–1830, 2007) where the second control is not localized and it allows to deduce the local controllability to trajectories with boundary controls. The proof relies on the Carleman inequality for the Stokes system of Imanuvilov et al. (Carleman estimates for second order nonhomogeneous parabolic equations, preprint) to deal with the velocity equations and on a new Carleman inequality for a Dynamo-type equation to deal with the magnetic field equations.     

Thin Shear Layer Structures in High Reynolds Number Turbulence

Three-dimensional tomographic time dependent PIV measurements of high Reynolds number (Re) laboratory turbulence are presented which show the existence of long-lived, highly sheared thin layer eddy structures with thickness of the order of the Taylor microscale and internal fluctuations. Highly sheared layer structures are also observed in direct numerical simulations of homogeneous turbulence at higher values of Re (Ishihara et al., Annu Rev Fluid Mech 41:165–180, 2009). But in the latter simulation, where the fluctuations are more intense, the layer thickness is greater. A rapid distortion model describes the structure and spectra for the velocity fluctuations outside and within ‘significant’ layers; their spectra are similar to the Kolmogorov (C R Acad Sci URSS 30:299–303, 1941) and Obukhov (Dokl Akad Nauk SSSR 32:22–24, 1941) statistical model (KO) for the whole flow. As larger-scale eddy motions are blocked by the shear layers, they distort smaller-scale eddies leading to local zones of down-scale and up-scale transfer of energy. Thence the energy spectrum for high wave number k is $E_X (k)\sim Bk^{-2p}$ . The exponent p depends on the forms of the large eddies. The non-linear interactions between the distorted inhomogeneous eddies produce a steady local structure, which implies that 2p?=?5/3 and a flux of energy into the thin-layers balancing the intense dissipation, which is much greater than the mean $\left<\epsilon\right>$ . Thence $B\sim\left<\epsilon\right>^{2/3}$ as in KO. Within the thin layers the inward flux energises extended vortices whose thickness and spacing are comparable with the viscous microscale. Although peak values of vorticity and velocity of these vortices greatly exceed those based on the KO scaling, the form of the viscous range spectrum is consistent with their model.     

On the Energy Dissipation Rate of Solutions to the Compressible Isentropic Euler System

In this paper we extend and complement the results in Chiodaroli et al. (Global ill-posedness of the isentropic system of gas dynamics, 2014) on the well-posedness issue for weak solutions of the compressible isentropic Euler system in 2 space dimensions with pressure law p(ρ) = ρ ^γ , γ ≥ 1. First we show that every Riemann problem whose one-dimensional self-similar solution consists of two shocks admits also infinitely many two-dimensional admissible bounded weak solutions (not containing vacuum) generated by the method of De Lellis and Székelyhidi (Ann Math 170:1417–1436, 2009), (Arch Ration Mech Anal 195:225–260, 2010). Moreover we prove that for some of these Riemann problems and for 1 ≤ γ < 3 such solutions have a greater energy dissipation rate than the self-similar solution emanating from the same Riemann data. We therefore show that the maximal dissipation criterion proposed by Dafermos in (J Diff Equ 14:202–212, 1973) does not favour the classical self-similar solutions.     

Estimation of the largest Lyapunov exponent-like (LLEL) stability measure parameter from the perturbation vector and its derivative dot product (part 2) experiment simulation

Controlling system dynamics with use of the Largest Lyapunov Exponent (LLE) is employed in many different areas of the scientific research. Thus, there is still need to elaborate fast and simple methods of LLE calculation. This article is the second part of the one presented in Dabrowski (Nonlinear Dyn 67:283–291, 2012). It develops method LLEDP of the LLE estimation and shows that from the time series of two identical systems, one can simply extract value of the stability parameter which value can be treated as largest LLE. Unlike the method presented in part, one developed method (LLEDPT) can be applied to the dynamical systems of any type, continuous, with discontinuities, with time delay and others. The theoretical improvement shows simplicity of the method and its obvious physical background. The proofs for the method effectiveness are based on results of the simulations of the experiments for Duffing and Van der Pole oscillators. These results were compared with ones obtained with use of the Stefanski method (Stefanski in Chaos Soliton Fract 11(15):2443–2451, 2000; Chaos Soliton Fract 15:233–244, 2003; Chaos Soliton Fract 23:1651–1659, 2005; J Theor Appl Mech 46(3):665–678, 2008) and LLEDP method. LLEDPT can be used also as the criterion of stability of the control system, where desired behavior of controlled system is explicitly known (Balcerzak et al. in Mech Mech Eng 17(4):325–339, 2013). The next step of development of the method can be considered in direction that allows estimation of LLE from the real time series, systems with discontinuities, with time delay and others.     

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.     

Kirchhoff plates with viscous boundary dissipation

In this paper, we confine our attention to Kirchhoff thin plates in presence of boundary viscoelastic dissipative mechanisms, in order to investigate the well-posedness and the asymptotic behavior within the minimal state approach, following the guidelines proposed in Deseri et al. (Arch Rational Mech Anal 181:43–96, 2006) [see also Fabrizio et al. (Arch Rational Mech Anal 198:189–232, 2010)].     

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.     

Existence Result for a Dynamical Equations of Generalized Marguerre-von Kármán Shallow Shells

In a recent work in the static case, Gratie (Appl. Anal. 81:1107–1126, 2002) has generalized the classical Marguerre-von Kármán equations studied by Ciarlet and Paumier in (Comput. Mech. 1:177–202, 1986), where only a portion of the lateral face of the shallow shell is subjected to boundary conditions of von Kármán type, while the remaining portion is subjected to boundary conditions of free edge. Then Ciarlet and Gratie (Math. Mech. Solids 11:83–100, 2006) have established an existence theorem for these equations. In Chacha et al. (Rev. ARIMA 13:63–76, 2010), we extended formally these studies to the dynamical case. More precisely, we considered a three-dimensional dynamical model for a nonlinearly elastic shallow shell with a specific class of boundary conditions of generalized Marguerre-von Kármán type. Using technics from formal asymptotic analysis, we showed that the scaled three-dimensional solution still leads to two-dimensional dynamical boundary value problem called the dynamical equations of generalized Marguerre-von Kármán shallow shells. In this paper, we establish the existence of solutions to these equations using a compactness method of Lions (Quelques Méthodes de Résolution des Problèmes aux Limites non Linéaires, Dunod, Paris, 1969).     

Well-Posedness for the Motion of Physical Vacuum of the Three-dimensional Compressible Euler Equations with or without Self-Gravitation

This paper concerns the well-posedness theory of the motion of a physical vacuum for the compressible Euler equations with or without self-gravitation. First, a general uniqueness theorem of classical solutions is proved for the three dimensional general motion. Second, for the spherically symmetric motions, without imposing the compatibility condition of the first derivative being zero at the center of symmetry, a new local-in-time existence theory is established in a functional space involving less derivatives than those constructed for three-dimensional motions in (Coutand et al., Commun Math Phys 296:559–587, 2010; Coutand and Shkoller, Arch Ration Mech Anal 206:515–616, 2012; Jang and Masmoudi, Well-posedness of compressible Euler equations in a physical vacuum, 2008) by constructing suitable weights and cutoff functions featuring the behavior of solutions near both the center of the symmetry and the moving vacuum boundary.