Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Stability and bifurcation behaviour of electrostatic torsional NEMS varactor influenced by dispersion forces

The influences of Casimir and van der Waals forces on the nano-electromechanical systems (NEMS) electrostatic torsional varactor are studied. A one degree of freedom, the torsional angle, is adopted, and the bifurcation behaviour of the NEMS torsional varactor is investigated. There are two bifurcation points, one of which is a Hopf bifurcation point and the other is an unstable saddle point. The phase portraits are also drawn, in which periodic orbits are around the Hopf bifurcation point, but the periodic orbit will break into a homoclinic orbit when meeting the unstable saddle point.

On the topological bifurcation of flows around a rotating circular cylinder

Flow fields around a rotating circular cylinder in a uniform stream are computed using a low dimensional Galerkin method. Results show that the formation of a Fopple vortex pair behind a stationary circular cylinder is caused by the structural instability in the vicinity of the saddle located at the rear of the cylinder. For rotating cylinder a bifurcation diagram with the consideration of two parameters, Reynolds number Re and rotation parameter a, is built by a kinematic analysis of the steady flow fields.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220

Numerical analyses for the material bifurcation in plane sheet under tension

The various patterns (shear banding, surface wrinkling and necking) of material bifurcation in plane sheet under tension are investigated in this paper by means of a numerical method. It is found that numerical analysis can provide better ground for searching for the lowest critical loads. The inhomogeneity caused by void damage and the nonuniformity in the stress distribution across sheet thickness are proved to have detrimental effects on the material bifurcation. Nevertheless, material stability can be promoted by any means of depressing void damage or alleviating stress, even locally across the thickness. Besides, the peculiar behaviour of material bifurcation under slight biaxiality state is demonstrated. Copyright (C) 1996 Elsevier Science Ltd

Hopf bifurcation in wakes behind a rotating and translating circular cylinder

A low-dimensional Galerkin method, initiated by Noack and Eckelmann [Physica D 56, 151 (1992)], for the prediction of the flow field around a stationary two-dimensional circular cylinder in a uniform stream at low Reynolds number is generalized to the case of a rotating and translating cylinder. The Hopf bifurcation describing the transition from steady to time-periodic solution is investigated. A curve indicating the transitional boundary is given in the two-dimensional parameter plane of Reynolds number Re and rotating parameter alpha. Our results show that rotation may delay the onset of vortex street and decrease the vortex-shedding frequency. (C) 1996 American Institute of Physics.

On dilatational plastic constitutive equation of ductile materials and plastic loading paths at bifurcation

The dilatational plastic constitutive equation presented in this paper is proved to be in a form of generality. Based on this equation, the constitutive behaviour of materials at the moment of bifurcation is demonstrated to follow a loading path with the response as "soft" as possible.

Two bifurcation processes for onset of oscillatory thermocapillary convection in a floating half zone

An overview on the onset of thermocapillary oscillatory convection in a floating half zone is provided, and it is a typical subject in the microgravity sciences related to the space materials science, especially the floating zone processing, and also to the microgravity fluid physics. The main interests are focused around the process for onset of oscillatory thermocapillary convection, which is known also as the bifurcation transition from quasi-steady convection to oscillatory convection. The onset of oscillation depends on a set of critical parameters, such as the Marangoni number, Prandtl number, geometrical parameters, and heat transfer parameters. Recent studies show that, there exists the bifurcation transition from steady and axial symmetric convection to the steady and axial non-symmetric convection before the onset of oscillation in cases of small Prandtl number fluids and in cases of larger Prandtl number fluids of fat liquid bridge with small aspect ratio. The transition process is a strong non-linear process because the velocity deviation has the same order of magnitude as that of an average flow after the onset of oscillation, and unsteady 3-D numerical simulation is suitable to do in depth analysis on strong non-linear process, and leads generally to a better comparison with the experimental results.

Bifurcation in a reaction diffusion system

The bifurcation and nonlinear stability properties of the Meinhardt-Gierer model for biochemical pattern formation are studied. Analyses are carried out in parameter ranges where the linearized system about a trivial solution loses stability through one to three eigenfunctions, yielding both time independent and periodic final states. Solution branches are obtained that exhibit secondary bifurcation and imperfection sensitivity and that appear, disappear, or detach themselves from other branches.

Topics in bifurcation theory

I. Existence and Structure of Bifurcation Branches

