Scaling for the Prandtl number of the natural convection boundary layer of an inclined flat plate under uniform surface heat flux

An improved scaling analysis and direct numerical simulations are performed for the unsteady natural convection boundary layer adjacent to a downward facing inclined plate with uniform heat flux. The development of the thermal or viscous boundary layers may be classified into three distinct stages: a start-up stage, a transitional stage and a steady stage, which can be clearly identified in the analytical as well as the numerical results. Previous scaling shows that the existing scaling laws of the boundary layer thickness, velocity and steady state time scale for the natural convection flow on a heated plate of uniform heat flux provide a very poor prediction of the Prandtl number dependency of the flow. However, those scalings perform very well with Rayleigh number and aspect ratio dependency. In this study, a modified Prandtl number scaling is developed using a triple layer integral approach for Pr > 1. It is seen that in comparison to the direct numerical simulations, the modified scaling performs considerably better than the previous scaling.

Modified beam theories for bending properties of nanowires considering surface/intrinsic effects and axial extension effect

Several studies of the surface effect on bending properties of a nanowire (NW) have been conducted. However, these analyses are mainly based on theoretical predictions, and there is seldom integration study in combination between theoretical predictions and simulation results. Thus, based on the molecular dynamics (MD) simulation and different modified beam theories, a comprehensive theoretical and numerical study for bending properties of nanowires considering surface/intrinsic stress effects and axial extension effect is conducted in this work. The discussion begins from the Euler-Bernoulli beam theory and Timoshenko beam theory augmented with surface effect. It is found that when the NW possesses a relatively small cross-sectional size, these two theories cannot accurately interpret the true surface effect. The incorporation of axial extension effect into Euler-Bernoulli beam theory provides a nonlinear solution that agrees with the nonlinear-elastic experimental and MD results. However, it is still found inaccurate when the NW cross-sectional size is relatively small. Such inaccuracy is also observed for the Euler-Bernoulli beam theory augmented with both contributions from surface effect and axial extension effect. A comprehensive model for completely considering influences from surface stress, intrinsic stress, and axial extension is then proposed, which leads to good agreement with MD simulation results. It is thus concluded that, for NWs with a relatively small cross-sectional size, a simple consideration of surface stress effect is inappropriate, and a comprehensive consideration of the intrinsic stress effect is required.

Natural convection in a triangular enclosure due to non-uniform cooling on top

Natural convection in a triangular enclosure subject to non-uniformly cooling at the inclined surfaces and uniformly heating at the base is investigated numerically. The numerical simulations of the unsteady flows over a range of Rayleigh numbers and aspect ratios are carried out using Finite Volume Method. Since the upper surface is cooled and the bottom surface is heated, the air flow in the enclosure is potentially unstable to Rayleigh Benard instability. It is revealed that the transient flow development in the enclosure can be classified into three distinct stages; an early stage, a transitional stage and a steady stage. It is also found that the flow inside the enclosure strongly depends on the governing parameters; Rayleigh number and aspect ratio. The asymmetric behaviour of the flow about the geometric centre line is discussed in detailed. The heat transfer through the roof and the ceiling as a form of Nusselt number is also reported in this study.

Free convection in a triangular enclosure with fluid-saturated porous medium and internal heat generation

Unsteady natural convection inside a triangular cavity has been studied in this study. The cavity is filled with a saturated porous medium with non-isothermal left inclined wall while the bottom surface is isothermally heated and the right inclined surface is isothermally cooled. An internal heat generation is also considered which is dependent on the fluid temperature. The governing equations are solved numerically by finite volume method. The Prandtl number, Pr of the fluid is considered as 0.7 (air) while the aspect ratio and the Rayleigh number, Ra are considered as 0.5 and 105 respectively. The effect of heat generation on the fluid flow and heat transfer have been presented as a form of streamlines and isotherms. The rate of heat transfer through three surfaces of the enclosure is also presented.

Double diffusive magneto-convection fluid flow in a strong cross magnetic field with uniform surface heat and mass flux

In this study, magnetohydrodynamic natural convection boundary layer flow of an electrically conducting and viscous incompressible fluid along a heated vertical flat plate with uniform heat and mass flux in the presence of strong cross magnetic field has been investigated. For smooth integrations the boundary layer equations are transformed in to a convenient dimensionless form by using stream function formulation as well as the free variable formulation. The nonsimilar parabolic partial differential equations are integrated numerically for Pr ≪1 that is appropriate for liquid metals against the local Hartmann parameter ξ . Further, asymptotic solutions are obtained near the leading edge using regular perturbation method for smaller values of ξ . Solutions for values of ξ ≫ 1 are also obtained by employing the matched asymptotic technique. The results obtained for small, large and all ξ regimes are examined in terms of shear stress, τw, rate of heat transfer, qw, and rate of mass transfer, mw, for important physical parameter. Attention has been given to the influence of Schmidt number, Sc, buoyancy ratio parameter, N and local Hartmann parameter, ξ on velocity, temperature and concentration distributions and noted that velocity and temperature of the fluid achieve their asymptotic profiles for Sc ≥ 10:0.

