Experimental investigation of the grid-generated turbulence features in a free surface flow

The present study aims to investigate the features of a grid-generated turbulence occurring in a current flow with a free surface flow. The interest is focused on the length and time scales of the turbulence. These are the macro, the micro and the Kolmogorov scales. To analyze the flow, a 2D LDV system has been used to measure , , u′ and w′. This non-intrusive and optical technique is really accurate (in terms of space and time resolution). Furthermore, it does not disturb the flow and provides a high data rate. Both horizontal and vertical velocities are recorded at the same time according to a coincidence window (τ_cw). Bias measurements are avoided by using a filtering technique during data processing. The improved homogeneity and isotropy of the turbulence downstream of the grid allows the use of the Taylor hypothesis. Thus, all length and time scales of the flow can be estimated. Results are discussed as well as the influence of the upcoming mean velocity on the turbulence properties.     

Variational formulation of ideal fluid flows according to gauge principle

On the basis of the gauge principle of field theory, a new variational formulation is presented for flows of an ideal fluid. The fluid is defined thermodynamically by mass density and entropy density, and its flow fields are characterized by symmetries of translation and rotation. The rotational transformations are regarded as gauge transformations as well as the translational ones. In addition to the Lagrangians representing the translation symmetry, a structure of rotation symmetry is equipped with a Lagrangian Λ_A including the vorticity and a vector potential bilinearly. Euler's equation of motion is derived from variations according to the action principle. In addition, the equations of continuity and entropy are derived from the variations. Equations of conserved currents are deduced as the Noether theorem in the space of Lagrangian coordinate a. Without Λ_A, the action principle results in the Clebsch solution with vanishing helicity. The Lagrangian Λ_A yields non-vanishing vorticity and provides a source term of non-vanishing helicity. The vorticity equation is derived as an equation of the gauge field, and the Λ_A characterizes topology of the field. The present formulation is comprehensive and provides a consistent basis for a unique transformation between the Lagrangian a space and the Eulerian x space. In contrast, with translation symmetry alone, there is an arbitrariness in the transformation between these spaces.     

A study of electromechanical switching in ferroelectric single crystals

This article documents both modeling and experimental studies developed to investigate the switching behavior of ferroelectric single crystals. The theoretical model makes a priori ansatz that switching follows the evolution of a particular domain pattern. The choice of this configuration is dictated by the requirement that domains remain compatible during evolution, giving rise to a low-energy path for the overall switching. The construction of this pattern is achieved using multirank laminates. It offers an advantage of specifying different types of domain wall movements, leading to a distinction for the switching types. A loading experiment is performed on a barium titanate (BaTiO₃) single crystal with a constant compressive stress and a cyclic electric field. Both 180 and 90 coercive fields are measured as input parameters required for the theoretical framework. The simulation results show good agreement with the observed strains measured by the present and other available experiments. It is found that depolarization has a non-trivial influence on attainable actuation strains.     

Analysis for a screw dislocation accelerating through the shear-wave speed barrier

An analysis is performed for an accelerating screw dislocation through the shear-wave speed barrier. At this instant, the function that determines the interval of the path of the dislocation motion that contributes to the wave front has roots that change from a pair of complex conjugate to a double real, which subsequently splits into two real ones. The analysis is performed at this transition to supersonic that occurs at the double root maximum of the function that defines the interval of the dislocation path that contributes to the field points. It is found that the stress has a log|ξ-ξ^*|/|ξ-ξ^*|^1/2 singularity in the coefficient of the delta function of the forming Mach front, implying that for this phenomenon the Volterra dislocation model has too strong a discontinuity (step-function) in the displacement to be meaningful. A ramp-core displacement dislocation model analysis, which removes the singularity in the stress, is presented. These results can be useful in a multiscale dislocation dynamics modeling with inertia effects.     

Heat transfer from a line source located in the periodic laminar near wake of a circular cylinder

