Nonlinear dynamic analysis of dragonfly inspired piezoelectrically driven flapping and pitching wing

The nonlinear equations for coupled elastic flapping-twisting motion of a dragonfly in- spired smart flapping wing are used for a flapping wing actuated from the root by a PZT unimorph in the piezofan configuration. Excitation by the piezoelectric harmonic force generates only the flap bending motion, which in turn, induces the elastic twist motion due to interaction between flexural and torsional vibrations modes. An unsteady aerodynamic model is used to obtain the aerodynamic forces. Numerical simulations are performed using a wing whose size is the same as the dragonfly Sympetrum Frequens wing. It is found that the value of average lift reaches to its maximum when the smart flapping wing is excited at a frequency closer to the natural frequency in torsion. Moreover, consideration of the elastic twisting of flapping wing leads to an increase in the lift force. It is also found that the flapping wing generates sufficient lift to support its own weight and carry a small pay- load. Therefore, the piezoelectrically actuated smart flapping wing based on the geometry of Sympetrum Frequens wing and undergoing flapping-twisting motions may be considered as a potential candidate for use in MAV applications.

Some problems in nonlinear elasticity

Two separate problems are discussed: axisymmetric equilibrium configurations of a circular membrane under pressure and subject to thrust along its edge, and the buckling of a circular cylindrical shell.

An ordinary differential equation governing the circular membrane is imbedded in a family of n-dimensional nonlinear equations. Phase plane methods are used to examine the number of solutions corresponding to a parameter which generalizes the thrust, as well as other parameters determining the shape of the nonlinearity and the undeformed shape of the membrane. It is found that in any number of dimensions there exists a value of the generalized thrust for which a countable infinity of solutions exist if some of the remaining parameters are made sufficiently large. Criteria describing the number of solutions in other cases are also given.

Donnell-type equations are used to model a circular cylindrical shell. The static problem of bifurcation of buckled modes from Poisson expansion is analyzed using an iteration scheme and pertubation methods. Analysis shows that although buckling loads are usually simple eigenvalues, they may have arbitrarily large but finite multiplicity when the ratio of the shell's length and circumference is rational. A numerical study of the critical buckling load for simple eigenvalues indicates that the number of waves along the axis of the deformed shell is roughly proportional to the length of the shell, suggesting the possibility of a "characteristic length." Further numerical work indicates that initial post-buckling curves are typically steep, although the load may increase or decrease. It is shown that either a sheet of solutions or two distinct branches bifurcate from a double eigenvalue. Furthermore, a shell may be subject to a uniform torque, even though one is not prescribed at the ends of the shell, through the interaction of two modes with the same number of circumferential waves. Finally, multiple time scale techniques are used to study the dynamic buckling of a rectangular plate as well as a circular cylindrical shell; transition to a new steady state amplitude determined by the nonlinearity is shown. The importance of damping in determining equilibrium configurations independent of initial conditions is illustrated.

I. Singular perturbation problems involving singular points and turning points. II. On the averaged Lagrangian technique for nonlinear dispersive waves

In Part I a class of linear boundary value problems is considered which is a simple model of boundary layer theory. The effect of zeros and singularities of the coefficients of the equations at the point where the boundary layer occurs is considered. The usual boundary layer techniques are still applicable in some cases and are used to derive uniform asymptotic expansions. In other cases it is shown that the inner and outer expansions do not overlap due to the presence of a turning point outside the boundary layer. The region near the turning point is described by a two-variable expansion. In these cases a related initial value problem is solved and then used to show formally that for the boundary value problem either a solution exists, except for a discrete set of eigenvalues, whose asymptotic behaviour is found, or the solution is non-unique. A proof is given of the validity of the two-variable expansion; in a special case this proof also demonstrates the validity of the inner and outer expansions.

Nonlinear dispersive wave equations which are governed by variational principles are considered in Part II. It is shown that the averaged Lagrangian variational principle is in fact exact. This result is used to construct perturbation schemes to enable higher order terms in the equations for the slowly varying quantities to be calculated. A simple scheme applicable to linear or near-linear equations is first derived. The specific form of the first order correction terms is derived for several examples. The stability of constant solutions to these equations is considered and it is shown that the correction terms lead to the instability cut-off found by Benjamin. A general stability criterion is given which explicitly demonstrates the conditions under which this cut-off occurs. The corrected set of equations are nonlinear dispersive equations and their stationary solutions are investigated. A more sophisticated scheme is developed for fully nonlinear equations by using an extension of the Hamiltonian formalism recently introduced by Whitham. Finally the averaged Lagrangian technique is extended to treat slowly varying multiply-periodic solutions. The adiabatic invariants for a separable mechanical system are derived by this method.

Analog methods of nonlinear vibration analysis

This thesis presents methods by which electrical analogies can be obtained for nonlinear systems. The accuracy of these methods is investigated and several specific types of nonlinear equations are studied in detail.

In Part I a general method is given for obtaining electrical analogs of nonlinear systems with one degree of freedom. Loop and node methods are compared and the stability of the loop analogy is briefly considered.

Parts II and III give a description of the equipment and a discussion of its accuracy. Comparisons are made between experimental and analytic solutions of linear systems.

Part IV is concerned with systems having a nonlinear restoring force. In particular, solutions of Duffing's equation are obtained, both by using the electrical analogy and also by approximate analytical methods.

Systems with nonlinear damping are considered in Part V. Two specific examples are chosen: (1) forced oscillations and (2) self-excited oscillations (van der Pol’s equation). Comparisons are made with approximate analytic solutions.

Part VI gives experimental data for a system obeying Mathieu's equation. Regions of stability are obtained. Examples of subharmonic, ultraharmonic, and ultrasubharmonic oscillat1ons are shown.

