Grain Floatation During Equiaxed Solidification of an Al-Cu Alloy in a Side-Cooled Cavity: Part II-Numerical Studies

In this article, a single-phase, one-domain macroscopic model is developed for studying binary alloy solidification with moving equiaxed solid phase, along with the associated transport phenomena. In this model, issues such as thermosolutal convection, motion of solid phase relative to liquid and viscosity variations of the solid-liquid mixture with solid fraction in the mobile zone are taken into account. Using the model, the associated transport phenomena during solidification of Al-Cu alloys in a rectangular cavity are predicted. The results for temperature variation, segregation patterns, and eutectic fraction distribution are compared with data from in-house experiments. The model predictions compare well with the experimental results. To highlight the influence of solid phase movement on convection and final macrosegregation, the results of the current model are also compared with those obtained from the conventional solidification model with stationary solid phase. By including the independent movement of the solid phase into the fluid transport model, better predictions of macrosegregation, microstructure, and even shrinkage locations were obtained. Mechanical property prediction models based on microstructure will benefit from the improved accuracy of this model.

Grain Floatation During Equiaxed Solidification of an Al-Cu Alloy in a Side-Cooled Cavity: Part I-Experimental Studies

Experimental studies were performed to investigate the role and influence of grain movement on macrosegregation and microstructure evolution during equiaxed solidification. Casting experiments were performed with a grain-refined Al-Cu alloy in a rectangular sand mold. For the aluminum alloy studied, the equiaxed grains are lighter than the bulk melt and thus float up. Experiments were designed to investigate floatation phenomena of equiaxed grains in the presence of thermosolutal convection. Cooling curves were recorded at key locations in both the casting and the chill. Quantitative image analysis and spatial chemical analysis were performed on the solidified casting to observe the chemical and microstructural inhomogeneity created by the melt convection and solid floatation. Several notable features that can be attributed to grain movement were observed in temperature histories, macrosegregation patterns, and microstructures. In our experiments, the floatation of grains influences the thermal conditions and the overall flow direction in the casting cavity. In some cases, the induced flow resulting from the grain movement caused a flow reversal. This in turn influences the solidification direction, microstructure evolution, and the overall macrosegregation behavior.

Thermal instability in a three-dimensional rigid container with prescribed heat flux at lower boundary

The Bénard–Marangoni convection is studied in a three-dimensional container with thermally insulated lateral walls and prescribed heat flux at lower boundary. The upper surface of the incompressible, viscous fluid is assumed to be flat with temperature dependent surface tension. A Galerkin–Tau method with odd and even trial functions satisfying all the essential boundary conditions except the natural boundary conditions at the free surface has been used to solve the problem. The critical Marangoni and Rayleigh numbers are determined for the onset of steady convection as a function of aspect ratios x0 and y0 for the cases of Bénard–Marangoni, pure Marangoni and pure Bénard convections. It is observed that critical parameters are decreasing with an increase in aspect ratios. The flow structures corresponding to the values of the critical parameters are presented in all the cases. It is observed that the critical parameters are higher for case with heat flux prescribed than those corresponding to the case with prescribed temperature. The critical Marangoni number for pure Marangoni convection is higher than critical Rayleigh number corresponding to pure Bénard convection for a given aspect ratio whereas the reverse was observed for two-dimensional infinite layer.

Analytical solutions for heat transfer during cyclic melting and freezing of a phase change material used in electronic or electrical packaging

In this paper we develop an analytical heat transfer model, which is capable of analyzing cyclic melting and solidification processes of a phase change material used in the context of electronics cooling systems. The model is essentially based on conduction heat transfer, with treatments for convection and radiation embedded inside. The whole solution domain is first divided into two main sub-domains, namely, the melting sub-domain and the solidification sub-domain. Each sub-domain is then analyzed for a number of temporal regimes. Accordingly, analytical solutions for temperature distribution within each subdomain are formulated either using a semi-infinity consideration, or employing a method of quasi-steady state, depending on the applicability. The solution modules are subsequently united, leading to a closed-form solution for the entire problem. The analytical solutions are then compared with experimental and numerical solutions for a benchmark problem quoted in the literature, and excellent agreements can be observed.

Dewetting on the Surface of Rutile