Mixed convection flow of micropolar fluid in an open ended arc-shape cavity

Numerical simulations for mixed convection of micropolar fluid in an open ended arc-shape cavity have been carried out in this study. Computation is performed using the Alternate Direct Implicit (ADI) method together with the Successive Over Relaxation (SOR) technique for the solution of governing partial differential equations. The flow phenomenon is examined for a range of values of Rayleigh number, 102 ≤ Ra ≤ 106, Prandtl number, 7 ≤ Pr ≤ 50, and Reynolds number, 10 ≤ Re ≤ 100. The study is mainly focused on how the micropolar fluid parameters affect the fluid properties in the flow domain. It was found that despite the reduction of flow in the core region, the heat transfer rate increases, whereas the skin friction and microrotation decrease with the increase in the vortex viscosity parameter, Δ.

A fundamental numerical and theoretical study for the vibrational properties of nanowires

Based on the molecular dynamics (MD) simulation and the classical Euler-Bernoulli beam theory, a fundamental study of the vibrational performance of the Ag nanowire (NW) is carried out. A comprehensive analysis of the quality (Q)-factor, natural frequency, beat vibration, as well as high vibration mode is presented. Two excitation approaches, i.e., velocity excitation and displacement excitation, have been successfully implemented to achieve the vibration of NWs. Upon these two kinds of excitations, consistent results are obtained, i.e., the increase of the initial excitation amplitude will lead to a decrease to the Q-factor, and moderate plastic deformation could increase the first natural frequency. Meanwhile, the beat vibration driven by a single relatively large excitation or two uniform excitations in both two lateral directions is observed. It is concluded that the nonlinear changing trend of external energy magnitude does not necessarily mean a nonconstant Q-factor. In particular, the first order natural frequency of the Ag NW is observed to decrease with the increase of temperature. Furthermore, comparing with the predictions by Euler- Bernoulli beam theory, the MD simulation provides a larger and smaller first vibration frequencies for the clamped-clamped and clamped-free thin Ag NWs, respectively. Additionally, for thin NWs, the first order natural frequency exhibits a parabolic relationship with the excitation magnitudes. The frequencies of the higher vibration modes tend to be low in comparison to Euler-Bernoulli beam theory predictions. A combined initial excitation is proposed which is capable to drive the NW under a multi-mode vibration and arrows the coexistence of all the following low vibration modes. This work sheds lights on the better understanding of the mechanical properties of NWs and benefits the increasing utilities of NWs in diverse nano-electronic devices.

Anomalous subdiffusion equations by the meshless collocation method

This paper aims to develop an implicit meshless collocation technique based on the moving least squares approximation for numerical simulation of the anomalous subdiffusion equation(ASDE). The discrete system of equations is obtained by using the MLS meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach related to the time discretization are investigated theoretically and numerically. The numerical examples with regular and irregular nodal distributions are used to the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling of ASDEs.

Prandtl number scaling of the unsteady natural convection boundary layer adjacent to a vertical flat plate for Pr > 1 subject to ramp surface heat flux

It is found in the literature that the existing scaling results for the boundary layer thickness, velocity and steady state time for the natural convection flow over an evenly heated plate provide a very poor prediction of the Prandtl number dependency of the flow. However, those scalings provide a good prediction of two other governing parameters’ dependency, the Rayleigh number and the aspect ratio. Therefore, an improved scaling analysis using a triple-layer integral approach and direct numerical simulations have been performed for the natural convection boundary layer along a semi-infinite flat plate with uniform surface heat flux. This heat flux is a ramp function of time, where the temperature gradient on the surface increases with time up to some specific time and then remains constant. The growth of the boundary layer strongly depends on the ramp time. If the ramp time is sufficiently long, the boundary layer reaches a quasi steady mode before the growth of the temperature gradient is completed. In this mode, the thermal boundary layer at first grows in thickness and then contracts with increasing time. However, if the ramp time is sufficiently short, the boundary layer develops differently, but after the wall temperature gradient growth is completed, the boundary layer develops as though the startup had been instantaneous.

A novel quadratic Edge-based Smoothed Conforming Point Interpolation Method (ES-CPIM) for elasticity problems

his paper formulates an edge-based smoothed conforming point interpolation method (ES-CPIM) for solid mechanics using the triangular background cells. In the ES-CPIM, a technique for obtaining conforming PIM shape functions (CPIM) is used to create a continuous and piecewise quadratic displacement field over the whole problem domain. The smoothed strain field is then obtained through smoothing operation over each smoothing domain associated with edges of the triangular background cells. The generalized smoothed Galerkin weak form is then used to create the discretized system equations. Numerical studies have demonstrated that the ES-CPIM possesses the following good properties: (1) ES-CPIM creates conforming quadratic PIM shape functions, and can always pass the standard patch test; (2) ES-CPIM produces a quadratic displacement field without introducing any additional degrees of freedom; (3) The results of ES-CPIM are generally of very high accuracy.

