Numerical simulations of unsteady wet steam flows with non-equilibrium condensation in the nozzle and the steam turbine

Numerical techniques for non-equilibrium condensing flows are presented. Conservation equations for homogeneous gas-liquid two-phase compressible flows are solved by using a finite volume method based on an approximate Riemann solver. The phase change consists of the homogeneous nucleation and growth of existing droplets. Nucleation is computed with the classical Volmer-Frenkel model, corrected for the influence of the droplet temperature being higher than the steam temperature due to latent heat release. For droplet growth, two types of heat transfer model between droplets and the surrounding steam are used: a free molecular flow model and a semi-empirical two-layer model which is deemed to be valid over a wide range of Knudsen number. The computed pressure distribution and Sauter mean droplet diameters in a convergent-divergent (Laval) nozzle are compared with experimental data. Both droplet growth models capture qualitatively the pressure increases due to sudden heat release by the non-equilibrium condensation. However the agreement between computed and experimental pressure distributions is better for the two-layer model. The droplet diameter calculated by this model also agrees well with the experimental value, whereas that predicted by the free molecular model is too small. Condensing flows in a steam turbine cascade are calculated at different Mach numbers and inlet superheat conditions and are compared with experiments. Static pressure traverses downstream from the blade and pressure distributions on the blade surface agree well with experimental results in all cases. Once again, droplet diameters computed with the two-layer model give best agreement with the experiments. Droplet sizes are found to vary across the blade pitch due to the significant variation in expansion rate. Flow patterns including oblique shock waves and condensation-induced pressure increases are also presented and are similar to those shown in the experimental Schlieren photographs. Finally, calculations are presented for periodically unsteady condensing flows in a low expansion rate, convergent-divergent (Laval) nozzle. Depending on the inlet stagnation subcooling, two types of self-excited oscillations appear: a symmetric mode at lower inlet subcooling and an asymmetric mode at higher subcooling. Plots of oscillation frequency versus inlet sub-cooling exhibit a hysteresis loop, in accord with observations made by other researchers for moist air flow. Copyright © 2006 by ASME.

Numerical simulation of the growth of ZnO nanostructures in a tube furnace by physical vapour deposition

Zinc oxide is a versatile II-VI naturally n-type semiconductor that exhibits piezoelectric properties. By controlling the growth kinetics during a simple carbothermal reduction process a wide range of 1D nanostructures such as nanowires, nanobelts, and nanotetrapods have been synthesized. The driving force: for the nanostructure growth is the Zn vapour supersaturation and supply rate which, if known, can be used to predict and explain the type of crystal structure that results. A model which attempts to determine the Zn vapour concentration as a function of position in the growth furnace is described. A numerical simulation package, COMSOL, was used to simultaneously model the effects of fluid flow, diffusion and heat transfer in a tube furnace made specifically for ZnO nanostructure growth. Parameters such as the temperature, pressure, and flow rate are used as inputs to the model to show the effect that each one has on the Zn concentration profile. An experimental parametric study of ZnO nanostructure growth was also conducted and compared to the model predictions for the Zn concentration in the tube. © 2008 Materials Research Society.

CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

Physics in fish processing technology

To design, develop and put into operation an equipment either to increase the productivity or to improve the existing technique to obtain a better quality of the product, the fishery engineer/scientist should have a comprehensive knowledge of fundamental principles involved in the process. Many a technique in fish processing technology, whether it applies to freezing, dehydration or canning, involves always a type of heat transfer, which is dependent to a certain extent on the external physical parameters like temperature. humidity, pressure, air flow etc. and also on the thermodynamic properties of fish muscle in the temperature ranges encountered. Similarly informations on other physical values like dielectric constant and dielectric loss in the design of quick trawlers and in quality assessment of frozen/iced fish, refractive index and viscosity in the measurement of the saturation and polymerisation of fish oils and shear strength in the judgement of textural qualities of cooked fish are also equally important.

Electricity Storage Using a Thermal Storage Scheme

The increasing use of renewable energy technologies for electricity generation, many of which have an unpredictably intermittent nature, will inevitably lead to a greater demand for large-scale electricity storage schemes. For example, the expanding fraction of electricity produced by wind turbines will require either backup or storage capacity to cover extended periods of wind lull. This paper describes a recently proposed storage scheme, referred to here as Pumped Thermal Storage (PTS), and which is based on "sensible heat" storage in large thermal reservoirs. During the charging phase, the system effectively operates as a high temperature-ratio heat pump, extracting heat from a cold reservoir and delivering heat to a hot one. In the discharge phase the processes are reversed and it operates as a heat engine. The round- trip efficiency is limited only by process irreversibilities (as opposed to Second Law limitations on the coefficient of performance and the thermal efficiency of the heat pump and heat engine respectively). PTS is currently being developed in both France and England. In both cases, the schemes operate on the Joule-Brayton (gas turbine) cycle, using argon as the working fluid. However, the French scheme proposes the use of turbomachinery for compression and expansion, whereas for that being developed in England reciprocating devices are proposed. The current paper focuses on the impact of the various process irreversibilities on the thermodynamic round-trip efficiency of the scheme. Consideration is given to compression and expansion losses and pressure losses (in pipe-work, valves and thermal reservoirs); heat transfer related irreversibility in the thermal reservoirs is discussed but not included in the analysis. Results are presented demonstrating how the various loss parameters and operating conditions influence the overall performance.

