Integrity bases for a symmetric tensor and a vector — The crystal classes

Integrity bases are derived for a symmetric second-order tensor and a vector for the transformation groups corresponding to each of the crystal classes. In each case the irreducibility of the integrity basis is proven.     

On a generalized nonlinearK-ɛ model for turbulence that models relaxation effects

In this paper, based on a similarity that exists between the constitutive relations for turbulent mean flow of a Newtonian fluid and that for the laminar flow of a non-Newtonian fluid, and making use of extended thermodynamics, we develop a generalized nonlinearK- model, whose approximate form includes the standardK- model and the nonlinearK- model of Speziale (1987) as special cases. Our nonlinearK- model, which is frame indifferent, can predict relaxation of the Reynolds stress, unlike most standardK- models. Also, our model is in keeping with that of Yakhotet al. (1992). Most interestingly, the linearized form of our model bears a striking resemblance to the model due to Yoshizawa and Nisizima (1993); however, it has been obtained from a totally different perspective.     

The compressible Coanda wall jet—an experimental study of jet structure and breakaway

An underexpanded jet issuing from a convergent slot and blowing over a surface of convex streamwise curvature was studied experimentally. The jet was confined between side walls, with the slot aspect ratio varying between 40 and 6, but tests showed that in the area of interest close to the slot the flow was effectively two-dimensional. The ratio of slot width to the radius of curvature of the downstream surface varied between 0.05 and 0.33. The main techniques used were Schlieren and shadowgraph to show the jet structure, and surface flow visualization which revealed areas of separation and reattachment. Surface static pressures were also measured on the curved surface. The curved jet proved to have a shock cell structure similar to that of a plane jet. However, the cell structure disappeared more rapidly as the outer shear layer grew more quickly due to the destabilizing effect of the curvature on the turbulence in the shear layer. Even at modest upstream jet pressures, a separated region on the Coanda surface became evident. This region was characterized by a stagnant constant pressure part followed by a region of strongly reversed flow before reattachment took place. The separation was caused by the compression at the end of the first shock cell, with reattachment taking place where expansion in the second cell started. The separated region grew rapidly as the upstream pressure was increased, until, finally, reattachment failed to occur and the jet suddenly broke away from the surface. This work is related to studies of the Coanda flare, where the jet is axisymmetric. The high level of turbulence causes rapid entrainment of air and so gives us clean combustion. However there should be more general application to devices that use the Coanda effect, varying from fluidic devices to blown jet flaps on wings.     

The characteristics of gas–solid flow and wall heat transfer in a fluidized bed reactor

Numerical study using computational fluid dynamics has been carried out to investigate the heat transfer characteristics of a laboratory fluidized bed reactor. The fluidized bed reactor of vTI (Johann Heinrich von Thünen-Institute)-Institute of Wood Technology and Wood Biology is modeled. For the simulation of multiphase flow and thermal fields, an Eulerian–Eulerian approach is applied. The flow and thermal characteristics of the reactor are fully investigated for the wide range of superficial gas velocities and two different particle diameters. In particular, the contributions of the gas bubble and emulsion phase flows on the wall heat transfer are scrutinized. From the predicted results, it is fully elucidated that particular near-wall bubble motions mainly govern the wall heat transfer.     

The Reissner–Sagoci Problem for a Non-homogeneous Half-space with a Surface Constraint

This paper deals with the problem of twisting a non-homogeneous, isotropic, half-space by rotating a circular part of its boundary surface (0 r < a, z = 0) through a given angle. A ring (a < r < b, z = 0) outside the circle is stress-free and the remaining part (r > b, z = 0) is rigidly clamped. The shear modulus is assumed to vary with the cylindrical coordinates, r, z by the relation (z) = ₁(c + z), c 0 where ₁, c and are real constants. Expressions for some quantities of physical importance, such as torque applied at the surface of the disk and stress intensity factors, are obtained. The effects of non-homogeneity on torque and stress intensity factor are illustrated graphically.     

The effect of sphere–wall interactions on particle motion in a viscoelastic suspension of FENE dumbbells