Prandtl number scaling of natural convection of the flow on an inclined at plate due to uniform surface heat flux

A new scaling analysis has been performed for the unsteady natural convection boundary layer under a downward facing inclined plate with uniform heat flux. The development of the thermal or viscous boundary layers may be classified into three distinct stages including an early stage, a transitional stage and a steady stage, which can be clearly identified in the analytical as well as numerical results. Earlier scaling shows that the existing scaling laws of the boundary layer thickness, velocity and steady state time scales for the natural convection flow on a heated plate of uniform heat flux provide a very poor prediction of the Prandtl number dependency. However, those scalings performed very well with Rayleigh number and aspect ratio dependency. In this study, a modifed Prandtl number scaling has been developed using a triple-layer integral approach for Pr > 1. It is seen that in comparison to the direct numerical simulations, the new scaling performs considerably better than the previous scaling.

Unsteady natural convection within a differentially heated enclosure of sinusoidal corrugated side walls

Numerically investigation of natural convection within a differentially heated modified square enclosure with sinusoidally corrugated side walls has been performed for different values of Rayleigh number. The fluid inside the enclosure considered is air and is quiescent, initially. The top and bottom surfaces are flat and considered as adiabatic. Results reveal three main stages: an initial stage, a transitory or oscillatory stage and a steady stage for the development of natural convection flow inside the corrugated cavity. The numerical scheme is based on the finite element method adapted to triangular non-uniform mesh element by a non-linear parametric solution algorithm. Investigation has been performed for the Rayleigh number, Ra ranging from 105 to 108 with variation of corrugation amplitude and frequency. Constant physical properties for the fluid medium have been assumed. Results have been presented in terms of the isotherms, streamlines, temperature plots, average Nusselt numbers, traveling waves and thermal boundary layer thickness plots, temperature and velocity profiles. The effects of sudden differential heating and its consequent transient behavior on fluid flow and heat transfer characteristics have been observed for the range of governing parameters. The present results show that the transient phenomena are greatly influenced by the variation of the Rayleigh Number with corrugation amplitude and frequency.

Orthopedic bone plates : evolution in structure implementation technique and biomaterial

With many important developments over the last century, nowadays orthopedic bone plate now excels over other types of internal fixators in bone fracture fixation. The developments involve the design, material and implementation techniques of the plates. This paper aims to review the evolution in implementation technique and biomaterial of the orthopedic bone plates. Plates were initially used to fix the underlying bones firmly. Accordingly, Compression plate (CP), Dynamic compression plate (DCP), Limited contact dynamic compression plate (LC-DCP) and Point contact fixator (PC-Fix) were developed. Later, the implementation approach was changed to locking, and the Less Invasive Stabilization System (LISS) plate was introduced as a result. Finally, a combination of both of these approaches has been used by introducing the Locking Compression Plate (LCP). Currently, precontoured LCPs are mainly used for bone fracture fixation. In parallel with structure and implementation techniques, numerous advances have occurred in biomaterials of the plates. Titanium and stainless steel alloys are now the most common biomaterials in production of orthopedic bone plates. However, regarding the biocompatibility, bioactivity and biodegradability characteristics of Mg alloys, Ta alloys, SMAs, carbon fiber composites and bioceramics, these materials are considered as potentially suitable for plates. However, due to poor mechanical properties, they have very limited applications. Therefore, further studies are required in future to solve these problems and make them feasible for heavy-duty bone plates.

On the oscillatory behavior of transient Rayleigh Benard convection of air for 2D channel flow at moderate Rayleigh numbers

Unsteady numerical simulation of Rayleigh Benard convection heat transfer from a 2D channel is performed. The oscillatory behavior is attributed to recirculation of ascending and descending flows towards the core of the channel producing organized rolled motions. Variation of the parameters such as Reynolds number, channel outlet flow area and inclination of the channel are considered. Increasing Reynolds number (for a fixed Rayleigh number), delays the generation of vortices. The reduction in the outflow area leads to the later and the less vortex generation. As the time progresses, more vortices are generated, but the reinforced mean velocity does not let the eddies to enter the core of the channel. Therefore, they attach to the wall and reduce the heat transfer area. The inclination of the channel (both positive and negative) induces the generated vortices to get closer to each other and make an enlarged vortex.

Atomistic investigations of single-crystal silicon with pre-existing defects

Molecular dynamics (MD) simulations have been employed to investigate the single-crystal Si properties with different pre-existing cavities under nanoindentation. Cavities with different radii and positions have been considered. It is found that pre-existing cavities in the Si substrate would obviously influence the mechanical properties of Si under nanoindentation. Furthermore, pre-existing cavities would absorb part of the strain energy during loading and then release during unloading. It would decrease plastic deformation to the substrate. Particularly, the larger of the cavity or the nearer of the cavity to the substrate’s top surface, the larger decrease of Young’s modulus and hardness is usually observed. Just as expected, the larger offset of the cavity in the lateral direction, the less influence is usually seen.