NUMERICAL PREDICTION OF PROCESS-INDUCED RESIDUAL STRESSES IN GLASS BULB PANEL

A numerical simulation model for predicting residual stresses which arise during the solidification process of pressed glass bulb panel was developed. The solidification of a molten layer of glass between cooled parallel plates was used to model the mechanics of the buildup of residual stresses in the forming process. A thermorheologically simple thermoviscoelastic model was assumed for the material. The finite element method employed was based on the theory of shells as an assembly of flat elements. This approach calculates residual stresses layer by layer like a truly three-dimensional calculation, which is well suited for thin pressed products of complex shape. An experimental comparison was employed to verify the proposed models and methods.     

THEORETICAL INVESTIGATION ON THE DYNAMIC PERFORMANCE OF CMUT FOR DESIGN OPTIMIZATION

Although many modeling approaches exist for analyzing the behavior of capacitive micro-machined ultrasonic transducers(CMUTs),the relation equation between the design parameters with input and output is still lacking.What there is can only be used to analyze the dynamic performance of CMUT indirectly and qualitatively,such as stiffness and sound pressure.A lumped-parameter theoretical model based on the dynamic theory is proposed in this paper.The relation equations between the design parameters with inputs and outputs are given.The results obtained by the proposed model agree well with those by finite element method(FEM)simulation.The dynamic and static behavior of CMUT can be clearly depicted,which is helpful for design and optimization iterations.This shows that the proposed model makes it easier to optimize the parameters of a CMUT with respect to output and bandwidth directly and to better understand the influence of each parameter.     

Rotating electroosmotic flow of an Eyring fluid

A perturbation analysis is presented in this paper for the electroosmotic(EO) flow of an Eyring fluid through a wide rectangular microchannel that rotates about an axis perpendicular to its own. Mildly shear-thinning rheology is assumed such that at the leading order the problem reduces to that of Newtonian EO flow in a rotating channel, while the shear thinning effect shows up in a higher-order problem.Using the relaxation time as the small ordering parameter,analytical solutions are deduced for the leading-as well as first-order problems in terms of the dimensionless Debye and rotation parameters. The velocity profiles of the Ekman–electric double layer(EDL) layer, which is the boundary layer that arises when the Ekman layer and the EDL are comparably thin, are also deduced for an Eyring fluid. It is shown that the present perturbation model can yield results that are close to the exact solutions even when the ordering parameter is as large as order unity. By this order of the relaxation time parameter, the enhancing effect on the rotating EO flow due to shear-thinning Eyring rheology can be significant.     

Denoising method of X-ray phase contrast DR image for TRISO-coated fuel particles

TRISO （tristructural-isotropic） fuel is a type of micro fuel particles used in high-temperature gas-cooled reactors （HTGRs）. Among the quality evaluation methods for such particles, inqine phase contrast imaging technique （PCI） is more feasible for nondestructive measurement. Due to imaging hardware limitations, high noise level is a distinct feature of PCI images, and as a result, the dimensional measurement accuracy of TRISO-coated fuel particles decreases. Therefore, we propose an improved denoising hybrid model named as NL P-M model which introduces non-local theory and retains the merits of the Perona-Malik （P-M） model. The improved model is applied to numerical simulation and practical PCI images. Quanti- tative analysis proves that this new anisotropic diffusion model can preserve edge or texture information effectively, while ruling out noise and distinctly decreasing staircasing artifacts. Especially during the process of coating layer thickness measurement, the NL P-M model makes it easy to obtain continuous contours without noisy points or fake contour segments, thus enhancing the measurement accuracy. To address calculation complexity, a graphic processing unit （GPU） is adopted to realize the acceleration of the NL P-M denoising.     

Modeling and simulation of circulating fluidized bed reactors applied to a carbonation/calcination loop

