Exact Solutions for Linearly Electroelastic Plate-Like Bodies

A three-dimensional, explicit and exact solution is derived for a transversely isotropic, linearly electroelastic body in the form of a right cylinder of arbitrary cross section, being simply supported and connected to ground over its lateral boundary, and subject to an arbitrary distribution of force and charge over its end faces. When electric phenomena are ignored, this solution reduces to the solution given in [9] for linearly elastic plate-like bodies.     

Experimental studies on heat and mass transfer performance of a coiled tube absorber for R134a-DMAC based absorption cooling system

Absorber is an important component in vapor absorption refrigeration system and its performance has greater influence in overall efficiency of absorption machines. Falling film heat and mass transfer in an absorber is greatly influenced by fluid properties, geometry of heat exchanger and its operating parameters. This paper presents on the results of experimental studies on the heat and mass transfer characteristics of a coiled tube falling film absorber, using 1,1,1,2-Tetrafluroethane(R-134a) and N-N Dimethyl Acetamide (DMAC) as working fluids. The effects of film Reynolds number, inlet solution temperature and cooling water temperature on absorber heat load, over all heat transfer coefficient and mass of refrigerant absorbed are presented and discussed. Normalized solution and coolant temperature profiles and refrigerant mass absorbed along the height of absorber are also observed from the experimental results. The optimum over all heat transfer coefficient for R-134a–DMAC solution found to be 726 W/m²K for a film Reynolds number of 350. The R-134a vapour absorption rate is maximum in the normalized coil height of 0.6 to 1.     

Computation of functionals over discontinuous sample paths of a stochastic van der pol oscillator

This paper deals with a random van der Pol oscillator. It is assumed that the oscillator is subjected to two different kinds of perturbation. The first kind of perturbation is represented by the standard Wiener process and the second kind by a homogeneous process with independent increments, finite second order moments, mean zero and no continuous sample functions. In order to measure quantitatively the stochastic stability of the oscillator, two functionals are defined over its phase plane sample paths. It is shown that each of these functionals is a solution to a corresponding partial integro-differential equation. A numerical procedure for the solution of these equations, is suggested, and its efficiency and applicability are demonstrated with examples.     

Two‐dimensional anisotropic Cartesian mesh adaptation for the compressible Euler equations

Simulating transient compressible flows involving shock waves presents challenges to the CFD practitioner in terms of the mesh quality required to resolve discontinuities and prevent smearing. This paper discusses a novel two‐dimensional Cartesian anisotropic mesh adaptation technique implemented for transient compressible flow. This technique, originally developed for laminar incompressible flow, is efficient because it refines and coarsens cells using criteria that consider the solution in each of the cardinal directions separately. In this paper, the method will be applied to compressible flow. The procedure shows promise in its ability to deliver good quality solutions while achieving computational savings. Transient shock wave diffraction over a backward step and shock reflection over a forward step are considered as test cases because they demonstrate that the quality of the solution can be maintained as the mesh is refined and coarsened in time. The data structure is explained in relation to the computational mesh, and the object‐oriented design and implementation of the code is presented. Refinement and coarsening algorithms are outlined. Computational savings over uniform and isotropic mesh approaches are shown to be significant. Copyright © 2004 John Wiley & Sons, Ltd.     

Dynamical behavior of a class of prey-predator system with impulsive state feedback control and Beddington–DeAngelis functional response

This paper studies systematically a Bedd-ington?CDeAngelis prey?Cpredator system with harvesting and impulsive state feedback control. Conditions for existence and stability of predator-free periodic solution are obtained. When the predator-free periodic solution loses its stability, the existence and stability of nontrivial period solution are also established. Furthermore, computer simulations show that this impulsive system displays a series of complex phenomena, including period-doubling bifurcation and cascade, period window, and chaotic bands. Through numerical simulation, it is also observed that capture capability can influence the amount of predator released and the interval of the stability for nontrivial period-1 solution. Moreover, the superiority of impulsive state feedback control strategy is also exhibited over the impulsive fixed-time control.     

A new exact solution for the flow of a Maxwell fluid past an infinite plate

A new exact solution corresponding to the flow of a Maxwell fluid over a suddenly moved flat plate is determined. This solution is in all accordance with a previous one and for λ→0 it goes to the well-known solution for Navier-Stokes fluids.     

Three-dimensional isothermal lava flows over a non-axisymmetric conical surface

