Averaging theorems for conservative systems and the weakly compressible Euler equations

A generic averaging theorem is proven for systems of ODEs with two-time scales that cannot be globally transformed into the usual action-angle variable normal form for such systems. This theorem is shown to apply to certain Fourier-space truncations of the non-isentropic slightly compressible Euler equations of fluid mechanics. For the full Euler equations, we derive formally the generic limit equations and analyze some of their properties. In the one-dimensional case, we prove a generic converic convergence result for the full Euler equations, analogous to the result for ODEs. By making use of special properties of the one-dimensional equations, we prove convergence to the solution of a more complicated set of averaged equations when the genericity assumptions fail.     

A note on geometric conditions for boundary control of wave equations with variable coefficients

The analytical condition given by Wyler for boundary stabilization of wave equations with variable coefficients is compared with the geometrical condition derived by Yao in terms of the Riemannian geometry method for exact controllability of wave equations with variable coefficients. It is shown that these two conditions are equivalent.     

Nonlinear dynamics of high-dimensional models of in-plane and out-of-plane vibration in an axially moving viscoelastic beam

Nonlinear dynamics of high-dimensional models of an axially moving viscoelastic beam with in-plane and out-of-plane vibration with combined parametric and forcing excitations are investigated by the incremental harmonic balance (IHB) method in this paper. Governing equations of transverse in-plane and out-of-plane and longitudinal vibration are obtained basing on the Hamilton's principle. The Galerkin method is used to separate time variable and spatial variable to obtain a set of multi-order differential equations. The IHB method with the fast Fourier transform (FFT) is used to solve periodic response of high-dimensional models of the beam for which convergent mode is reached. Stability of the steady-state periodic solutions is analyzed using the multivariable Floquet theory. Particular attention is paid to in-plane and out-of-plane vibration on convergent mode of the beam with combined parametric and forcing excitations. Multiple solutions are observed, and jump phenomena between in-plane and out-of-plane vibration with different transverse cross sections are discovered.     

Nonlinear Wave–Current Interactions in Shallow Water

We study here the propagation of long waves in the presence of vorticity. In the irrotational framework, the Green–Naghdi equations (also called Serre or fully nonlinear Boussinesq equations) are the standard model for the propagation of such waves. These equations couple the surface elevation to the vertically averaged horizontal velocity and are therefore independent of the vertical variable. In the presence of vorticity, the dependence on the vertical variable cannot be removed from the vorticity equation but it was however shown in 1 that the motion of the waves could be described using an extended Green–Naghdi system. In this paper, we propose an analysis of these equations, and show that they can be used to get some new insight into wave–current interactions. We show in particular that solitary waves may have a drastically different behavior in the presence of vorticity and show the existence of solitary waves of maximal amplitude with a peak at their crest, whose angle depends on the vorticity. We also show some simple numerical validations. Finally, we give some examples of wave–current interactions with a nontrivial vorticity field and topography effects.     

Analytical solution for the fully coupled static response of variable stiffness composite beams

An analytical solution is presented for the 3D static response of variable stiffness non-uniform composite beams. Based on Euler-Bernoulli theory, a set of governing differential equations are obtained, in which four degrees of freedom are fully coupled. For the variable stiffness beam, the governing field equations have variable coefficients reflecting the stiffness variation along the beam. Using the direct integration technique, the general analytical solution is derived in the integral form and the closed-form expressions of the obtained solutions are presented employing a series expansion approximation. The series expansion representation enables the proposed approach to be applicable for variable stiffness composite beams with arbitrary span-wise variation of properties. As an alternative solution, the Chebyshev collocation method is applied to the proposed formulation to verify the results obtained from the analytical solution. A number of variable stiffness composite beams made by fibre steering with various boundary conditions and stacking sequences are considered as the test cases. The static response are presented based on the analytical solution and Chebyshev collocation method and excellent agreement is observed for all test cases. The proposed model presents a reliable and efficient approach for capturing the complicated behaviour of variable stiffness non-uniform composite beams.     