After annealing a continuous SiO2 film on the (001) surface of TiO2, the film dewets and then spreads to form a complex pattern. The final droplet morphology displays a densely branching morphology similar to those seen in computer-simulated models. It is proposed that Bénard-Marangoni convection cells form within the film before dewetting occurs. The formation of Bénard-Marangoni convection cells prior to dewetting results in the uniform size and spacing of the droplets on the surface. These convection cells form at temperature when the TiO2 substrate dissolves into the SiO2 thin film. The change in composition results in regions of differing surface tensions and therefore leads to the formation of the convection cells.

The origin of the solar magnetic cycle

After summarizing the relevant observational data, we discuss how a study of flux tube dynamics in the solar convection zone helps us to understand the formation of sunspots. Then we introduce the flux transport dynamo model and assess its success in modelling both the solar cycle and its departures from strictly periodic behaviour.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter=9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 108. Under these conditions the convection is turbulent, and the time-averaged velocity at any point is ‘zero’. The Reynolds number based on the Taylor microscale, Reλ, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6–7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as Ra1/2Sc1/2, and the Reynolds number would scale as Ra1/2Sc−1/2. The velocity and the flux measurements appear to be consistent with the Ra1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are much higher compared to what would be obtained in Rayleigh–Bénard (R–B) convection for similar density differences.

On breaks of the Indian monsoon

For over a century, the term break has been used for spells in which the rainfall over the Indian monsoon zone is interrupted. The phenomenon of 'break monsoon' is of great interest because long intense breaks are often associated with poor monsoon seasons. Such breaks have distinct circulation characteristics (heat trough type circulation) and have a large impact on rainfed agriculture. Although interruption of the monsoon rainfall is considered to be the most important feature of the break monsoon, traditionally breaks have been identified on the basis of the surface pressure and wind patterns over the Indian region. We have defined breaks (and active spells) on the basis of rainfall over the monsoon zone. The rainfall criteria are chosen so as to ensure a large overlap with the traditional breaks documented by Ramamurthy (1969) and De et al (1998). We have identified these rainbreaks for 1901-89. We have also identified active spells on the basis of rainfall over the Indian monsoon zone. We have shown that the all-India summer monsoon rainfall is significantly negatively correlated with the number of rainbreak days (correlation coefficient -0.56) and significantly positively correlated with the number of active days (correlation coefficient 0.47). Thus the interannual variation of the all-India summer monsoon rainfall is shown to be related to the number of days of rainbreaks and active spells identified here. There have been several studies of breaks (and also active spells in several cases) identified on the basis of different criteria over regions differing in spatial scales (e.g., Webster et al 1998; Krishnan et al 2000; Goswami and Mohan 2000; and Annamalai and Slingo 2001). We find that there is considerable overlap between the rainbreaks we have identified and breaks based on the traditional definition. There is some overlap with the breaks identified by Krishnan et al (2000) but little overlap with breaks identified by Webster et al (1998). Further, there are three or four active-break cycles in a season according to Webster et al (1998) which implies a time scale of about 40 days for which Goswami and Mohan (2000), and Annamalai and Slingo'(2001) have studied breaks and active minus break fluctuations. On the other hand, neither the traditional breaks (Ramamurthy 1969; and De et al 1998) nor the rainbreaks occur every year. This suggests that the 'breaks' in these studies axe weak spells of the intraseasonal variation of the monsoon, which occur every year. We have derived the OLR and circulation patterns associated with rainbreaks and active spells and compared them with the patterns associated with breaks/active minus break spells from these studies. Inspite of differences in the patterns over the Indian region, there is one feature which is seen in the OLR anomaly patterns of breaks identified on the basis of different criteria as well as the rainbreaks identified in this paper viz., a quadrapole over the Asia-west Pacific region arising from anomalies opposite (same) in sign to those over the Indian region occurring over the equatorial Indian Ocean and northern tropical (equatorial) parts of the west Pacific. Thus it appears that this quadrapole is a basic feature of weak spells of the intraseasonal,variation over the Asia-west Pacific region. Since the rainbreaks are intense weak spells, this basic feature is also seen in the composite patterns of these breaks. We find that rainbreaks (active spells) are also associated with negative (positive) anomalies over a part of the cast Pacific suggesting that the convection over the Indian region is linked to that over the east Pacific not only on the interannual scale (as evinced by the link between the Indian summer monsoon rainfall and ENSO) but on the intraseasonal scale as well.

