On the Chaotic Instability of a Nonsliding Liquid-Filled Top with a Small Spheroidal Base via Melnikov-Holmes-Marsden Integrals

Chaotic orientations of a top containing a fluid filled cavity are investigated analytically and numerically under small perturbations. The top spins and rolls in nonsliding contact with a rough horizontal plane and the fluid in the ellipsoidal shaped cavity is considered to be ideal and describable by finite degrees of freedom. A Hamiltonian structure is established to facilitate the application of Melnikov-Holmes-Marsden (MHM) integrals. In particular, chaotic motion of the liquid-filled top is identified to be arisen from the transversal intersections between the stable and unstable manifolds of an approximated, disturbed flow of the liquid-filled top via the MHM integrals. The developed analytical criteria are crosschecked with numerical simulations via the 4th Runge-Kutta algorithms with adaptive time steps.     

Chaotic attitude motion of gyrostat statellites in a central force field

The nonlinear attitude motion of gyrostat satellites in a central force field is investigated, with particular emphasis on their long-time dynamic behavior for a wide range of parameters. The numbers of equilibrium solutions, as well as their stability, vary with the rotor speed, and bifurcation diagrams have been obtained. Various dynamic behaviors of gyrostat satcllites, e.g. periodic, quasiperiodic, and chaotic, are studied via the Poincaré map technique. It is shown that the rotor speed has a significant effect on the dynamic behavior of gyrostat satellites.     

Parametric identification of a chaotic base-excited double pendulum experiment

The parametric identification of a chaotic system was investigated for a double pendulum. From recorded experimental response data, the unstable periodic orbits (UPOs) were extracted and then used in a harmonic balance identification process. By applying digital filtering, digital differentiation and linear regression techniques for optimization, the results were improved. Verification of the related simulation system and linearized system also corroborated the success of the identification algorithm.     

On the existence of intermediate magnetohydrodynamic shock waves

We consider the questions related to the structure of shock waves for a system of magnetohydrodynamic equations. Using Conley's connection matrix, we recover and extend earlier results due to C. Conley and J. Smoller. In particular, we give a simpler proof of the existence of fast and slow shocks with structure. We also demonstrate that for some viscosity parameters intermediate shocks occur. Furthermore, under an assumption of transversality, we show that there exist multi-parameter families of these intermediate shocks.This research was done while both authors were visiting the Lefschetz Center for Dynamical Systems at Brown University.Supported in part by the NSF under Grant DMS-8507056.Supported in part by AFOSR 87-0347.     

About the basic integral variants of holonomic nonconservative dynamical systems

In this paper, we prove that for holonomic nonconservative dynamical system the Poincaré and Poincaré-Cartan integral invariants do not exist. Instead of them, we introduce the integral variants of Poincaré-Cartan's type and of Poincare's type for holonomic nonconservative dynamical systems, and use these variants to solve the problem of nonlinear vibration. We also prove that the integral invariants introduced in references [1] and [2] are merely the basic integral variants given by this paper. Project supported by National Natural Science Foundation of China     

Numerical modeling of the rotation of a solid body with a liquid-filled cavity having radial ribs

The boundary-value problem of unsteady vortex flow of a viscous incompressible fluid in a cylindrical vessel with radial ribs rotating at a variable angular velocity is solved using a finite-difference method. The results of the solution are used to calculate the motion of a system of a solid body and a cavity filled with a liquid. The results are compared with available experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 2, pp. 135–139, March–April, 2007.     

On the rotational motion of a gyrostat about a fixed point with mass distribution

The perturbed rotational motion of a gyrostat about a fixed point with mass distribution near to Lagrange’s case is investigated. The gyrostat is subjected under the influence of a variable restoring moment vector, a perturbing moment vector, and a third component of a gyrostatic moment vector. It is assumed that the angular velocity of the gyrostat is sufficiently large, its direction is close to the axis of dynamic symmetry, and the perturbing moments are small as compared to the restoring ones. These assumptions permit us to introduce a small parameter. Averaged systems of the equations of motion in the first and second approximations are obtained. Also, the evolution of the precession angle up to the second approximation is determined. The graphical representations of the nutation and precession angles are presented to describe the motion at any time.     

On the computation of the near-tip stress intensities for three-dimensional wedges via two-state M-integral

Singular stress fields around three-dimensional wedges are examined, and the near-tip intensity is calculated via the two-state M-integral with the aid of the domain integral representation. Two numerical examples demonstrate the effectiveness and accuracy of the present scheme for computing the stress intensities of singular stresses near the generic three-dimensional wedges.     

Solution of the laminar duct flow problem by a boundary integral equation method