Well-posedness,stability and numerical results for the thermoelastic behavior of a coupled joint-beam PDE-ODE system modeling the transverse motions of the antennas of a space structure

A mathematical model for both axial and transverse motions of two beams with cylindrical cross-sections coupled through a joint is presented and analyzed. The motivation for this problem comes from the need to accurately model damping and joint dynamics for the next generation of inflatable/rigidizable space structures. Thermo-elastic damping is included in the two beams and the motions are coupled through a joint which includes an internal moment. Thermal response in each beam is modeled by two temperature fields. The first field describes the circumferentially averaged temperature along the beam, and is linked to the axial deformation of the beam. The second describes the circumferential variation and is coupled to transverse bending. The resulting equations of motion consist of four, second-order in time, partial differential equations, four, first-order in time, partial differential equations, four second order ordinary differential equations, and certain compatibility boundary conditions. The system is written as an abstract differential equation in an appropriate Hilbert space, consisting of function spaces describing the distributed beam deflections and temperature fields, and a finite-dimensional space that projects important features at the joint boundary. Semigroup theory is used to prove that the system is well-posed, and that with positive damping parameters the resulting semigroup is exponentially stable. Steady states are characterized and several numerical approximation results are presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Self-Similarity and the Evaluation of Long-Time Nonhydrostatic Effects in Compositionally Driven Gravity Flows in Deep Surroundings

In the study of compositionally driven gravity currents involving one or more homogeneous fluid layers, it has been customary to adopt the hydrostatic assumption for the pressure field in each layer which, in turn, leads to a depth‐independent horizontal velocity field in each of these layers and significant simplifications to the governing equations. Under this hydrostatic paradigm, each layer will then have its motion governed by the well‐known reduced dimension shallow‐water equations. For the so‐called ‐layer or reduced gravity shallow‐water equations, similarity solutions for fixed volume gravity currents released in rectangular geometry have been found. Very few attempts have been made to evaluate contributions arising from the possible loss of hydrostatic balance in the context of the problems treated using the classic shallow‐water approach. Where such attempts have been pursued, they have usually been carried out in a time‐independent context or using layer‐averaged equations and very small amplitude disturbances. The vast majority of these studies into nonhydrostatic effects do not include any relevant numerical work to assess these effects. In this paper, we develop an approach for evaluating nonhydrostatic contributions to the flow field for bottom gravity currents in deep surroundings and rectangular geometry. Our approach makes no assumptions on the amplitudes of the disturbances and does not depend on layer‐averaging in the governing equations. We seek asymptotic expansions of the solutions to the Euler equations for a shallow fluid by using the small parameter δ², where δ is the aspect ratio of the flow regime. At leading order the equations enforce hydrostatic balance while those obtained at first order retain certain nonhydrostatic effects which we evaluate. Our method for evaluation of these first‐order contributions employs the self‐similar nature of the solution to the leading‐order equations in the new first‐order equations without any vertical averaging procedures being employed.     

Symmetry breaking, coupling management, and localized modes in dual-core discrete nonlinear-Schrödinger lattices

We introduce a system of two linearly coupled discrete nonlinear Schrödinger equations (DNLSEs), with the coupling constant subject to a rapid temporal modulation. The model can be realized in bimodal Bose-Einstein condensates (BEC). Using an averaging procedure based on the multiscale method, we derive a system of averaged (autonomous) equations, which take the form of coupled DNLSEs with additional nonlinear coupling terms of the four-wave-mixing type. We identify stability regions for fundamental onsite discrete symmetric solitons (single-site modes with equal norms in both components), as well as for two-site in-phase and twisted modes, the in-phase ones being completely unstable. The symmetry-breaking bifurcation, which destabilizes the fundamental symmetric solitons and gives rise to their asymmetric counterparts, is investigated too. It is demonstrated that the averaged equations provide a good approximation in all the cases. In particular, the symmetry-breaking bifurcation, which is of the pitchfork type in the framework of the averaged equations, corresponds to a Hopf bifurcation in terms of the original system.     

