Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Jounal of Mechanics and Physics of Solids 42, 813-843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis er al. (Journal of the Mechanics and Physics of Solids 43, 319-341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice ei al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity. (C) 1997 Elsevier Science Ltd.

Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a = 10(-4) to 10(-2) where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions. (C) 2002 Published by Elsevier Science Ltd.

Self-similarity limits of genomic signatures

It is shown that metric representation of DNA sequences is one-to-one. By using the metric representation method, suppression of nucleotide strings in the DNA sequences is determined. For a DNA sequence, an optimal string length to display genomic signature in chaos game representation is obtained by eliminating effects of the finite sequence. The optimal string length is further shown as a self-similarity limit in computing information dimension. By using the method, self-similarity limits of bacteria complete genomic signatures are further determined.

Two-phase flow patterns in a 90 degrees bend at microgravity

Bends are widely used in pipelines carrying single- and two-phase fluids in both ground and space applications. In particular, they play more important role in space applications due to the extreme spatial constraints. In the present study, a set of experimental data of two-phase flow patterns and their transitions in a 90degrees bend with inner diameter of 12.7 mm. and curvature radius of 76.5 mm at microgravity conditions are reported. Gas and liquid superficial velocities are found to range from (1.0 similar to 23.6) m/s for gas and (0.09 similar to 0.5) m/s for liquid, respectively. Three major flow patterns, namely slug, slug-annular transitional, and annular flows, are observed in this study. Focusing on the differences between flow patterns in bends and their counterparts in straight pipes, detailed analyses of their characteristics are made. The transitions between adjoining flow patterns are found to be more or less the same as those in straight pipes, and can be predicted using Weber number models satisfactorily. The reasons for such agreement are carefully examined.

On crack path stability in a layered material

Crack paths in an elastic layer on top of a substrate are considered. Crack growth is initiated from an edge crack in the layer. The plane of the initially straight crack forms an angle to the free surface. The load consists of a pair of forces applied at the crack mouth and parallel to the interface. Crack paths are calculated using a boundary element method. Crack growth is assumed to proceed along a path for which the mode II stress intensity factor vanishes. The inclination and the length of the initial crack are varied. The effect of two different substrates on the crack path evolution is demonstrated. A crack path initially leading perpendicularly to the interface is shown to be directionally unstable for a rigid substrate. Irrespective of its initial angle, the crack does not reach the interface, but reaches the free surface if the layer is infinitely long. At finite layer length the crack reaches the upper free surface if the initial crack inclination to the surface is small enough. For an inextendable flexible substrate, on the other hand, the crack reaches the interface if its initial inclination is large enough. For the flexible substrate an unstable path parallel with the sides of an infinitely long layer is identified. The results are compared with experimental results and discussed in view of characterisation of directionally unstable crack paths. The energy release rate for an inclined edge crack is determined analytically.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

Computation of turbulent-flow in a square duct - aspects of the secondary flow and its origin

A numerical study of turbulent flow in a straight duct of square cross-section is made. An order-of-magnitude analysis of the 3-D, time-averaged Navier-Stokes equations resulted in a parabolic form of the Navier-Stokes equations. The governing equations, expressed in terms of a new vector-potential formulation, are expanded as a multi-deck structure with each deck characterized by its dominant physical forces. The resulting equations are solved using a finite-element approach with a bicubic element representation on each cross-sectional plane. The numerical integration along the streamwise direction is carried out with finite-difference approximations until a fully-developed state is reached. The computed results agree well with other numerical studies and compare very favorably with the available experimental data. One important outcome of the current investigation is the interpretation analytically that the driving force of the secondary flow in a square duct comes mainly from the second-order terms of the difference in the gradients of the normal and transverse Reynolds stresses in the axial vorticity equation.

A Finite Element Analysis of a Crack Penetrating or Deflecting into an Interface in a Thin Laminate

The crack tip driving force of a crack growing from a pre-crack that is perpendicular to and terminating at an interface between two materials is investigated using a linear fracture mechanics theory. The analysis is performed both for a crack penetrating the interface, growing straight ahead, and for a crack deflecting into the interface. The results from finite element calculations are compared with asymptotic solutions for infinitesimally small crack extensions. The solution is found to be accurate even for fairly large amounts of crack growth. Further, by comparing the crack tip driving force of the deflected crack with that of the penetrating crack, it is shown how to control the path of the crack by choosing the adhesion of the interface relative to the material toughness.

Two-phase Flow in Fuel Cells in Short-term Microgravity Condition

A visual observation of liquid-gas two-phase flow in anode channels of a direct methanol proton exchange membrane fuel cells in microgravity has been carried out in a drop tower. The anode flow bed consisted of 2 manifolds and 11 parallel straight channels. The length, width and depth of single channel with rectangular cross section was 48.0 mm, 2.5 mm and 2.0 mm, respectively. The experimental results indicated that the size of bubbles in microgravity condition is bigger than that in normal gravity. The longer the time, the bigger the bubbles. The velocity of bubbles rising is slower than that in normal gravity because buoyancy lift is very weak in microgravity. The flow pattern in anode channels could change from bubbly flow in normal gravity to slug flow in microgravity. The gas slugs blocked supply of reactants from channels to anode catalyst layer through gas diffusion layer. When the weakened mass transfer causes concentration polarization, the output performance of fuel cells declines.

郭永怀文集

Two-Phase Flow In Anode Flow Field Of A Small Direct Methanol Fuel Cell In Different Gravities

An in-situ visualization of two-phase flow inside anode flow bed of a small liquid fed direct methanol fuel cells in normal and reduced gravity has been conducted in a drop tower. The anode flow bed consists of 11 parallel straight channels. The length, width and depth of single channel, which had rectangular cross section, are 48.0, 2.5 and 2.0 mm, respectively. The rib width was 2.0 mm. The experimental results indicated that when the fuel cell orientation is vertical, two-phase flow pattern in anode channels can evolve from bubbly flow in normal gravity into slug flow in microgravity. The size of bubbles in the reduced gravity is also bigger. In microgravity, the bubbles rising speed in vertical channels is obviously slower than that in normal gravity. When the fuel cell orientation is horizontal, the slug flow in the reduced gravity has almost the same characteristic with that in normal gravity. It implies that the effect of gravity on two-phase flow is small and the bubbles removal is governed by viscous drag. When the gas slugs or gas columns occupy channels, the performance of liquid fed direct methanol fuel cells is failing rapidly. It infers that in long-term microgravity, flow bed and operating condition should be optimized to avoid concentration polarization of fuel cells.

A numerical investigation on the flow physics of koi carp (Cyprinus carpio koi) burst-and-coast swimming

Abstract—Burst-and-coast is the most common locomotion type in freely routine swimming of koi carps (Cyprinus carpio koi), which consists of a burst phase and a coast phase in each cycle and mostly leads to a straight-line trajectory. Combining with the tracking experiment, the flow physics of koi carp’s burst-andcoast swimming is investigated using a novel integrated CFD method solving the body-fluid interaction problem. The dynamical equations of a deforming body are formulated. Following that, the loose-coupled equations of the body dynamics and the fluid dynamics are numerically solved with the integrated method. The two burst modes, MT (Multiple Tail-beat) and HT (Half Tail-beat), which have been reported by the experiments, are investigated by numerical simulations in this paper. The body kinematics is predicted and the flow physics is visualized, which are in good agreement with the corresponding experiments. Furthermore, the optimization on the energy cost and several critical control mechanisms in burst-and-coast swimming of koi carps are explored, by varying the parameters in its selfpropelled swimming. In this paper, energetics is measured by the two mechanical quantities, total output power CP and Froude efficiency Fr. Results and discussion show that from the standpoint of mechanical energy, burst-and-coast swimming does not actually save energy comparing with steady swimming at the same average speed, in that frequently changing of speed leads to decrease of efficiency.

平面弯曲波导耦合器的特性分析

根据耦合模理论和弯曲波导耦合器的结构特点,对平面弯曲波导耦合器的特性进行了分析,结果表明弯曲波导耦合器的弯曲半径和最小间距两个可调变量,增加了波导器件设计的灵活性;同时由于等效耦合长度的调制作用使得弯曲波导耦合器在波分复用/解复用中比平行直波导耦合器具有更大的复用带宽;分析了弯曲半径和最小间距对弯曲波导耦合器复用带宽的影响,为实际波导器件的设计制作提供了一定的理论依据.

Phase-step calibration technique based on a two-run-times-two-frame phase-shift method

A novel phase-step calibration technique is presented on the basis of a two-run-times-two-frame phase-shift method. First the symmetry factor M is defined to describe the distribution property of the distorted phase due to phase-shifter miscalibration; then the phase-step calibration technique, in which two sets of two interferograms with a straight fringe pattern are recorded and the phase step is obtained by calculating M of the wrapped phase map, is developed. With this technique, a good mirror is required, but no uniform illumination is needed and no complex mathematical operation is involved. This technique can be carried out in situ and is applicable to any phase shifter, whether linear or nonlinear. (c) 2006 Optical Society of America.

Characteristics of 980/1550 nm WDM coupler based on planar curved waveguides

A wavelength division multiplexer (WDM) for 980/1550 nm based on planar curved waveguide coupler (CWC) is proposed. Compared with conventional parallel straight waveguide coupler (SWC), this structure has more flexibility with two variable parameters of bending radius R and minimum edge-to-edge spacing do, which are the two main parameters for the splitting ratio of coupler and decrease the complexity of device design and fabrication. Based on coupled mode theory (CMT) and waveguide theory, R and do of the WDM CWC are designed to be R = 13.28 m and d(0) = 4.39 mu m. The contrast ratio (CR) and insertion loss (IL) for 980 and 1550 nm are CR1 = 24.62 dB, CR2 = 24.56 dB and IL1 = 0.014 dB, IL2 = 0.015 dB, respectively. The 3D beam propagation method (BPM) is used to verify the validity of the design result. The influence of R and d(0) variations on the device performance is analyzed. For CR > 20 dB, the variation ranges of R and d(0) should be within -0.10 to +0.44 m and -0.05 to + 0.02 mu m, respectively. (c) 2006 Elsevier GmbH. All rights reserved.