Modeling of the non-isothermal liquid droplet impact on a heated solid substrate with heterogeneous wettability

A comprehensive numerical investigation on the impingement and spreading of a non-isothermal liquid droplet on a solid substrate with heterogeneous wettability is presented in this work. The time-dependent incompressible Navier-Stokes equations are used to describe the fluid flow in the liquid droplet, whereas the heat transfer in the moving droplet and in the solid substrate is described by the energy equation. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the time-dependent incompressible Navier-Stokes equation and the energy equation in the time-dependent moving domain. Moreover, the Marangoni convection is included in the variational form of the Navier-Stokes equations without calculating the partial derivatives of the temperature on the free surface. The heterogeneous wettability is incorporated into the numerical model by defining a space-dependent contact angle. An array of simulations for droplet impingement on a heated solid substrate with circular patterned heterogeneous wettability are presented. The numerical study includes the influence of wettability contrast, pattern diameter, Reynolds number and Weber number on the confinement of the spreading droplet within the inner region, which is more wettable than the outer region. Also, the influence of these parameters on the total heat transfer from the solid substrate to the liquid droplet is examined. We observe that the equilibrium position depends on the wettability contrast and the diameter of the inner surface. Consequently. the heat transfer is more when the wettability contrast is small and/or the diameter of inner region is large. The influence of the Weber number on the total heat transfer is more compared to the Reynolds number, and the total heat transfer increases when the Weber number increases. (C) 2015 Elsevier Ltd. All rights reserved.

Isolated singularities of polyharmonic operator in even dimension

We consider the equation Delta(2)u = g(x, u) >= 0 in the sense of distribution in Omega' = Omega\textbackslash {0} where u and -Delta u >= 0. Then it is known that u solves Delta(2)u = g(x, u) + alpha delta(0) - beta Delta delta(0), for some nonnegative constants alpha and beta. In this paper, we study the existence of singular solutions to Delta(2)u = a(x) f (u) + alpha delta(0) - beta Delta delta(0) in a domain Omega subset of R-4, a is a nonnegative measurable function in some Lebesgue space. If Delta(2)u = a(x) f (u) in Omega', then we find the growth of the nonlinearity f that determines alpha and beta to be 0. In case when alpha = beta = 0, we will establish regularity results when f (t) <= Ce-gamma t, for some C, gamma > 0. This paper extends the work of Soranzo (1997) where the author finds the barrier function in higher dimensions (N >= 5) with a specific weight function a(x) = |x|(sigma). Later, we discuss its analogous generalization for the polyharmonic operator.

Phase boundary at x=0.03 and its anomalous influence on the structure and properties in the lead-free piezoelectric (1-x)Na1/2Bi1/2TiO3-(x)BaTiO3

The complexity associated with local structures continues to pose challenges with regard to the understanding of the structure-property relationship in Na1/2Bi1/2TiO3-based lead-free piezoceramics. (1-x)Na1/2Bi1/2TiO3-(x)BaTiO3 is an extensively studied system because of its interesting piezoelectric properties. Recently, a room temperature phase boundary was reported at x = 0.03 in this system Ma et al., Adv. Funct. Mater. 23, 5261 (2013)]. In the present work we have examined this subtle phase boundary using x-ray diffraction, neutron diffraction, dielectric measurements as a function of composition (x < 0.06), temperature, and electric field. Our results show that this boundary separates an R3c + Cc-like structural state for x < 0.03 from an R3c+ cubiclike structural state for 0.03 <= x <= 0.05 in the unpoled specimens. This phase boundary is characterized by an anomalous reduction in the depolarization temperature, and a suppression of the tetragonal distortion of the high temperature P4bm phase. Our results also provide the clue to understand the pathway leading to the cubiclike structure of the critical composition x = 0.06, known for its highest piezoelectric response.

On the Equivalence of Acoustic Impedance and Squeeze Film Impedance in Micromechanical Resonators