The Tartar Equation for Homogenization of a Model of the Dynamics of Fine-Dispersion Mixtures

We consider a mathematical model describing a nonstationary Stokes flow in a fine-dispersion mixture of viscous incompressible fluids with rapidly oscillating initial data. We perform homogenization of the model as the frequency of oscillations tends to infinity; this leads to the problem of finding effective coefficients of the averaged equations. To solve this problem, we propose and implement a method which bases on supplementing the averaged system with the Cauchy problem for the kinetic Tartar equation whose unique solution is the Tartar H-measure. Thereby we construct a correct closed model for describing the motion of a homogeneous mixture, because the effective coefficients of the averaged equations are uniquely expressed in terms of the H-measure.     

Random Fields with Multifractional Regularity Order on Heterogenous Fractal Domains

In this article, we study the effect of the geometry of a domain with variable local dimension on the regularity/singularity of the restriction of a multifractional random field on such a domain. The theories of reproducing kernel Hilbert spaces (RKHS) and generalized random fields are applied. Fractional Sobolev spaces of variable order are considered as RKHSs of random fields satisfying certain elliptic multifractional pseudodifferential equations. The multifractal spectra of these random fields are trivial due to the regularity assumptions on the variable order of the fractional derivatives. In this article, we introduce a family of RKHSs defined by isomorphic identification with the trace on a compact heterogeneous fractal domain of a fractional Sobolev space of variable order. The local regularity/singularity order of functions in these spaces, which depends on the variable order of the fractional Sobolev space considered and on the local dimension of the domain, is derived. We also study the spectral properties of the family of models introduced in the mean-square sense. In the Gaussian case, random fields with sample paths having multifractional local Hölder exponent are covered in this framework.     

Two-dimensional analysis of simply supported piezoelectric beams with variable thickness

A two-dimensional analysis is presented for piezoelectric beam with variable thickness which is simply supported and grounded along its two ends. According to the governing equations of plane stress problems, the displacement solutions, which exactly satisfy the governing differential equations and the simply-supported boundary conditions at two ends of the beam, are derived. The unknown coefficients in the solution are then determined by using the Fourier sinusoidal series expansion to the boundary equations on the upper and lower surfaces of the beams. The present solutions show a good convergence and the numerical results are presented and compared with those available in the literature. The method could be applied to control engineering and other projects with highly accurate demand on stress and displacement analysis such as the design of micro-mechanical apparatuses.     

Asymptotics for the solution of averaged equations for a system of coupled oscillators

The equations obtained by averaging a system of three weakly nonlinear oscillators are considered. Distinctive properties of asymptotics at infinity for the solutions of averaged equations for both constant and variable proper frequencies are studied. __________ Translated from Fundamentalnaya i Prikladnaya Matematika, Vol. 12, No. 6, pp. 99–113, 2006.     

Deformations of algebraic curves and integrable nonlinear evolution equations

A brief exposition of applications of the methods of algebraic geometry to systems integrable by the IST method with variable spectral parameters is presented. Usually, theta-functional solutions for these systems are generated by some deformations of algebraic curves. The deformations of algebraic curves are also related with theta-functional solutions of Yang-Mills self-duality equations which contain special systems with a variable spectral parameter as a special reduction. Another important situation in which the deformations of algebraic curves naturally occur is the KdV equation with string-like boundary conditions. Most important concrete examples of systems integrable by the IST method with variable spectral parameter having different properties within a framework of the behavior of moduli of underlying curves, analytic properties of the Baker-Akhiezer functions, and the qualitative behavior of the solutions, are vacuum axially symmetric Einstein equations, the Heisenberg cylindrical magnet equation, the deformed Maxwell-Bloch system, and the cylindrical KP equation.Dedicated to the memory of J.-L. Verdier     

Cramér's asymptotics in systems with fast and slow motions

