Local numerical modelling of magnetoconvection and turbulence - implications for mean-field theories

During the last decades mean-field models, in which large-scale magnetic fields and differential rotation arise due to the interaction of rotation and small-scale turbulence, have been enormously successful in reproducing many of the observed features of the Sun. In the meantime, new observational techniques, most prominently helioseismology, have yielded invaluable information about the interior of the Sun. This new information, however, imposes strict conditions on mean-field models. Moreover, most of the present mean-field models depend on knowledge of the small-scale turbulent effects that give rise to the large-scale phenomena. In many mean-field models these effects are prescribed in ad hoc fashion due to the lack of this knowledge. With large enough computers it would be possible to solve the MHD equations numerically under stellar conditions. However, the problem is too large by several orders of magnitude for the present day and any foreseeable computers. In our view, a combination of mean-field modelling and local 3D calculations is a more fruitful approach. The large-scale structures are well described by global mean-field models, provided that the small-scale turbulent effects are adequately parameterized. The latter can be achieved by performing local calculations which allow a much higher spatial resolution than what can be achieved in direct global calculations. In the present dissertation three aspects of mean-field theories and models of stars are studied. Firstly, the basic assumptions of different mean-field theories are tested with calculations of isotropic turbulence and hydrodynamic, as well as magnetohydrodynamic, convection. Secondly, even if the mean-field theory is unable to give the required transport coefficients from first principles, it is in some cases possible to compute these coefficients from 3D numerical models in a parameter range that can be considered to describe the main physical effects in an adequately realistic manner. In the present study, the Reynolds stresses and turbulent heat transport, responsible for the generation of differential rotation, were determined along the mixing length relations describing convection in stellar structure models. Furthermore, the alpha-effect and magnetic pumping due to turbulent convection in the rapid rotation regime were studied. The third area of the present study is to apply the local results in mean-field models, which task we start to undertake by applying the results concerning the alpha-effect and turbulent pumping in mean-field models describing the solar dynamo.

Smoldering combustion of "incense" sticks - Experiments and modeling

This paper is concerned with the experimental and modeling studies on the smoldering rates of incense sticks as a function of ambient oxygen fraction in air, the flow velocity and size. The experimental results are obtained both for forward and reverse smolder conditions. The results are explained on the basis of surface combustion due to diffusion of oxygen to the surface by both free and forced convection supporting the heat transfer into the solid by conduction, into the stream by convection and the radiant heat transfer from the surface. The heat release at the surface is controlled by the convective transport of the oxidizer to the surface. To obtain the diffusion rates particularly for the reverse smolder, CFD calculations of fluid flow with along with a passive scalar are needed; these calculations have been made both for forward and reverse smolder. The interesting aspect of the CFD calculations is that while the Nusselt umber for forward smolder shows a clear root( Re-u) dependence ( Re-u = Flow Reynolds Number), the result for reverse smolder shows a peak in the variation with Reynolds number with the values lower than for forward smolder and unsteadiness in the flow beyond a certain flow rate. The results of flow behavior and Nusselt number are used in a simple model for the heat transfer at the smoldering surface to obtain the dependence of the smoldering rate on the diameter of the incense stick, the flow rate of air and the oxygen fraction. The results are presented in terms of a correlation for the non-dimensional smoldering rate with radiant flux from the surface and heat generation rate at the surface. The correlations appear reasonable for both forward and reverse smolder cases.

Raman spectra of ferro-electric crystals Part IV. Lithium hydrazinium sulphate (LiN2H5SO4)

The Raman spectrum of lithium hydrazinium sulphate has been recorded both in the single crystal form and in aqueous solutions. The crystal exhibits thirty-eight Raman lines having the frequency shifts 52, 70, 104, 146, 174, 220, 260, 302, 350, 454, 470, 610, 630, 715, 977, 1094, 1115, 1132, 1177, 1191, 1260, 1444, 1493, 1577, 1630, 1670, 2205, 2484, 2553, 2655, 2734, 2848, 2894, 2939, 3028, 3132, 3290 and 3330 cm.−1 The aqueous solution gave rise to six Raman lines at 452, 980, 1050–1200, 1260, 1425 and 1570 cm.−1 apart from a maximum at 180 cm.−1 in the ‘wing’ accompanying the Rayleigh line. The observed Raman lines have been assigned as arising from the vibrations of the SO4 ion, N2H5+ ion, Li-O4 group, hydrogen bond and the lattice. The influence of the hydrogen bond on the N-H stretching vibrations has been pointed out. The various features of the observed spectrum strongly support the hypothesis that the NH3 group in the crystal is rotating around the N-N axis at room temperature.

Open and closed boundaries in large-scale convective dynamos

Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. To test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation level of large-scale dynamos. Open (vertical field) and closed (perfect conductor) boundary conditions are used for the magnetic field. The contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate in this regime increases proportional to Rm^(1/2). In the nonlinear regime, the saturation level of the energy of the mean magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

The buoyant plume above a heat source

