Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

SOLUTIONS OF A SYSTEM OF FORCED BURGERS EQUATION IN TERMS OF GENERALIZED LAGUERRE POLYNOMIALS

In this article, we obtain explicit solutions of a linear PDE subject to a class of radial square integrable functions with a monotonically increasing weight function |x|(n-1)e(beta vertical bar x vertical bar 2)/2, beta >= 0, x is an element of R-n. This linear PDE is obtained from a system of forced Burgers equation via the Cole-Hopf transformation. For any spatial dimension n > 1, the solution is expressed in terms of a family of weighted generalized Laguerre polynomials. We also discuss the large time behaviour of the solution of the system of forced Burgers equation.

Using heat treatment effects and EBSD analysis to tailor microstructure of hybrid Mg nanocomposite for enhanced overall mechanical response

In this study, a detailed investigation on the effect of heat treatment on the microstructural characteristics, texture evolution and mechanical properties of Mg-(5.6Ti+2.5B(4)C)(BM) hybrid nanocomposite is presented. Optimised heat treatment parameters, namely, heat treatment temperature and heat treatment time, were first identified through grain size and microhardness measurements. Initially, heat treatment of composites was conducted at temperature range between 100 and 300 degrees C for 1 h. Based on optical microscopic analysis and microhardness measurements, it was evident that significant grain growth and reduction in microhardness occurred for temperatures > 200 degrees C. The cutoff temperature that caused significant grain growth/matrix softening was thus identified. Second, at constant temperature (200 degrees C), the effect of variation of heat treatment time was carried out (ranging between 1 and 5 h) so as to identify the range wherein increase in average grain size and reduction in microhardness occurred. Furthering the study, the effect of optimised heat treatment parameters (200 degrees C, 5 h) on the microstructural texture evolution and hence, on the tensile and compressive properties of the Mg-(5.6Ti+2.5B(4)C)(BM) hybrid nanocomposite was carried out. From electron backscattered diffraction (EBSD) analysis, it was identified that the optimised heat treatment resulted in recrystallisation and residual stress relaxation, as evident from the presence of similar to 87% strain free grains, when compared to that observed in the non-heat treated/as extruded condition (i.e. 2.2 times greater than in the as extruded condition). For the heat treated composite, under both tensile and compressive loads, a significant improvement in fracture strain values (similar to 60% increase) was observed when compared to that of the non-heat treated counterpart, with similar to 20% reduction in yield strength. Based on structure-property correlation, the change in mechanical characteristics is identified to be due to: (1) the presence of less stressed matrix/reinforcement interface due to the relief of residual stresses and (2) texture weakening due to matrix recrystallisation effects, both arising due to heat treatment.

An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale

In this paper we present a massively parallel open source solver for Richards equation, named the RichardsFOAM solver. This solver has been developed in the framework of the open source generalist computational fluid dynamics tool box OpenFOAM (R) and is capable to deal with large scale problems in both space and time. The source code for RichardsFOAM may be downloaded from the CPC program library website. It exhibits good parallel performances (up to similar to 90% parallel efficiency with 1024 processors both in strong and weak scaling), and the conditions required for obtaining such performances are analysed and discussed. These performances enable the mechanistic modelling of water fluxes at the scale of experimental watersheds (up to few square kilometres of surface area), and on time scales of decades to a century. Such a solver can be useful in various applications, such as environmental engineering for long term transport of pollutants in soils, water engineering for assessing the impact of land settlement on water resources, or in the study of weathering processes on the watersheds. (C) 2014 Elsevier B.V. All rights reserved.

Modelling of transient heat conduction with diffuse interface methods

We present a survey on different numerical interpolation schemes used for two-phase transient heat conduction problems in the context of interface capturing phase-field methods. Examples are general transport problems in the context of diffuse interface methods with a non-equal heat conductivity in normal and tangential directions to the interface. We extend the tonsorial approach recently published by Nicoli M et al (2011 Phys. Rev. E 84 1-6) to the general three-dimensional (3D) transient evolution equations. Validations for one-dimensional, two-dimensional and 3D transient test cases are provided, and the results are in good agreement with analytical and numerical reference solutions.

Discerning heat transfer in building materials

