Non-equilibirum condensation effects on the flow field and the performance of a low pressure steam turbine

The influence of non-equilibrium condensation on the flow field and performance of a three stage low pressure model steam turbine is examined using modern three dimensional CFD techniques. An equilibrium steam model and a non-equilibrium steam model, which accounts for both subcooling and condensation effects, are used, and have been verified by comparison with test data in an earlier publication [1]. The differences in the calculated flow field and turbine performance with these models show that the latent heat released during condensation influences both the thermodynamic and the aerodynamic performance of the turbine, leading to a change in inlet flow angles of about 5°. The calculated three dimensional flowfield is used to investigate the magnitude and distribution of the additional thermo-dynamic wetness loss arising from steam condensation under non-equilibrium flow conditions. Three simple methods are described to calculate this, and all show that this amounts to around 6.5% of the total losses at the design condition. At other load conditions the wetness losses change in magnitude and axial distribution in the turbine. © 2010 by ASME.

Analysis and simulation of combustion oscillations in a laboratory spray combustor

The unstable combustion that can occur in combustion chambers is a major problem for aeroengines and ground-based industrial gas turbines. Nowadays, CFD provides a flexible, low cost tool to supplement direct measurement. This paper presents simulations of combustion oscillations in a liquid-fuelled experimental rig at the University of Cambridge. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone and to develop inlet and outlet boundary conditions just upstream and downstream of the combustion region enabling the CFD calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in good agreement with experimental measurements. More details about the unstable combustion can be obtained from the simulation results. The approach developed here is expected to provide a powerful tool for the design and operation of stable combustion systems. Copyright © 2009 by ASME.

Microstructural optimization of cellular acoustic materials

A microstructure based acoustic model is introduced, which can be used to optimize the microstructure of cellular materials and thus to obtain their optimal acoustic property. This acoustic model is an unsteady one which is appropriate in the limit of low Reynolds numbers. The model involves three elements. This first involves the propagation of acoustic waves passing the cylinders whose axes are aligned parallel to the direction of propagation. The second model relates to the propagation of acoustic waves passing the cylinders whose axes are aligned perpendicular to the direction of propagation. In both cases the interaction between adjacent cylinders is taken into account by considering the effect of polygonal periodic boundary conditions. As these two models are linear they are combined to give the characteristics of propagation at arbitrary incidence. The third model involves propagation passing spheres in order to represent the joints. Heat transfer is also included. These three models are then used to expand the design space and calculate the optimum cell structure for desired acoustic performance in a number of different applications. Moreover, the application fields are also analyzed.

Investigation of the losses in the stator duct spacers of large A.C. motors using the A -ip 3-D eddy current formulation

This paper describes how the A -if) formulation may be applied to determine the losses in the stator duct spacers of large a.c. motors. The model is described in terms of its geometry and boundary conditions. The novel aspects of the application of the formulation to this problem are explained. These include the modelling of fixed currents sources (the stator windings), the location of the necessary cut surfaces and the determination of their magnetic scalar potential differences, and the implementation of periodic boundary conditions for vector variables. Results are presented showing how the duct spacer losses vary with load, and with the relative permeability of the spacer material. The effects of modelling iron nonlinearity, of both the spacer and the steel laminations, are also illustrated. © 1996 IEEE.

On the domain decomposition and transmission line modelling finite element method for time-domain induction motor analysis

This paper demonstrates how a finite element model which exploits domain decomposition is applied to the analysis of three-phase induction motors. It is shown that a significant gain in cpu time results when compared with standard finite element analysis. Aspects of the application of the method which are particular to induction motors are considered: the means of improving the convergence of the nonlinear finite element equations; the choice of symmetrical sub-domains; the modelling of relative movement; and the inclusion of periodic boundary conditions. © 1999 IEEE.

Adjoint linearised Euler solver in the frequency domain for jet noise modelling

The paper is devoted to extending the new efficient frequency-domain method of adjoint Green's function calculation to curvilinear multi-block RANS domains for middle and farfield sound computations. Numerical details of the method such as grids, boundary conditions and convergence acceleration are discussed. Two acoustic source models are considered in conjunction with the method and acoustic modelling results are presented for a benchmark low-Reynolds-number jet case.

TOPS (Task optimization in the presence of signal-dependent noise) model

