Nonlinear static and dynamic damped response of a composite rectangular plate resting on a two-parameter elastic foundation

Nonlinear static and dynamic response analyses of a clamped. rectangular composite plate resting on a two-parameter elastic foundation have been studied using von Karman's relations. Incorporating the material damping, the governing coupled, nonlinear partial differential equations are obtained for the plate under step pressure pulse load excitation. These equations have been solved by a one-term solution and by applying Galerkin's technique to the deflection equation. This yields an ordinary nonlinear differential equation in time. The nonlinear static solution is obtained by neglecting the time-dependent variables. Thc nonlinear dynamic damped response is obtained by applying the ultraspherical polynomial approximation (UPA) technique. The influences of foundation modulus, shear modulus, orthotropy, etc. upon the nonlinear static and dynamic responses have been presented.

Unsteady compressible flow due to relative rotation and translation of a viscous fluid on a disk

The modification of the axisymmetric viscous flow due to relative rotation of the disk or fluid by a translation of the boundary is studied. The fluid is taken to be compressible, and the relative rotation and translation velocity of the disk or fluid are time-dependent. The nonlinear partial differential equations governing the motion are solved numerically using an implicit finite difference scheme and Newton's linearisation technique. Numerical solutions are obtained at various non-dimensional times and disk temperatures. The non-symmetric part of the flow (secondary flow) describing the translation effect generates a velocity field at each plane parallel to the disk. The cartesian components of velocity due to secondary flow exhibit oscillations when the motion is due to rotation of the fluid on a translating disk. Increase in translation velocity produces an increment in the radial skin friction but reduces the tangential skin friction.

Unsteady MUD rotating flow over a rotating sphere near the equator

The unsteady rotating flow of a laminar incompressible viscous electrically conducting fluid over a rotating sphere in the vicinity of the equator has been studied. The fluid and the body rotate either in the same direction or in opposite directions. The effects of surface suction and magnetic field have been included in the analysis. There is an initial steady state that is perturbed by a sudden change in the rotational velocity of the sphere, and this causes unsteadiness in the flow field. The nonlinear coupled parabolic partial differential equations governing the boundary-layer flow have been solved numerically by using an implicit finite-difference scheme. For large suction or magnetic field, analytical solutions have also been obtained. The magnitude of the radial, meridional and rotational velocity components is found to be higher when the fluid and the body rotate in opposite directions than when they rotate in the same direction. The surface shear stresses in the meridional and rotational directions change sign when the ratio of the angular velocities of the sphere and the fluid lambda greater than or equal to lambda(0). The final (new) steady state is reached rather quickly which implies that the spin-up time is small. The magnetic field and surface suction reduce the meridional shear stress, but increase the surface shear stress in the rotational direction.

Longitudinal Magnetic Field Effect on Nonlocal Ultrasonic Vibration Analysis of Single-Walled Carbon Nanotubes Based on Wave Propagation Approach

This paper studies the effect of longitudinal magnetic field on ultrasonic vibration in single walled carbon nanotubes (CNTs) based on nonlocal continuum medium theory. Governing partial differential equations of CNTs are derived by considering the Lorentz magnetic forces applied on CNTs induced by a longitudinal magnetic field through Maxwell equations. The vibration characteristics of CNTs under a longitudinal magnetic field are obtained by solving the governing equations via wave propagation approach. The effects of longitudinal magnetic field on vibration of CNTs are discussed through numerical experiments. The present analysis show that vibration frequencies of CNTs drops dramatically in the presence of the magnetic field for various circumferential wavenumbers. Such effect is also observed for various boundary conditions of the CNT. New features for the effect of longitudinal magnetic field on ultrasonic vibration of CNTs, presented in this paper are useful in the design of nano-drive device, nano-oscillator and actuators and nano-electron technology, where carbon nanotubes act as basic elements.

A Finite Element Method Based Approach for Evaluating the Sensitivity of Faraday-Type Electromagnetic Flow Meter for Liquid Sodium Flow

Faraday-type electromagnetic flow meters are employed for measuring the flow rate of liquid sodium in fast breeder reactors. The calibration of such flow meters, owing to the required elaborative arrangements is rather difficult. On the other hand, theoretical approach requires solution of two coupled electromagnetic partial differential equation with profile of the flow and applied magnetic field as the inputs. This is also quite involved due to the 3D nature of the problem. Alternatively, Galerkin finite element method based numerical solution is suggested in the literature as an attractive option for the required calibration. Based on the same, a computer code in Matlab platform has been developed in this work with both 20 and 27 node brick elements. The boundary conditions are correctly defined and several intermediate validation exercises are carried out. Finally it is shown that the sensitivities predicted by the code for flow meters of four different dimensions agrees well with the results given by analytical expression, thereby providing strong validation. Sensitivity for higher flow rates, for which analytical approach does not exist, is shown to decrease with increase in flow velocity.

