Coupling time-stepping numerical methods and standard aerodynamics codes for instability analysis of flows in complex geometries

The development of a global instability analysis code coupling a time-stepping approach, as applied to the solution of BiGlobal and TriGlobal instability analysis 1, 2 and finite-volume-based spatial discretization, as used in standard aerodynamics codes is presented. The key advantage of the time-stepping method over matrix-formulation approaches is that the former provides a solution to the computer-storage issues associated with the latter methodology. To-date both approaches are successfully in use to analyze instability in complex geometries, although their relative advantages have never been quantified. The ultimate goal of the present work is to address this issue in the context of spatial discretization schemes typically used in industry. The time-stepping approach of Chiba 3 has been implemented in conjunction with two direct numerical simulation algorithms, one based on the typically-used in this context high-order method and another based on low-order methods representative of those in common use in industry. The two codes have been validated with solutions of the BiGlobal EVP and it has been showed that small errors in the base flow do not have affect significantly the results. As a result, a three-dimensional compressible unsteady second-order code for global linear stability has been successfully developed based on finite-volume spatial discretization and time-stepping method with the ability to study complex geometries by means of unstructured and hybrid meshes

Wave-like Disturbances on the Downstream Wall of an Open Cavity

This contribution presents results of an incompressible two-dimensional flow over an open cavity of fixed aspect ratio (length/depth) L/D = 2 and the coupling between the three dimensional low frequency oscillation mode confined in the cavity and the wave-like disturbances evolving on the downstream wall of the cavity in the form of Tollmien-Schlichting waves. BiGlobal instability analysis is conducted to search the global disturbances superimposed upon a two-dimensional steady basic flow. The base solution is computed by the integration of the laminar Navier-Stokes equations in primitive variable formulation, while the eigenvalue problem (EVP) derived from the discretization of the linearized equations of motion in the BiGlobal framework is solved using an iterative procedure. The formulation of the BiGlobal EVP for the unbounded flow in the open cavity problem introduces additional difficulties regarding the flow-through boundaries. Local analysis has been utilized for the determination of the proper boundary conditions in the upper limit of the downstream region

Global Linear Instability at the Dawn of its 4th Decade: A List of Challenges (A Practical Guide on how to Contain the Euphoria and Avoid the Oversell)

Global linear instability theory is concerned with the temporal or spatial development of small-amplitude perturbations superposed upon laminar steady or time-periodic threedimensional flows, which are inhomogeneous in two (and periodic in one) or all three spatial directions.1 The theory addresses flows developing in complex geometries, in which the parallel or weakly nonparallel basic flow approximation invoked by classic linear stability theory does not hold. As such, global linear theory is called to fill the gap in research into stability and transition in flows over or through complex geometries. Historically, global linear instability has been (and still is) concerned with solution of multi-dimensional eigenvalue problems; the maturing of non-modal linear instability ideas in simple parallel flows during the last decade of last century2–4 has given rise to investigation of transient growth scenarios in an ever increasing variety of complex flows. After a brief exposition of the theory, connections are sought with established approaches for structure identification in flows, such as the proper orthogonal decomposition and topology theory in the laminar regime and the open areas for future research, mainly concerning turbulent and three-dimensional flows, are highlighted. Recent results obtained in our group are reported in both the time-stepping and the matrix-forming approaches to global linear theory. In the first context, progress has been made in implementing a Jacobian-Free Newton Krylov method into a standard finite-volume aerodynamic code, such that global linear instability results may now be obtained in compressible flows of aeronautical interest. In the second context a new stable very high-order finite difference method is implemented for the spatial discretization of the operators describing the spatial BiGlobal EVP, PSE-3D and the TriGlobal EVP; combined with sparse matrix treatment, all these problems may now be solved on standard desktop computers.

An efficient nonlinear transformation for the numerical computation of the singular integrals appearing in the 2-D Boundary Element Method

We discuss several methods, based on coordinate transformations, for the evaluation of singular and quasisingular integrals in the direct Boundary Element Method. An intrinsec error of some of these methods is detected. Two new transformations are suggested which improve on those currently available.

Physics and applications of electrodynamic space tethers

Basic effects and dynamical and electrical contact issues in the physics of (electrodynamic space) bare tethers are discussed. Scientific experiments and powerpropulsion applications, including a paradoxical use of bare tethers in outer-planet exploration,are considered.

