First steps towards modeling a multi-scale Earth system

Recent advances in computational geodynamics are applied to explore the link between Earth’s heat, its chemistry and its mechanical behavior. Computational thermal-mechanical solutions are now allowing us to understand Earth patterns by solving the basic physics of heat transfer. This approach is currently used to solve basic convection patterns of terrestrial planets. Applying the same methodology to smaller scales delivers promising similarities between observed and predicted structures which are often the site of mineral deposits. The new approach involves a fully coupled solution to the energy, momentum and continuity equations of the system at all scales, allowing the prediction of fractures, shear zones and other typical geological patterns out of a randomly perturbed initial state. The results of this approach are linking a global geodynamic mechanical framework over regional-scale mineral deposits down to the underlying micro-scale processes. Ongoing work includes the challenge of incorporating chemistry into the formulation.

Steady natural convection of non-Newtonian power-law fluid in a trapezoidal enclosure

Numerically investigation of free convection heat transfer in a differentially heated trapezoidal cavity filled with non-Newtonian Power-law fluid has been performed in this study. The left inclined surface is uniformly heated whereas the right inclined surface is maintained as uniformly cooled. The top and bottom surfaces are kept adiabatic with initially quiescent fluid inside the enclosure. Finite volume based commercial software FLUENT 14.5 is used to solve the governing equations. Dependency of various flow parameters of fluid flow and heat transfer is analyzed including Rayleigh number, Ra ranging from 10^5 to 10^7, Prandtl number, Pr of 100 to 10,000 and power index, n of 0.6 to 1.4. Outcomes have been reported in terms of isotherms, streamline, and local Nusselt number for various Ra, Pr, n and inclined angles. Grid sensitivity analysis is performed and numerically obtained results have been compared with those results available in the literature and found good agreement.

Liquid-vapour phase change : nucleate boiling of pure fluid and nanofluid under different gravity levels

This research was a step towards the comprehension of the nano-particles interaction with bubbles created during boiling. It was aimed at solving the controversies of whether the heat transfer is enhanced or deteriorated during the boiling of the nanofluid. Experiments were conducted in normal gravity and reduced gravity environments on-board the European Space Agency Parabolic Flight Program. The local modification of the thermo-physical properties of the fluid and moreover the modification experienced in the liquid microlayer under the growing vapour bubble were the dominant factors in explaining the mechanisms of the boiling behaviour of the nanofluid.

Development of empirical equations for irradiance profile of a standard parabolic trough collector using Monte Carlo ray tracing technique

Irradiance profile around the receiver tube (RT) of a parabolic trough collector (PTC) is a key effect of optical performance that affects the overall energy performance of the collector. Thermal performance evaluation of the RT relies on the appropriate determination of the irradiance profile. This article explains a technique in which empirical equations were developed to calculate the local irradiance as a function of angular location of the RT of a standard PTC using a vigorously verified Monte Carlo ray tracing model. A large range of test conditions including daily normal insolation, spectral selective coatings and glass envelop conditions were selected from the published data by Dudley et al. [1] for the job. The R2 values of the equations are excellent that vary in between 0.9857 and 0.9999. Therefore, these equations can be used confidently to produce realistic non-uniform boundary heat flux profile around the RT at normal incidence for conjugate heat transfer analyses of the collector. Required values in the equations are daily normal insolation, and the spectral selective properties of the collector components. Since the equations are polynomial functions, data processing software can be employed to calculate the flux profile very easily and quickly. The ultimate goal of this research is to make the concentrating solar power technology cost competitive with conventional energy technology facilitating its ongoing research.

Thin single-walled carbon nanotubes with narrow chirality distribution : constructive interplay of plasma and Gibbs–Thomson effects

Multiscale, multiphase numerical modeling is used to explain the mechanisms of effective control of chirality distributions of single-walled carbon nanotubes in direct plasma growth and suggest effective approaches to further improvement. The model includes an unprecedented combination of the plasma sheath, ion/radical transport, species creation/loss, plasma–surface interaction, heat transfer, surface/bulk diffusion, graphene layer nucleation, and bending/lift-off modules. It is shown that the constructive interplay between the plasma and the Gibbs–Thomson effect can lead to the effective nucleation and lift-off of small graphene layers on small metal catalyst nanoparticles. As a result, much thinner nanotubes with narrower chirality distributions can nucleate at much lower process temperatures and pressures compared to thermal CVD. This approach is validated by a host of experimental results, substantially reduces the amounts of energy and atomic matter required for the nanotube growth, and can be extended to other nanoscale structures and materials systems, thereby nearing the ultimate goal of energy- and matter-efficient nanotechnology.