A fluid dynamic model for a gas-solid circulating fluidized bed （CFB） designed using two coupled riser reactors is developed and implemented numerically with code programmed in Matlab. The fluid dynamic model contains heat and species mass balances to calculate temperatures and compositions for a carbonation/calcination loop process. Because of the high computational costs required to resolve the three-dimensional phenomena, a model representing a trade-offbetween computational time requirements and accuracy is developed. For dynamic processes with a solid flux between the two reactor units that depends on the fluid dynamics of both risers, a dynamic one-dimensional two-fluid model is sufficient. A two-fluid model using the constant particle viscosity closure for the stress term is used for the solid phase, and an algebraic turbulence model is applied to the gas phase. The numerical model implementa- tion is based on the finite volume method with a staggered grid scheme. The exchange of solids between the reactor units constituting the circulating fluidized bed （solid flux） is implemented through additional mass source/sink terms in the continuity equations of the two phases, For model validation, a relevant experimental analysis provided in the literature is reproduced by the numerical simulations, The numerical analysis indicates that sufficient heat integration between the two reactor units is important for the performance of the circulating fluidized bed system, The two-fluid model performs fairly well for this chemical process operated in a CFB designed as two coupled riser reactors. Further analysis and optimization of the solution algorithms and the reactor coupling strategy is warranted.     

Numerical simulation of pedestrian flow past a circular obstruction

In this paper, a revisiting Hughes' dynamic continuum model is used to investigate and predict the essential macroscopic characteristics of pedestrian flow, such as flow, density and average speed, in a two dimensional continuous walking facility scattered with a circular obstruction. It is assumed that pedestrians prefer to walk a path with the lowest instantaneous travel cost from origin to destination, under the consideration of the current traffic conditions and the tendency to avoid a high-density region and an obstruction. An algorithm for the pedestrian flow model is based on a cellcentered finite volume method for a scalar conservation law equation, a fast sweeping method for an Eikonal-type equation and a second-order TVD Runge-Kutta method for the time integration on unstructured meshes. Numerical results demonstrate the effectiveness of the algorithm. It is verified that density distribution of pedestrian flow is influenced by the position of the obstruction and the path-choice behavior of pedestrians.     

A novel type of transverse surface wave propagating in a layered structure consisting of a piezoelectric layer attached to an elastic half-space

The existence and propagation of transverse surface waves in piezoelectric coupled solids is investigated, in which perfect bonding between a metal/dielectric substrate and a piezoelectric layer of finite-thickness is assumed. Dis- persion equations relating phase velocity to material con- stants for the existence of various modes are obtained in a simple mathematical form for a piezoelectric material of class 6mm. It is discovered and proved by numerical examples in this paper that a novel Bleustein-Gulyaev （B-G） type of transverse surface wave can exist in such piezoelectric cou- pled solid media when the bulk-shear-wave velocity in the substrate is less than that in the piezoelectric layer but greater than the corresponding B-G wave velocity in the same pie- zoelectric material with an electroded surface. Such a wave does not exist in such layered structures in the absence of pie- zoelectricity. The mode shapes for displacement and electric potential in the piezoelectric layer are obtained and discussed theoretically. The study extends the regime of transverse sur- face waves and may lead to potential applications to surface acoustic wave devices.     

Experiment-based identification of time delays in linear systems

This paper presents an identification approach to time delays in single-degree-of-freedom (SDOF) and multiple-degree-of-freedom (MDOF) systems. In an SDOF system, the impedance function of the delayed system is expressed by the system parameters, the feedback gain, and the time delay. The time delay can be treated as the "fre-quency"of the difference between the impedance function of the delayed system and that of the corresponding uncon-trolled system. Thus, it can be identified from the Fourier transform of the difference between the two impedance func-tions. In an MDOF system, the pseudo-impedance functions are defined. The relationships between the time delay and the pseudo-impedance functions of the delayed system and uncontrolled system are deduced. Similarly, the time delay can be identified from the Fourier transform of the difference between the two pseudo-impedance functions. The results of numerical examples and experimental tests show that the identification approach to keeps a relatively high accuracy.     

Precracking and interfacial delamination in a bi-material structure: Static and dynamic loadings