This experimental study is devoted to the diffusion of a passive scalar downstream a line source located in a Bénard–von Kármán street. Measurements of velocity and temperature have been performed using LDA and cold wire thermometer with a phase reference. Information on the initial evolution of mean, fluctuating velocity and temperature and associated shear-stresses and heat fluxes fields are presented for two locations of the source: ( and 1). The results show that the velocity field in the wake is strongly related to the geometric structure of vortices while the temperature field is controlled by both the time scale of rotation of the vortices and the location of the heated fluid within the vortex street.     

Separation of variables in one problem of motion of the generalized Kowalevski top

In the problem of motion of the Kowalevski top in a double force field the new case of reduction to a Hamiltonian system with two degrees of freedom was pointed out by Kharlamov [Kharlamov, M.P., 2004. Mekh. Tverd. Tela 34, 47–58]. We show that the equations of motion in this case can be separated by the appropriate change of variables, the new variables U,V being hyperelliptic functions of time. The natural phase variables (components of the angular velocity and the direction vectors of the forces with respect to the movable basis) are expressed via U,V explicitly in elementary algebraic functions.     

Numerical simulation of the noise generated by a low Mach number, low Reynolds number jet

In this paper we study the sound field produced by a turbulent round jet with a Mach number of 0.6 based on the centerline velocity and the ambient speed of sound c_∞. The turbulent flow field is found by solving the fully compressible Navier–Stokes equations with help of high-order compact finite difference schemes. It is shown that the simulated flow field is in good agreement with experiments. The corresponding sound field has been obtained with help of the Lighthill equation using two different formulations for the Lighthill stress tensor T_ij. In the first formulation of T_ij the fluctuating density is taken into account. In the second formulation the density is assumed to be constant. As an additional check we have also performed an acoustic calculation using a formulation in which a homogeneous wave equation is solved. The boundary conditions for this homogeneous wave equation are obtained from the numerical simulation of the Navier–Stokes equation. The results obtained with both formulations of the Lighthill stress tensor are nearly identical. This implies that an incompressible formulation of the conservations laws could be used to predict jet noise at low Mach numbers.     

Experimental study of liquid sloshing dynamics in a barge carrying tank

An experimental work has been carried out to study the phenomena of sloshing of liquid in partially filled tanks mounted on a barge exposed to regular beam waves. Three liquid fill levels with liquid depth, h_s to length of tank, l ratio (h_s/l) of 0.163, 0.325 and 0.488, are studied. The time histories of sloshing oscillation are measured along the length of container at predefined locations. The nonlinear behaviour of sloshing oscillation is observed for the regular wave excitation. The spectra of the sloshing oscillation and their qualitative assessment are reported. Attempts are made to evaluate the harmonics present in the sloshing oscillation and compare with the results of earlier studies. The effects of wave excitation frequency and wave height on the sloshing oscillation as well as on the response of the barge are studied.     

Perturbation theory for the double sine-Gordon equation

This paper presents a perturbation theory for the double sine-Gordon equation. We obtain a system of differential equations that shows the soliton parameters modification under the influence of the perturbation. For the particular case λ=0 the results transform into the well-known perturbation theory for the sine-Gordon equation.     

On the flow in a uniformly porous pipe

This article considers fully laminar flow of an incompressible viscous fluid in a uniformly porous pipe with suction and injection. An exact solution of the Navier–Stokes equations is given. The velocity filed can be expressed in a series form in terms of the modified Bessel function of the first kind of order n. The volume flux across a plane normal to the flow, the vorticity and the stress on the boundary are presented. The flow properties depend on the cross-Reynolds number, Ua/ν, where U is the suction velocity, a is the radius of the pipe and ν is the kinematic viscosity of the fluid. It is found that for large values of the cross-Reynolds number, the flow near the region of the suction shows a boundary layer character. In this region the velocity and the vorticity vary sharply. Outside the boundary layer, the velocity and the vorticity do not show an appreciable change.     

Direct numerical simulations of two-layer viscosity-stratified flow

