A novel transient wall heat transfer approach for the start-up of SI engines with gasoline direct injection

The introduction of CO₂-reduction technologies like Start–Stop or the Hybrid-Powertrain and the worldwide stringent emission legislation require a detailed optimization of the engine start-up. The combustion concept development as well as the calibration of the engine control unit makes an explicit thermodynamic analysis of the combustion process during the start-up necessary. Initially, the well-known thermodynamic analysis of in-cylinder pressure at stationary condition was transmitted to the highly non-stationary engine start-up. For this running mode of the engine the current models for calculation of the transient wall heat fluxes were found to be misleading. With a fraction of nearly 45% of the burned fuel energy, the wall heat is very important for the calculation of energy balance and for the combustion process analysis. Based on the measurements of transient wall heat transfer densities during the start-up presented in a former work (Lejsek and Kulzer in Investigations on the transient wall heat transfer at start-up for SI engines with gasoline direct injection. SAE Paper), the paper describes the development of adaptations to the known correlations by Woschni (MTZ 31:491, 1970), Hohenberg (Experimentelle Erfassung der Wandwärme von Kolbenmotoren. TU Graz, Habil., 1980) and Bargende (Ein Gleichungsansatz zur Berechnung der instationären Wandwärmeverluste im Hochdruckteil von Ottomotoren. TH Darmstadt, PhD-Thesis, 1991) for the application during engine start-up. To demonstrate the high accuracy of the model, the results of the cyclic resolved thermodynamic analysis using the presented novel approaches were compared with the results of the measurements. It is shown, that the novel heat flux models for the engine start-up process gives a cyclic resolved thermodynamic analysis to optimize the engine start-up pretty efficient.     

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

Centreline velocity decay characterisation in low-velocity jets downstream from an extended conical diffuser

The centreline velocity decay of round airflow jets issuing from extended conical diffusers with length-to-diameter ratio 1.2≤L _t /d≤20 is studied for moderate bulk Reynolds numbers 1131≤Re _b ≤9054. The centreline velocity decay varies as a function of the initial conditions. The functional correlation between the centreline velocity decay coefficient and the initial centreline turbulence level observed on convergent nozzles (Malmström et al. in J. Fluid Mech. 246:363–377, 1997) breaks down as the initial centreline turbulence level exceeds 20 %. In addition, the centreline velocity decay coefficient expressed as function of the bulk velocity U _b decreases for U _b <3 m/s instead of initial mean velocity U ₀<6 m/s as reported for convergent nozzles (Malmström et al. in J. Fluid Mech. 246:363–377, 1997). The asymptotic values of the decay coefficient for U _b >3 m/s decrease linearly when expressed as function of the initial centreline turbulence intensity u ₀/U ₀. Studied flow and geometrical conditions are relevant to flow through the human upper airways.     

PIV experiments in rough-wall, laminar-to-turbulent, oscillatory boundary-layer flows

Exploratory measurements of oscillatory boundary layers were conducted over a smooth and two different rough beds spanning the laminar, transitional and turbulent flow regimes using a multi-camera 2D-PIV system in a small oscillatory-flow tunnel (Admiraal et al. in J Hydraul Res 44(4):437–450, 2006). Results show how the phase lag between bed shear stress and free-stream velocity is better defined when the integral of the momentum equation is used to estimate the bed shear stress. Observed differences in bed shear stress and phase lag between bed shear stress and free-stream velocity are highly sensitive to the definition of the bed position (y = b). The underestimation of turbulent stresses close to the wall is found to explain such differences when using the addition of Reynolds and viscous stresses to define both the bed shear stress and the phase lag. Regardless of the flow regime, in all experiments, boundary-layer thickness reached its maximum value at a phase near the flow reversal at the wall. Friction factors in smooth walls are better estimated using a theoretical equation first proposed by Batchelor (An introduction to fluid dynamics. Cambridge University Press, Cambridge, 1967) while the more recent empirical predictor of Pedocchi and Garcia (J Hydraul Res 47(4):438–444, 2009a) was found to be appropriate for estimating friction coefficients in the laminar-to-turbulent transition regime.     