The behavior of a precracked bi-material structure interface under given static and dynamic axial loading is an interest object in the present paper.Firstly,it is shown that the shear-lag model is a proper tool to analyze a delamination process in a precracked bi-material structure undergoing static loading.Secondly,the"shear-lag model"is applied to the structure under dynamic loading.To solve the problem for an interface delamination of the structure and to determine the debond length along the interface,our own 2D boundary element method(BEM)code is proposed in the case of static loading,and the shear-lag model together with the Laplace transforms and half-analytical calculations are used in the case of dynamic loading.The interface layer is assumed as a very thin plate compared with the other two.The parametric(geometric and elastic)analysis of the debond length and interface shear stress is done. The results from the 2D BEM code proved the validity of analytical solutions to the shear-lag model.In the dynamic case,the influence of loading characteristics,i.e.,frequencies and amplitude fluctuations on the shear stress and the value of debond length for an interval of time,is discussed. The analysis of the obtained results is illustrated by an example of the modern ceramic-metal composite,namely cermet, and depicted in figures.     

The effect of the couple stress parameter and Prandtl number on the transient natural convection flow over a vertical cylinder

An analysis is performed to study transient free convective boundary layer flow of a couple stress fluid over a vertical cylinder, in the absence of body couples. The solution of the time-dependent non-linear and coupled governing equations is carried out with the aid of an unconditionally stable Crank-Nicolson type of numerical scheme. Numerical results for the steady-state velocity, temperature as well as the time histories of the skin-friction coefficient and Nus- selt number are presented graphically and discussed. It is seen that for all flow variables as the couple stress control parameter, Co, is amplified, the time required for reaching the temporal maximum increases but the steady-state decreases.     

Experimental and Theoretical Analysis of Lamb Wave Generation by Piezoceramic Actuators for Structural Health Monitoring

Piezoceramic transducers, acting as actuators and sensors, are attractive for generation and reception of Lamb waves in Structural Health Monitoring (SHM) systems. To get insight into the source-mechanisms of Lamb waves, the vibrations of piezoceramic actuators are analyzed for the free and bonded state of the piezoceramic by analytical and finite element (FEM) calculations. Mode shapes and spectra of piezoceramic actuators and Lamb wave fields are experimentally recorded by scanning laser vibrometry. The analytical solutions for bending modes are shown to be valid for large diameter-to-thickness-relations of a free piezoactuator (D/H > 10) only. For thicker piezoceramics, a FEM-solution gives better results. Calculated frequencies for radial modes of vibration are confirmed by 3-D-laser-vibrometry and measurements of electrical impedance. The bonded case of a piezoactuator exhibits a broad resonance peak resulting from the strong coupling between radial and bending modes. The assumption that optimal excitation of Lamb modes occurs for a matching of the wavelengths to the diameter of the piezoceramic holds only for thin ceramics. Otherwise the distinct modes of out-of-plane and in-plane vibrations control the excitation of the Lamb modes more than the wavelength matching.     

Load sharing in bioinspired fibrillar adhesives with backing layer interactions and interfacial misalignment

Bio-inspired fibrillar adhesives rely on the utilization of short-range intermolecular forces harnessed by intimate contact at fibril tips. The combined adhesive strength of multiple fibrils can only be utilized if equal load sharing (ELS) is obtained at detachment. Previous investigations have highlighted that mechanical coupling of fibrils through a compliant backing layer gives rise to load concentration and the nucleation and propagation of interfacial flaws. However, misalignment of the adhesive and contacting surface has not been considered in theoretical treatments of load sharing with backing layer interactions. Alignment imperfections are difficult to avoid for a flat-on-flat interfacial configuration. In this work we demonstrate that interfacial misalignment can significantly alter load sharing and the kinematics of detachment in a model adhesive system. Load sharing regimes dominated by backing layer interactions and misalignment are revealed, the transition between which is controlled by the misalignment angle, fibril separation, and fibril compliance. In the regime dominated by misalignment, backing layer deformation can counteract misalignment giving rise to improved load sharing when compared to an identical fibrillar array with a rigid backing layer. This result challenges the conventional belief that stiffer (and thinner) backing layers consistently reduce load concentration among fibrils. Finally, we obtain analytically the fibril compliance distribution required to harness backing layer interactions to obtain ELS. Through fibril compliance optimization, ELS can be obtained even with misalignment. However, since misalignment is typically not deterministic, it is of greater practical significance that the array optimized for perfect alignment exhibits load sharing superior to that of a homogeneous array subject to misalignment. These results inform the design of fibrillar arrays with graded compliance capable of exhibiting improved load sharing over large areas.     

