Invariant solutions of the Westervelt model of nonlinear hydroacoustics without dissipation

We study three-dimensional Westervelt model of a nonlinear hydroacoustics without dissipation. We received all of its invariant submodels. We studied all invariant submodels described by the invariant solutions of rank 0 and 1. All invariant solutions of rank 0 and 1 are found either explicitly, or their search is reduced to the solution of the nonlinear integral equations. With a help of these invariant solutions we researched: (1) a propagation of the intensive acoustic waves (self-similar, axisymmetric, planar and one-dimensional) for which the acoustic pressure and a speed of its change, or the acoustic pressure and its derivative in the direction of one of the axes are specified at the initial moment of the time at a fixed point , (2) a spherically symmetric ultrasonic field for which the acoustic pressure and a speed of its change, or the acoustic pressure and its radial derivative are specified at the initial moment of the time at a fixed point. Solving of the boundary value problems describing these processes is reduced to the solving of nonlinear integral equations. We are established the existence and uniqueness of solutions of these boundary value problems under some additional conditions. Mechanical relevance of the obtained solutions is as follows: (1) these solutions describe nonlinear and diffraction effects in ultrasonic fields of a special kind, (2) these solutions can be used as a test solutions in the numerical calculations performed in studies of ultrasonic fields generated by powerful emitters. We found all the conservation laws of the first order for the Westerveld equation written in dimensionless variables.     

Invariant submodels and exact solutions of Khokhlov–Zabolotskaya–Kuznetsov model of nonlinear hydroacoustics with dissipation

We study three-dimensional Khokhlov–Zabolotskaya–Kuznetsov (KZK) model of the nonlinear hydroacoustics with dissipation. This model is described by third order quasilinear partial differential equation of the (KZK). We obtained that the (KZK) equation admits an infinite Lie group of the transformations, depending on the three arbitrary functions. This is due to the fact that in the (KZK) model the main direction of the wave’s propagation is singled out. The submodels of the (KZK) model.are described by the invariant solutions of the (KZK) equation. We studied essentially distinct, not linked by means of the point transformations, invariant solutions of rank 0 and 1 of this equation. Also we considered the invariant solutions of rank 2 and 3. The invariant solutions of rank 0 and 1 are found either explicitly, or their search is reduced to the solution of the nonlinear integro-differential equations. For example, we obtained the invariant solutions that we called by “Ultrasonic knife” and “Ultrasonic destroyer”. The submodel “Ultrasonic knife” have the following property: at each fixed moment of the time in the field of the existence of the solution near a some plane the pressure increases indefinitely and becomes infinite on this plane. The submodel “Ultrasonic destroyer” contains a countable number of “Ultrasonic knives”. The presence of the arbitrary constants in the integro-differential equations, that determine invariant solutions of rank 1 provides a new opportunities for analytical and numerical study of the boundary value problems for the received submodels, and, thus, for the original (KZK) model. With a help of these invariant solutions we researched a propagation of the intensive acoustic waves (one-dimensional, axisymmetric and planar) for which the acoustic pressure, speed and acceleration of its change, or the acoustic pressure , speed and acceleration of its change in the radial direction, or the acoustic pressure, speed and acceleration of its change in the direction of one of the axes are specified at the initial moment of the time at a fixed point. Under the certain additional conditions, we established the existence and the uniqueness of the solutions of boundary value problems, describing these wave processes. Mechanical relevance of the obtained solutions is as follows: (1) these solutions describe nonlinear and diffraction effects in ultrasonic fields of a special kind, (2) these solutions can be used as a test solutions in the numerical calculations performed in studies of ultrasonic fields generated by powerful emitters. Application of the obtained formula generating the new solutions for the found solutions gives families of the solutions containing three arbitrary functions.     

Invariant submodels of the model of thermal motion of gas in a rarefied space

A model describing the thermal motion of a gas in a rarefied space is investigated. This model can be used in the study of the motion of gas in outer space, and the processes occurring inside the tornado, and the state of the medium behind the shock front of the wave after a very intense explosion. For a given initial pressure distribution, a special choice of mass Lagrange variables leads to a reduced system of differential equations describing this motion, in which the number of independent variables is one less than the original system. This means that there is a stratification of a highly rarefied gas with respect to pressure. Namely, in a strongly rarefied space for each given initial pressure distribution, at each instant of time all gas particles are localized on a two-dimensional surface moving in this space. At each point of this surface, the acceleration vector is collinear with its normal vector. The resulting system admits an infinite Lie transformation group. All significantly various submodels that are invariant with respect to the subgroups of its eight-parameter subgroup generated by the transfer, extension, rotation, and hyperbolic rotation operators (the Lorentz operator) are found. For invariant submodels of rank 1, the basic mechanical characteristics of the gas flow described by them are obtained. Conditions for the existence of these submodels are given. For invariant submodels of rank 2, integral equations describing these submodels are obtained. For some submodels, the problem of describing the gas flow from the initial location of its particles and the distribution of their velocities has been investigated.     