In this work, we address the issue of modeling squeeze film damping in nontrivial geometries that are not amenable to analytical solutions. The design and analysis of microelectromechanical systems (MEMS) resonators, especially those that use platelike two-dimensional structures, require structural dynamic response over the entire range of frequencies of interest. This response calculation typically involves the analysis of squeeze film effects and acoustic radiation losses. The acoustic analysis of vibrating plates is a very well understood problem that is routinely carried out using the equivalent electrical circuits that employ lumped parameters (LP) for acoustic impedance. Here, we present a method to use the same circuit with the same elements to account for the squeeze film effects as well by establishing an equivalence between the parameters of the two domains through a rescaled equivalent relationship between the acoustic impedance and the squeeze film impedance. Our analysis is based on a simple observation that the squeeze film impedance rescaled by a factor of jx, where x is the frequency of oscillation, qualitatively mimics the acoustic impedance over a large frequency range. We present a method to curvefit the numerically simulated stiffness and damping coefficients which are obtained using finite element analysis (FEA) analysis. A significant advantage of the proposed method is that it is applicable to any trivial/nontrivial geometry. It requires very limited finite element method (FEM) runs within the frequency range of interest, hence reducing the computational cost, yet modeling the behavior in the entire range accurately. We demonstrate the method using one trivial and one nontrivial geometry.

On the Structural Stability of Melt Spun Ribbons of Fe95-x Zr (x) B4Cu1 (x=7 and 9) Alloys and Correlation with Their Magnetic Properties

Melt spun ribbons of Fe95-x Zr (x) B4Cu1 with x = 7 (Z7B4) and 9 (Z9B4) alloys have been prepared, and their structure and magnetic properties have been evaluated using XRD, DSC, TEM, VSM, and Mossbauer spectroscopy. The glass forming ability (GFA) of both alloys has been calculated theoretically using thermodynamical parameters, and Z9B4 alloy is found to possess higher GFA than that of Z7B4 alloy which is validated by XRD results. On annealing, the amorphous Z7B4 ribbon crystallizes into nanocrystalline alpha-Fe, whereas amorphous Z9B4 ribbon shows two-stage crystallization process, first partially to bcc solid solution which is then transformed to nanocrystalline alpha-Fe and Fe2Zr phases exhibiting bimodal distribution. A detailed phase analysis using Mossbauer spectroscopy through hyperfine field distribution of phases has been carried out to understand the crystallization behavior of Z7B4 and Z9B4 alloy ribbons. In order to understand the phase transformation behavior of Z7B4 and Z9B4 ribbons, molar Gibbs free energies of amorphous, alpha-Fe, and Fe2Zr phases have been evaluated. It is found that in case of Z7B4, alpha-Fe is always a stable phase, whereas Fe2Zr is stable at higher temperature for Z9B4. (C) The Minerals, Metals & Materials Society and ASM International 2015

Macropores generated during shrinkage in two paddy soils using X-ray micro-computed tomography

Soil shrinkage curve represents a decrease of total porosity or an increase of bulk density with water loss. However, our knowledge of the dynamics of pores and their geometry during soil shrinkage is scarce, partially due to lack of reliable methods for determining soil pores in relation to change in soil water. This study aimed to investigate the dynamics of macropores (>30 mu m) of paddy soils during shrinkage. Two, paddy soils, which were sampled from one paddy field cultivated for 20 years (YPF) and the other one for over 100 years (OPF), represented difference in crack geometry in the field. Macropore parameters (volume, connectivity, and orientation of pores) and soil shrinkage parameters were determined on the same undisturbed soil cores by X-ray microtomography and shrinkage curve, respectively. Macroporosity was on average four times larger in the YPF than in the OPF whereas the shrinkage capacity was lower in the YPF as compared to the OPF (0.09 vs. 0.15 COLE). Soil shrinkage increased the volume of pores by 3.7% in the YPF and by 1.6% in the OPF as well as their connectivity. The formation of macropores occurred mostly in the proportional shrinkage phase. As a result, the slope of the proportional shrinkage phase was smaller in the YPF (0.65) than in the OPF (0.89). New macropores were cracks and extended pre-existing pores in the range of 225-1215 pm size without any preferential orientation. This work provides image evidences that in paddy soils with high shrinkage capacity more macropores are generated in the soil presenting a smaller proportional shrinkage slope. (C) 2015 Elsevier B.V. All rights reserved.