Dynamics of coaxial torsional resonators consisting or segments with different radii,material properties and surrounding media

High-Q torsional resonators constitute the most sensitive transducers for high frequency dynamic viscoelastic measurements of dilute polymer solutions. Most such resonators described in the literature are segmented. Because of the need for torque and torsional displacement transducers the Q-value of the individual segments most often differ, but normally all segments have the same radius.A detailed analysis of the dynamics of such resonators when both the radii, material properties and surrounding media may be different for each segment, is presented. For resonators where all segment lengths equal an integer multiple of a quarter of the torsional wavelength, we find that the Q-value of the resonator as a whole is mainly determined by the Q-value of the segment with the smallest radius. We further find that reduction of the radius of the segment surrounded by polymer solution results in a stronger mechanical coupling between the resonator as a whole and the polymer solution. These findings suggest that the segment radii are important optimization parameters of segmented torsional resonators used to measure the high frequency dynamic viscoelastic properties of e.g. polymer solutions.     

植被层湍流的大涡模拟

研究植被层湍流的大涡模拟,发展了一个TSF（transientstructurefunction）亚格于模式,尽可能真实地处理植被湍流这种既有强剪切,又有热对流的流动．我们建立了植被湍流数据库,并进行了较为详细的分析研究．湍流统计量如平均风速剖面、雷诺应力、湍流脉动能等等,与有关观测结果作了对比,符合较好．大涡模拟计算同样发现已由现场观测到的、在强对流情况时出现的温度场斜坡型有组织结构．     

Interaction of a high-speed combustion front with a closely packed bed of spheres

The head-on collision of a combustion front with a closely packed bed of ceramic-oxide spheres was investigated in a vertical 76.2 mm diameter tube containing a nitrogen diluted stoichiometric ethylene–oxygen mixture. A layer of spherical beads in the diameter range of 3–12.7 mm was placed at the bottom of the tube and a flame was ignited at the top endplate. Four orifice plates spaced at one tube diameter were placed at the ignition end of the tube in order to accelerate the flame to either a “fast-flame” or a detonation wave before the bead layer face. The mixture reactivity was adjusted by varying the initial mixture pressure between 10 and 100 kPa absolute. The pressure before and within the bead layer was measured by flush wall-mounted pressure transducers. For initial pressures where a fast-flame interacts with the bead layer peak pressures recorded at the bead layer face were as high as five times the reflected Chapman–Jouget detonation pressure. The explosion resulting from the interaction developed by two distinct mechanisms; one due to the shock reflection off the bead layer face, and the other due to shock transmission and mixing of burned and unburned gas inside the bead layer. The measured explosion delay time (time after shock reflection from the bead layer face) was found to be independent of the incident shock velocity. As a result, the explosion initiation is not the direct result of the shock reflection process but instead is more likely due to the interaction of the reflected shock wave and the trailing flame. The bead layer was found to be very effective in attenuating the explosion front transmitted through the bead layer and thus isolating the tube endplate. This paper is based on work that was presented at the 21th International Colloquium on the Dynamics of Explosions and Reactive Systems, Poitiers, France, July 23–27, 2007.     

Temporal evolution and instability in a viscoelastic dielectric elastomer