Capital cost model for Robert evaporators

ROBERT EVAPORATORS in Australian sugar factories are traditionally constructed with 44.45 mm outside diameter stainless steel tubes of ~2 m length for all stages of evaporation. There are a few vessels with longer tubes (up to 2.8 m) and smaller and larger diameters (38.1 and 50.8 mm). Queensland University of Technology is undertaking a study to investigate the heat transfer performance of tubes of different lengths and diameters for the whole range of process conditions typically encountered in the evaporator set. Incorporation of these results into practical evaporator designs requires an understanding of the cost implications for constructing evaporator vessels with calandrias having tubes of different dimensions. Cost savings are expected for tubes of smaller diameter and longer length in terms of material, labour and installation costs in the factory. However these savings must be considered in terms of the heat transfer area requirements for the evaporation duty, which will likely be a function of the tube dimensions. In this paper a capital cost model is described which provides a relative cost of constructing and installing Robert evaporators of the same heating surface area but with different tube dimensions. Evaporators of 2000, 3000, 4000 and 5000 m2 are investigated. This model will be used in conjunction with the heat transfer efficiency data (when available) to determine the optimum tube dimensions for a new evaporator at a specified evaporation duty. Consideration is also given to other factors such as juice residence time (and implications for sucrose degradation and control) and droplet de-entrainment in evaporators of different tube dimensions.

CFD convective flow simulation of the varying properties of CO2-H20 mixtures in geothermal systems

Numerical simulation of a geothermal reservoir, modelled as a bottom-heated square box, filled with water-CO2 mixture is presented in this work. Furthermore, results for two limiting cases of a reservoir filled with either pure water or CO2 are presented. Effects of different parameters including CO2 concentration as well as reservoir pressure and temperature on the overall performance of the system are investigated. It has been noted that, with a fixed reservoir pressure and temperature, any increase in CO2concentration leads to better performance, i.e. stronger convection and higher heat transfer rates. With a fixed CO2 concentration, however, the reservoir pressure and temperature can significantly affect the overall heat transfer and flow rate from the reservoir. Details of such variations are documented and discussed in the present paper.

Computational fluid dynamics simulation and turbomachinery code validation of a high pressure ratio radial-inflow turbine

The present study explores reproducing the closest geometry of a high pressure ratio single stage radial-inflow turbine applied in the Sundstrans Power Systems T-100 Multipurpose Small Power Unit. The commercial software ANSYS-Vista RTD along with a built in module, BladeGen, is used to conduct a meanline design and create 3D geometry of one flow passage. Carefully examining the proposed design against the geometrical and experimental data, ANSYS-TurboGrid is applied to generate computational mesh. CFD simulations are performed with ANSYS-CFX in which three-dimensional Reynolds-Averaged Navier-Stokes equations are solved subject to appropriate boundary conditions. Results are compared with numerical and experimental data published in the literature in order to generate the exact geometry of the existing turbine and validate the numerical results against the experimental ones.

Solar evaporator-collectors : analyses and applications

Taking into consideration of growing energy needs and concern for environmental degradation, clean and inexhaustible energy source, such as solar energy, is receiving greater attention for various applications. The use of solar energy system reduces pollution, waste and has little or no harmful effects on the environment. It is appreciated that this source of energy can be complementary rather than being competitive to conventional energy sources. In order to collect and harness energy from the sun, a solar collector is essential. A solar collector is basically a heat exchanger that transforms solar radiant energy into heat or thermal energy. Improvement of performance is essential for commercial acceptance of their use in such applications. Many studies have been undertaken on the enhancement of thermal performance of solar collectors, using diverse materials of various shapes, dimensions and layouts. In the literature, various collector designs have been proposed and tested with the objective of meeting these requirements [1-8]. Omer et al. [1] found the efficiency of a solar collector of about 70% in a solar assisted heat pump system. Traditional solar collectors are single phase collectors, in which the working fluid is either air or water. Different modifications are suggested and applied to improve the heat transfer between the absorber and working fluid in a collector. These modifications include the use of absorber with fins attached [2,3], corrugated absorber [4,5], matrix type absorber [6], V-groove solar air collector [7]. Karim et al. [8] approached a review of design and construction of three types (flat, vee-grooved, and finned) of air collectors. Two-phase collectors, on the other hand, have significant potential for continuous operation round the clock, when used in conjunction with a compressor, as found in a solar assisted heat-pump cycle.