Two-dimensional simulations of flow instability at the interface of a two-layer, density-matched, viscosity-stratified Poiseuille flow are performed using a front-tracking/finite difference method. We present results for the small-amplitude (linear) growth rate of the instability at small to medium Reynolds number for varying thickness ratio n, viscosity ratio m, and wavenumber. We also present results for large-amplitude non-linear evolution of the interface for varying viscosity ratio and interfacial tension. For the linear case, the interfacial mode is neutrally stable for as predicted by analysis. The growth rate is proportional to Reynolds number for small Re, and increases with viscosity ratio. The growth rate also increases when the thickness of the more viscous layer is reduced. Strong non-linear behavior is observed for relatively large initial perturbation amplitude. The higher viscosity fluid is drawn out as a finger that penetrates into the lower viscosity layer. The simulated interface shape compares well with previously reported experiments. Increasing interfacial tension retards the growth rate of the interface as expected, whereas increasing the viscosity ratio enhances it. Drop formation at the small Reynolds number considered in this study is precluded by the two-dimensional nature of the calculations.     

Autoparametric quasiperiodic excitation

The dynamics of an autonomous conservative three degree of freedom system which exhibits autoparametric quasiperiodic excitation is investigated. The system is a generalization of a classical system known as the “particle in the plane”. The system exhibits a motion, the z=0 mode, whose stability is governed by a linear second order ODE with quasiperiodic coefficients. The behavior of the latter ODE is studied by using three different methods: numerical integration, harmonic balance and perturbation methods.     

Shroud heat transfer measurements inside a heated multiple rotating cavity with axial throughflow

Experimental measurements of heat transfer are made from the inner peripheral surface of a rotating test rig designed to be similar to a gas turbine high pressure compressor internal air system. The test rig comprises a number of annular discs sealed at their periphery by a shroud. An axial throughflow of cooling air enters the test rig and flows through the annular section between the disc bores and a central shaft. Tests were carried out for the following range of rotational speeds and axial throughflow rates: 540 < N_R < 10,800 rev/min and (corresponding to the range of rotational and axial Reynolds numbers 4 × 10⁵ < Re < 7.7 × 10⁶ and 3.3 × 10⁴ < Re_z < 2.2 × 10⁵).

The shroud Nusselt numbers are found to depend on the shroud Grashof number. They are relatively insensitive to changes in axial Reynolds number and two geometrically similar cavities give similar values of Nusselt number. The heat transfer from the shroud is governed by the mechanism of free convection. It is recommended that a modified form of a correlation for Rayleigh–Bénard convection in a gravitational force field be used, with appropriate modification, to predict shroud heat transfer.     

A complete and irreducible dynamic SGS heat-flux modelling based on the strain rate tensor for large-eddy simulation of thermal convection

In this paper, a general family of explicit algebraic tensor diffusivity functions based on the resolved temperature gradient vector and strain rate tensor is studied and applied to the construction of new constitutive relations for modelling the subgrid-scale (SGS) heat flux (HF). Based on Noll’s formulation, dynamic linear and nonlinear tensor diffusivity models are proposed for large-eddy simulation of thermal convection. The constitutive relations for these two proposed models are complete and irreducible. These two new models include several existing dynamic SGS HF models as special cases. It is shown that in contrast to the conventional modelling approach, the proposed models embody more degrees of freedom, permit non-alignment between the SGS HF and resolved temperature gradient vectors, reflect near-wall flow physics at the subgrid scale, and therefore, allow for a more realistic geometrical representation of the SGS heat flux for large-eddy simulation of thermal convection. Numerical simulations have been performed using a benchmark test case of a combined forced and natural convective flow in a vertical channel with a Reynolds number of and a Grashof number of Gr = 9.6 × 10⁵. The results obtained using the two proposed SGS HF models are compared with reported direct numerical simulation (DNS) data as well as predictions obtained using several conventional dynamic SGS HF models.     

Linear and nonlinear sloshing in a circular conical tank

Linear and nonlinear fluid sloshing problems in a circular conical tank are studied in a curvilinear coordinate system. The linear sloshing modes are approximated by a series of the solid spheric harmonics. These modes are used to derive a new nonlinear modal theory based on the Moiseyev asymptotics. The theory makes it possible to both classify steady-state waves occurring due to horizontal resonant excitation and visualise nonlinear wave patterns. Secondary (internal) resonances and shallow fluid sloshing (predicted for the semi-apex angles >60) are extensively discussed.     