Dielectric elastomer transducers are being developed for applications in stretchable electronics, tunable optics, biomedical devices, and soft machines. These transducers exhibit highly nonlinear electromechanical behavior: a dielectric membrane under voltage can form wrinkles, undergo snap-through instability, and suffer electrical breakdown. We investigate temporal evolution and instability by conducting a large set of experiments under various prestretches and loading rates, and by developing a model that allows viscoelastic instability. We use the model to classify types of instability, and map the experimental observations according to prestretches and loading rates. The model describes the entire set of experimental observations. A new type of instability is discovered, which we call wrinkle-to-wrinkle transition. A flat membrane at a critical voltage forms wrinkles and then, at a second critical voltage, snaps into another state of winkles of a shorter wavelength. This study demonstrates that viscoelasticity is essential to the understanding of temporal evolution and instability of dielectric elastomers.     

Dynamic response of tubular dielectric elastomer transducers

In this paper, a numerical model for the dynamic response of tubular dielectric elastomer transducers is presented and validated with experimental results for the first time. Dielectric elastomers (DE) are soft polymer based smart materials that can be potentially employed in applications such as actuation, sensing and energy harvesting (Kornbluh, 2004, Carpi et al., 2005, Waki et al., 2008). In our previous work, the quasi-static response of tubular DE transducers was studied (Goulbourne et al., 2007, Son and Goulbourne, 2009). Here, a numerical model is developed to predict the dynamic response of tubular DE transducers. Inertia effects are included in our previous static model which yields a system of partial differential equations. The results of the dynamic response of the tubular DE transducers are obtained by numerically solving the simplified partial different equations using a finite difference scheme. The capacitance change induced by the dynamic deformation of the tubular DE is also calculated by a simple electrostatic model, illustrating dynamic passive sensing.Several tubular DE transducer samples (VHB 4905 and silicone) were fabricated and an experimental setup was developed to investigate the dynamic response by measuring capacitance and radial deformation. In the sensing experiments, a sweep of dynamic pressure profiles (0–5 Hz) are applied. It is observed that silicone transducers have a larger dynamic sensing range. In the actuation experiments, the deformation of the silicone actuator is monitored while a voltage signal (4.5 kV) is applied from 0 to 30 Hz. The silicone actuator shows a good actuation response. The comparison between numerical and experimental results for the DE transducers shows an overall error of 3%.     

氢燃料缓燃向爆震转捩过程中波与火焰的匹配特性

为了解氢燃料爆震过程中压力波与火焰之间相互匹配的特性,在60mm60mm2000mm 方 爆震管内,用氢气和空气混合物进行了单爆震性能研究。在爆震转捩区内布置压力传感器与离子探针,用来 监控压力波和火焰的信号,同时利用高速摄影仪集中拍摄转捩区域。根据压力波和火焰面在爆震管不同时刻 的强度特性、速度特性及位置特性来分析爆震过程中波与火焰匹配的规律。结果表明:压力波和火焰的强度 呈现为相互正反馈匹配性质;缓燃向爆震转捩(DDT)过程中,压力波和火焰的速度表现为相互交替的变化过 程,且缓燃阶段中火焰速度的增幅大于压力波速度的增幅;当火焰面追赶上激波时,产生过爆,火焰面会临时 位于激波前面;在过爆衰减为正常爆震波的过程中,激波在火焰前面。     

Mimicking Cortex Convolutions Through the Wrinkling of Growing Soft Bilayers

Brain convolutions are a specificity of mammals. Varying in intensity according to the animal species, it is measured by an index called the gyrification index, ratio between the effective surface of the cortex compared to its apparent surface. Its value is close to 1 for rodents (smooth brain), 2.6 for newborns and 5 for dolphins. For humans, any significant deviation is a signature of a pathology occurring in fetal life, which can be detected now by magnetic resonance imaging (MRI). We propose a simple model of growth for a bilayer made of the grey and white matter, the grey matter being in cortical position. We analytically solved the neo-Hookean approximation in the short and large wavelength limits. When the upper layer is softer than the bottom layer, the selection mechanism is shown to be dominated by the physical properties of the upper layer. When the anisotropy favors the growth tangentially as for the human brain, it decreases the threshold value for gyri formation. The gyrification index is predicted by a post-buckling analysis and compared with experimental data. We also discuss some pathologies in the model framework.