Wind tunnel analysis on the influence of cantilever parapets on the wind loads on curved roofs

Different methods to reduce the high suction caused by conical vortices have been reported in the literature: vertical parapets, either solids or porous, placed at the roof edges being the most analysed configuration. Another method for alleviating the high suction peaks due to conical vortices is to round the roof edges. Very recently, the use of some non-standard parapet configurations, like cantilever parapets, has been suggested. In this paper, its efficiency to reduce suction loads on curved roofs is experimentally checked by testing the pressure distribution on the curved roof of a low-rise building model in a wind tunnel. Very high suction loads have been measured on this model, the magnitude of these high suction loads being significantly decreased when cantilever...

Experimental analysis of a photovoltaic concentrator based on a single reflective stage immersed in an optical fluid

This article reviews all the experimental tests carried out to analyze the performance of a FluidReflex photovoltaic concentrator. This novel concentrator concept consists of a single reflective stage immersed in an optical fluid. The presence of the fluid entails significant advantages. It not only allows a high system optical efficiency and increases the attainable concentration but also enhances the heat dissipation from the cell. In addition, the electrical insulation is improved, and the problem of water vapor condensation inside the module is avoided. A complete characterization is addressed in this paper. Among the experimental results, a measured optical efficiency of 83.5% for a concentration of 1035× stands out

Dynamics of thermal ignition of spray flames in mixing layers

Conditions are identified under which analyses of laminar mixing layers can shed light on aspects of turbulent spray combustion. With this in mind, laminar spray-combustion models are formulated for both non-premixed and partially premixed systems. The laminar mixing layer separating a hot-air stream from a monodisperse spray carried by either an inert gas or air is investigated numerically and analytically in an effort to increase understanding of the ignition process leading to stabilization of high-speed spray combustion. The problem is formulated in an Eulerian framework, with the conservation equations written in the boundary-layer approximation and with a one-step Arrhenius model adopted for the chemistry description. The numerical integrations unveil two different types of ignition behaviour depending on the fuel availability in the reaction kernel, which in turn depends on the rates of droplet vaporization and fuel-vapour diffusion. When sufficient fuel is available near the hot boundary, as occurs when the thermochemical properties of heptane are employed for the fuel in the integrations, combustion is established through a precipitous temperature increase at a well-defined thermal-runaway location, a phenomenon that is amenable to a theoretical analysis based on activation-energy asymptotics, presented here, following earlier ideas developed in describing unsteady gaseous ignition in mixing layers. By way of contrast, when the amount of fuel vapour reaching the hot boundary is small, as is observed in the computations employing the thermochemical properties of methanol, the incipient chemical reaction gives rise to a slowly developing lean deflagration that consumes the available fuel as it propagates across the mixing layer towards the spray. The flame structure that develops downstream from the ignition point depends on the fuel considered and also on the spray carrier gas, with fuel sprays carried by air displaying either a lean deflagration bounding a region of distributed reaction or a distinct double-flame structure with a rich premixed flame on the spray side and a diffusion flame on the air side. Results are calculated for the distributions of mixture fraction and scalar dissipation rate across the mixing layer that reveal complexities that serve to identify differences between spray-flamelet and gaseous-flamelet problems.

The breaking of axisymmetric slender liquid bridges

Liquids held by surface tension forces can bridge the gap between two solid bodies placed not too far apart from each other. The equilibrium conditions and stability criteria for static, cylindrical liquid bridges are well known. However, the behaviour of an unstable liquid bridge, regarding both its transition toward breaking and the resulting configuration, is a matter for discussion. The dynamical problem of axisymmetric rupture of a long liquid bridge anchored at two equal coaxial disks is treated in this paper through the adoption of one-dimensional theories which are widely used in capillary jet problems

One-dimensional self-similar solution of the dynamics of axisymmetric slender liquid bridges

In this paper the dynamics of axisymmetric, slender, viscous liquid bridges having volume close to the cylindrical one, and subjected to a small gravitational field parallel to the axis of the liquid bridge, is considered within the context of one-dimensional theories. Although the dynamics of liquid bridges has been treated through a numerical analysis in the inviscid case, numerical methods become inappropriate to study configurations close to the static stability limit because the evolution time, and thence the computing time, increases excessively. To avoid this difficulty, the problem of the evolution of these liquid bridges has been attacked through a nonlinear analysis based on the singular perturbation method and, whenever possible, the results obtained are compared with the numerical ones.

