Energy density and release rate study of an elliptic-arc interface crack in piezoelectric solid under anti-plane shear

The anti-plane problem of an elliptical inhomogeneity with an interfacial crack in piezoelectric materials is investigated. The system is subjected to arbitrary singularity loads (point charge and anti-plane concentrated force) and remote anti-plane mechanical and in-plane electrical loads. Using the complex variable method, the explicit series form solutions for the complex potentials in the matrix and the inclusion regions are derived. The electroelastic field intensity factors, the corresponding energy release rates and the generalized strain energy density at the cracks tips are then provided. The influence of the aspect ratio of the ellipse, the crack geometry and the electromechanical coupling coefficient on the energy release rate and the strain energy density is discussed and shown in graphs. The results indicate that the energy release rate increases with increment of the aspect ratio of the ellipse and the influence of electromechanical coupling coefficient on the energy release rate is significant. The strain energy density decreases with increment of the aspect radio of the ellipse and it is always positive for the cases discussed. The energy release rate, however, can be negative when both mechanical and fields are applied.     

The effect of air-water interface on the vortex shedding from a vertical circular cylinder

The flow past an interface piercing circular cylinder at the Reynolds number Re=2.7×10⁴ and the Froude numbers Fr=0.2 and 0.8 is investigated using large-eddy simulation. A Lagrangian dynamic subgrid-scale model and a level set based sharp interface method are used for the spatially filtered turbulence closure and the air-water interface treatment, respectively. The mean interface elevation and the rms of interface fluctuations from the simulation are in excellent agreement with the available experimental data. The organized periodic vortex shedding observed in the deep flow is attenuated and replaced by small-scale vortices at the interface. The streamwise vorticity and the outward transverse velocity generated near the edge of the separated region, which enforces the separated shear layers to deviate from each other and restrains their interaction, are primarily responsible for the devitalization of the periodic vortex shedding at the interface. The lateral gradient of the difference between the vertical and transverse Reynolds normal stresses, increasing with the Froude number, is the main source of the streamwise vorticity and the outward transverse velocity at the interface.     

THE COMPACTION OF TIME-DEPENDENT VISCOUS GRANULAR MATERIALS CONSIDERING INERTIAL FORCES

In the present paper, compactions of time-dependent viscous granular materials are simulated step by step using the automatic adaptive mesh generation schemes. Inertial forces of the viscous incompressible aggregates are taken into account. The corresponding conservation equations, the weighted-integral formulations, and penalty finite element model are investigated. The fully discrete finite element equations for the simulation are derived. Polygonal particles of aggregates are simplified as mixed three-node and four-node elements. The automatic adaptive mesh generation schemes include contact detection algorithms, and mesh upgrade schemes. Solutions of the numerical simulation are in good agreement with some results from literatures. With minor modification, the proposed numerical model can be applied in several industries, including the pharmaceutical, ceramic, food, and household product manufacturing.     

THREE-DIMENSIONAL THERMOELASTIC ANALYSIS OF FUNCTIONALLY GRADED PLATE

Three-dimensional thermoelastic analysis is presented for an orthotropic functionally graded rectangular plate,which is simply supported and isothermal on its four lateral edges.With the assumption that material properties have arbitrary dependence on the thickness-coordinate,a Peano-Baker series solution is obtained for the thermoelastic fields of the functionally graded plate subjected to mechanical and thermal loads on its upper and lower surfaces by means of state space method.The correctness of the obtained series solution is validated through numerical examples.The influence of different material properties distributions on the structural response of the plate is also studied.     

平面滑动型岩质边坡稳定性极限分析上限法

在极限分析理论框架体系的基础上,提出了平面滑动型岩质边坡极限分析上限法。该方法按照外力所做的功率等于内力所消耗的功率,即滑体处于极限状态时两功率相等的条件——虚功率方程,综合考虑作用在岩质边坡上的后缘裂缝静水压力、沿滑面扬压力、重力、水平地震作用力、锚固力等外力,按照强度折减法,推导得出了岩质边坡稳定性评价的极限分析上限解。锦屏一级水电站右岸泄洪洞引渠内侧工程边坡属于典型的多级台阶状岩质高边坡,该边坡岩体结构为层状结构,变形破坏模式为平面滑动。对引渠内侧岩质边坡稳定性分析结果表明,极限分析上限法对平面滑动型岩质边坡有较好的适用性。     

Rolling/sliding of a particle on a flat wall in a linear shear flow at finite Re