This paper proposes a movement trajectory planning model, which is a maximum task achievement model in which signal-dependent noise is added to the movement command. In the proposed model, two optimization criteria are combined, maximum task achievement and minimum energy consumption. The proposed model has the feature that the end-point boundary conditions for position, velocity, and acceleration need not be prespecified. Consequently, the method can be applied not only to the simple point-to-point movement, but to any task. In the method in this paper, the hand trajectory is derived by a psychophysical experiment and a numerical experiment for the case in which the target is not stationary, but is a moving region. It is shown that the trajectory predicted from the minimum jerk model or the minimum torque change model differs considerably from the results of the psychophysical experiment. But the trajectory predicted from the maximum task achievement model shows good qualitative agreement with the hand trajectory obtained from the psychophysical experiment. © 2004 Wiley Periodicals, Inc.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Modelling of long-term ground response to tunnelling under St James's Park, London

Following a tunnel excavation in low-permeability soil, it is commonly observed that the ground surface continues to settle and ground loading on the tunnel lining changes, as the pore pressures in the ground approach a new equilibrium condition. The monitored ground response following the tunnelling under St James's Park, London, shows that the mechanism of subsurface deformation is composed of three different zones: swelling, consolidation and rigid body movement. The swelling took place in a confined zone above the tunnel crown, extending vertically to approximately 5 m above it. On the sides of the tunnel, the consolidation of the soil occurred in the zone primarily within the tunnel horizon, from the shoulder to just beneath the invert, and extending laterally to a large offset from the tunnel centreline. Above these swelling and consolidation zones the soil moved downward as a rigid body. In this study, soil-fluid coupled three-dimensional finite element analyses were performed to simulate the mechanism of long-term ground response monitored at St James's Park. An advanced critical state soil model, which can simulate the behaviour of London Clay in both drained and undrained conditions, was adopted for the analyses. The analysis results are discussed and compared with the field monitoring data. It is found that the observed mechanism of long-term subsurface ground and tunnel lining response at St James's Park can be simulated accurately only when stiffness anisotropy, the variation of permeability between different units within the London Clay and non-uniform drainage conditions for the tunnel lining are considered. This has important implications for future prediction of the long-term behaviour of tunnels in clays.

Vibration of fluid loaded conical shells.

An analytical model is presented to describe the vibration of a truncated conical shell with fluid loading in the low frequency range. The solution for the dynamic response of the shell is presented in the form of a power series. Fluid loading is taken into account by dividing the shell into narrow strips which are considered to be locally cylindrical. Analytical results are presented for different boundary conditions and have been compared with the computational results from a boundary element model. Limitations of the model to the low frequency range are discussed.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Synchronization of modulated traveling baroclinic waves in a periodically forced, rotating fluid annulus.

Frequency entrainment and nonlinear synchronization are commonly observed between simple oscillatory systems, but their occurrence and behavior in continuum fluid systems are much less well understood. Motivated by possible applications to geophysical fluid systems, such as in atmospheric circulation and climate dynamics, we have carried out an experimental study of the interaction of fully developed baroclinic instability in a differentially heated, rotating fluid annulus with an externally imposed periodic modulation of the thermal boundary conditions. In quasiperiodic and chaotic amplitude-modulated traveling wave regimes, the results demonstrate a strong interaction between the natural periodic modulation of the wave amplitude and the externally imposed forcing. This leads to partial or complete phase synchronization. Synchronization effects were observed even with very weak amplitudes of forcing, and were found with both 1:1 and 1:2 frequency ratios between forcing and natural oscillations.

QM/MM simulation of liquid water with an adaptive quantum region.

The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 Å) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl(-) ion.

Corner effect and asymmetry in transonic channel flows

An experimental and numerical investigation into transonic shock/boundary-layer interactions in rectangular ducts has been performed. Experiments have shown that flow development in the corners of transonic shock/boundary-layer interactions in confined channels can have a significant impact on the entire flowfield. As shock strength is increased from M∞ = 1:3 to 1.5, the flowfield becomes very slightly asymmetrical. The interaction of corner flows with one another is thought to be a potential cause of this asymmetry. Thus, factors that govern the size of corner interactions (such as interaction strength) and their proximity to one another (such as tunnel aspect ratio) can affect flow symmetry. The results of the computational study show reasonable agreement with experiments, although simulations with particular turbulence models predict highly asymmetrical solutions for flows that were predominantly symmetrical in experiments. These discrepancies are attributed to the tendency of numerical schemes to overprediction corner-interaction size, and this also accounts for why computational fluid dynamics predicts the onset of asymmetry at lower shock strengths than in experiments. The findings of this study highlight the importance of making informed decisions about imposing artificial constraints on symmetry and boundary conditions for internal transonic flows. Future effort into modeling corner flows accurately is required. Copyright © 2011 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.