Submodels of model of nonlinear diffusion with non-stationary absorption

We study the model, describing a nonlinear diffusion process (or a heat propagation process) in an inhomogeneous medium with non-stationary absorption (or source). We found tree submodels of the original model of the nonlinear diffusion process (or the heat propagation process), having different symmetry properties. We found all invariant submodels. All essentially distinct invariant solutions describing these invariant submodels are found either explicitly, or their search is reduced to the solution of the nonlinear integral equations. For example, we obtained the invariant solution describing the nonlinear diffusion process (or the heat distribution process) with two fixed "black holes", and the invariant solution describing the nonlinear diffusion process (or the heat distribution process) with the fixed "black hole" and the moving "black hole". The presence of the arbitrary constants in the integral equations, that determine these solutions provides a new opportunities for analytical and numerical study of the boundary value problems for the received submodels, and, thus, for the original model of the nonlinear diffusion process (or the heat distribution process). For the received invariant submodels we are studied diffusion processes (or heat distribution process) for which at the initial moment of the time at a fixed point are specified or a concentration (a temperature) and its gradient, or a concentration (a temperature) and its rate of change. Solving of boundary value problems describing these processes are reduced to the solving of nonlinear integral equations. We are established the existence and uniqueness of solutions of these boundary value problems under some additional conditions. The obtained results can be used to study the diffusion of substances, diffusion of conduction electrons and other particles, diffusion of physical fields, propagation of heat in inhomogeneous medium.     

Invariant formulation of the electromechanical enthalpy function of transversely isotropic piezoelectric materials

Summary Theoretical and numerical aspects of the formulation of electromechanically coupled, transversely isotropic solids are discussed within the framework of the invariant theory. The main goal is the representation of the governing constitutive equations for reversible material behaviour based on an anisotropic electromechanical enthalpy function, which automatically fulfills the requirements of material symmetry. The introduction of a preferred direction in the argument list of the enthalpy function allows the construction of isotropic tensor functions, which reflect the inherent geometrical and physical symmetries of the polarized medium. After presenting the general framework, we consider two important model problems within this setting: i) the linear piezoelectric solid; and ii) the nonlinear electrostriction. A parameter identification of the invariant- and the common coordinate-dependent formulation is performed for both cases. The tensor generators for the stresses, electric displacements and the moduli are derived in detail, and some representative numerical examples are presented.We thank Dipl.-Ing. H. Romanowski for his support and helpful remarks.     

A numerical model for the nonlinear interaction of elastic waves with cracks

Microstructures such as cracks and microfractures play a significant role in the nonlinear interactions of elastic waves, but the precise mechanism of why and how this works is less clear. Here we simulate wave propagation to understand these mechanisms, complementing existing theoretical and experimental works.We implement two models, one of homogeneous nonlinear elasticity and one of perturbations to cracks, and then use these models to improve our understanding of the relative importance of cracks and intrinsic nonlinearity. We find, by modeling the perturbations in the speed of a low-amplitude P-wave caused by the propagation of a large-amplitude S-wave that the nonlinear interactions of P- and S-waves with cracks are significant when the particle motion is aligned with the normal to the crack face, resulting in a larger magnitude crack dilation. This improves our understanding of the relationship between microstructure orientations and nonlinear wave interactions to allow for better characterization of fractures for analyzing processes including earthquake response, reservoir properties, and non-destructive testing.     

地震作用下摩擦耗能减震结构优化

提出了地震作用下摩擦耗能支撑参数优化的一种新的数学模型。它的特点是在给定的几条地震波作用下,在满足框架的规范层间位移角限值要求下,框架各层安装的耗能支撑刚度之和最小,从而实现安装较少的耗能装置而能达到相同的抗震要求。本文利用遗传算法求解该优化问题,并编制了基于遗传算法的几条地震波分别作用下的摩擦耗能框架结构支撑参数弹塑性动力优化分析程序GAOFF。最后给出了一6层框架在三条地震波分别输入下的优化算例,分析结果表明本文方法是有效的。     

Dynamic response of a deepwater riser subjected to combined axial and transverse excitation by the nonlinear coupled model

