All-regime Lagrangian-Remap numerical schemes for the gas dynamics equations. Applications to the large friction and low Mach coefficients

In this talk, we propose an all regime Lagrange-Projection like numerical scheme for the gas dynamics equations. By all regime, we mean that the numerical scheme is able to compute accurate approximate solutions with an under-resolved discretization with respect to the Mach number M, i.e. such that the ratio between the Mach number M and the mesh size or the time step is small with respect to 1. The key idea is to decouple acoustic and transport phenomenon and then alter the numerical flux in the acoustic approximation to obtain a uniform truncation error in term of M. This modified scheme is conservative and endowed with good stability properties with respect to the positivity of the density and the internal energy. A discrete entropy inequality under a condition on the modification is obtained thanks to a reinterpretation of the modified scheme in the Harten Lax and van Leer formalism. A natural extension to multi-dimensional problems discretized over unstructured mesh is proposed. Then a simple and efficient semi implicit scheme is also proposed. The resulting scheme is stable under a CFL condition driven by the (slow) material waves and not by the (fast) acoustic waves and so verifies the all regime property. Numerical evidences are proposed and show the ability of the scheme to deal with tests where the flow regime may vary from low to high Mach values.

Low cost 3D global instability analysis and flow sensitivity based on dynamic mode decomposition and high-order numerical tools

We explore the recently developed snapshot-based dynamic mode decomposition (DMD) technique, a matrix-free Arnoldi type method, to predict 3D linear global flow instabilities. We apply the DMD technique to flows confined in an L-shaped cavity and compare the resulting modes to their counterparts issued from classic, matrix forming, linear instability analysis (i.e. BiGlobal approach) and direct numerical simulations. Results show that the DMD technique, which uses snapshots generated by a 3D non-linear incompressible discontinuous Galerkin Navier?Stokes solver, provides very similar results to classical linear instability analysis techniques. In addition, we compare DMD results issued from non-linear and linearised Navier?Stokes solvers, showing that linearisation is not necessary (i.e. base flow not required) to obtain linear modes, as long as the analysis is restricted to the exponential growth regime, that is, flow regime governed by the linearised Navier?Stokes equations, and showing the potential of this type of analysis based on snapshots to general purpose CFD codes, without need of modifications. Finally, this work shows that the DMD technique can provide three-dimensional direct and adjoint modes through snapshots provided by the linearised and adjoint linearised Navier?Stokes equations advanced in time. Subsequently, these modes are used to provide structural sensitivity maps and sensitivity to base flow modification information for 3D flows and complex geometries, at an affordable computational cost. The information provided by the sensitivity study is used to modify the L-shaped geometry and control the most unstable 3D mode.

The importance of a triple bottom line approach for safeguarding urban water quality

Water environments are greatly valued in urban areas as ecological and aesthetic assets. However, it is the water environment that is most adversely affected by urbanisation. Urban land use coupled with anthropogenic activities alters the stream flow regime and degrade water quality with urban stormwater being a significant source of pollutants. Unfortunately, urban water pollution is difficult to evaluate in terms of conventional monetary measures. True costs extend beyond immediate human or the physical boundaries of the urban area and affect the function of surrounding ecosystems. Current approaches for handling stormwater pollution and water quality issues in urban landscapes are limited as these are primarily focused on ‘end-of-pipe’ solutions. The approaches are commonly based either on, insufficient design knowledge, faulty value judgements or inadequate consideration of full life cycle costs. It is in this context that the adoption of a triple bottom line approach is advocated to safeguard urban water quality. The problem of degradation of urban water environments can only be remedied through innovative planning, water sensitive engineering design and the foresight to implement sustainable practices. Sustainable urban landscapes must be designed to match the triple bottom line needs of the community, starting with ecosystem services first such as the water cycle, then addressing the social and immediate ecosystem health needs, and finally the economic performance of the catchment. This calls for a cultural change towards urban water resources rather than the current piecemeal and single issue focus approach. This paper discusses the challenges in safeguarding urban water environments and the limitations of current approaches. It then explores the opportunities offered by integrating innovative planning practices with water engineering concepts into a single cohesive framework to protect valuable urban ecosystem assets. Finally, a series of recommendations are proposed for protecting urban water resources within the context of a triple bottom line approach.

A finite volume method for solving the two-sided time-space fractional advection-dispersion equation

The field of fractional differential equations provides a means for modelling transport processes within complex media which are governed by anomalous transport. Indeed, the application to anomalous transport has been a significant driving force behind the rapid growth and expansion of the literature in the field of fractional calculus. In this paper, we present a finite volume method to solve the time-space two-sided fractional advection dispersion equation on a one-dimensional domain. Such an equation allows modelling different flow regime impacts from either side. The finite volume formulation provides a natural way to handle fractional advection-dispersion equations written in conservative form. The novel spatial discretisation employs fractionally-shifted Gr¨unwald formulas to discretise the Riemann-Liouville fractional derivatives at control volume faces in terms of function values at the nodes, while the L1-algorithm is used to discretise the Caputo time fractional derivative. Results of numerical experiments are presented to demonstrate the effectiveness of the approach.