The boundary-layer type conservation equations of mass, momentum and energy for the steady free turbulent flow in gravitational convection over heat sources are set up for both two-dimensional and axisymmetric cases. These are reduced to ordinary differential equations in a similarity parameter by suitable transformations. The three classical hypotheses of turbulent diffusion-the Constant Exchange Coefficient hypothesis, Prandtl's Momentum Transfer theory and Taylor's Vorticity Transfer theory-are then incorporated into these equations in succession. The resulting equations are solved numerically and the results compared with some experimental results on gravitational convection over heat sources reported by Rouse et al.

Large-scale structures in turbulent and reverse-transitional sink flow boundary layers

Aspects of large-scale organized structures in sink flow turbulent and reverse-transitional boundary layers are studied experimentally using hot-wire anemometry. Each of the present sink flow boundary layers is in a state of 'perfect equilibrium' or 'exact self-preservation' in the sense of Townsend (The Structure of Turbulent Shear Flow, 1st and 2nd edns, 1956, 1976, Cambridge University Press) and Rotta (Progr. Aeronaut. Sci., vol. 2, 1962, pp. 1-220) and conforms to the notion of 'pure wall-flow' (Coles, J. Aerosp. Sci., vol. 24, 1957, pp. 495-506), at least for the turbulent cases. It is found that the characteristic inclination angle of the structure undergoes a systematic decrease with the increase in strength of the streamwise favourable pressure gradient. Detectable wall-normal extent of the structure is found to be typically half of the boundary layer thickness. Streamwise extent of the structure shows marked increase as the favourable pressure gradient is made progressively severe. Proposals for the typical eddy forms in sink flow turbulent and reverse-transitional flows are presented, and the possibility of structural self-organization (i.e. individual hairpin vortices forming streamwise coherent hairpin packets) in these flows is also discussed. It is further indicated that these structural ideas may be used to explain, from a structural viewpoint, the phenomenon of soft relaminarization or reverse transition of turbulent boundary layers when subjected to strong streamwise favourable pressure gradients. Taylor's 'frozen turbulence' hypothesis is experimentally shown to be valid for flows in the present study even though large streamwise accelerations are involved, the flow being even reverse transitional in some cases. Possible conditions, which are required to be satisfied for the safe use of Taylor's hypothesis in pressure-gradient-driven flows, are also outlined. Measured convection velocities are found to be fairly close to the local mean velocities (typically 90% or more) suggesting that the structure gets convected downstream almost along with the mean flow.

On the use of a coordinate transformation for analysis of axisymmetric vibration of polar orthotropic annular plates

Estimates of natural frequencies corresponding to axisymmetric modes of flexural vibration of polar orthotropic annular plates have been obtained for various combinations of clamped, simply supported and free edge conditions. A coordinate transformation in the radial direction has been used to obtain effective solutions by the classical Rayleigh-Ritz method. The analysis with this transformation has been found to be advantageous in computations, particularly for large hole sizes, over direct analysis. Numerical results have been obtained for various values of hole sizes and rigidity ratio. The eigenvalue parameter has been found to vary more or less linearly with the rigidity ratio. A comparison with the results for isotropic plates has brought out some interesting features.

An Integrated Design Procedure for High Frequency AC Inductors for Inverters

Inductors are important energy storage elements that are used as filters in switching power converters. The operating efficiency of power inductors depend on the initial design choices and they remain as one of the most inefficient elements in a power converter. The focus of this paper is to explore the inductor design procedure from the point of efficiency and operating temperature. A modified form of the area product approach is used as starting point for the inductor design. The equations which estimate the power loss in core and copper winding are described. The surface temperature of the inductor is modelled using heat transfer equations for radiation and natural convection. All design assumptions are verified by actual experimental data and results show a good match with the analysis.

Robust/optimal temperature profile control using neural networks

An approximate dynamic programming (ADP) based neurocontroller is developed for a heat transfer application. Heat transfer problem for a fin in a car's electronic module is modeled as a nonlinear distributed parameter (infinite-dimensional) system by taking into account heat loss and generation due to conduction, convection and radiation. A low-order, finite-dimensional lumped parameter model for this problem is obtained by using Galerkin projection and basis functions designed through the 'Proper Orthogonal Decomposition' technique (POD) and the 'snap-shot' solutions. A suboptimal neurocontroller is obtained with a single-network-adaptive-critic (SNAC). Further contribution of this paper is to develop an online robust controller to account for unmodeled dynamics and parametric uncertainties. A weight update rule is presented that guarantees boundedness of the weights and eliminates the need for persistence of excitation (PE) condition to be satisfied. Since, the ADP and neural network based controllers are of fairly general structure, they appear to have the potential to be controller synthesis tools for nonlinear distributed parameter systems especially where it is difficult to obtain an accurate model.

Optimal Weighted Antenna Selection For Imperfect Channel Knowledge From Training