Within the thin-layer approximation for a highly-viscous heavy incompressible fluid, a hydrodynamicmodel of a 3D isothermal lava flow over a non-axisymmetric conical surface is constructed. Using analytical methods, a self-similar solution for the law of leading-edge propagation is obtained in the case of a flow from a non-axisymmetric source located at the apex of a conical surface with smoothly varying properties. In the case of a flow over a substantially non-axisymmetric surface, it is shown that there exists a self-similar solution for the free-surface shape and the law of leading-edge motion. This solution is studied numerically for particular examples of the substrate surface and the source. In the general case of a non-self-similar flow over a substantially non-axisymmetric conical surface, a local analytical solution is obtained for the free-surface shape and the velocity field near the leading flow front.     

On time independent Schr?dinger equations in quantum mechanics by the homotopy analysis method

A general analytic approach, namely the homotopy analysis method(HAM), is applied to solve the time independent Schr?dinger equations. Unlike perturbation method, the HAM-based approach does not depend on any small physical parameters, corresponding to small disturbances.Especially, it provides a convenient way to gain the convergent series solution of quantum mechanics. This study illustrates the advantages of this HAM-based approach over the traditional perturbative approach, and its general validity for the Schr?dinger equations. Note that perturbation methods are widely used in quantum mechanics, but perturbation results are hardly convergent. This study suggests that the HAM might provide us a new, powerful alternative to gain convergent series solution for some complicated problems in quantum mechanics, including many-body problems, which can be directly compared with the experiment data to improve the accuracy of the experimental findings and/or physical theories.     

基于运动约束解过约束并联机构变形协调方程

提出利用运动约束关系来间接求解过约束并联机构变形协调方程．首先介绍了该方法的原理,接着分别针对平面和空间过约束并联机构,详述该方法的解决步骤,结果验证了该方法的正确性,从中还可看出该方法在求解复杂过约束并联机构时非常简洁,最后介绍了采用该方法解决多度过约束问题．     

Micropolar fluid flow over a shrinking sheet

An analysis is carried out for the steady two-dimensional flow of a micropolar fluid over a shrinking sheet in its own plane. The shrinking velocity is assumed to vary linearly with the distance from a fixed point on the sheet. The features of the flow and heat transfer characteristics are analyzed and discussed. It is found that the solution exists only if adequate suction through the permeable sheet is introduced. Moreover, stronger suction is necessary for the solution to exist for a micropolar fluid compared to a classical Newtonian fluid. Dual solutions are obtained for certain suction and material parameters.     

The transient response of second order non-linear equations

A method is presented for calculating the transient response of second order non-linear differential equations. The increase in accuracy over such methods as the Kryloff-Bogoliuboff technique, physically speaking, is due to the inclusion of the phase-lag effects caused by the linear damping term in the equation. Mathematically the increase in accuracy is obtained by an appropriate choice of the zero-th order terms for the dependent variable and its first derivative. This inclusion requires little more calculation effort than that required in the Kryloff-Bogoliuoff technique. Additionally, this method may be iterated to obtain solutions of increased accuracy ; the next higher order solution being carried out in detail in the text. Examples are presented to show the increase in accuracy of the present method over the Kryloff-Bogoliuboff technique.     

Buckling of Columns with Movable Boundaries∗

The Euler buckling problem of a slender tubular column subject to its own weight, tension or compression exerted at its top, and internal and external variable static fluid pressure is studied. This problem finds many applications in drilling and production risers, mining risers, hydraulic columns, and legs of Tension Leg Platforms. The supports at the upper end of the column are considered movable to properly simulate drill ships or platforms that support such columns. The corresponding eigenvalue problem is comprised of a fourth-order differential equation with variable coefficients and four homogeneous boundary conditions. Two methods of solution derived in previous work are implemented numerically. The first solution is expressed in terms of Airy functions of the first and second kind and the second in terms of power series. The combined results of the two methods yield the critical buckling curves over the entire domain of practical interest. The critical curves are plotted in the plane of the two loading variables for the first six buckling modes for four different sets of boundary conditions. The results reveal the dependence of the asymptotic behavior of the critical curves for long columns on the boundary conditions.     

A solution for the vertical variation of stress,rather than velocity,in a three-dimensional circulation model