Visualisation of complex flows using texture-based techniques

Complex flow datasets are often difficult to represent in detail using traditional vector visualisation techniques such as arrow plots and streamlines. This is particularly true when the flow regime changes in time. Texture-based techniques, which are based on the advection of dense textures, are novel techniques for visualising such flows (i.e., complex dynamics and time-dependent). In this paper, we review two popular texture-based techniques and their application to flow datasets sourced from real research projects. The texture-based techniques investigated were Line Integral Convolution (LIC), and Image-Based Flow Visualisation (IBFV). We evaluated these techniques and in this paper report on their visualisation effectiveness (when compared with traditional techniques), their ease of implementation, and their computational overhead.

Visualisation of complex flows using texture-based techniques

Detailed representations of complex flow datasets are often difficult to generate using traditional vector visualisation techniques such as arrow plots and streamlines. This is particularly true when the flow regime changes in time. Texture-based techniques, which are based on the advection of dense textures, are novel techniques for visualising such flows. We review two popular texture based techniques and their application to flow datasets sourced from active research projects. The techniques investigated were Line integral convolution (LIC) [1], and Image based flow visualisation (IBFV) [18]. We evaluated these and report on their effectiveness from a visualisation perspective. We also report on their ease of implementation and computational overheads.

Partitioning the variation in stream fish assemblages within a spatio-temporal hierarchy

This paper describes the relative influence of: (i) landscape scale environmental and hydrological factors; (ii) local scale environmental conditions including recent flow history, and; (iii) spatial effects (proximity of sites to one another) on the spatial and temporal variation in local freshwater fish assemblages in the Mary River, south-eastern Queensland, Australia. Using canonical correspondence analysis, each of the three sets of variables explained similar amounts of variation in fish assemblages (ranging from 44 to 52%). Variation in fish assemblages was partitioned into eight unique components: pure environmental, pure spatial, pure temporal, spatially structured environmental variation, temporally structured environmental variation, spatially structured temporal variation, the combined spatial/temporal component of environmental variation and unexplained variation. The total variation explained by these components was 65%. The combined spatial/temporal/environmental component explained the largest component (30%) of the total variation in fish assemblages, whereas pure environmental (6%), temporal (9%) and spatial (2%) effects were relatively unimportant. The high degree of intercorrelation between the three different groups of explanatory variables indicates that our understanding of the importance to fish assemblages of hydrological variation (often highlighted as the major structuring force in river systems) is dependent on the environmental context in which this role is examined.

Terrestrial invertebrates of dry river beds are not simply subsets of riparian assemblages

Dry river beds are common worldwide and are rapidly increasing in extent due to the effects of water management and prolonged drought periods due to climate change. While attention has been given to the responses of aquatic invertebrates to drying rivers, few studies exist on the terrestrial invertebrates colonizing dry river beds. Dry river beds are physically harsh and they often differ substantially in substrate, topography, microclimate and inundation frequency from adjacent riparian zones. Given these differences, we predicted that dry river beds provide a unique habitat for terrestrial invertebrates, and that their assemblage composition differs from that in adjacent riparian zones. Dry river beds and riparian zones in Australia and Italy were sampled for terrestrial invertebrates with pitfall traps. Sites differed in substrate type, climate and flow regime. Dry river beds contained diverse invertebrate assemblages and their composition was consistently different from adjacent riparian zones, irrespective of substrate, climate or hydrology. Although some taxa were shared between dry river beds and riparian zones, 66 of 320 taxa occurred only in dry river beds. Differences were due to species turnover, rather than shifts in abundance, indicating that dry river bed assemblages are not simply subsets of riparian assemblages. Some spatial patterns in invertebrate assemblages were associated with environmental variables (irrespective of habitat type), but these associations were statistically weak. We suggest that dry river beds are unique habitats in their own right. We discuss potential human stressors and management issues regarding dry river beds and provide recommendations for future research.

Performance assessment in a heat exchanger tube fitted with double counter twisted tape inserts

The present study explored the effects of the double counter twisted tapes on heat transfer and fluid friction characteristics in a heat exchanger tube. The double counter twisted tapes were used as counter-swirl flow generators in the test section. The experiments were performed with double counter twisted tapes of four different twist ratios (y = 1.95, 3.85, 5.92 and 7.75) using air as the testing fluid in a circular tube turbulent flow regime where the Reynolds number was varied from 6950 to 50,050. The experimental results demonstrated that the Nusselt number, friction factor and thermal enhancement efficiency were increased with decreasing twist ratio. The results also revealed that the heat transfer rate in the tube fitted with double counter twisted tape was significantly increased with corresponding increase in pressure drop. In the range of the present work, heat transfer rate and friction factor were obtained to be around 60 to 240% and 91 to 286% higher than those of the plain tube values, respectively. The maximum thermal enhancement efficiency of 1.34 was achieved by the use of double counter twisted tapes at constant blower power. In addition, the empirical correlations for the Nusselt number, friction factor and thermal enhancement efficiency were also developed, based on the experimental data.

Drop impact on a solid surface comprising micro groove structure