In offshore engineering long slender risers are simultaneously subjected to both axial and transverse excitations. The axial load is the fluctuating top tension which is induced by the floater’s heave motion, while the transverse excitation comes from environmental loads such as waves. As the time-varying axial load may trigger classical parametric resonance, dynamic analysis of a deepwater riser with combined axial and transverse excitations becomes more complex. In this study, to fully capture the coupling effect between the planar axial and transverse vibrations, the nonlinear coupled equations of a riser’s dynamic motion are formulated and then solved by the central difference method in the time domain. For comparison, numerical simulations are carried out for both linear and nonlinear models. The results show that the transverse displacements predicted by both models are similar to each other when only the random transverse excitation is applied. However, when the combined axial dynamic tension and transverse wave forces are both considered, the linear model underestimates the response because it ignores the coupling effect. Thus the coupled model is more appropriate for deep water. It is also found that the axial excitation can significantly increase the riser’s transverse response and hence the bending stress, especially for cases when the time-varying tension is located at the classical parametric resonance region. Such time-varying effects should be taken into account in fatigue safety assessment.     

On the response of rubbers at high strain rates—I. Simple waves

In this series of papers, we examine the propagation of waves of finite deformation in rubbers through experiments and analysis; in the present paper, Part I, attention is focused on the propagation of one-dimensional waves in strips of natural, latex and synthetic, nitrile rubber. Tensile wave propagation experiments were conducted at high strain rates by holding one end fixed and displacing the other end at a constant velocity. A high-speed video camera was used to monitor the motion and to determine the evolution of strain and particle velocity in rubber strips. Analysis of the response through the theory of finite waves indicated a need for an appropriate constitutive model for rubber; by quantitative matching between the experimental observations and analytical predictions, an appropriate instantaneous elastic response for the rubbers was obtained. This matching process suggested that a simple power-law constitutive model was capable of representing the high strain-rate response for both rubbers used.     

Theoretical and experimental study of energy dissipation in the course of strain localization in iron

Plastic strain localization is studied with the use of a high-sensitivity infrared camera. Plastic strain localization in iron is demonstrated to be accompanied by the emergence of heat waves and their propagation over the sample surface. Constitutive equations that describe the energy balance in the material under plastic strains are derived, and the plastic flow of iron is analyzed. The results of research are compared with the data obtained by infrared scanning. The proposed model of strain localization in the form of soliton-like waves (phase triggering waves) is demonstrated to agree with the kinetics of temperature waves characterizing dissipation inherent in the development of plastic deformation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 153–164, January–February, 2009.     

Group classification of one-dimensional equations of fluids with internal energy depending on the density and the gradient of the density

Group analysis provides a regular procedure for mathematical modeling by classifying differential equations with respect to arbitrary elements. This article presents the group classification of one-dimensional equations of fluids, where the internal energy is a function of the density and the gradient of the density. The equivalence Lie group and the admitted Lie group are provided. The group classification separates all models into 21 different classes according to the admitted Lie group. Invariant solutions of one particular model are obtained.      

Effects of viscous dissipation on thermally developing forced convection in a porous saturated circular tube with an isoflux wall

The viscous dissipation effect on forced convection in a porous saturated circular tube with an isoflux wall is investigated on the basis of the Brinkman flow model. For the thermally developing region, a numerical study is reported while a perturbation analysis is presented to find expressions for the temperature profile and the Nusselt number for the fully developed region. The fully developed Nusselt number found by numerical solution for the developing region is compared with that of asymptotic analysis and a good degree of agreement is observed.     

Solution to the Riemann problem for drift-flux model with modified Chaplygin two-phase flows

In this paper, we concern about the Riemann problem for compressible no-slip drift-flux model which represents a system of quasi-linear partial differential equations derived by averaging the mass and momentum conservation laws with modified Chaplygin two-phase flows. We obtain the exact solution of Riemann problem by elaborately analyzing characteristic fields and discuss the elementary waves namely, shock wave, rarefaction wave and contact discontinuity wave. By employing the equality of pressure and velocity across the middle characteristic field, two nonlinear algebraic equations with two unknowns as gas density ahead and behind the middle wave are formed. The Newton–Raphson method of two variables is applied to find the unknowns with a series of initial data from the literature. Finally, the exact solution for the physical quantities such as gas density, liquid density, velocity, and pressure are illustrated graphically.     

Homogenization of the Prager model in one-dimensional plasticity

We propose a new method for the homogenization of hysteresis models of plasticity. For the one-dimensional wave equation with an elasto-plastic stress-strain relation we derive averaged equations and perform the homogenization limit for stochastic material parameters. This generalizes results of the seminal paper by Franců and Krejčí. Our approach rests on energy methods for partial differential equations and provides short proofs without recurrence to hysteresis operator theory.      