The mechanical integrity of fuel pin cladding in a pulsed-beam accelerator driven subcritical reactor

The Accelerator Driven Subcritical Reactor (ADSR) is one of the reactor designs proposed for future nuclear energy production. Interest in the ADSR arises from its enhanced and intrinsic safety characteristics, as well as its potential ability to utilize the large global reserves of thorium and to burn legacy actinide waste from other reactors and decommissioned nuclear weapons. The ADSR concept is based on the coupling of a particle accelerator and a subcritical core by means of a neutron spallation target interface. One of the candidate accelerator technologies receiving increasing attention, the Fixed Field Alternating Gradient (FFAG) accelerator, generates a pulsed proton beam. This paper investigates the impact of pulsed proton beam operation on the mechanical integrity of the fuel pin cladding. A pulsed beam induces repetitive temperature changes in the reactor core which lead to cyclic thermal stresses in the cladding. To perform the thermal analysis aspects of this study a code that couples the neutron kinetics of a subcritical core to a cylindrical geometry heat transfer model was developed. This code, named PTS-ADS, enables temperature variations in the cladding to be calculated. These results are then used to perform thermal fatigue analysis and to predict the stress-life behaviour of the cladding. © 2011 Elsevier Ltd. All rights reserved.

Effects of flow geometry and Damkohler number on turbulent premixed flame

DNS of planar turbulent flame and turbulent V-flame has been conducted to investigate turbulence-scalar interaction in relatively practical turbulent combustion. Several turbulence quantities are examined for the understandings of fundamental characteristics of flow field in V-flame. Due to the additional turbulence production by the hot-rod, turbulence does not simply decay in V-flame. Turbulence-scalar interaction, scalar alignments with the principal strain rate in other words, is then clarified. The competition of turbulence and dilatation can be found in the conditional PDF of flame normal alignment. The results suggests that the alignment characteristics in high Da flames are applicable to low Da flames in the region of intense heat release.

Increased gas production from hydrates by combining depressurization with heating of the wellbore

To extract gas from hydrate reservoirs, it has to be dissociated in situ. This endothermic dissociation process absorbs heat energy from the formation and pore fluid. The heat transfer governs the dissociation rate, which is proportional to the difference between the actual temperature and the equilibrium temperature. This study compares three potential gas production schemes from hydrate-bearing soil, where the radial heat transfer is governing. Cylindrical samples with 40% pore-filling hydrate saturation were tested. The production tests were carried out over 90 min by dissociating the hydrate from a centered miniature wellbore, by either lowering the pressure to 6, 4, or 6 MPa with simultaneous heating of the wellbore to 288 K. All tests were replicated by a numerical simulation. With additional heating at the same wellbore pressure, the gas production from hydrates could, on average, be increased by 1.8 and 3.6 times in the simulation and experiments, respectively. If the heat influx from the outer boundary is limited, a simulation showed that the specific heat of the formation is rapidly used up when the wellbore is only depressurized and not heated. © 2012 American Chemical Society.

Strategies for modeling turbulent flows in electronics

Hybrid methods based on the Reynolds Averaged Navier Stokes (RANS) equations and the Large Eddy Simulation (LES) formulation are investigated to try and improve the accuracy of heat transfer and surface temperature predictions for electronics systems and components. Two relatively low Reynolds number flows are studied using hybrid RANS-LES, RANS-Implicit-LES (RANS-ILES) and non-linear LES models. Predictions using these methods are in good agreement with each other, even using different grid resolutions. © 2008 IEEE.

Predictive techniques for turbulent oscillatory flows

The prediction of turbulent oscillatory flow at around transitional Reynolds numbers is considered for an idealized electronics system. To assess the accuracy of turbulence models, comparison is made with measurements. A stochastic procedure is used to recover instantaneous velocity time traces from predictions. This procedure enables more direct comparison with turbulence intensity measurements which have not been filtered to remove the oscillatory flow component. Normal wall distances, required in some turbulence models, are evaluated using a modified Poisson equation based technique. A range of zero, one and two equation turbulence models are tested, including zonal and a non-linear eddy viscosity models. The non-linear and zonal models showed potential for accuracy improvements.

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.