Spreading and receding processes of water drops impacting on a stainless steel surface comprising rectangular shaped parallel grooves are studied experimentally. The study was confined to the impact of drops in inertia dominated flow regime with Weber number in the range 15 - 257. Measurements of spreading drop diameter and drop height were obtained during the impact process as function of time. Experimental measurements of spreading drop diameter and drop height obtained for the grooved surface were compared with those obtained for a smooth surface to elucidate the influence of surface grooves on the impact process. The grooves definitely influence both spreading and receding processes of impacting liquid drops. A more striking observation from this study is that the receding process of impacting liquid drops is dramatically changed by the groove structure for all droplet Weber number.

Ion irradiation enhances the mechanical performance of metallic glasses

We demonstrate that irradiation may enhance the plasticity in metallic glasses by increasing the free-volume content via micropillar compression experiments on an ion-irradiated bulk metallic glass (BMG). Results show that irradiation decreases the flow stress and enhances the shear band formation by lowering the magnitude of stress serrations in plastic flow regime. These results highlight that amorphous alloys can mitigate the deleterious affects of severe ion irradiation as compared to their crystalline counterparts. (C) 2009 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

The collective dynamics of self-propelled particles

We propose a method for the dynamic simulation of a collection of self-propelled particles in a viscous Newtonian fluid. We restrict attention to particles whose size and velocity are small enough that the fluid motion is in the creeping flow regime. We propose a simple model for a self-propelled particle, and extended the Stokesian Dynamics method to conduct dynamic simulations of a collection of such particles. In our description, each particle is treated as a sphere with an orientation vector p, whose locomotion is driven by the action of a force dipole Sp of constant magnitude S0 at a point slightly displaced from its centre. To simplify the calculation, we place the dipole at the centre of the particle, and introduce a virtual propulsion force Fp to effect propulsion. The magnitude F0 of this force is proportional to S0. The directions of Sp and Fp are determined by p. In isolation, a self-propelled particle moves at a constant velocity u0 p, with the speed u0 determined by S0. When it coexists with many such particles, its hydrodynamic interaction with the other particles alters its velocity and, more importantly, its orientation. As a result, the motion of the particle is chaotic. Our simulations are not restricted to low particle concentration, as we implement the full hydrodynamic interactions between the particles, but we restrict the motion of particles to two dimensions to reduce computation. We have studied the statistical properties of a suspension of self-propelled particles for a range of the particle concentration, quantified by the area fraction φa. We find several interesting features in the microstructure and statistics. We find that particles tend to swim in clusters wherein they are in close proximity. Consequently, incorporating the finite size of the particles and the near-field hydrodynamic interactions is of the essence. There is a continuous process of breakage and formation of the clusters. We find that the distributions of particle velocity at low and high φa are qualitatively different; it is close to the normal distribution at high φa, in agreement with experimental measurements. The motion of the particles is diffusive at long time, and the self-diffusivity decreases with increasing φa. The pair correlation function shows a large anisotropic build-up near contact, which decays rapidly with separation. There is also an anisotropic orientation correlation near contact, which decays more slowly with separation. Movies are available with the online version of the paper.

Model for Predicting the Mean Drop Size in Effervescent Sprays

This paper is focused on the development of a model for predicting the mean drop size in effervescent sprays. A combinatorial approach is followed in this modeling scheme, which is based on energy and entropy principles. The model is implemented in cascade in order to take primary breakup (due to exploding gas bubbles) and secondary breakup (due to shearing action of surrounding medium) into account. The approach in this methodology is to obtain the most probable drop size distribution by maximizing the entropy while satisfying the constraints of mass and energy balance. The comparison of the model predictions with the past experimental data is presented for validation. A careful experimental study is conducted over a wide range of gas-to-liquid ratios, which shows a good agreement with the model predictions: It is observed that the model gives accurate results in bubbly and annular flow regimes. However, discrepancies are observed in the transitional slug flow regime of the atomizer.

Mass Distribution Studies in Effervescent Sprays

An experimental study on liquid mass distribution in effervescent sprays using water and air as working fluids is presented in this paper. Optical patternation and techniques of image processing are employed for analyzing the spray. The flow regime inside the effervescent atomizer largely dictates the mass distribution patterns. The patterns are seen to vary from concentrated, poorly atomized liquid lumps to uniformly distributed, fine droplets as the flow regime changes from bubbly flow to annular flow. A large variety of instantaneous spray patterns are observed in bubbly flow regime indicating a highly unsteady atomization process. However, relatively better consistency in spray patterns is observed at higher gas flow rates. Thus, the degree of unsteadiness gradually diminishes as gas flow rate is increased. The axial evolution of the spray in annular mode shows good mixing of liquid and gas across the interface.

Anomalous Stress Profile in a Sheared Granular Column

We present measurements of the stress as a function of vertical position in a column of granular material sheared in a cylindrical Couette device. All three components of the stress tensor on the outer cylinder were measured as a function of distance from the free surface at shear rates low enough that the material was in the dense, slow flow regime. We find that the stress profile differs fundamentally from that of fluids, from the predictions of plasticity theories, and from intuitive expectation. We argue that the anomalous stress profile is due to an anisotropic fabric caused by the combined action of gravity and shear.