Optimal Operations Management and Network Planning of a District Heating System with a Combined Heat and Power Plant

District heating plants are becoming more common in European cities. These systems make it possible to furnish users with warm water while locating the production plants in the outskirts having the double benefit of lowering the impact of pollution on the center of the city and achieving better conversion performances. In order to amortize the costs throughout the year, the system often includes a combined heat and power (CHP) plant, to exploit the energy during the summer as well, when the demand for warm water decreases. A linear programming model for the optimal resource management of such a plant is presented and some results for a real case are reported. A distribution network design problem is also addressed and solved by means of mixed integer linear programming.     

Real-time coloured visualization of phase flows by the schlieren method

A coloured real-time visualizer of optical inhomogeneities comprising a bichromatic schlieren system, video camera and colour monitor has been developed. The schlieren system represents a function Foucault-Hilbert transformation provided with an amplitude spatial frequency filter, or a quadrant Foucault knife edge. Two colour-coded complementary Toepler-grams are obtained in the exit plane of this schlieren system. Their summed image is then recorded by the video camera and displayed on the screen of the colour monitor. The schlieren photograph of internal gravity waves, generated by the cylindrical body motion in the reservoir filled with the stratified liquid, is presented.     

A proof of Howard’s conjecture in homogeneous parallel shear flows

A rigorous mathematical proof of Howard's conjecture which states that the growth rate of an arbitrary unstable wave must approach zero, as the wave length decreases to zero, in the linear instability of nonviscous homogeneous parallel shear flows, is presented here for the first time under the restriction of the boundedness of the second derivative of the basic velocity field with respect to the vertical coordinate in the concerned flow domain.     

Artificial compressibility, characteristics-based schemes for variable density, incompressible, multi-species flows. Part I. Derivation of different formulations and constant density limit

The paper presents various formulations of characteristics-based schemes in the framework of the artificial-compressibility method for variable-density incompressible flows. In contrast to constant-density incompressible flows, where the characteristics-based variables reconstruction leads to a single formulation, in the case of variable density flows three different schemes can be obtained henceforth labeled as: transport, conservative and hybrid schemes. The conservative scheme results in pseudo-compressibility terms in the (multi-species) density reconstruction. It is shown that in the limit of constant density, the transport scheme becomes the (original) characteristics-based scheme for incompressible flows, but the conservative and hybrid schemes lead to a new characteristics-based variant for constant density flows. The characteristics-based schemes are combined with second and third-order interpolation for increasing the computational accuracy locally at the cell faces of the control volume. Numerical experiments for constant density flows reveal that all the characteristics-based schemes result in the same flow solution, but they exhibit different convergence behavior. The multigrid implementation and numerical studies for variable density flows are presented in Part II of this study.     

Confined ferrofluid droplet in crossed magnetic fields

When a ferrofluid drop is trapped in a horizontal Hele-Shaw cell and subjected to a vertical magnetic field, a fingering instability results in the droplet evolving into a complex branched structure. This fingering instability depends on the magnetic field ramp rate but also depends critically on the initial state of the droplet. Small perturbations in the initial droplet can have a large influence on the resulting final pattern. By simultaneously applying a stabilizing (horizontal) azimuthal magnetic field, we gain more control over the mode selection mechanism. We perform a linear stability analysis that shows that any single mode can be selected by appropriately adjusting the strengths of the applied fields. This offers a unique and accurate mode selection mechanism for this confined magnetic fluid system. We present the results of numerical simulations that demonstrate that this mode selection mechanism is quite robust and “overpowers” any initial perturbations on the droplet. This provides a predictable way to obtain patterns with any desired number of fingers.     

Ultrasound velocimetry of ferrofluid spin-up flow measurements using a spherical coil assembly to impose a uniform rotating magnetic field