Summary A boundary integral equation method is proposed for approximate numerical and exact analytical solutions to fully developed incompressible laminar flow in straight ducts of multiply or simply connected cross-section. It is based on a direct reduction of the problem to the solution of a singular integral equation for the vorticity field in the cross section of the duct. For the numerical solution of the singular integral equation, a simple discretization of it along the cross-section boundary is used. It leads to satisfactory rapid convergency and to accurate results. The concept of hydrodynamic moment of inertia is introduced in order to easily calculate the flow rate, the main velocity, and the fRe-factor. As an example, the exact analytical and, comparatively, the approximate numerical solution of the problem of a circular pipe with two circular rods are presented. In the literature, this is the first non-trivial exact analytical solution of the problem for triply connected cross section domains. The solution to the Saint-Venant torsion problem, as a special case of the laminar duct-flow problem, is herein entirely incorporated.     

On the Chaotic Behavior of a Third-Order Log-Domain Filter

Operation of log-domain filters revolves around the large-signal forward-active mode exponential current-voltage relationship of the bipolar junction transistor, which is used to map the input currents to the logarithmic domain, where the analog processing takes place, and to convert the filtered voltage waveforms back to the linear domain. Sometimes, the inherent internal nonlinearity of such circuits is the cause of unwanted externally-nonlinear effects. In this paper a differential third-order log-domain Chebyshev low-pass filter is designed by applying the method of operational simulation of doubly-terminated LC ladders. When the six grounded shunt capacitors of this filter are properly replaced with three half-sized floating capacitors, as is common practice in internally-linear fully-differential capacitively-loaded circuits, under particular conditions the resulting system exhibits limit-cycle oscillations, period-doubling bifurcations and chaotic behavior for zero-input. The presence of chaos is confirmed by extracting the spectrum of Lyapunov exponents of the system from the time series of the output voltage.     

On the dynamics of free and excited oscillations of a simple portal frame foundation

In this paper we search for the dynamics of a simple portal structure in the free and in the periodic excitation cases. By using the Center Manifold approach and Averaging Method, we obtain results on both stability and bifurcation of equilibrium points and periodic orbits.     

A domain independent integral expression for the crack extension force of a curved crack in three dimensions

An integral expression that is domain independent in curvilinear coordinates and compatible with zero divergence of Eshelby's (Phil. Trans. Roy. Soc. (London) 244 (1951) 87.) energy momentum tensor was obtained from the principle of virtual work. By applying Eshelby's definition of the force on a material defect a general expression of the crack extension force for a curved crack in three dimensions, here called the F-integral, was derived from the domain independent integral expression. The F-integral is given explicitly for a number of curved cracks and found to be in agreement with previously known solutions, when available. The influence of crack surface and crack front curvature upon the various forms of the F-integral is discussed. The F-integral presented in this work is a generalisation of the J-integral (Rice, J. Appl. Mech. 35 (1968) 379.) to curved cracks in orthogonal curvilinear coordinates.     

On partial sliding along a planar crack: the case of a circular sliding zone

Summary A problem of partial sliding along a planar crack with a local drop in frictional resistance is investigated. A sliding zone initiates in the area of reduced friction, and then propagates as the applied shear load is monotonously increased. The problem is formulated in general terms, and then solved for the case when sliding spreads as a penny-shaped zone. Conditions under which the front of the zone stays circular during sliding are analyzed. It is observed that the axisymmetry of the profile of frictional resistance does not necessarily guarantee uniform propagation of sliding in the radial direction. The circular shape becomes the most favorable growth condition only if the shear modes are related in a certain way. The problem is studied based on the criterion of propagation that stress intensity factors(SIFs) for II and III modes vanish on the boundary of the sliding zone. The singular integrals in expressions for the SIFs are reduced to non-singular ones. Analytical solutions are derived for a number of special cases where the radius of the sliding zone is related to the applied shear load.This work is supported by the National Science Foundation through grant DGE-0209543 to the University of New Mexico.     

On anti-plane shear behavior of a Griffith permeable crack in piezoelectric materials by use of the non-local theory

In this paper, the non-local theory of elasticity is applied to obtain the behavior of a Griffith crack in the piezoelectric materials under anti-plane shear loading for permeable crack surface conditions. By means of the Fourier transform the problem can be solved with the help of a pair of dual integral equations with the unknown variable being the jump of the displacement across the crack surfaces. These equations are solved by the Schmidt method. Numerical examples are provided. Unlike the classical elasticity solutions, it is found that no stress and electric displacement singularity is present at the crack tip. The non-local elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum stress hypothesis. The finite hoop stress at the crack tip depends on the crack length and the lattice parameter of the materials, respectively. The project supported by the National Natural Science Foundation of China (50232030 and 10172030)     

On the Form of the Contact Reaction in a Solid Circular Plate Simply Supported Along Two Antipodal Edge Arcs and Deflected by a Central Transverse Concentrated Force