Filling the gap between transient and steady shear rheology of aqueous graphene oxide dispersions

Even though the rheological behavior of aqueous graphene oxide (G-O) dispersions has been shown to be strongly time-dependent, only few transient measurements have been reported in the literature. In this work, we attempt to fill the gap between transient and steady shear rheological characterizations of aqueous G-O dispersions in the concentration range of 0.004 < ? <?3.5 wt%, by conducting comprehensive rheological measurements, including oscillatory shear flow, transient shear flow, and steady shear flow. Steady shear measurements have been performed after the evaluation of transient properties of the G-O dispersions, to assure steady-state conditions. We identify the critical concentration ? _c =?0.08 wt% (where G-O sheets start to interact) from oscillatory shear experiments. We find that the rheology of G-O dispersions strongly depends on the G-O concentration ?. Transient measurements of shear viscosity and first normal stress difference suggest that G-O dispersions behave like nematic polymeric liquid crystals at ?/? _c =?25, in agreement with other work reported in the literature. G-O dispersions also display a transition from negative to positive values of the first normal stress difference with increasing shear rates. Experimental findings of aqueous graphene oxide dispersions are compared and discussed with models and experiments reported for nematic polymeric liquid crystals, laponite, and organoclay dispersions.     

Dynamics of shear aligning of nematic liquid crystal monodomains

The equations of linear and angular momentum for nematic liquid crystals have been described with Ericksen's transversely isotropic fluid [TIF] model and solved for start-up of shear flow at constant rate and varying initial alignment conditions. An analytical solution for the rotation provides predictions of the nematic director which closely agree with experimental results of Boudreau et al. (1999), supporting the validity of Ericksen's TIF model. The solution is limited to flows where the effects of director gradients are negligible. Received: 13 September 1999/Accepted: 24 January 2000     

Jeffrey-Hamel flow of Leslie-Ericksen nematic liquids

A similarity solution of the Leslie-Ericksen equations for nematic liquid crystals is obtained for flow between converging and diverging planar walls (Jeffrey-Hamel flow). There are three regions in the flow: extensional or compressional flow near the centerline, shear near the wall, and a wall boundary layer in which elastic stresses control the transition from the wall-induced orientation to the bulk behavior. The boundary layer thickness is obtained in closed form; the scaling with the Ericksen number depends on whether or not the boundary layer extends into the region of extensional flow. Imposition of a magnetic field with an azimuthal component in a converging flow can result in a Freedericksz-like transition from radial to transverse orientation at the center line at a critical field strength. This transition provides a new means to measure the irrotational viscosity λ₂.     

Effects of relaxation time on start-up time for starting flow of Maxwell fluid in a pipe

Although the analytical solution of the starting flow of Maxwell fluid in a pipe has been derived for a long time, the effect of relaxation time λ on start-up time ts of this flow is still not well understood. Especially, there exist a series of jumps on the ts-λ. curve. In this paper we introduce a normalized mechanical energy by mode decomposition and mathematical analogy to describe the start-up process. An improved definition of start-up time is presented based on the normalized mechanical energy. It is proved that the ts-λ. curve contains a series of jumps if λ is larger than a critical value. The exact positions of the jumps are determined and the physical reason of the jumps is discussed.     

Scaling relations for elongational flow of polystyrene melts and concentrated solutions of polystyrene in oligomeric styrene