Wave scattering from a rough surface: Solution by an iterative method

Two iterative solutions of the Helmholtz equation for a scalar field in above a rough surface that admits the Dirichlet boundary condition are derived. The bases for the two iterative methods are two different boundary integral equations that represent the solution. The first integral equation is classified as a Fredholm integral equation of the first kind. The second is classified as a Fredholm integral equation of the second kind. This classification suggests that it is easier to find stable solution methods to the second equation. In both methods, the boundary integral was separated into a major part which is easy to calculate and a local residual part. The major part is a convolution and thus can be calculated using FFT in complexity O(N log N), where N is the number of surface points in which the surface height and its first derivatives together with the incoming wave and its normal derivative are all known. The residual element of the equations can be approximated efficiently only for surfaces where their amplitude is less than the wavelength of the incoming wave. The iterative schemes were tested numerically against a reference solution in order to examine the applicability range, the error estimation and the stability of the schemes. All tests supported the superiority of the second method. In particular the error estimation and stability tests indicated good performance for surfaces with slope up to 1. Yet, being an equation in the scattered field alone, makes the first method useful as a benchmark solution in its domain of applicability.     

A secant-viscosity composite model for the strain-rate sensitivity of nanocrystalline materials

In order to address the strain-rate sensitivity of nanocrystalline solids, a secant-viscosity composite model is developed in this article. The microgeometry of the composite is taken to consist of the grain-interior phase and the grain-boundary affected zone (GBAZ) as suggested by Schwaiger et al. [Schwaiger, R., Moser, B., Dao, M., Chollacoop, N., Suresh, S., 2003. Some critical experiments on the strain-rate sensitivity of nanocrystalline nickel. Acta Mater. 51, 5159–5172], while the constituent properties are modeled by a unified viscoplastic constitutive law. The drag stress of the grain interior is assumed to follow the Hall–Petch relation, but that of the GBAZ is independent of grain size, d. Then in terms of the secant viscosity of the constituent phases, the strain-rate sensitivity of the nanocrystalline solid is determined with the help of a linear viscous comparison composite and a field-fluctuation approach. To test the applicability of the developed model, it is applied to predict the strain-rate effect of a nanocrystalline Ni, and the grain-size dependence of its stress–strain relations. Our theoretical calculations indicate that the tensile strength of a nanocrystalline Ni with d = 40 nm is about five times that of a microcrystalline one with d = 10 μm under the same strain rate of , and that the nanocrystalline Ni exhibits a much stronger strain-rate effect. These predictions are found to be consistent with the experimental data of Schwaiger et al. Possible grain-size softening with further grain-size reduction such as reported in molecular dynamic simulations is also demonstrated.     

Predicting transition in two- and three-dimensional separated flows

This paper is concerned with the numerical prediction of two- and three-dimensional transitional separated flows of turbomachinery interest. The recently proposed single-point transition model based on the use of a laminar kinetic energy transport equation is considered, insofar as it does not require to evaluate any integral parameter, such as boundary-layer thickness, and is thus directly applicable to three-dimensional flows. A well established model, combining a transition-onset correlation with an intermittency transport equation, is also used for comparison. Both models are implemented within a Reynolds-averaged Navier–Stokes solver employing a low-Reynolds-number k–ω turbulence model. The performance of the transition models have been evaluated and tested versus well-documented incompressible flows past a flat plate with semi-circular leading edge, namely: tests T3L2, T3L3, T3L5, and T3LA1 of ERCOFTAC, with different Reynolds numbers and free-stream conditions, the last one being characterized by a non-zero pressure gradient. In all computations, the first model has proven as adequate as or superior to the second one and has been then applied with success to two more complex test cases, for which detailed experimental data are available in the literature, namely: the two- and three-dimensional flows through the T106 linear turbine cascade.     

On the stability of solutions of a non-linear, non-autonomous equation

Sufficient conditions for the stability of solutions of the equation

Full-size image

are deduced. These conditions depend on coefficient ratios a/c and b/c as well as on initial conditions. As an example the pendulum of variable length is discussed.