A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil. (C) 2011 Elsevier Ltd. All rights reserved.

Effect of side confining walls on the growth rate of compressible mixing layers

Model free simulations are performed to study the effect of the presence of side wall in compressible mixing of two parallel dissimilar gaseous streams with significant temperature difference. The turbulence statistics shows the three dimensional nature of the flow with and without the presence of side walls. The presence of side wall neither makes the flow field two dimensional, nor suppresses three dimensional disturbances. However, the comparison of shear layer growth rate and wall pressures reveal a better match with the two dimensional simulation results. This better match is explained on the basis of formation of oblique structures due to the presence of side walls which also suppress the distribution of momentum in third direction making the pressures to be higher as compared with the case without side walls. (C) 2013 Elsevier Ltd. All rights reserved.

Computational fluid-structure interaction of DGB parachutes in compressible fluid flow

This paper describes large-scale simulations of compressible flows over a supersonic disk-gap-band parachute system. An adaptive mesh refinement method is used to resolve the coupled fluid-structure model. The fluid model employs large-eddy simulation to describe the turbulent wakes appearing upstream and downstream of the parachute canopy and the structural model employed a thin-shell finite element solver that allows large canopy deformations by using subdivision finite elements. The fluid-structure interaction is described by a variant of the Ghost-Fluid method. The simulation was carried out at Mach number 1.96 where strong nonlinear coupling between the system of bow shocks, turbulent wake and canopy is observed. It was found that the canopy oscillations were characterized by a breathing type motion due to the strong interaction of the turbulent wake and bow shock upstream of the flexible canopy. Copyright © 2010 by ASME.

An efficient shock capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.

A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids

In addition to the Hamiltonian functional itself, non-canonical Hamiltonian dynamical systems generally possess integral invariants known as ‘Casimir functionals’. In the case of the Euler equations for a perfect fluid, the Casimir functionals correspond to the vortex topology, whose invariance derives from the particle-relabelling symmetry of the underlying Lagrangian equations of motion. In a recent paper, Vallis, Carnevale & Young (1989) have presented algorithms for finding steady states of the Euler equations that represent extrema of energy subject to given vortex topology, and are therefore stable. The purpose of this note is to point out a very general method for modifying any Hamiltonian dynamical system into an algorithm that is analogous to those of Vallis etal. in that it will systematically increase or decrease the energy of the system while preserving all of the Casimir invariants. By incorporating momentum into the extremization procedure, the algorithm is able to find steadily-translating as well as steady stable states. The method is applied to a variety of perfect-fluid systems, including Euler flow as well as compressible and incompressible stratified flow.

On the validity of single-parcel energetics to assess the importance of internal energy and compressibility effects in stratified fluids

It is often assumed on the basis of single-parcel energetics that compressible effects and conversions with internal energy are negligible whenever typical displacements of fluid parcels are small relative to the scale height of the fluid (defined as the ratio of the squared speed of sound over gravitational acceleration). This paper shows that the above approach is flawed, however, and that a correct assessment of compressible effects and internal energy conversions requires considering the energetics of at least two parcels, or more generally, of mass conserving parcel re-arrangements. As a consequence, it is shown that it is the adiabatic lapse rate and its derivative with respect to pressure, rather than the scale height, which controls the relative importance of compressible effects and internal energy conversions when considering the global energy budget of a stratied fluid. Only when mass conservation is properly accounted for is it possible to explain why available internal energy can account for up to 40 percent of the total available potential energy in the oceans. This is considerably larger than the prediction of single-parcel energetics, according to which this number should be no more than about 2 percent.

Contribution of thickened drinks, food and enteral and parenteral fluids to fluid intake in hospitalised patients with dysphagia