Review of computational aeroacoustics for application in electronics cooler noise

A review of computational aeroacoustics (CCA) was made for application in electronics cooler noise. Computational aeroacoustics encompasses all numerical methods where the purposes is to predict the noise emissions from a simulated flow. Numerical simulation of the flow inside and around heat sinks and fans can lead to a prediction of the emitted noise while they are still in the design phase. Direct CCA is theoretically the best way to predict flow-based acoustic phenomena numerically. It is typically used only for low-frequency sound prediction. The boundary element method offers low computational cost and does not use a computational grid, but instead use vortex-surface calculations to determine tonal noise. Axial fans are commonly used to increase the airflow and thus the heat transfer over the heat sinks within the computer cases. Very detailed source simulations in the fan and heat sink region coupled with the use of analogy methods could result in excellent simulation results with a reasonable computational effort.

Winglets for Improved Aerothermal Performance of High Pressure Turbines

This paper investigates the design of winglet tips for unshrouded high pressure turbine rotors, considering aerodynamic and thermal performance simultaneously. A novel parameterization method has been developed to alter the tip geometry of a rotor blade. A design survey of un-cooled, flat-tipped winglets is performed using RANS calculations for a single rotor at engine representative operating conditions. Compared to a plain tip, large efficiency gains can be realized by employing an overhang around the full perimeter of the blade, but the overall heat load rises significantly. By employing an overhang on only the early suction surface, significant efficiency improvements can be obtained without increasing the overall heat transfer to the blade. The flow physics are explored in detail to explain the results. For a plain tip, the leakage and passage vortices interact to create a three-dimensional impingement onto the blade suction surface, causing high heat transfer. The addition of an overhang on the early suction surface displaces the tip leakage vortex away from the blade, weakening the impingement effect and reducing the heat transfer on the blade. The winglets reduce the aerodynamic losses by unloading the tip section, reducing the leakage flow rate, turning the leakage flow in a more streamwise direction and reducing the interaction between the leakage fluid and endwall flows. Generally these effects are most effective close to the leading edge of the tip, where the leakage flow is subsonic.

Measurement of the unburnt gas temperature in an IC engine by means of a pressure transducer

A novel method of measuring cylinder gas temperature in an internal combustion engine cylinder is introduced. The physical basis for the technique is that the flow rate through an orifice is a function of the temperature of the gas flowing through the orifice. Using a pressure transducer in the cylinder, and another in a chamber connected to the cylinder via an orifice, it is shown how the cylinder temperature can be determined with useful sensitivity. In this paper the governing equations are derived, which show that the heat transfer characteristics of the chamber are critical to the performance of the system, and that isothermal or adiabatic conditions give the optimum performance. For a typical internal combustion engine, it is found that the pre-compression cylinder temperature is related to the chamber pressure late in the compression process with sensitivity of the order of 0.005 bar/K. Copyright © 2010 SAE International.

Investigation of the dilution process for measurement of particulate matter from spark-ignition engines

Measurements of particulate matter (PM) from spark ignition (SI) engine exhaust using dilution tunnels will become more prevalent as emission standards are tightened. Hence, a study of the dilution process was undertaken in order to understand how various dilution related parameters affect the accuracy with which PM sizes and concentrations can be determined. A SI and a compression ignition (CI) engine were separately used to examine parameters of the dilution process; the present work discusses the results in the context of SI exhaust dilution. A Scanning Mobility Particle Sizer (SMPS) was used to measure the size distribution, number density, and volume fraction of PM. Temperature measurements in the exhaust pipe and dilution tunnel reveal the degree of mixing between exhaust and dilution air, the effect of flowrate on heat transfer from undiluted and diluted exhaust to the environment, and the minimum permissible dilution ratio for a maximum sample temperature of 52°C. Measurements of PM concentrations as a function of dilution ratio show the competing effects of temperature and particle/vapor concentrations on particle growth dynamics, which result in a range of dilution ratios-from 13 to 18-where the effect of dilution ratio, independent of flowrate, is kept to a minimum. This range of dilution ratios is therefore optimal in order to achieve repeatable PM concentration measurements. Particle dynamics during transit through the tunnel operating at the optimal dilution ratio was found statistically insignificant compared to data scatter. Such small differences in number concentration may be qualitatively representative of particle losses for SI exhaust, but small increases in PM volume fraction during transit through the tunnel may significantly underestimate accretion of mass due to unburned hydrocarbons (HCs) emitted by SI engines. The fraction of SI-derived PM mass due to adsorbed/absorbed vapor, estimated from these data, is consistent with previous chemical analyses of PM. © 1998 Society of Automotive Engineers, Inc.