The problem of bifurcation is formulated as an operator equation in a Banach space, depending on relevant control parameters, say of the form G(u,λ) = 0. If dimN(G_u(u_O,λ_O)) = m the method of Lyapunov-Schmidt reduces the problem to the solution of m algebraic equations. The possible structure of these equations and the various types of solution behaviour are discussed. The equations are normally derived under the assumption that G^O_λεR(G^O_u). It is shown, however, that if G^O_λεR(G^O_u) then bifurcation still may occur and the local structure of such branches is determined. A new and compact proof of the existence of multiple bifurcation is derived. The linearized stability near simple bifurcation and "normal" limit points is then indicated.

II. Constructive Techniques for the Generation of Solution Branches

A method is described in which the dependence of the solution arc on a naturally occurring parameter is replaced by the dependence on a form of pseudo-arclength. This results in continuation procedures through regular and "normal" limit points. In the neighborhood of bifurcation points, however, the associated linear operator is nearly singular causing difficulty in the convergence of continuation methods. A study of the approach to singularity of this operator yields convergence proofs for an iterative method for determining the solution arc in the neighborhood of a simple bifurcation point. As a result of these considerations, a new constructive proof of bifurcation is determined.

Some modified bifurcation problems with application to imperfection sensitivity in buckling

The branching theory of solutions of certain nonlinear elliptic partial differential equations is developed, when the nonlinear term is perturbed from unforced to forced. We find families of branching points and the associated nonisolated solutions which emanate from a bifurcation point of the unforced problem. Nontrivial solution branches are constructed which contain the nonisolated solutions, and the branching is exhibited. An iteration procedure is used to establish the existence of these solutions, and a formal perturbation theory is shown to give asymptotically valid results. The stability of the solutions is examined and certain solution branches are shown to consist of minimal positive solutions. Other solution branches which do not contain branching points are also found in a neighborhood of the bifurcation point.

The qualitative features of branching points and their associated nonisolated solutions are used to obtain useful information about buckling of columns and arches. Global stability characteristics for the buckled equilibrium states of imperfect columns and arches are discussed. Asymptotic expansions for the imperfection sensitive buckling load of a column on a nonlinearly elastic foundation are found and rigorously justified.

Broadband high-gain optical parametric chirped-pulse amplification near 780 nm with LBO-I near-degenerate near-collinear phase match

Near-degenerative near-collinear phase-match geometry for broadband optical parametric chirped-pulse amplification (OPCPA) at approximate to 780 nm is calculated in comparison with nondegenerate noncollinear phase-match geometry. In an experiment on LBO-I near-degenerate near-collinear OPCPA, high gain with broad gain bandwidth (approximate to 71 nm, FWHM) at approximate to 780 nm is achieved by using an approximate to 390-nm pumping pulse. The stretched broadband chirped signal pulse near 780 nm is amplified to approximate to 412 mu J with a pumping energy of approximate to 15 mJ, and the total gain is > 3.7 X 10(6), which agrees well with the calculation. For a broadband (covering approximate to 100 nm) chirped signal pulse, the theoretical gain bandwidth has been attained experimentally for the first time. (c) 2005 Society of Photo-Optical Instrumentation Engineers.

Broadband high-gain optical parametric chirped-pulse amplification near 780 nm with LBO-I near-degenerate near-collinear phase match

Near-degenerative near-collinear phase-match geometry for broadband optical parametric chirped-pulse amplification (OPCPA) at approximate to 780 nm is calculated in comparison with nondegenerate noncollinear phase-match geometry. In an experiment on LBO-I near-degenerate near-collinear OPCPA, high gain with broad gain bandwidth (approximate to 71 nm, FWHM) at approximate to 780 nm is achieved by using an approximate to 390-nm pumping pulse. The stretched broadband chirped signal pulse near 780 nm is amplified to approximate to 412 mu J with a pumping energy of approximate to 15 mJ, and the total gain is > 3.7 X 10(6), which agrees well with the calculation. For a broadband (covering approximate to 100 nm) chirped signal pulse, the theoretical gain bandwidth has been attained experimentally for the first time. (c) 2005 Society of Photo-Optical Instrumentation Engineers.

Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

Polarization evolution of a few-cycle ultrashort laser pulse propagating in a degenerate three-level medium

The propagation of an arbitrary polarized few-cycle ultrashort laser pulse in a degenerate three-level medium is investigated by using an iterative predictor-corrector finite-difference time-domain method. It is found that the polarization evolution of the ultrashort laser pulse is dependent not only on the initial atomic coherence of the medium but also on the polarization condition of the incident laser pulse. When the initial effective area is equal to 2 pi, complete linear-to-circular and circular-to-linear polarization conversion of few-cycle ultrashort laser pulses can be achieved due to the quantum interference effects between the two different transition paths.