A simple technique is presented that allows a numerical solution to be sought for the vertical variation of shear stress as a substitute for the vertical variation of velocity in a three-dimensional hydrodynamic model. In its most general form the direct stress solution (DSS) method depends only upon the validity of an eddy viscosity relation between the shear stress and the vertical gradient of velocity. The rationale for preferring a numerical solution for shear stress to one for velocity is that shear stress tends to vary more slowly over the vertical than velocity, particularly near boundaries. Consequently, a numerical solution can be obtained much more efficiently for shear stress than for velocity. When needed, the velocity profile can be recovered from the stress profile by solving a one-dimensional integral equation over the vertical. For most practical problems this equation can be solved in closed form. Comparisons are presented between the DSS technique, the standard velocity solution technique and analytical solutions for wind-driven circulation in an unstratified, closed, rectangular channel governed by the linear equations of motion. In no case was the computational effort required by the velocity solution competitive with the DSS when a physically realistic boundary layer was included. The DSS technique should be particularly beneficial in numerical models of relatively shallow water bodies in which the bottom and surface boundary layers occupy a significant portion of the water column.     

高超声速尾迹流场稳定性数值研究

通过数值模拟, 对高超声速尾迹流场进行了研究, 对其尾迹流动的失稳过程进行了分析.选取计算模型为圆球,Ma= 6.0, Re = 1.71\times 10^6(Re以球头半径为参考长度). 通过数值模拟,首先得到的流动是稳定解,在底部发展出一个主分离区和一个二次分离区,流动是轴对称状态. 不添加任何扰动继续进行计算,发现底部流场缓慢发展出微弱的非定常流动. 随后,该现象继续发展,出现明显的结构失稳,得到了无量纲周期为12.0的周期解. 给出了高超声速圆球绕流尾迹结构的周期性演化过程,对其涡系结构的演化及奇点特征进行了分析. 研究表明该数值模拟方法可用于底部流动稳定性问题的研究,同时证实了高超声速底部流动也存在流动不稳定性.      

Dynamics of a Pipeline under the Action of Internal Shock Pressure

The static and dynamic bending of a pipeline in the vertical plane under the action of its own weight is considered with regard to the interaction of the internal pressure with the curvature of the axial line and the axisymmetric deformation. The pressure consists of a constant and timevarying parts and is assumed to be uniformly distributed over the entire span between the supports. The pipeline reaction to the stepwise increase in the pressure is analyzed in the case where it is possible to determine the exact solution of the problem. The initial stage of bending determined by the smallness of elastic forces as compared to the inertial forces is introduced into the consideration. At this stage, the solution is sought in the form of power series and the law of pressure variation can be arbitrary. This solution provides initial conditions for determining the further process. The duration of the inertial stage is compared with the times of sharp changes of the pressure and the shock waves in fluids. The structure parameters are determined in the case where the shock pressure is accepted only by the inertial forces in the pipeline.     

A discontinuous Galerkin method with plane waves and Lagrange multipliers for the solution of short wave exterior Helmholtz problems on unstructured meshes

Recently, a discontinuous Galerkin method with plane wave basis functions and Lagrange multiplier degrees of freedom was proposed for the efficient solution of the Helmholtz equation in the mid-frequency regime. This method was fully developed however only for regular meshes, and demonstrated only for interior Helmholtz problems. In this paper, we extend it to irregular meshes and exterior Helmholtz problems in order to expand its scope to practical acoustic scattering problems. We report preliminary results for two-dimensional short wave problems that highlight the superior performance of this discontinuous Galerkin method over the standard finite element method.     

Wave propagation in a fluid flowing through a curved thin-walled elastic tube

The pulsatile flow in a curved elastic pipe of circular cross section is investigated. The unsteady flow of a viscous fluid and the wall motion equations are written in a toroidal coordinate system, superimposed and linearized over a steady state solution. Being the main application relative to the vascular system, the radius of the pipe is assumed small compared with the radius of curvature. This allows an asymptotic analysis over the curvature parameter. The model results an extension of the Womersley's model for the straight elastic tube. A numerical solution is found for the first order approximation and computational results are finally presented, demonstrating the role of curvature in the wave propagation and in the development of a secondary flow.     

Computation of the stabilization parameter for the finite element solution of advective–diffusive problems

In a previous paper a general procedure for deriving stabilized finite element schemes for advective type problems based on invoking higher order balance laws over finite size domains was presented. This provides an expression for the element stabilization parameter in terms of the solution residual and its first derivatives in a kind of iterative or adaptative manner. Details of the application of this procedure to 1D and 2D advective–diffusive problems are given. Some examples of applications showing the potential of the new approach are presented. © 1997 John Wiley & Sons, Ltd.