Ferrofluid spin-up flow is studied within a sphere subjected to a uniform rotating magnetic field from two surrounding spherical coils carrying sinusoidally varying currents at right angles and 90° phase difference. Ultrasound velocimetry measurements in a full sphere of ferrofluid shows no measureable flow. There is significant bulk flow in a partially filled sphere (1-14 mm/s) of ferrofluid or a finite height cylinder of ferrofluid with no cover (1-4 mm/s) placed in the spherical coil apparatus. The flow is due to free surface effects and the non-uniform magnetic field associated with the shape demagnetizing effects. Flow is also observed in the fully filled ferrofluid sphere (1-20 mm/s) when the field is made non-uniform by adding a permanent magnet or a DC or AC excited small solenoidal coil. This confirms that a non-uniform magnetic field or a non-uniform distribution of magnetization due to a non-uniform magnetic field are causes of spin-up flow in ferrofluids with no free surface, while tangential magnetic surface stress contributes to flow in the presence of a free surface.Recent work has fitted velocity flow measurements of ferrofluid filled finite height cylinders with no free surface, subjected to uniform rotating magnetic fields, neglecting the container shape effects which cause non-uniform demagnetizing fields, and resulting in much larger non-physical effective values of spin viscosity η′∼10⁻⁸−10⁻¹² N s than those obtained from theoretical spin diffusion analysis where η′≤10⁻¹⁸ N s. COMSOL Multiphysics finite element computer simulations of spherical geometry in a uniform rotating magnetic field using non-physically large experimental fit values of spin viscosity η′∼10⁻⁸−10⁻¹² N s with a zero spin-velocity boundary condition at the outer wall predicts measureable flow, while simulations setting spin viscosity to zero (η′=0) results in negligible flow, in agreement with the ultrasound velocimetry measurements. COMSOL simulations also confirm that a non-uniform rotating magnetic field or a uniform rotating magnetic field with a non-uniform distribution of magnetization due to an external magnet or a current carrying coil can drive a measureable flow in an infinitely long ferrofluid cylinder with zero spin viscosity (η′=0).     

Modeling nongray gas-phase and soot radiation in luminous turbulent nonpremixed jet flames

Much progress has been made in radiative heat transfer modeling with respect to treatment of nongray radiation from both gas-phase species and soot particles, while radiation modeling in turbulent flame simulations is still in its infancy. Aiming at reducing this gap, this paper introduces state-of-the-art models of gas-phase and soot radiation to turbulent flame simulations. The full-spectrum k-distribution method (Modest, M.F., 2003, Journal of Quantitative Spectroscopy & Radiative Transfer, 76, 69–83) is implemented into a three-dimensional unstructured CFD code for nongray radiation modeling. The mixture full-spectrum k-distributions including nongray absorbing soot particles are constructed from a narrow-band k-distribution database created for individual gas-phase species, and an efficient scheme is employed for their construction in CFD simulations. A detailed reaction mechanism including NO _x and soot kinetics is used to predict flame structure, and a detailed soot model using a method of moments is employed to determine soot particle size distributions. A spherical-harmonic P₁ approximation is invoked to solve the radiative transfer equation. An oxygen-enriched, turbulent, nonpremixed jet flame is simulated, which features large concentrations of gas-phase radiating species and soot particles. Nongray soot modeling is shown to be of greater importance than nongray gas modeling in sooty flame simulations, with gray soot models producing large errors. The nongray treatment of soot strongly influences flame temperatures in the upstream and the flame-tip region and is essential for accurate predictions of NO. The nongray treatment of gases, however, weakly influences upstream flame temperatures and, therefore, has only a small effect on NO predictions. The effect of nongray soot radiation on flame temperature is also substantial in downstream regions where the soot concentration is small. Limitations of the P₁ approximation are discussed for the jet flame configuration; the P₁ approximation yields large errors in the spatial distribution of the computed radiative heat flux for highly anisotropic radiation fields such as those in flames with localized, near-opaque soot regions.     

Momentum transport and torque scaling in Taylor-Couette flow from an analogy with turbulent convection

We generalize an analogy between rotating and stratified shear flows. This analogy is summarized in Table 1. We use this analogy in the unstable case (centrifugally unstable flow vs. convection) to compute the torque in Taylor-Couette configuration, as a function of the Reynolds number. At low Reynolds numbers, when most of the dissipation comes from the mean flow, we predict that the non-dimensional torque G = T/ν² L, where L is the cylinder length, scales with Reynolds number R and gap width η, G = 1.46η^3/2(1 - η)^-7/4 R ^3/2. At larger Reynolds number, velocity fluctuations become non-negligible in the dissipation. In these regimes, there is no exact power law dependence the torque versus Reynolds. Instead, we obtain logarithmic corrections to the classical ultra-hard (exponent 2) regimes: G = 0.50 . These predictions are found to be in excellent agreement with avail-able experimental data. Predictions for scaling of velocity fluctuations are also provided. Received 7 June 2001 and Received in final form 7 December 2001