Scaling of thermo-magnetic convection

The scaling to characterize unsteady boundary layer development for thermo-magnetic convection of paramagnetic fluids with the Prandtl number greater than one is developed. Under the consideration is a square cavity with initially quiescent isothermal fluid placed in microgravity condition (g = 0) and subject to a uniform, vertical gradient magnetic field. A distinct magnetic thermal-boundary layer is produced by sudden imposing of a higher temperature on the vertical sidewall and as an effect of magnetic body force generated on paramagnetic fluid. The transient flow behavior of the resulting boundary layer is shown to be described by three stages: the start-up stage, the transitional stage and the steady state. The scaling is verified by numerical simulations with the magnetic momentum parameter m variation and the parameter γRa variation.

An improved boundary layer scaling with ramp heating on a sloping plate

A scaling analysis for the natural convection boundary layer adjacent to an inclined semi-infinite plate subject to a non-instantaneous heating in the form of an imposed wall temperature which increases linearly up to a prescribed steady value over a prescribed time is reported. The development of the boundary layer flow from start-up to a steady-state has been described based on scaling analyses and verified by numerical simulations. The analysis reveals that, if the period of temperature growth on the wall is sufficiently long, the boundary layer reaches a quasi-steady mode before the growth of the temperature is completed. In this mode the thermal boundary layer at first grows in thickness and then contracts with increasing time. However, if the imposed wall temperature growth period is sufficiently short, the boundary layer develops differently, but after the wall temperature growth is completed, the boundary layer develops as though the startup had been instantaneous. The steady state values of the boundary layer for both cases are ultimately the same.

Natural convection in the attic-shaped spaces

Heat transfer through an attic space into or out of buildings is an important issue for attic-shaped houses in both hot and cold climates. One of the important objectives for design and construction of houses is to provide thermal comfort for occupants. In the present energy-conscious society, it is also a requirement for houses to be energy efficient, i.e. the energy consumption for heating or air-conditioning houses must be minimized. Relevant to these objectives, research into heat transfer in attics has been conducted for about three decades. The transient behaviour of an attic space is directly relevant to our daily life. Certain periods of the day or night may be considered as having a constant ambient temperature (e.g. during 11am - 2pm or 11pm - 2am). However, at other times during the day or night the ambient temperature changes with time (e.g. between 5am - 9am or 5pm - 9pm). Therefore, the analysis of steady state solution is not sufficient to describe the fluid flow and heat transfer in the attic space. The discussion of the transient development of the boundary is required. A theoretical understanding of the transient behaviour of the flow in the enclosure is performed through scaling analysis for sudden and ramp heating conditions. A proper identification of the timescales, the velocity and the thickness relevant to the flow that develops inside the cavity makes it possible to predict theoretically the basic flow features that will survive once the thermal flow in the enclosure reaches a steady state. Those scaling predictions have been verified by a series of numerical simulations.

The influence of the motion of a rotor grid on a stator grid

In this paper we have used the method of characteristics developed for two dimensional unsteady flow problems to study a simplified axial turbine problem. The system consists of two sets of blades —the guiding vanes which are fixed and the rotor blades which move perpendicular to these vanes. The initial undisturbed constant flow in the system is perturbed by introducing a small velocity normal to the rotor blades to simulate a slight constant inclination. The resulting perturbed flow is periodic after the first three cycles. We have studied the perturbed density distribution throughout the system during a period.

Effect of wave parameters on flood wave subsidence

A detailed analy~is on the propagation of a sinusoidal flood wave in a wide prismatic open channel b.as hen made by numc? ii.~ll~integrating We govemins nondimenional equations of unsteady flow in an open chamei. EmpE:dsis has been laid on the effect of wave parmefen on th propagation of 6.8 sinusoidal wave. Results show that the amount of subsidence is more in the case of small wave anplltude and wave duration cases. Further, wave duration has been noticed to have a relatively Vier influence on subsidence than wave amplitude. The speed at which the peak of the wave moves is observed to be a function of only the wave amplitude.

Parametric Study of Flood Wave Propagation

A parametric study of the flood wave propagation problem is made, based on numerical solution of the nondimensionalized unsteady flow equations of open channels. The propagation of a sinusoidal flood wave in a prismatic channel is studied for uniform initial flow. The governing parameters (initial uniform flow Froude number, wave amplitude, wave duration, channel width parameter and side slope) are varied over a wide range. In all, 49 cases are studied. Effects of these governing parameters on the subsidence of stage and discharge and the speed of the wave peak are described in detail. The relative wave amplitude is found to vary linearly with F0, the initial uniform flow froude number, for lower F0 values. Wave duration has a very pronounced effect on subsidence with greater subsidence at lower wave duration values.

Laminar and turbulent natural convection in rectangular cavities

