High temperature thermoelectric properties of Zr and Hf based transition metal dichalcogenides: A first principles study

We investigate the electronic and thermal transport properties of bulk MX2 compounds (M = Zr, Hf and X = S, Se) by first-principles calculations and semi-classical Boltzmann transport theory. The band structure shows the confinement of heavy and light bands along the out of plane and in-plane directions, respectively. This results in high electrical conductivity (sigma) and large thermopower leading to a high power factor (S-2 sigma) for moderate n-type doping. The phonon dispersion demonstrates low frequency flat acoustical modes, which results in low group velocities (v(g)). Consequently, lowering the lattice thermal conductivity (kappa(latt)) below 2 W/m K. Low kappa(latt) combined with high power factor results in ZT > 0.8 for all the bulk MX2 compounds at high temperature of 1200 K. In particular, the ZT(max) of HfSe2 exceeds 1 at 1400 K. Our results show that Hf/Zr based dichalcogenides are very promising for high temperature thermoelectric application. (C) 2015 AIP Publishing LLC.

Assessment of driver vision functions in relation to their crash involvement in India

Among the human factors that influence safe driving, visual skills of the driver can be considered fundamental. This study mainly focuses on investigating the effect of visual functions of drivers in India on their road crash involvement. Experiments were conducted to assess vision functions of Indian licensed drivers belonging to various organizations, age groups and driving experience. The test results were further related to the crash involvement histories of drivers through statistical tools. A generalized linear model was developed to ascertain the influence of these traits on propensity of crash involvement. Among the sampled drivers, colour vision, vertical field of vision, depth perception, contrast sensitivity, acuity and phoria were found to influence their crash involvement rates. In India, there are no efficient standards and testing methods to assess the visual capabilities of drivers during their licensing process and this study highlights the need for the same.

Increasing Critical Sensitivity of the Load/Unload Response Ratio before Large Earthquakes with Identified Stress Acclumulation Pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.

Study of Gas Combustion in a Vented Cylindrical Vessel

Complicated interaction of a flame front with a turbulent flow induced by venting is studied during combustion of the stoichiometric propane/air mixture in a relatively large vented cylindrical vessel. Flame position, its shape, and combustion pressure were measured as a function of time and vent parameters. The experimental data were used to verify numerical simulation of the combustion process. The proposed numerical model satisfactorily simulates the main features of combustion in a closed and vented vessel such as flame configuration, flow and temperature fields, and pressure variation pattern. Simulated velocity and temperature distribution are very useful pieces of information because they are not available from experiments.

Modeling Study on the Characteristics of Laminar and Turbulent Argon Plasma Jets Impinging Normally upon a Flat Plate in Ambient Air

Modeling study is performed to compare the flow and heat transfer characteristics of laminar and turbulent argon thermal-plasma jets impinging normally upon a flat plate in ambient air. The combined-diffusion-coefficient method and the turbulence-enhanced combined-diffusion-coefficient method are employed to treat the diffusion of argon in the argon-air mixture for the laminar and the turbulent cases, respectively. Modeling results presented include the flow, temperature and argon concentration fields, the air mass flow-rates entrained into the impinging plasma jets, and the distributions of the heat flux density on the plate surface. It is found that the formation of a radial wall jet on the plate surface appreciably enhances the mass flow rate of the ambient air entrained into the laminar or turbulent plasma impinging-jet. When the plate standoff distance is comparatively small, there exists a significant difference between the laminar and turbulent plasma impinging-jets in their flow fields due to the occurrence of a large closed recirculation vortex in the turbulent plasma impinging-jet, and no appreciable difference is found between the two types of jets in their maximum values and distributions of the heat flux density at the plate surface. At larger plate standoff distances, the effect of the plate on the jet flow fields only appears in the region near the plate, and the axial decaying-rates of the plasma temperature, axial velocity and argon mass fraction along the axis of the laminar plasma impinging-jet become appreciably less than their turbulent counterparts.

Laboratory Study on Sediment Diffusion and Deposition into Blind Channels

Any waterway with one end closed and the other open is generally called a blind channel. The main flow tends to expand, separate, and cause circulation at the mouth of blind channels. The main flow continuously transfers momentum and sediment into the circulation region through the turbulent mixing region (TMR) between them, thus leading to a large amount of sediment deposition in the blind channels. This paper experimentally investigated the properties of the water flow and sediment diffusion in TMR, demonstrating that both water flow and sediment motion in TMR approximately coincide with a similar structure as in the free mixing layer induced by a jet. The similarity functions of flow velocity and sediment concentration are then assumed, based on observation, and the resulting calculation of these functions is substantially facilitated. For the kind of low velocity flow system of blind channels with a finite width, a simple formula for the sediment deposition rate in blind channels is established by analyzing the gradient of crosswise velocity and sediment concentration in TMR.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Pull-in Instability Study of Carbon Nanotube Tweezers Under the Influence of Van der Waals Forces

The pull-in instability of two nanotubes under van der Waals force is studied. The cantilever beam with large deformation model is used. The influence of nanotube parameters such as the interior radius, the gap distance between the two nanotubes, etc, on the pull-in instability is studied. The critical nanotube length is determined for each specific set of nanotube parameters. The Galerkin method is applied to discretize the governing equations, and it shows good convergence.

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

Experimental Study of the Earthquake Recurrence Period and the Trend of Post-Seismic Development

In order to study the earthquake recurrence and the characteristics of earthquake series, rupture tests of rock samples and plexiglass samples were made. On rock samples, a number of acoustic emission (AE) and strain measuring points were deployed; the load was one side direct shear. The variation characteristics of AE and strain at different detecting points around the extra large fracture were observed and studied. On plexiglass samples, a series of inclined cracks were prefabricated by a small-scale compressive testing machine. The samples were then loaded on a shockproof platen, when the samples were loaded, the stress intensity factor (SIF) was determined by the laser interferometric technique and shadow optical method of caustics. The fracture conditions such as material toughness around the extra large fracture were also studied. From those experimental results and the theory of fracture mechanics, the earthquake recurrence period and the trend of post-seismic development were studied.

Laboratory study on cracks in saturated sands

It has been reported([1]) that when a loosely packed column of saturated sand in a vertical cylindrical container is shock loaded axially by dropping to the floor, large horizontal cracks initiate, grow and eventually fade away in the sand as it settles under gravity. This paper shows that a similar phenomenon can also be observed when shock loading is replaced by forcing water to percolate upward through the sand column. It is believed that our result sheds further light on the physics of formation of these cracks.

Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES

Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Experimental study on two-phase gas-liquid flow patterns at normal and reduced gravity conditions

Experimental studies have been performed for horizontal two-phase air-water flows at normal and reduced gravity conditions in a square cross-section channel. The experiments at reduced gravity are conducted on board the Russian IL-76 reduced gravity airplane. Four flow patterns, namely bubble, slug, slug-annular transition and annular flows, are observed depending on the liquid and gas superficial velocities at both conditions. Semi-theoretical Weber number model is developed to include the shape influence on the slug-annular transition. It is shown that its prediction is in reasonable agreement with the experimental slug-annular transition under both conditions. For the case of two-phase gas-liquid flow with large value of the Froude number, the drift-flux model can predict well the observed boundary between bubble and slug flows.

A Numerical Study for Turbulent Flow and Thermal Influence Over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier-Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052-2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349-377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369-398)'s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243-258)'s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.