A consistent model of the rheology of polymer melts and concentrated solutions is presented, based on the idea that the pressures exerted by a polymer chain on the walls of an anisotropic confinement are anisotropic (Doi and Edwards. The Theory of Polymer Dynamics, Oxford University Press, 1986). In a tube model with variable tube diameter, chain stretch and tube diameter reduction are related, and at deformation rates larger than the inverse Rouse time τ _R, the chain is stretched and its confining tube becomes increasingly anisotropic. Tube diameter reduction leads to an interchain pressure in the lateral direction of the tube (Marrucci and Ianniruberto. Macromolecules 37:3934-3942, 2004). Chain stretch is balanced by interchain tube pressure in the lateral direction, which is proportional to the third power of stretch, and by a spring force in the longitudinal direction of the tube, which is linear in stretch. Analyzing elongational viscosity data of Huang et al. (Macromolecules 46:5026-5035, 2013a; ACS Macro Letters 2:741-744, 2013b) shows that dilution of polystyrene by oligomeric styrene does not change the relative interchain tube pressure. Based on this extended interchain pressure concept, scaling relations for linear viscoelasticity and elongational viscosity of polystyrene melts and concentrated solutions of polystyrene in oligomeric styrene are presented based exclusively on the relaxation modulus of a reference polymer melt, the volume fraction of polymer in the solution, and the time-molar-mass shift as well as the time-temperature shift caused by the reduction of the glass transition temperature T _g of the polymer in a solution relative to T _g of the melt.     

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.     

Normal Mode Analysis of the High Speed Channel Flow in a Bidisperse Porous Medium

The linear Darcy–Brinkman model of the high speed flow in a bidisperse porous medium proposed by Nield and Kuznetsov (Transport Phenomena in Porous Media, 2005) is revisited in this paper. For the steady unidirectional flow in a parallel plane channel the exact analytical solutions for the fluid velocities are worked out by the normal-mode reduction of the governing equations. The limiting cases of the weak and strong momentum transfer between the flows in the fracture and porous phases are discussed in detail. A comparison to the nonlinear Forchheimer extension of the model proposed recently by Nield and Kuznetsov (Transport Porous Media, 2013) shows that, in the considered parameter range, the nonlinear effect of the Forchheimer drag is negligibly small. Even the simplest zero-momentum transfer solution yields an acceptable approximation.     

Passive and reactive scalar measurements in a transient high-Schmidt-number Rayleigh–Taylor mixing layer

An experimental study of mixing induced by Rayleigh?CTaylor (RT) instability at an Atwood number (A _t ) ~7.5?×?10^?4 and Schmidt number (Sc) ~1,000 has been performed. A new transient experimental facility developed on the working principles of the draw-tank facility at Cambridge (Dalziel et al. in J Fluid Mech, 399:1?C48, 1999) has been established and enhanced to observe a higher (2×) Reynolds number regime. Water and brine were used to produce the RT density stratification. The evolution of the instability was studied using passive and reactive scalar techniques and quantified using optical diagnostic methods. The data were combined to estimate local and global mixing metrics representative of the mixing mechanism across the mixing layer. In comparison with parameters reported from analogous experiments, the mixing phenomenon at a high Sc shows a strong dependency on the initial conditions prevailing at the onset of the instability and the evidence of a delay in the mixing transition. Values of global and integral mixing parameters did not reach late-time asymptotic values that have been reported previously from steady-state experiments (Texas A&M Water Channel) and may be attributed to the effect of the barrier pull and the overturning mechanism that is thought to hinder the progress of the mixing layer.     

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

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.     

Transient director patterns upon flow start-up of nematic liquid crystals (an explanation for stress oscillation damping)

In this work we attempt to determine the origin of damped stress oscillations upon flow start-up of a nematic liquid crystalline monodomain. These damped stress oscillations were first observed by Gu et al. (1993) in the cone-plate flow cell and have since also been observed by Mather et al. (1997) in the parallel disk cell. Although Mather's work explained the cause of the stress oscillation damping in the torsional flow cell, the origin of the damping in the cone-plate device remains a mystery. Here we report finding similar damped stress oscillations in the cylindrical Couette cell and combined with the optical experiments reported earlier by Cladis and Torza (1975, 1976) we are able to propose an explanation for the damping in this geometry. We also report new optical experiments using the cone-plate cell in hopes of determining a cause to the damping in the cone-plate cell. Received: 11 August 2000 Accepted: 31 October 2000