A computer program has been developed for the prediction of buoyancy-driven laminar and turbulent flow in rectangular air-filled two-dimensional cavities with differentially heated side walls. Laminar flow predictions for a square cavity and Rayleigh numbers from Ra = 10^3 up to the onset of unsteady flow have been obtained. Accurate solutions for Ra = 5 x 10^6, 10^7, 5 x 10^7 and 10^8 are presented and an estimate for the critical Rayleigh number at which the steady laminar flow becomes unsteady is given for this geometry. Numerical predictions of turbulent flow have been obtained for RaH~0(10^9 -10^11 ) and compared with existing experimental data. A previously developed second moment closure model (Behnia et al. 1987) has been used to model the turbulence. Results indicate that a second moment closure model is capable of predicting the observed flow features.

Time-dependent analysis of an N2O gasdynamic laser

A simplified two-temperature model is presented for the vibrational energy levels of the N2O and N2 molecules of an N2O-N2-He gasdynamic laser (GDL), and the governing equations for the unsteady flow of the gas mixture in a convergent-divergent contour nozzle are solved using a time-dependent numerical technique. Final steady-state distributions are obtained for vibrational temperatures, population inversion, and the small-signal laser gain along the nozzle. It is demonstrated that, for plenum temperatures lower than 1200 K, an N2O GDL such as the present is more efficient than a CO2 GDL in identical operating conditions

An efficient marching algorithm for waterhammer analysis by the method of characteristics

A new fast and efficient marching algorithm is introduced to solve the basic quasilinear, hyperbolic partial differential equations describing unsteady, flow in conduits by the method of characteristics. The details of the marching method are presented with an illustration of the waterhammer problem in a simple piping system both for friction and frictionless cases. It is shown that for the same accuracy the new marching method requires fewer computational steps, less computer memory and time.

Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

Time-domain-finite-wave analysis of the engine exhaust system by means of the stationary-frame method of characteristics. Part I. Theory

Time-domain-finite-wave analysis of the engine exhaust system is usually done using the method of characteristics. This makes use of either the moving frame method, or the stationary frame method. The stationary frame method is more convenient than its counterpart inasmuch as it avoids the tedium of graphical computations. In this paper (part I), the stationary-frame computational scheme along with the boundary conditions has been implemented. The analysis of a uniform tube, cavity-pipe junction including the engine and the radiation ends, and also the simple area discontinuities has been presented. The analysis has been done accounting for wall friction and heat-transfer for a one-dimensional unsteady flow. In the process, a few inconsistencies in the formulations reported in the literature have been pointed out and corrected. In the accompanying paper (part II) results obtained from the simulation are shown to be in good agreement with the experimental observations.

A bicharacteristic formulation of the ideal MHD equations

On a characteristic surface Omega of a hyperbolic system of first-order equations in multi-dimensions (x, t), there exits a compatibility condition which is in the form of a transport equation along a bicharacteristic on Omega. This result can be interpreted also as a transport equation along rays of the wavefront Omega(t) in x-space associated with Omega. For a system of quasi-linear equations, the ray equations (which has two distinct parts) and the transport equation form a coupled system of underdetermined equations. As an example of this bicharacteristic formulation, we consider two-dimensional unsteady flow of an ideal magnetohydrodynamics gas with a plane aligned magnetic field. For any mode of propagation in this two-dimensional flow, there are three ray equations: two for the spatial coordinates x and y and one for the ray diffraction. In spite of little longer calculations, the final four equations (three ray equations and one transport equation) for the fast magneto-acoustic wave are simple and elegant and cannot be derived in these simple forms by use of a computer program like REDUCE.

Abel Inversion of Axially-Symmetric Shock Wave Flows

Finite-fringe interferograms produced for axisymmetric shock wave flows are analyzed by Fourier transform fringe analysis and an Abel inversion method to produce density field data for the validation of numerical models. For the Abel inversion process, we use basis functions to model phase data from axially-symmetric shock wave structure. Steady and unsteady flow problems are studied, and compared with numerical simulations. Good agreement between theoretical and experimental results is obtained when one set of basis functions is used during the inversion process, but the shock front is smeared when another is used. This is because each function in the second set of basis functions is infinitely differentiable, making them poorly-suited to the modelling of a step function as is required in the representation of a shock wave.

Grid-refined LES predictions for Fan-OGV broadband noise

This work forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. In this paper, we focus on LES calculations of noise sources on and close to a blade trailing edge. We consider two test cases; one an isolated NACA0012 airfoil in flow, and the other an industry-standard rotating fan. In the first case, turbulent mean and RMS velocities and energy spectra at different locations are compared with those from experiment. 1,2The sound generated by the unsteady pressure fluctuations on the airfoil surface and by the flow turbulence will be predicted using a Ffowcs Williams Hawkings (FW-H) surface. In the second case, unsteady flow and acoustic fields around the blade passage 3 are presented for a refined mesh, and the rotor-stator tonal noise will be predicted by using the rotor-wake mean velocity profile and the methodology described in Lloyd & Peake 4. Copyright © 2009 by Qinling Li, Nigel Peake & Mark Savill.

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.