A new method for simulating the motion of particles in viscoelastic Boger fluids is extended to problems with bounded geometries. Viscoelasticity is incorporated into the Stokesian dynamics method by modeling a viscoelastic fluid as a suspension of finite-extension nonlinear-elastic (FENE) dumbbells. Wall–particle and wall–bead interactions are included by using the image system method of Blake; particle–particle and particle–bead interactions are also modified by the presence of the wall. The method of incorporating sphere–wall interactions is verified by doing calculations for several problems involving particle–wall interactions in Newtonian fluids. The method is then used to study particle–wall interactions in viscoelastic dumbbell suspensions by examining several problems of interest: the sedimentation of a spherical particle near vertical and tilted walls; the sedimentation of a nonspherical particle between two flat plates; and the migration of a neutrally buoyant sphere in plane Poiseuille flow. We find that a single sphere falling near a wall moves toward the wall and exhibits anomalous rotation. When the wall is tilted by an amount less than a few degrees, the sphere still moves toward the wall, but tilting the wall greater than an angle of approximately 1.5° results in the sphere falling away from the wall. A nonspherical particle settling in a channel exhibits an oscillatory motion, but ultimately becomes centered in the channel with its long axis parallel to gravity. Finally, it is shown that a neutrally buoyant sphere in plane Poiseuille flow migrates to the channel center in wide channels, but migrates to the walls when the sphere is sufficiently large relative to the channel width.     

The onset of Görtler vortices in laminar boundary layer flow over a slightly concave wall

The onset of convective instability in the laminar boundary layer over the slightly curved wall is analyzed theoretically and compared with the existing experimental data. A new set of stability equations are derived by the propagation theory considering the relative instability under the linear stability theory. In this analysis the disturbances are assumed to have the form of longitudinal vortices and also to grow themselves in streamwise direction. The critical position to mark the onset of Görtler instability is obtained as a function of the Görtler number, where disturbances at the critical state are mainly confined to the hydrodynamic boundary layer. Comparing the theoretical predictions with available experimental and other theoretical results, the present predictions follow experimental trends fairly well with slightly higher critical Görtler numbers than those from the local stability theory. The propagation theory commanding the local eigenvalue analysis is successful to obtain stability conditions reasonably in Görtler vortex problems, relaxing the limitations by the conventional analyses.     

The Modulational Instability for a Generalized Korteweg–de Vries Equation

We study the spectral stability of a family of periodic standing wave solutions to the generalized Korteweg–de Vries in a neighborhood of the origin in the spectral plane using what amounts to a rigorous Whitham modulation theory calculation. In particular we are interested in understanding the role played by the null directions of the linearized operator in the stability of the traveling wave to perturbations of long wavelength. A study of the normal form of the characteristic polynomial of the monodromy map (the periodic Evans function) in a neighborhood of the origin in the spectral plane leads to two different instability indices. The first, an orientation index, counts modulo 2 the total number of periodic eigenvalues on the real axis. This index is a generalization of the one which governs the stability of the solitary wave. The second, a modulational instability index, provides a necessary and sufficient condition for the existence of a long-wavelength instability. This index is essentially the quantity calculated by Hǎrǎguş and Kapitula in the small amplitude limit. Both of these quantities can be expressed in terms of the map between the constants of integration for the ordinary differential equation defining the traveling waves and the conserved quantities of the partial differential equation. These two indices together provide a good deal of information about the spectrum of the linearized operator. We sketch the connection of this calculation to a study of the linearized operator—in particular we perform a perturbation calculation in terms of the Floquet parameter. This suggests a geometric interpretation attached to the vanishing of the orientation index previously mentioned.     

A solution of the E-ε model for nearly homogeneous turbulence with a mean shear

An analytical solution of the E- model for the downstream evolution of a stationary and nearly homogeneous turbulent shear flow is presented. In case that the turbulent time scale has adjusted itself to the time scale imposed by the shear, an asymptotic solution can be derived from the full solution, which shows that both E and increase downstream exponentially. By comparing this asymptotic solution with experimental data a value for the unknown constant c _l , in the -equation, is derived. Moreover, we find an expression for the downstream development of the variance of a scalar, which is also compared with experimental data. The analytical solution shows that a homogeneous shear flow with a uniform velocity gradient can only be obtained if the shear is sufficiently small. In the experiments this condition is not always satisfied. A discussion is given on how a nearly homogeneous shear flow can be obtained over a limited downstream interval by changing the initial conditions in E and , and a comparison is made with experimental data. Finally it is shown that better transverse homogeneity can be obtained by taking an exponential velocity profile instead of a linear profile.     

An adaptive finite element method for viscoplastic flows in a square pipe with stick–slip at the wall

This paper presents the numerical solution of non-linear yield stress phenomena by using a new mixed anisotropic auto-adaptive finite element method. The Poiseuille flow of a Bingham fluid with slip yield boundary condition at the wall is considered. Despite its practical interest, for instance for pipeline flows of yield-stress fluids such as concrete and cements, this problem has not been addressed yet to our knowledge. The case of a pipe with a square section has been investigated in detail. The computations cover the full range of the two main dimensionless numbers and exhibit complex flow patterns: all the different flow regimes are completely identified.     