The function of a building is to ensure safety and thermal comfort for healthy living conditions. Buildings primarily comprise an envelope, which acts as an interface separating the external environment from the indoors environment. The building envelope is primarily responsible for regulating indoor thermal comfort in response to external climatic conditions. It usually comprises a configuration of building materials to thus far provide requisite structural performance. However, studies into building-envelope configurations to provide a particular thermal performance are limited. As the building envelope is exposed to the external environment there will be heat and moisture transfer to the indoor environment through it. The overall phenomenon of heat and moisture transfer depends on the microstructure and configuration within the building material. Further, thermal property of a material is generally dependent on its microstructure, which comprises a network of pores and particles arranged in a definite structure. Thermal behaviour of a building material thus depends on the thermal conductivities of the solid particles, pore micro-structure and its constituent fluid (air and/or moisture). The thermal response of a building envelope is determined by the thermal characteristics of the individual building materials and its configuration. Understanding the heat transfer influenced by the complex networks of pores and particles is a relatively new study in the area of building climatic-response. The current study reviews the heat-transfer mechanisms that determine the thermal performance of a building material attributed to its micro-structure. A theoretical basis for the same is being evolved and its relevance in regulating heat-transfer through building envelopes, walls in particular, is reviewed in this paper. (C) 2014 N.C. Balaji. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).

Heat kernels on the AdS(2) cone and logarithmic corrections to extremal black hole entropy

We develop new techniques to efficiently evaluate heat kernel coefficients for the Laplacian in the short-time expansion on spheres and hyperboloids with conical singularities. We then apply these techniques to explicitly compute the logarithmic contribution to black hole entropy from an N = 4 vector multiplet about a Z(N) orbifold of the near-horizon geometry of quarter-BPS black holes in N = 4 supergravity. We find that this vanishes, matching perfectly with the prediction from the microstate counting. We also discuss possible generalisations of our heat kernel results to higher-spin fields over ZN orbifolds of higher-dimensional spheres and hyperboloids.

Large carbon cluster thin film gauges for measuring aerodynamic heat transfer rates in hypersonic shock tunnels

Different types of Large Carbon Cluster (LCC) layers are synthesized by a single-step pyrolysis technique at various ratios of precursor mixture. The aim is to develop a fast responsive and stable thermal gauge based on a LCC layer which has relatively good electrical conduction in order to use it in the hypersonic flow field. The thermoelectric property of the LCC layer has been studied. It is found that these carbon clusters are sensitive to temperature changes. Therefore suitable thermal gauges were developed for blunt cone bodies and were tested in hypersonic shock tunnels at a flow Mach number of 6.8 to measure aerodynamic heating. The LCC layer of this thermal gauge encounters high shear forces and a hostile environment for test duration in the range of a millisecond. The results are favorable to use large carbon clusters as a better sensor than a conventional platinum thin film gauge in view of fast responsiveness and stability.

Assessment of land surface temperature and heat fluxes over Delhi using remote sensing data

Surface energy processes has an essential role in urban weather, climate and hydrosphere cycles, as well in urban heat redistribution. The research was undertaken to analyze the potential of Landsat and MODIS data in retrieving biophysical parameters in estimating land surface temperature & heat fluxes diurnally in summer and winter seasons of years 2000 and 2010 and understanding its effect on anthropogenic heat disturbance over Delhi and surrounding region. Results show that during years 2000-2010, settlement and industrial area increased from 5.66 to 11.74% and 4.92 to 11.87% respectively which in turn has direct effect on land surface temperature (LST) and heat fluxes including anthropogenic heat flux. Based on the energy balance model for land surface, a method to estimate the increase in anthropogenic heat flux (Has) has been proposed. The settlement and industrial areas has higher amounts of energy consumed and has high values of Has in all seasons. The comparison of satellite derived LST with that of field measured values show that Landsat estimated values are in close agreement within error of 2 degrees C than MODIS with an error of 3 degrees C. It was observed that, during 2000 and 2010, the average change in surface temperature using Landsat over settlement & industrial areas of both seasons is 1.4 degrees C & for MODIS data is 3.7 degrees C. The seasonal average change in anthropogenic heat flux (Has) estimated using Landsat & MODIS is up by around 38 W/m(2) and 62 W/m(2) respectively while higher change is observed over settlement and concrete structures. The study reveals that the dynamic range of Has values has increased in the 10 year period due to the strong anthropogenic influence over the area. The study showed that anthropogenic heat flux is an indicator of the strength of urban heat island effect, and can be used to quantify the magnitude of the urban heat island effect. (C) 2013 Elsevier Ltd. All rights reserved.

Fragility correlates thermodynamic and kinetic properties of glass forming liquids