Discretization and solution of the inhomogeneous Bogoliubov-de Gennes equations

In the presence of inhomogeneities, defects and currents, the equations describing a Bose-condensed ensemble of alkali atoms have to be solved numerically. By combining both linear and nonlinear equations within a Discrete Variable Representation framework, we describe a computational scheme for the solution of the coupled Bogoliubov-de Gennes (BdG) and nonlinear Schrodinger (NLS) equations for fields in a 3D spheroidal potential. We use the method to calculate the collective excitation spectrum and quasiparticle mode densities for excitations of a Bose condensed gas in a spheroidal trap. The method is compared against finite-difference and spectral methods, and we find the DVR computational scheme to be superior in accuracy and efficiency for the cases we consider. (C) 2004 Elsevier B.V. All rights reserved.

Multiple Roots of Systems of Equations by Repulsion Merit Functions

In this paper we address the problem of computing multiple roots of a system of nonlinear equations through the global optimization of an appropriate merit function. The search procedure for a global minimizer of the merit function is carried out by a metaheuristic, known as harmony search, which does not require any derivative information. The multiple roots of the system are sequentially determined along several iterations of a single run, where the merit function is accordingly modified by penalty terms that aim to create repulsion areas around previously computed minimizers. A repulsion algorithm based on a multiplicative kind penalty function is proposed. Preliminary numerical experiments with a benchmark set of problems show the effectiveness of the proposed method.

Initial boundary value problems for the nonlinear Schrödinger equation

A new spectral method for solving initial boundary value problems for linear and integrable nonlinear partial differential equations in two independent variables is applied to the nonlinear Schrödinger equation and to its linearized version in the domain {x≥l(t), t≥0}. We show that there exist two cases: (a) if l″(t)<0, then the solution of the linear or nonlinear equations can be obtained by solving the respective scalar or matrix Riemann-Hilbert problem, which is defined on a time-dependent contour; (b) if l″(t)>0, then the Riemann-Hilbert problem is replaced by a respective scalar or matrix problem on a time-independent domain. In both cases, the solution is expressed in a spectrally decomposed form.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Interaction of pulses in the nonlinear Schrödinger model

A study was conducted on the interaction of two pulses in the nonlinear Schrodinger (NLS) model. The presence of different scenarios of the behavior depending on the initial parameters of the pulses, such as the pulse areas, the relative phase shift, the spatial and frequency separations were shown. It was observed that a pure real initial condition of the NLS equation can result in additional moving solitons.

Asymptotic approach for the nonlinear equatorial long wave interactions

In the present work we use an asymptotic approach to obtain the long wave equations. The shallow water equation is put as a function of an external parameter that is a measure of both the spatial scales anisotropy and the fast to slow time ratio. The values given to the external parameters are consistent with those computed using typical values of the perturbations in tropical dynamics. Asymptotically, the model converge toward the long wave model. Thus, it is possible to go toward the long wave approximation through intermediate realizable states. With this approach, the resonant nonlinear wave interactions are studied. To simplify, the reduced dynamics of a single resonant triad is used for some selected equatorial trios. It was verified by both theoretical and numerical results that the nonlinear energy exchange period increases smoothly as we move toward the long wave approach. The magnitude of the energy exchanges is also modified, but in this case depends on the particular triad used and also on the initial energy partition among the triad components. Some implications of the results for the tropical dynamics are disccussed. In particular, we discuss the implications of the results for El Nĩo and the Madden-Julian in connection with other scales of time and spatial variability. © Published under licence by IOP Publishing Ltd.

Bistable elliptic equations with fractional diffusion

This work concerns the study of bounded solutions to elliptic nonlinear equations with fractional diffusion. More precisely, the aim of this thesis is to investigate some open questions related to a conjecture of De Giorgi about the one-dimensional symmetry of bounded monotone solutions in all space, at least up to dimension 8. This property on 1-D symmetry of monotone solutions for fractional equations was known in dimension n=2. The question remained open for n>2. In this work we establish new sharp energy estimates and one-dimensional symmetry property in dimension 3 for certain solutions of fractional equations. Moreover we study a particular type of solutions, called saddle-shaped solutions, which are the candidates to be global minimizers not one-dimensional in dimensions bigger or equal than 8. This is an open problem and it is expected to be true from the classical theory of minimal surfaces.

Intermediate asymptotics in nonlinear optical systems

Using a fiber laser system as a specific illustrative example, we introduce the concept of intermediate asymptotic states in finite nonlinear optical systems. We show that intermediate asymptotics of nonlinear equations (e.g., coherent structures with a finite lifetime or distance) can be used in applications similar to those of truly stable asymptotic solutions, such as, e.g., solitons and dissipative nonlinear waves. Applying this general idea to a particular, albeit practically important, physical system, we demonstrate a novel type of nonlinear pulse-shaping regime in a mode-locked fiber laser leading to the generation of linearly chirped pulses with a triangular distribution of the intensity.

Intermediate asymptotics in nonlinear optical systems

Using a fiber laser system as a specific illustrative example, we introduce the concept of intermediate asymptotic states in finite nonlinear optical systems. We show that intermediate asymptotics of nonlinear equations (e.g., coherent structures with a finite lifetime or distance) can be used in applications similar to those of truly stable asymptotic solutions, such as, e.g., solitons and dissipative nonlinear waves. Applying this general idea to a particular, albeit practically important, physical system, we demonstrate a novel type of nonlinear pulse-shaping regime in a mode-locked fiber laser leading to the generation of linearly chirped pulses with a triangular distribution of the intensity.

Steffensen Methods for Solving Generalized Equations

2000 Mathematics Subject Classification: 65G99, 65K10, 47H04.