Closed-form solutions for axially functionally graded Timoshenko beams having uniform cross-section and fixed-fixed boundary condition

In this paper, we study the free vibration of axially functionally graded (AFG) Timoshenko beams, with uniform cross-section and having fixed-fixed boundary condition. For certain polynomial variations of the material mass density, elastic modulus and shear modulus, along the length of the beam, there exists a fundamental closed form solution to the coupled second order governing differential equations with variable coefficients. It is found that there are an infinite number of non-homogeneous Timoshenko beams, with various material mass density, elastic modulus and shear modulus distributions having simple polynomial variations, which share the same fundamental frequency. The derived results can be used as benchmark solutions for testing approximate or numerical methods used for the vibration analysis of non-homogeneous Timoshenko beams. They can also be useful for designing fixed-fixed non-homogeneous Timoshenko beams which may be required to vibrate with a particular frequency. (C) 2013 Elsevier Ltd. All rights reserved.

Analytical test functions for free vibration analysis of rotating non-homogeneous Timoshenko beams

In this paper, the governing equations for free vibration of a non-homogeneous rotating Timoshenko beam, having uniform cross-section, is studied using an inverse problem approach, for both cantilever and pinned-free boundary conditions. The bending displacement and the rotation due to bending are assumed to be simple polynomials which satisfy all four boundary conditions. It is found that for certain polynomial variations of the material mass density, elastic modulus and shear modulus, along the length of the beam, the assumed polynomials serve as simple closed form solutions to the coupled second order governing differential equations with variable coefficients. It is found that there are an infinite number of analytical polynomial functions possible for material mass density, shear modulus and elastic modulus distributions, which share the same frequency and mode shape for a particular mode. The derived results are intended to serve as benchmark solutions for testing approximate or numerical methods used for the vibration analysis of rotating non-homogeneous Timoshenko beams.

A C-0 Interior Penalty Method for a Fourth-order Variational Inequality of the Second Kind

In this article, we propose a C-0 interior penalty ((CIP)-I-0) method for the frictional plate contact problem and derive both a priori and a posteriori error estimates. We derive an abstract error estimate in the energy norm without additional regularity assumption on the exact solution. The a priori error estimate is of optimal order whenever the solution is regular. Further, we derive a reliable and efficient a posteriori error estimator. Numerical experiments are presented to illustrate the theoretical results. (c) 2015Wiley Periodicals, Inc.

Magnetostatic equilibrium of isolated magnetic-flux tube

In the present paper, an isolated axisymmetric flux tube is discussed for slender magnetic configuration. The magnetostatic model and the stratified atmospheric model are applied, respectively, to the regions inside and outside the flux tube. The problem is described mathematically by the nonlinear partial differential equations under the nonlinear boundary condition at the free boundary of flux tube. According to the approximation of a small expansive angle, the solutions of series expressions are obtained formally. The model of polytropic plasma is discussed in detail especially. The results show the distributions of thermodynamic quantities and magnetic field extending from the high β region to the low β region, and the flux tube may be either divergent or convergent according to the pressure difference outside and inside the flux tube.

The thermal and mechanical structure of a two-dimensional plume in the earths mantle

On the condition that the distribution of velocity and temperature at the mid-plane of a mantle plume has been obtained (pages 213–218, this issue), the problem of determining the lateral structure of the plume at a given depth is reduced to solving an eigenvalue problem of a set of ordinary differential equations with five unknown functions, with an eigenvalue being related to the thermal thickness of the plume at this depth. The lateral profiles of upward velocity, temperature and viscosity in the plume and the thickness of the plume at various depths are calculated for two sets of Newtonian rheological parameters. The calculations show that the precondition for the existence of the plume, δT/L 1 (L = the height of the plume, δT = lateral distance from the mid-plane), can be satisfied, except for the starting region of the plume or near the base of the lithosphere. At the lateral distance, δT, the upward velocity decreases to 0.1 – 50% of its maximum value at different depths. It is believed that this model may provide an approach for a quantitative description of the detailed structure of a mantle plume.

强外加电场与大调制度下光折变动力学光栅形成研究