Some salient features of the atmosphere observed over the north Bay of Bengal during BOBMEX

This paper describes the near surface characteristics and vertical variations based on the observations made at 17.5degreesN and 89degreesE from ORV Sagar Kanya in the north Bay of Bengal during the Bay of Bengal Monsoon Experiment (BOBMEX) carried out in July-August 1999. BOBMEX captured both the active and weak phases of convection. SST remained above the convection threshold throughout the BOBMEX. While the response of the SST to atmospheric forcing was clearly observed, the response of the atmosphere to SST changes was not clear. SST decreased during periods of large scale precipitation, and increased during a weak phase of convection. It is shown that the latent heat flux at comparable wind speeds was about 25-50% lower over the Bay during BOBMEX compared to that over the Indian Ocean during other seasons and tropical west Pacific. On the other hand, the largest variations in the surface daily net heat flux are observed over the Bay during BOBMEX. SST predicted using observed surface fluxes showed that 1-D heat balance model works sometime but not always, and horizontal advection is important. The high resolution Vaisala radiosondes launched during BOBMEX could clearly bring out the changes in the vertical structure of the atmosphere between active and weak phases of convection. Convective Available Potential Energy of the surface air decreased,by 2-3 kJ kg(-1) following convection, and recovered in a time period of one or two days. The mid tropospheric relative humidity and water vapor content, and wind direction show the major changes between the active and weak phases of convection.

The Indian monsoon and its variability

For over 300 years, the monsoon has been viewed as a gigantic land-sea breeze. It is shown in this paper that satellite and conventional observations support an alternative hypothesis, which considers the monsoon as a manifestation of seasonal migration of the intertropical convergence zone (ITCZ). With the focus on the Indian monsoon, the mean seasonal pattern is described, and why it is difficult to simulate it is discussed. Some facets of the intraseasonal variation, such as active-weak cycles; break monsoon; and a special feature of intraseasonal variation over the region, namely, poleward propagations of the ITCZ at intervals of 2-6 weeks, are considered. Vertical moist stability is shown to be a key parameter in the variation of monthly convection over ocean and land as well as poleward propagations. Special features of the Bay of Bengal and the monsoon brought out by observations during a national observational experiment in 1999 are briefly described.

Laboratory simulations show diabatic heating drives cumulus-cloud evolution and entrainment

Clouds are the largest source of uncertainty in climate science, and remain a weak link in modeling tropical circulation. A major challenge is to establish connections between particulate microphysics and macroscale turbulent dynamics in cumulus clouds. Here we address the issue from the latter standpoint. First we show how to create bench-scale flows that reproduce a variety of cumulus-cloud forms (including two genera and three species), and track complete cloud life cycles-e.g., from a ``cauliflower'' congestus to a dissipating fractus. The flow model used is a transient plume with volumetric diabatic heating scaled dynamically to simulate latent-heat release from phase changes in clouds. Laser-based diagnostics of steady plumes reveal Riehl-Malkus type protected cores. They also show that, unlike the constancy implied by early self-similar plume models, the diabatic heating raises the Taylor entrainment coefficient just above cloud base, depressing it at higher levels. This behavior is consistent with cloud-dilution rates found in recent numerical simulations of steady deep convection, and with aircraft-based observations of homogeneous mixing in clouds. In-cloud diabatic heating thus emerges as the key driver in cloud development, and could well provide a major link between microphysics and cloud- scale dynamics.

Filtered density function for large eddy simulation of turbulent reacting flows