This paper is devoted to the averaging principle for stochastic systems with slow and intermixing fast motions. Here we (i) prove the existence of the Cramér type asymptotics for the probabilities of large deviations from an averaged motion, which implies the central limit theorem, and (ii) develop an analytic procedure for computation of this asymptotics. We use general apparatus of superregular perturbations of fiber ergodic semigroups to investigate two systems in the same way. The first of them is a cascade in which slow motions are determined by a vector field depending both on slow and fast variables, and fast motions compose a Markov chain depending on the slow variable. The second is a process defined by a system of two stochastic differential equations.     

An improved modeling approach for circular flexure hinges with application to geometry optimization

Within this paper, a modeling approach for flexure hinges based on the Euler-Bernoulli beam theory for beams of variable cross section is investigated in a static analysis. The proposed approach is implemented in a finite beam element routine, for which two different discretizations are discussed. The results are compared to a full scale three dimensional model. It is shown that a circular flexure hinge cannot be modeled accurately with one element. An improved model with three elements across the flexure hinge length is presented which shows excellent accordance with the reference model. A geometry optimization is realized based on the improved, low-DOF model. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical solution of the elasticity problem of settled of the wronkler ground with variable coefficients

In this paper, we have given numerical solution of the elasticity problem of settled on the wronkler ground with variable coefficient. The approximation solution of boundary value problem which is pertinent to this has been converted to integral equations, and then by using the successive approximation methods, has been reached. In addition to this, the approximation solution of the problem was put into Padé series form. We applied these methods to an example which is the elasticity problem of unit length homogeny beam, which is a special form of boundary value problem. First we calculate the successive approximation of the given boundary value problem then transform it into Padé series form, which give an arbitrary order for solving differential equation numerically.     

变系数波方程振动传递的边界镇定

本文讨论了变系数波方程振动传递的边界镇定,应用黎曼几何方法和迹的正则性得到了所讨论问题的能量一致衰减率.     

Propagation of socillations in the Solutions of 1-D Compressible Fluid Equations

We study the propagation of initial osillations in the solutions of one-dimensional inviscid gas dynamic equations and the compressible Navier-Stokes equations. Using Multiple scale analysis, we derivbe the homogenized equations which take the form of n averaged system coupled with a dynamic cell-problem. We prove rigorous error estimates to justify the validity of these equations. We also show that the weak limits of the osicllatory solytions satisfy gas dynamic equations with an equation of state depeding on the microstructurer of the inital data     

Internal and Boundary Observability Estimates for the Heterogeneous Maxwell's System

Observability estimates for Maxwell's system with variable coefficients are established using the differential geometry method recently developed for scalar wave equations. The main tool is that Maxwell's system is reducible to a perturbed vectorial wave equation with a decoupled principal part.     

A Schr(o)idinger formulation research for light beam propagation

The wave equation of light beam propagation was written in the form of an axial-coordinate-dependent Schrdinger equation, and the expectation value of a dynamical variable, the trial function of variational approach and the ABCD law were discussed by use of quantum mechanics approach. In view of the evolution equations of expectation values of dynamical variables in the framework of quantum mechanics, the definition of a potential function representing the beam propagation stability and its universal formula with the quality factor, the universal formula of beam width and curvature radius for a paraxial beam and cylindrically symmetric non-paraxial beam, the general formula of second derivative of beam width with respect to the axial coordinate of beam for a paraxial beam, and the general criteria of the conservation of beam quality factor and the existence of a potential well of a potential function for a paraxial beam, were given or derived, respectively. Starting with the same trial function, the comparative research of our formulation with variational approach was done, which gave some further insight into the physical nature of a beam propagation parameters. The ABCD law of non-paraxial beam was discussed in terms of the definition of the non-paraxial expectation value of a dynamical variable for the first time. The applications to the media of constant second derivative of beam width with respect to the axial coordinate of a beam, square law media and the media of constant refractive index in the momentum representation were discussed, respectively.