强外加电场与大调制度在光折变效应的研究中已经得到了广泛应用。采用PDECOL算法, 严格求解光折变带输运方程, 得到外加电场时不同调制度下光折变晶体中随时间变化的空间电荷场、载流子浓度, 并讨论了外加电场对它们的影响。通过将物质方程与耦合波方程联立数值求解, 可得到光折变光栅形成过程中两波耦合增益系数以及光束条纹相位的变化。模拟结果表明, 在强外加电场作用下, 两束记录光之间的光强与相位耦合都得到了增强, 而原有的解析式忽视了强外加电场与大调制度对空间电荷场相位耦合的影响, 此时不再适用。同时发现折射率光

On the Cesari fixed point method in a Banach space

In a paper published in 1961, L. Cesari [1] introduces a method which extends certain earlier existence theorems of Cesari and Hale ([2] to [6]) for perturbation problems to strictly nonlinear problems. Various authors ([1], [7] to [15]) have now applied this method to nonlinear ordinary and partial differential equations. The basic idea of the method is to use the contraction principle to reduce an infinite-dimensional fixed point problem to a finite-dimensional problem which may be attacked using the methods of fixed point indexes.

The following is my formulation of the Cesari fixed point method:

Let B be a Banach space and let S be a finite-dimensional linear subspace of B. Let P be a projection of B onto S and suppose Г≤B such that pГ is compact and such that for every x in PГ, P^-1x∩Г is closed. Let W be a continuous mapping from Г into B. The Cesari method gives sufficient conditions for the existence of a fixed point of W in Г.

Let I denote the identity mapping in B. Clearly y = Wy for some y in Г if and only if both of the following conditions hold:

(i) Py = PWy.

(ii) y = (P + (I - P)W)y.

Definition. The Cesari fixed paint method applies to (Г, W, P) if and only if the following three conditions are satisfied:

(1) For each x in PГ, P + (I - P)W is a contraction from P^-1x∩Г into itself. Let y(x) be that element (uniqueness follows from the contraction principle) of P^-1x∩Г which satisfies the equation y(x) = Py(x) + (I-P)Wy(x).

(2) The function y just defined is continuous from PГ into B.

(3) There are no fixed points of PWy on the boundary of PГ, so that the (finite- dimensional) fixed point index i(PWy, int PГ) is defined.

Definition. If the Cesari fixed point method applies to (Г, W, P) then define i(Г, W, P) to be the index i(PWy, int PГ).

The three theorems of this thesis can now be easily stated.

Theorem 1 (Cesari). If i(Г, W, P) is defined and i(Г, W, P) ≠0, then there is a fixed point of W in Г.

Theorem 2. Let the Cesari fixed point method apply to both (Г, W, P₁) and (Г, W, P₂). Assume that P₂P₁=P₁P₂=P₁ and assume that either of the following two conditions holds:

(1) For every b in B and every z in the range of P₂, we have that ‖b=P₂b‖ ≤ ‖b-z‖

(2)P₂Г is convex.

Then i(Г, W, P₁) = i(Г, W, P₂).

Theorem 3. If Ω is a bounded open set and W is a compact operator defined on Ω so that the (infinite-dimensional) Leray-Schauder index i_LS(W, Ω) is defined, and if the Cesari fixed point method applies to (Ω, W, P), then i(Ω, W, P) = i_LS(W, Ω).

Theorems 2 and 3 are proved using mainly a homotopy theorem and a reduction theorem for the finite-dimensional and the Leray-Schauder indexes. These and other properties of indexes will be listed before the theorem in which they are used.

SBLOCK: A framework for efficient stencil-based PDE solvers on multi-core platforms

We present a new software framework for the implementation of applications that use stencil computations on block-structured grids to solve partial differential equations. A key feature of the framework is the extensive use of automatic source code generation which is used to achieve high performance on a range of leading multi-core processors. Results are presented for a simple model stencil running on Intel and AMD CPUs as well as the NVIDIA GT200 GPU. The generality of the framework is demonstrated through the implementation of a complete application consisting of many different stencil computations, taken from the field of computational fluid dynamics. © 2010 IEEE.

A semi-analytical theory of coastal upwelling and jet

The general forms of the conservation of momentum, temperature and potential vorticity of coastal ocean are obtained in the x-z plane for the nonlinear ocean circulation of Boussinesq fluid, and a elliptic type partial differential equations of second order are derived. Solution of the partial differential equations are obtained under the conditions that the fluid moves along the topography. The numerical results show that there exist both upwelling and downwelling along coastline that mainly depends on the large scale ocean condition. Numerically results of the upwelling (downwelling), coastal jet and temperature front zone are favorable to the observations.