On the asymptotic behavior of a phase-field model for elastic phase transitions

We study the asymptotic behavior of a one-dimensional, dynamical model of solid-solid elastic transitions in which the phase is determined by an order parameter. The system is composed of two coupled evolution equations, the mechanical equation of elasticity which is hyperbolic and a parabolic equation in the order parameter. Due to the strong coupling and the lack of smoothing in the hyperbolic equation, the asymptotic behavior of solutions is difficult to determine using standard methods of gradient-like systems. However, we show that under suitable assumptions all solutions approach the equilibrium set weakly, while the phase field stabilizes strongly.     

刻槽弹侵彻混凝土受力模型研究

利用混凝土材料的动态球形空腔膨胀理论,建立了针对刻槽弹体的低速花瓣形受力模型和高速圆孔形受力模型,并采用这两种模型计算了刻槽弹体侵彻混凝土的侵深。结果表明:当初速低于1 000 m/s时,运用低速花瓣形受力模型计算得出的侵深和实验值的误差小于11%;当初速高于1 000 m/s时,运用高速圆孔形受力模型计算出的侵深和实验值的误差约为20%。综合实验过程和实验误差分析可知,建立的刻槽弹侵彻混凝土受力模型可用于刻槽弹对混凝土的侵彻能力分析。     

On the propagation waves in the theory of thermoelasticity with microtemperatures

The present paper studies the propagation of plane time harmonic waves in an infinite space filled by a thermoelastic material with microtemperatures. It is found that there are seven basic waves traveling with distinct speeds: (a) two transverse elastic waves uncoupled, undamped in time and traveling independently with the speed that is unaffected by the thermal effects; (b) two transverse thermal standing waves decaying exponentially to zero when time tends to infinity and they are unaffected by the elastic deformations; (c) three dilatational waves that are coupled due to the presence of thermal properties of the material. The set of dilatational waves consists of a quasi-elastic longitudinal wave and two quasi-thermal standing waves. The two transverse elastic waves are not subjected to the dispersion, while the other two transverse thermal standing waves and the dilatational waves present the dispersive character. Explicit expressions for all these seven waves are presented. The Rayleigh surface wave propagation problem is addressed and the secular equation is obtained in an explicit form. Numerical computations are performed for a specific model, and the results obtained are depicted graphically.     

Reduced model of a rotor-shaft system using modified SEREP

A modified System Equivalent Reduction Expansion Process (SEREP) is reported to reduce large linear system-equations of motion, which also include gyroscopic effect, internal and external damping. The process is employed on equations of motion written in state-space and utilizes both right and left eigenvectors of the original system unlike only the former in conventional SEREP approach. Numerical example with a finite element model representing multi-disc rotor-shaft-bearing-dynamics shows effectiveness of the process. Nice agreement of Campbell diagrams, FRF plots and unbalance response amplitude between reduced and original systems show usefulness of the reduction process.     

Experimental study and electromechanical model analysis of the nonlinear deformation behavior of IPMC actuators

This paper develops analytical electromechanical formulas to predict the mechanical deformation of ionic polymer–metal composite(IPMC) cantilever actuators under DC excitation voltages. In this research, IPMC samples with Pt and Ag electrodes were manufactured, and the large nonlinear deformation and the effect of curvature on surface electrode resistance of the IPMC samples were investigated experimentally and theoretically. A distributed electrical model was modified for calculating the distribution of voltage along the bending actuator. Then an irreversible thermodynamic model that could predict the curvature of a unit part of an IPMC actuator is combined with the electrical model so that an analytical electromechanical model is developed. The electromechanical model is then validated against the experimental results obtained from Pt-and Ag-IPMC actuators under various excitation voltages. The good agreement between the electromechanical model and the actuators shows that the analytical electromechanical model can accurately describe the large nonlinear quasi-static deflection behavior of IPMC actuators.     

水饱和花岗岩本构模型及地下爆炸波的数值研究

采用广义自洽方法描述介质的细观几何结构,基于混合物理论建立了水饱和岩石含损伤弹塑性本构模型,提出了一个适用于爆炸波数值计算的饱和水体积分数演化方程,也即给出了本构关系中的闭合方程.使用Biot有效应力描述水压对岩石强度的影响,在有效应力空间里建立了水饱和岩石剪切屈服和破坏的强度准则,并引入了剪胀效应对强度的影响因子.对地下封闭核爆炸的力学效应开展了数值模拟,系统地研究了含水花岗岩中地下爆炸波的传播特征以及含水量对爆炸应力波传播规律的影响.