One-dimensional linear analysis of the compound jet

The stability of an infinitely long compound liquid column is analysed by using a one-dimensional inviscid slice model. Results obtained from this one-dimensional linear analysis are applicable to the study of compound capillary jets, which are used in the ink-jet printing technique. Stability limits and the breaking regimes of such fluid configurations are established, and, whenever possible, theoretical results are compared with experimental ones.

Axisymmetric long liquid bridges in a time-dependent microgravity field

This paper deals with the dynamics of liquid bridges when subjected to an oscillatory microgravity field. The analysis has been performed by using a one-dimensional slice model, already used in liquid bridge problems, which allows to calculate not only the resonance frequencies of a wide range of such fluid configurations but also the dependence of the dynamic response of the liquid bridge on the frequency on the imposed perturbations. Theoretical results are compared with experimental ones obtained aboard Spacelab-Dl, the agreement between theoretical and experimental results being satisfactory

Influence of a liquid on the natural frequencies of almost circular plates

In this work, the influence of the surrounding fluid on the dynamic characteristics of almost circular plates is investigated. First the natural frequencies and normal modes for the plates in vacuum are calculated by a perturbation procedure. The method is applied for the case of elliptical plates with a low value of eccentricity. The results are compared with other available methods for this type of plates with good agreement. Next, the effect of the fluid is considered. The normal modes of the plate in vacuum are used as a base to express the vibration mode of the coupled plate-fluid system. By applying the Hankel transformation the nondimensional added virtual mass 2 increment (NAVMI) are calculated for elliptical plates. Results of the NAVMI factors and the effect of the fluid on the natural frequencies are given and it is shown that when the eccentricity of the plate is reduced to zero (circular plate) the known results of the natural frequencies for circular plates surrounded by liquid are recovered.

Dynamics of gas bubbles encapsulated by a viscoelastic fluid shell under acoustic fields

The dynamics of a gas-filled microbubble encapsulated by a viscoelastic fluid shell immersed in a Newtonian liquid and subject to an external pressure field is theoretically studied. The problem is formulated by considering a nonlinear Oldroyd type constitutive equation to model the rheological behavior of the fluid shell. Heat and mass transfer across the surface bubble have been neglected but radiation losses due to the compressibility of the surrounding liquid have been taken into account. Bubble collapse under sudden increase of the external pressure as well as nonlinear radial oscillations under ultrasound fields are investigated. The numerical results obtained show that the elasticity of the fluid coating intensifies oscillatory collapse and produces a strong increase of the amplitudes of radial oscillations which may become chaotic even for moderate driving pressure amplitudes. The role played by the elongational viscosity has also been analyzed and its influence on both, bubble collapse and radial oscillations, has been recognized. According to the theoretical predictions provided in the present work, a microbubble coated by a viscoelastic fluid shell is an oscillating system that, under acoustic driving, may experience volume oscillations of large amplitude, being, however, more stable than a free bubble. Thus, it could be expected that such a system may have a suitable behavior as an echogenic agent.

Numerical issues in Lagrangian tracking and topological evolution of fluid particles in wallbounded turbulent flows

The determination of the local Lagrangian evolution of the flow topology in wall-bounded turbulence, and of the Lagrangian evolution associated with entrainment across the turbulent / non-turbulent interface into a turbulent boundary layer, require accurate tracking of a fluid particle and its local velocity gradients. This paper addresses the implementation of fluid-particle tracking in both a turbulent boundary layer direct numerical simulation and in a fully developed channel flow simulation. Determination of the sub-grid particle velocity is performed using both cubic B-spline, four-point Hermite spline and higher-order Hermite spline interpolation. Both wall-bounded flows show similar oscillations in the Lagrangian tracers of both velocity and velocity gradients, corresponding to the movement of particles across the boundaries of computational cells. While these oscillation in the particle velocity are relatively small and have negligible effect on the particle trajectories for time-steps of the order of CFL = 0.1, they appear to be the cause of significant oscillations in the evolution of the invariants of the velocity gradient tensor.