A static, purely flexural mechanical analysis is presented for a Kirchhoff solid circular plate, deflected by a transverse central force, and bilaterally supported along two antipodal periphery arcs, the remaining part of the boundary being free. Two kinds of contact reactions are considered, namely the case of distributed reaction force alone, and the situation in which the distributed force is added to a distributed couple of properly selected profile. For both cases this plate problem is formulated in terms of an integral equation of the Prandtl type, coupled with two constraint conditions. The existence of solutions in an appropriate scaled weighted Sobolev space is discussed, and the behaviour of the solution at the endpoints of the support is exhibited. This revised version was published online in June 2006 with corrections to the Cover Date.     

Antiplane internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane

This paper deals with the antiplane magnetoelectroelastic problem of an internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane. The properties of the material such as elastic modulus, piezoelectric constant, dielectric constant, piezomagnetic coefficient, magnetoelectric coefficient and magnetic permeability are assumed in exponential forms and vary along the crack direction. Fourier transforms are used to reduce the impermeable and permeable crack problems to a system of singular integral equations, which is solved numerically by using the Gauss-Chebyshev integration technique. The stress, electric displacement and magnetic induction intensity factors at the crack tips are determined numerically. The energy density theory is applied to study the effects of nonhomogeneous material parameter β, edge conditions, location of the crack and load ratios on the fracture behavior of the internal crack.     

Analysis of iterative methods for the viscous/inviscid coupled problem via a spectral element approximation

Based on a new global variational formulation, a spectral element approximation of the incompressible Navier–Stokes/Euler coupled problem gives rise to a global discrete saddle problem. The classical Uzawa algorithm decouples the original saddle problem into two positive definite symmetric systems. Iterative solutions of such systems are feasible and attractive for large problems. It is shown that, provided an appropriate pre‐conditioner is chosen for the pressure system, the nested conjugate gradient methods can be applied to obtain rapid convergence rates. Detailed numerical examples are given to prove the quality of the pre‐conditioner. Thanks to the rapid iterative convergence, the global Uzawa algorithm takes advantage of this as compared with the classical iteration by sub‐domain procedures. Furthermore, a generalization of the pre‐conditioned iterative algorithm to flow simulation is carried out. Comparisons of computational complexity between the Navier–Stokes/Euler coupled solution and the full Navier–Stokes solution are made. It is shown that the gain obtained by using the Navier–Stokes/Euler coupled solution is generally considerable. Copyright © 2000 John Wiley & Sons, Ltd.     

On the interference of a weak and a strong shock wave in a dissociating nitrogen flow

The influence of the nitrogen dissociation on the interactions due to the interference of two planar shock waves in a hypersonic high enthalpy flow is theoretically investigated for infinite reaction rates. The two limiting cases of infinitely slow and infinitely fast reactions are modelled as a perfect gas and an ideal dissociating gas in chemical equilibrium.To investigate the influence of finite reaction rates on the interactions of shock waves, experiments are performed in the high enthalpy shock tunnel Göttingen (HEG) with a wind tunnel model consisting of a wedge type shock generator and a transversally mounted cylinder. The pressure and heat transfer loads resulting from the shock wave interferences are measured and the flow field is visualized by means of interferograms. The experimental results are compared with the results of a numerical simulation for a dissociating nitrogen flow and with the experimental results for a perfect gas flow.     

On the Displacement of NAPL Lenses and Plumes in a Phreatic Aquifer

The movement of an LNAPL lens above a sloping or horizontal water table, and a DNAPL lens above an impermeable surface is discussed. The governing equations are derived, using the vertical equilibrium approach and assuming the water mobility to be much greater than that of the NAPL. Analytical solutions are obtained for onedimensional movement of a lens along a sloping water table. They describe the lens movement with the formation of a jump at the leading front (large-scale approximation), and the distribution of NAPL in the transition zone near the jump (smallscale approximation). A model,describing the movement of a lens, taking into account NAPL retention,is proposed. Approximate onedimensional solutions for the movement of a NAPL lens along sloping or horizontal surfaces under such conditions are presented for this model. Some approximate analytical solutions for twodimensional lens (plume) formation and movement are obtained for the case of a point source at a sloping surface.     

On the case when steady converging/diverging flow of a non-Newtonian fluid in a round cone permits an exact solution

This communication considers the steady converging/diverging flow of a non-Newtonian viscous power-law fluid in a round cone. The motion is driven by a sink/source of mass at the origin. It is shown that the problem permits exact similarity solution for a particular value (n=4/3) of the fluid index. In this case a complete set of governing equations can be reduced to an ordinary differential equation, which is solved numerically for different values of the main non-dimensional parameters (the cone angle and the dimensionless sink/source intensity).