In our earlier communication we proposed a simple fragility determining function, (NBO]/(VmTg)-T-3), which we have now used to analyze several glass systems using available thermal data. A comparison with similar fragility determining function, Delta C-p/C-p(1), introduced by Chryssikos et al. in their investigation of lithium borate glasses has also been performed and found to be more convenient quantity for discussing fragilities. We now propose a new function which uses both Delta C-p and Delta T and which gives a numerical fragility parameter, F whose value lies between 0 and 1 for glass forming liquids. F can be calculated through the use of measured thermal parameters Delta C-p, C-p(1), T-g and T-m. Use of the new fragility values in reduced viscosity equation reproduces the whole range of viscosity curves of the Angell plot. The reduced viscosity equation can be directly compared with the Adam-Gibbs viscosity equation and a heat capacity function can be formulated which reproduces satisfactorily the Delta C-p versus In(T-r) curves and hence the configurational entropy. (C) 2014 Elsevier Ltd. All rights reserved.

Fragility correlates thermodynamic and kinetic properties of glass forming liquids

In our earlier communication we proposed a simple fragility determining function, (NBO]/(VmTg)-T-3), which we have now used to analyze several glass systems using available thermal data. A comparison with similar fragility determining function, Delta C-p/C-p(1), introduced by Chryssikos et al. in their investigation of lithium borate glasses has also been performed and found to be more convenient quantity for discussing fragilities. We now propose a new function which uses both Delta C-p and Delta T and which gives a numerical fragility parameter, F whose value lies between 0 and 1 for glass forming liquids. F can be calculated through the use of measured thermal parameters Delta C-p, C-p(1), T-g and T-m. Use of the new fragility values in reduced viscosity equation reproduces the whole range of viscosity curves of the Angell plot. The reduced viscosity equation can be directly compared with the Adam-Gibbs viscosity equation and a heat capacity function can be formulated which reproduces satisfactorily the Delta C-p versus In(T-r) curves and hence the configurational entropy. (C) 2014 Elsevier Ltd. All rights reserved.

Analysis of Entropy Generation During Mixed Convective Heat Transfer of Nanofluids Past a Rotating Circular Cylinder

The entropy generation due to mixed convective heat transfer of nanofluids past a rotating circular cylinder placed in a uniform cross stream is investigated via streamline upwind Petrov-Galerkin based finite element method. Nanosized copper (Cu) particles suspended in water are used with Prandtl number (Pr)=6.9. The computations are carried out at a representative Reynolds number (Re) of 100. The dimensionless cylinder rotation rate, a, is varied between 0 and 2. The range of nanoparticle volume fractions (phi) considered is 0 <= phi <= 5%. Effect of aiding buoyancy is brought about by considering two fixed values of the Richardson number (Ri) as 0.5 and 1.0. A new model for predicting the effective viscosity and thermal conductivity of dilute suspensions of nanoscale colloidal particles is presented. The model addresses the details of the agglomeration-deagglomeration in tune with the pertinent variations in the effective particulate dimensions, volume fractions, as well as the aggregate structure of the particulate system. The total entropy generation is found to decrease sharply with cylinder rotation rates and nanoparticle volume fractions. Increase in nanoparticle agglomeration shows decrease in heat transfer irreversibility. The Bejan number falls sharply with increase in alpha and phi.

REDUCED ORDER MODELLING OF COMBUSTION INSTABILITY IN A BACKWARD FACING STEP COMBUSTOR

This paper develops a fully coupled time domain Reduced Order Modelling (ROM) approach to model unsteady combustion dynamics in a backward facing step combustor The acoustic field equations are projected onto the canonical acoustic eigenmodes of the systems to obtain a coupled system of modal evolution equations. The heat release response of the flame is modelled using the G-equation approach. Vortical velocity fluctuations that arise due to shear layer rollup downstream of the step are modelled using a simplified 1D-advection equation whose phase speed is determined from a linear, local, temporal stability analysis of the shear layer just downstream of the step. The hydrodynamic stability analysis reveals a abrupt change in the value of disturbance phase speed from unity for Re < Re-crit to 0.5 for Re > Re-crit, where Remit for the present geometry was found to be approximate to 10425. The results for self-excited flame response show highly wrinkled flame shapes that are qualitatively similar to those seen in prior experiments of acoustically forced flames. The effect of constructive and destructive interference between the two contributions to flame surface wrinkling results in high amplitude wrinkles for the case when K-c -> 1.

L-p estimates for the wave equation associated to the Grushin operator

We prove that the solution of the wave equation associated to the Grushin operator G = -Delta -vertical bar x vertical bar(2)partial derivative(2)(t) is bounded on L-P (Rn+1), with 1 < p < infinity, when vertical bar 1/p - 1/2 vertical bar < 1/n+2.