Recently Lee and Balachandar proposed analytically-based expressions for drag and lift coefficients for a spherical particle moving on a flat wall in a linear shear flow at finite Reynolds number. In order to evaluate the accuracy of these expressions, we have conducted direct numerical simulations of a rolling particle for shear Reynolds number up to 100. We assume that the particle rolls on a horizontal flat wall with a small gap separating the particle from the wall (L = 0.505) and thus avoiding the logarithmic singularity. The influence of the shear Reynolds number and the translational velocity of the particle on the hydrodynamic forces of the particle was investigated under both transient and the final drag-free and torque-free steady state. It is observed that the quasi-steady drag and lift expressions of Lee and Balachandar provide good approximation for the terminal state of the particle motion ranging from perfect sliding to perfect rolling. With regards to transient particle motion in a wall-bounded shear flow it is observed that the above validated quasi-steady drag and lift forces must be supplemented with appropriate wall-corrected added-mass and history forces in order to accurately predict the time-dependent approach to the terminal steady state. Quantitative comparison with the actual particle motion computed in the numerical simulations shows that the theoretical models quite effective in predicting rolling/sliding motion of a particle in a wall-bounded shear flow at moderate Re.     

The characteristic‐based split (CBS) meshfree method for free surface flow problems in ALE formulation

A semi‐implicit characteristic‐based split (CBS) meshfree algorithm in the arbitrary Lagrangian Eulerian (ALE) framework is proposed for the numerical solution of incompressible free surface flow problem in the paper. The algorithm is the extension of general CBS method which was initially introduced in finite element framework, this is due to the fact that CBS method not only can enhance the stability, but also avoid LBB condition when equal order basis function is used to approximate velocity and pressure variables. Meanwhile, a simple way for node update and node speed calculation is developed which is used to capture the free surface exactly. The numerical solutions are compared with available analytical and numerical solutions, which shows that the proposed method has better ability to simulate the free surface incompressible flow problem. Copyright © 2009 John Wiley & Sons, Ltd.     

Mobile and immobile fluid transport: Coupling framework

We introduce a solver method for mobile and immobile transport regions. The motivation is driven by transport processes in porous media (e.g. waste disposal, chemical deposition processes). We analyze the coupled transport‐reaction equation with mobile and immobile areas. We apply analytical methods, such as Laplace‐transformation, and for the numerical methods we apply Godunov's scheme, see (Mat. Sb. 1959; 47 :271–306; Finite Volume Methods for Hyperbolic Problems. Cambridge University Press: Cambridge, 2002). The method is based numerically on flux‐based characteristic methods and is an attractive alternative to the classical higher‐order TVD methods, see (J. Comput. Phys. 1993; 49 :357–393). In this paper, we will focus on the derivation of analytical solutions for general and special solutions of the characteristic methods that are embedded in a finite‐volume method. At the end of the paper, we illustrate the higher‐order method for different benchmark problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Internal flow numerical simulation of a cross‐flow fan in refrigerant operating condition using user‐defined functions

In the paper, a cross‐flow fan in refrigerant operating condition is systematically simulated using user‐defined functions. Three‐dimensional simulations are acquired with Navier–Stokes equations coupled with k–ε turbulence model, and internal flow characteristics of an indoor split‐type air conditioner are obtained, which is mainly composed of cross‐flow fan and heat exchanger. It has systematically been simulated in the isothermal flow condition that the performance of cross‐flow fan may be reduced easily with dry or humid air, and in the refrigerant operating condition in which user‐defined functions are applied to the humid air, considered as a mixture of dry air and vapor. A density‐modulated function is adopted to deal with the condensation of the vapor at the heat‐transfer region approximately. The results show flow mechanism of the two gas‐phase flow, including phase‐vary process. The distribution of the parameters is not uniform at the inlet of the machine, the intensity and position of pressure and velocity vary along the axial direction of the fan, the distribution of vapor volume fraction and turbulent intensity in heat‐transfer region is obtained, and the external characteristic data of the indoor machine are obtained and analyzed. Compared with the experimental data, the calculated characteristic curves and designed parameters are on target. © British Crown Copyright 2010/MOD. Reproduced with permission. Published by John Wiley & Sons, Ltd.     

A local domain‐free discretization method for simulation of incompressible flows over moving bodies

This paper presents a local domain‐free discretization (DFD) method for the simulation of unsteady flows over moving bodies governed by the incompressible Navier–Stokes equations. The discretization strategy of DFD is that the discrete form of partial differential equations at an interior point may involve some points outside the solution domain. All the mesh points are classified as interior points, exterior dependent points and exterior independent points. The functional values at the exterior dependent points are updated at each time step by the approximate form of solution near the boundary. When the body is moving, only the status of points is changed and the mesh can stay fixed. The issue of ‘freshly cleared nodes/cells’ encountered in usual sharp interface methods does not pose any particular difficulty in the presented method. The Galerkin finite‐element approximation is used for spatial discretization, and the discrete equations are integrated in time via a dual‐time‐stepping scheme based on artificial compressibility. In order to validate the present method for moving‐boundary flow problems, two groups of flow phenomena have been simulated: (1) flows over a fixed circular cylinder, a harmonic in‐line oscillating cylinder in fluid at rest and a transversely oscillating cylinder in uniform flow; (2) flows over a pure pitching airfoil, a heaving–pitching airfoil and a deforming airfoil. The predictions show good agreement with the published numerical results or experimental data. Copyright © 2010 John Wiley & Sons, Ltd.