A methodology termed the “filtered density function” (FDF) is developed and implemented for large eddy simulation (LES) of chemically reacting turbulent flows. In this methodology, the effects of the unresolved scalar fluctuations are taken into account by considering the probability density function (PDF) of subgrid scale (SGS) scalar quantities. A transport equation is derived for the FDF in which the effect of chemical reactions appears in a closed form. The influences of scalar mixing and convection within the subgrid are modeled. The FDF transport equation is solved numerically via a Lagrangian Monte Carlo scheme in which the solutions of the equivalent stochastic differential equations (SDEs) are obtained. These solutions preserve the Itô-Gikhman nature of the SDEs. The consistency of the FDF approach, the convergence of its Monte Carlo solution and the performance of the closures employed in the FDF transport equation are assessed by comparisons with results obtained by direct numerical simulation (DNS) and by conventional LES procedures in which the first two SGS scalar moments are obtained by a finite difference method (LES-FD). These comparative assessments are conducted by implementations of all three schemes (FDF, DNS and LES-FD) in a temporally developing mixing layer and a spatially developing planar jet under both non-reacting and reacting conditions. In non-reacting flows, the Monte Carlo solution of the FDF yields results similar to those via LES-FD. The advantage of the FDF is demonstrated by its use in reacting flows. In the absence of a closure for the SGS scalar fluctuations, the LES-FD results are significantly different from those based on DNS. The FDF results show a much closer agreement with filtered DNS results. © 1998 American Institute of Physics.

Command Tracking in High Performance Aircrafts: A New Dynamic Inversion Design

This paper proposes a new straight forward technique based on dynamic inversion, which is applied for tracking the pilot commands in high performance aircrafts.Pilot commands assumed in longitudinal mode are normal acceleration and total velocity(while roll angle and lateral acceleration are maintained at zero). In lateral mode, roll rate and total velocity are used as pilot commands (while climb rate and lateral acceleration are maintained at zero). Ensuring zero lateral acceleration leads to a better turn co-ordination. A six degree-of-freedom model of F-16 aircraft is used for both control design as well as simulation studies. Promising results are obtained which are found to be superior as compared to an existing approach (which is also based on dynamic inversion). The new approach has two potential benefits, namely reduced oscillatory response and reduced control magnitude. Another advantage of this approach is that it leads to a significant reduction of tuning parameters in the control design process.

Model-following Neuro-adaptive approach for robust control of high performance aircrafts

Based on dynamic inversion, a relatively straightforward approach is presented in this paper for nonlinear flight control design of high performance aircrafts, which does not require the normal and lateral acceleration commands to be first transferred to body rates before computing the required control inputs. This leads to substantial improvement of the tracking response. Promising results are obtained from six degree-offreedom simulation studies of F-16 aircraft, which are found to be superior as compared to an existing approach (which is also based on dynamic inversion). The new approach has two potential benefits, namely reduced oscillatory response (including elimination of non-minimum phase behavior) and reduced control magnitude. Next, a model-following neuron-adaptive design is augmented the nominal design in order to assure robust performance in the presence of parameter inaccuracies in the model. Note that in the approach the model update takes place adaptively online and hence it is philosophically similar to indirect adaptive control. However, unlike a typical indirect adaptive control approach, there is no need to update the individual parameters explicitly. Instead the inaccuracy in the system output dynamics is captured directly and then used in modifying the control. This leads to faster adaptation, which helps in stabilizing the unstable plant quicker. The robustness study from a large number of simulations shows that the adaptive design has good amount of robustness with respect to the expected parameter inaccuracies in the model.

Single Particle and Packed Bed Combustion in Modern Gasifier Stoves—Density Effects

This article is concerned with a study of an unusual effect due to density of biomass pellets in modern stoves based on close-coupled gasification-combustion process. The two processes, namely, flaming with volatiles and glowing of the char show different effects. The mass flux of the fuel bears a constant ratio with the air flow rate of gasification during the flaming process and is independent of particle density; char glowing process shows a distinct effect of density. The bed temperatures also have similar features: during flaming, they are identical, but distinct in the char burn (gasification) regime. For the cases, wood char and pellet char, the densities are 350, 990 kg/m(3), and the burn rates are 2.5 and 3.5 g/min with the bed temperatures being 1380 and 1502 K, respectively. A number of experiments on practical stoves showed wood char combustion rates of 2.5 +/- 0.5 g/min and pellet char burn rates of 3.5 +/- 0.5 g/min. In pursuit of the resolution of the differences, experimental data on single particle combustion for forced convection and ambient temperatures effects have been obtained. Single particle char combustion rate with air show a near-d(2) law and surface and core temperatures are identical for both wood and pellet char. A model based on diffusion controlled heat release-radiation-convection balance is set up. Explanation of the observed results needs to include the ash build-up over the char. This model is then used to explain observed behavior in the packed bed; the different packing densities of the biomass chars leading to different heat release rates per unit bed volume are deduced as the cause of the differences in burn rate and bed temperatures.