Particle size distribution effects on preferential deposition areas in metal foam wrapped tube bundle

This paper presents a numerical model for understanding particle transport and deposition in metal foam heat exchangers. Two-dimensional steady and unsteady numerical simulations of a standard single row metal foam-wrapped tube bundle are performed for different particle size distributions, i.e. uniform and normal distributions. Effects of different particle sizes and fluid inlet velocities on the overall particle transport inside and outside the foam layer are also investigated. It was noted that the simplification made in the previously-published numerical works in the literature, e.g. uniform particle deposition in the foam, is not necessarily accurate at least for the cases considered here. The results highlight the preferential particle deposition areas both along the tube walls and inside the foam using a developed particle deposition likelihood matrix. This likelihood matrix is developed based on three criteria being particle local velocity, time spent in the foam, and volume fraction. It was noted that the particles tend to deposit near both front and rear stagnation points. The former is explained by the higher momentum and direct exposure of the particles to the foam while the latter only accommodate small particles which can be entrained in the recirculation region formed behind the foam-wrapped tubes.

Onset of transition in mixed convection of a lid-driven trapezoidal enclosure filled with water-Al2O3 nanofluid

A numerical study is carried out to investigate the transition from laminar to chaos in mixed convection heat transfer inside a lid-driven trapezoidal enclosure. In this study, the top wall is considered as isothermal cold surface, which is moving in its own plane at a constant speed, and a constant high temperature is provided at the bottom surface. The enclosure is assumed to be filled with water-Al2O3 nanofluid. The governing Navier–Stokes and thermal energy equations are expressed in non-dimensional forms and are solved using Galerkin finite element method. Attention is paid in the present study on the pure mixed convection regime at Richandson number, Ri = 1. The numerical simulations are carried out over a wide range of Reynolds (0.1 ≤ Re ≤ 103) and Grashof (0.01 ≤ Gr ≤ 106) numbers. Effects of the presence of nanofluid on the characteristics of mixed convection heat transfer are also explored. The average Nusselt numbers of the heated wall are computed to demonstrate the influence of flow parameter variations on heat transfer. The corresponding change of flow and thermal fields is visualized from the streamline and the isotherm contour plots.

Energetic and thermal properties of tilt grain boundaries in graphene/hexagonal boron nitride heterostructures

Graphene/hexagonal boron nitride (G/h-BN) heterostructure has attracted tremendous research efforts owing to its great potential for applications in nano-scale electronic devices. In such hybrid materials, tilt grain boundaries (GBs) between graphene and h-BN grains may have unique physical properties, which have not been well understood. Here we have conducted non-equilibrium molecular dynamics simulations to study the energetic and thermal properties of tilt GBs in G/h-BN heterostructures. The effect of misorientation angles of tilt GBs on both GB energy and interfacial thermal conductance are investigated.

Enhanced interfacial thermal transport across graphene–polymer interfaces by grafting polymer chains

Thermal transport in graphene-polymer nanocomposite is complicated and has not been well understood. The interfacial thermal transport between graphene nanofiller and polymer matrix is expected to play a key role in controlling the overall thermal performance of graphene-polymer nanocomposite. In this work, we investigated the thermal transport across graphene-polymer interfaces functionalized with end-grafted polymer chains using molecular dynamics simulations. The effects of grafting density, chain length and initial morphology on the interfacial thermal transport were systematically investigated. It was found that end-grafted polymer chains could significantly enhance interfacial thermal transport and the underlying mechanism was considered to be the enhanced vibration coupling between graphene and polymer. In addition, a theoretical model based on effective medium theory was established to predict the thermal conductivity in graphene-polymer nanocomposites.

Double diffusive Marangoni convection flow of electrically conducting fluid in a square cavity with chemical reaction

Double diffusive Marangoni convection flow of viscous incompressible electrically conducting fluid in a square cavity is studied in this paper by taking into consideration of the effect of applied magnetic field in arbitrary direction and the chemical reaction. The governing equations are solved numerically by using alternate direct implicit (ADI) method together with the successive over relaxation (SOR) technique. The flow pattern with the effect of governing parameters, namely the buoyancy ratio W, diffusocapillary ratio w, and the Hartmann number Ha, is investigated. It is revealed from the numerical simulations that the average Nusselt number decreases; whereas the average Sherwood number increases as the orientation of magnetic field is shifted from horizontal to vertical. Moreover, the effect of buoyancy due to species concentration on the flow is stronger than the one due to thermal buoyancy. The increase in diffusocapillary parameter, w caus