Receive antenna selection (AS) reduces the hardware complexity of multi-antenna receivers by dynamically connecting an instantaneously best antenna element to the available radio frequency (RF) chain. Due to the hardware constraints, the channels at various antenna elements have to be sounded sequentially to obtain estimates that are required for selecting the ``best'' antenna and for coherently demodulating data. Consequently, the channel state information at different antennas is outdated by different amounts. We show that, for this reason, simply selecting the antenna with the highest estimated channel gain is not optimum. Rather, the channel estimates of different antennas should be weighted differently, depending on the training scheme. We derive closed-form expressions for the symbol error probability (SEP) of AS for MPSK and MQAM in time-varying Rayleigh fading channels for arbitrary selection weights, and validate them with simulations. We then derive an explicit formula for the optimal selection weights that minimize the SEP. We find that when selection weights are not used, the SEP need not improve as the number of antenna elements increases, which is in contrast to the ideal channel estimation case. However, the optimal selection weights remedy this situation and significantly improve performance.

Aerodynamic noise from wave-turbulence interaction

Abstract. In order to estimate the acoustic energy scattered when a unit volume of free turbulence, such as in free jets, interacts with a plane steady sound wave, theoretical expressions are derived for two simple models of turbulence: eddy model and isotropic model. The effect of convection by mean motion of the energy-bearing eddies on the incident sound wave and on the sound generated from wave-turbulence interaction is taken into account. Finally, by means of a representative calculation,the directionality pattern and Mach number dependence of the noise so generated is discussed.

Towards understanding the unusual Indian monsoon in 2009

The Indian summer monsoon season of 2009 commenced with a massive deficit in all-India rainfall of 48% of the average rainfall in June. The all-India rainfall in July was close to the normal but that in August was deficit by 27%. In this paper, we first focus on June 2009, elucidating the special features and attempting to identify the factors that could have led to the large deficit in rainfall. In June 2009, the phase of the two important modes, viz., El Nino and Southern Oscillation (ENSO) and the equatorial Indian Ocean Oscillation (EQUINOO) was unfavourable. Also, the eastern equatorial Indian Ocean (EEIO) was warmer than in other years and much warmer than the Bay. In almost all the years, the opposite is true, i.e., the Bay is warmer than EEIO in June. It appears that this SST gradient gave an edge to the tropical convergence zone over the eastern equatorial Indian Ocean, in competition with the organized convection over the Bay. Thus, convection was not sustained for more than three or four days over the Bay and no northward propagations occurred. We suggest that the reversal of the sea surface temperature (SST) gradient between the Bay of Bengal and EEIO, played a critical role in the rainfall deficit over the Bay and hence the Indian region. We also suggest that suppression of convection over EEIO in association with the El Nino led to a positive phase of EQUINOO in July and hence revival of the monsoon despite the El Nino. It appears that the transition to a negative phase of EQUINOO in August and the associated large deficit in monsoon rainfall can also be attributed to the El Nino.

Circulation and salinity distribution in coastal inlets

The mathematical model developed by Hansen and Rattray based on Pritchard's equations for a coastal-plain estuary has been analysed to study the circulation and salinity distributions in coastal inlets with constant width and depth. Numerical solutions of the basic equations have been obtained without placing any restriction on Rayleigh numbers. A noteworthy contribution of the present analysis is that solutions of equations have been obtained for higher Rayleigh numbers, which was not possible in the earlier model. It is found that the effect of higher Rayleigh numbers is to increase the vertical advection, making the salinities in the upper and lower layers more uniform with a distinct halocline near the mid-depths. Solutions are discussed for some special cases of practical interest.

Optimal Receive Antenna Selection in Time-Varying Fading Channels with Practical Training Constraints

Hardware constraints, which motivate receive antenna selection, also require that various antenna elements at the receiver be sounded sequentially to obtain estimates required for selecting the `best' antenna and for coherently demodulating data thereafter. Consequently, the channel state information at different antennas is outdated by different amounts and corrupted by noise. We show that, for this reason, simply selecting the antenna with the highest estimated channel gain is not optimum. Rather, a preferable strategy is to linearly weight the channel estimates of different antennas differently, depending on the training scheme. We derive closed-form expressions for the symbol error probability (SEP) of AS for MPSK and MQAM in time-varying Rayleigh fading channels for arbitrary selection weights, and validate them with simulations. We then characterize explicitly the optimal selection weights that minimize the SEP. We also consider packet reception, in which multiple symbols of a packet are received by the same antenna. New suboptimal, but computationally efficient weighted selection schemes are proposed for reducing the packet error rate. The benefits of weighted selection are also demonstrated using a practical channel code used in third generation cellular systems. Our results show that optimal weighted selection yields a significant performance gain over conventional unweighted selection.

Origin of the plateau in the low-temperature thermal conductivity of silica

Thermal conductivities of glasses at low temperatures show strikingly similar behavior irrespective of their chemical composition. While for T<1 K the thermal conductivity can be understood in the phenomenological tunneling model; the ‘‘universal plateau’’ in the temperature interval 15>T>2 K is totally unexplained. While Rayleigh scattering of phonons by structural disorder should be the natural cause for limiting the mean free path of phonons in this temperature range, it has been concluded before that in glasses a strong enough source of such scattering does not exist. In this study we show by a proper structural analysis in at least one material (namely, silica) that a strong enough source of Rayleigh scattering of phonons in glasses does exist so that the ‘‘universal plateau’’ can be explained without invoking any new mechanism. This may be for the first time that the low-temperature property of a structural glass has been correlated to its structure.