Joint modeling of lithosphere and mantle dynamics: Evaluation of constraints from global tomography models

With the advances in technology, seismological theory, and data acquisition, a number of high-resolution seismic tomography models have been published. However, discrepancies between tomography models often arise from different theoretical treatments of seismic wave propagation, different inversion strategies, and different data sets. Using a fixed velocity-to-density scaling and a fixed radial viscosity profile, we compute global mantle flow models associated with the different tomography models and test the impact of these for explaining surface geophysical observations (geoid, dynamic topography, stress, and strain rates). We use the joint modeling of lithosphere and mantle dynamics approach of Ghosh and Holt (2012) to compute the full lithosphere stresses, except that we use HC for the mantle circulation model, which accounts for the primary flow-coupling features associated with density-driven mantle flow. Our results show that the seismic tomography models of S40RTS and SAW642AN provide a better match with surface observables on a global scale than other models tested. Both of these tomography models have important similarities, including upwellings located in Pacific, Eastern Africa, Iceland, and mid-ocean ridges in the Atlantic and Indian Ocean and downwelling flows mainly located beneath the Andes, the Middle East, and central and Southeast Asia.

Reconciling results of MOCVD of a CNT composite with equilibrium thermodynamics

Composition and microstructure of the composite films can be tailored by controlling the CVD process parameters if an appropriate model can be suggested for quantitative prediction of growth. This is possible by applying equilibrium thermodynamics. A modification of such standard modeling procedure was required to account for the deposition of a hybrid film comprised of carbon nanotubes (CNTs), metallic iron (Fe), and magnetite (Fe3O4), a composite useful for energy storage. In contrast with such composite nature of the deposits obtained by inert-ambient CVD using Fe(acac)3 as precursor, equilibrium thermodynamic modeling with standard procedure predicts the deposition of only Fe3C and carbon, without any co-deposition of Fe and Fe3O4. A modification of the procedure comprising chemical reasoning is therefore proposed herein, which predicts simultaneous deposition of FeO1-x, Fe3C, Fe3O4 and C. At high temperatures and in a carbon-rich atmosphere, these convert to Fe3O4, Fe and C, in agreement with experimental CVD. Close quantitative agreement between the modified thermodynamic modeling and experiment validates the reliability of the modified procedure. Understanding of the chemical process through thermodynamic modeling provides potential for control of CVD process parameters to achieve desired hybrid growth. (C) 2016 Elsevier B.V. All rights reserved.

COUPLED OPTICAL-THERMAL-FLUID MODELING OF A DIRECTLY HEATED TUBULAR SOLAR RECEIVER FOR SUPERCRITICAL CO2 BRAYTON CYCLE

Recent studies have evaluated closed-loop supercritical carbon dioxide (s-CO2) Brayton cycles to be a higher energy density system in comparison to conventional superheated steam Rankine systems. At turbine inlet conditions of 923K and 25 MPa, high thermal efficiency (similar to 50%) can be achieved. Achieving these high efficiencies will make concentrating solar power (CSP) technologies a competitive alternative to current power generation methods. To incorporate a s-CO2 Brayton power cycle in a solar power tower system, the development of a solar receiver capable of providing an outlet temperature of 923 K (at 25 MPa) is necessary. The s-CO2 will need to increase in temperature by similar to 200 K as it passes through the solar receiver to satisfy the temperature requirements of a s-CO2 Brayton cycle with recuperation and recompression. In this study, an optical-thermal-fluid model was developed to design and evaluate a tubular receiver that will receive a heat input similar to 2 MWth from a heliostat field. The ray-tracing tool SolTrace was used to obtain the heat-flux distribution on the surfaces of the receiver. Computational fluid dynamics (CFD) modeling using the Discrete Ordinates (DO) radiation model was used to predict the temperature distribution and the resulting receiver efficiency. The effect of flow parameters, receiver geometry and radiation absorption by s-CO2 were studied. The receiver surface temperatures were found to be within the safe operational limit while exhibiting a receiver efficiency of similar to 85%.