Background Studies amongst older people with acute dysphagic stroke requiring thickened fluids have assessed fluid intakes from combinations of beverage, food, enteral and parenteral sources, but not all sources simultaneously. The study aimed to comprehensively assess total water intake from food, beverages, enteral and parenteral sources amongst dysphagic adult in-patients receiving thickened fluids. Methods Patients requiring thickened fluid following dysphagia diagnosis were recruited consecutively from a tertiary teaching hospital’s medical and neurosurgical wards. Fluid intake from food and beverages was assessed by wastage, direct observation and quantified from enteral and parenteral sources through clinical medical records. Results No patients achieved their calculated fluid requirements unless enteral or parenteral fluids were received. The mean daily fluid intake from food was greater than from beverages whether receiving diet alone (food 807±363mL, food and beverages 370±179mL, p<0.001) or diet with enteral or parenteral fluid support (food 455±408mL, food and beverages 263±232mL, p<0.001). Greater daily fluid intakes occurred when receiving enteral and parenteral fluid in addition to oral dietary intake, irrespective of age group, whether assistance was required, diagnosis and whether stage 3 or stage 2 thickened fluids were required (p<0.05). After enteral and parenteral sources, food provided the most important contribution to daily fluid intakes. Conclusions The greatest contribution to oral fluid intake was from food, not beverages. Designing menus and food services which promote and encourage the enjoyment of fluid dense foods, in contrast to thickened beverages, may present an important way to improve fluid intakes of those with dysphagia. Supplemental enteral or parenteral fluid may be necessary to achieve minimum calculated fluid requirements.

Candidate radial-inflow turbines and high-density working fluids for geothermal power systems

Optimisation of Organic Rankine Cycle (ORCs) for binary-cycle geothermal applications could play a major role in determining the competitiveness of low to moderate temperature geothermal resources. Part of this optimisation process is matching cycles to a given resource such that power output can be maximised. Two major and largely interrelated components of the cycle are the working fluid and the turbine. Both components need careful consideration: the selection of working fluid and appropriate operating conditions as well as optimisation of the turbine design for those conditions will determine the amount of power that can be extracted from a resource. In this paper, we present the rationale for the use of radial-inflow turbines for ORC applications and the preliminary design of several radial-inflow machines based on a number of promising ORC systems that use five different working fluids: R134a, R143a, R236fa, R245fa and n-Pentane. Preliminary meanline analysis lead to the generation of turbine designs for the various cycles with similar efficiencies (77%) but large differences in dimensions (139–289 mm rotor diameter). The highest performing cycle, based on R134a, was found to produce 33% more net power from a 150 °C resource flowing at 10 kg/s than the lowest performing cycle, based on n-Pentane.

Candidate radial-inflow turbines and high-density working fluids for geothermal power systems

Optimisation of Organic Rankine Cycles (ORCs) for binary-cycle geothermal applications could play a major role in the competitiveness of low to moderate temperature geothermal resources. Part of this optimisation process is matching cycles to a given resource such that power output can be maximised. Two major and largely interrelated components of the cycle are the working fluid and the turbine. Both components need careful consideration. Due to the temperature differences in geothermal resources a one-size-fits-all approach to surface power infrastructure is not appropriate. Furthermore, the traditional use of steam as a working fluid does not seem practical due to the low temperatures of many resources. A variety of organic fluids with low boiling points may be utilised as ORC working fluids in binary power cycle loops. Due to differences in thermodynamic properties, certain fluids are able to extract more heat from a given resource than others over certain temperature and pressure ranges. This enables the tailoring of power cycle infrastructure to best match the geothermal resource through careful selection of the working fluid and turbine design optimisation to yield the optimum overall cycle performance. This paper presents the rationale for the use of radial-inflow turbines for ORC applications and the preliminary design of several radial-inflow turbines based on a selection of promising ORC cycles using five different high-density working fluids: R134a, R143a, R236fa, R245fa and n-Pentane at sub- or trans-critical conditions. Numerous studies published compare a variety of working fluids for various ORC configurations. However, there is little information specifically pertaining to the design and implementation of ORCs using realistic radial turbine designs in terms of pressure ratios, inlet pressure, rotor size and rotational speed. Preliminary 1D analysis leads to the generation of turbine designs for the various cycles with similar efficiencies (77%) but large differences in dimensions (139289 mm rotor diameter). The highest performing cycle (R134a) was found to produce 33% more net power from a 150°C resource flowing at 10 kg/s than the lowest performing cycle (n-Pentane).