The cross sectional area difference method — a new technique for determination of particle concentration by laser doppler anemometry

A new technique for the determination of particle concentration from the signals of a laser Doppler anemometer (LDA) is described. It is based on a statistical relation between the number of Doppler periods, or the amplitude of the Doppler signals, and the particle concentration. The technique allows the mass flux of the dispersed phase of a two-phase flow to be obtained from the data set of a conventional one-dimensional (ID) LDA. The technique has been called the cross sectional area difference method. Simulations and first experimental results are presented and discussed.List of symbols a, b, c half-axes of measurement control volume (mcv) - a ₁, b ₁, c ₁ half-axes of detection volume - c _L velocity of light - d _m beam waist diameter - d _p particle diameter - d _pc diameter of the calibration particle - d _pmin minimum detectable particle diameter - e elementary charge - h Planck's constant - i number of particle size classes - k wavenumber - m visibility - m refractive index - n(d _p ) particle concentration - n(d _pi ) concentration of ith particle class - n vector of n(d _pi ) - q exponent of size dependence of G(d _p ) - v _x x-velocity component - x fringe spacing - y ₀, z ₀ coordinates of particle trajectory and cross sectional area - A cross sectional area of mcv - A matrix of A ₁ - a ₁ cross sectional area of detection volume - A ₁ difference of neighbouring cross sectional areas - C _A normalisation constant for linear graduation of amplitude - C _N normalisation constant for Doppler periods - C _scat non-size-dependent factor of G(d _p ) - C _x normalisation constant for nonlinear graduation of amplitude - F() power spectral density - G(d _p ) integral scattering function - H number of accumulated counts - H _max maximum number of accumulated counts - I amplitude of Doppler signal - I _max I for a particle passing through the origin of the mcv - I _s trigger level - K logarithmic amplitude ratio - K _max logarithmic amplitude ratio for I _s - K _x degree of linear class width of amplitude - K _A degree of nonlinear class width of amplitude - N number of Doppler periods - N _m number of Doppler periods required by signal validation - N _max N for a particle passing through the origin of the mcv - N ₀ fringe number inside mcv along x-axis - P _L laser power - S ₀ particle arrival rate - S ₁ trigger rate - S ₁ contribution to trigger rate coming from A ₁ - S ₁ vector of S _1i - S _1i contribution to trigger rate coming from ith class of distribution - _Q quantum efficiency - wavelength of laser light - off-axis angle - elevation angle - angular frequency - beam intersection angle - phase difference     

A reformulated synthetic turbulence generation method for a zonal RANS–LES method and its application to zero-pressure gradient boundary layers

A synthetic turbulence generation (STG) method for subsonic and supersonic flows at low and moderate Reynolds numbers to provide inflow distributions of zonal Reynolds-averaged Navier–Stokes (RANS) – large-eddy simulation (LES) methods is presented. The STG method splits the LES inflow region into three planes where a local velocity signal is decomposed from the turbulent flow properties of the upstream RANS solution. Based on the wall-normal position and the local flow Reynolds number, specific length and velocity scales with different vorticity content are imposed at the inlet plane of the boundary layer. The quality of the STG method for incompressible and compressible zero-pressure gradient boundary layers is shown by comparing the zonal RANS–LES data with pure LES, pure RANS, and direct numerical simulation (DNS) solutions. The distributions of the time and spanwise wall-shear stress, Reynolds stress distributions, and two point correlations of the zonal RANS–LES simulations are smooth in the transition region and in good agreement with the pure LES and reference DNS findings. The STG approach reduces the RANS-to-LES transition length to less than four boundary-layer thicknesses.     

Analytical solution for a three-dimensional non-homogeneous bivariate population balance equation—a special case

There has been a dramatic increase in the number of research publications using the population balance equation (PBE). The PBE allows the prediction of the spatial distribution of the dispersed phase size for an accurate estimation of the flow fields in multiphase flows. A few recent studies have proposed new efficient numerical methods to solve non-homogeneous multivariate PBE and implemented the same in computational fluid dynamics (CFD) codes. However, these codes are generally benchmarked against other numerical methods and applied without verification. To address this gap, an analytical solution for a three-dimensional non-homogeneous bivariate PBE is presented here for the first time. The method of manufactured solutions (MMS) has been used to construct a solution of the non-homogeneous PBE containing breakage and coalescence terms, and an additional source term appearing as a result of this method. The analytical solution presented in this work can be used for the rigorous verification of computer codes written to solve the non-homogeneous bivariate PBE. Quantification of the errors due to different numerical schemes will also become possible with the availability of this analytical solution for the PBE.