On the stability of some exact solutions to the generalized convection-reaction-diffusion equation

Stability of a set of traveling wave solutions to the convection-reaction-diffusion equation taking into account the effects of memory is studied by means of the qualitative methods and numerical simulation.     

Numerical modelling of surface water waves generated by moving underwater landslide on irregular bottom

The numerical modelling of surface water waves generated by a moving underwater landslide on irregular bottom is considered. The non-linear shallow water model is used with taking into account bottom mobility. The equations are obtained for an underwater landslide movement under the action of gravity force, buoyancy force, friction force and water resistance force. The predictor-corrector scheme [5], preserving the monotonicity of the numerical solution profiles in a linear case, is used on adaptive grids. The scheme is validated on the problem with a known analytical solution. The analysis is done for the dependencies of wave regimes on bottom slope, initial landslide depth, its length and width. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Scattering of a SH-wave by an elastic fiber of nonclassical cross section with an interface crack

The problem of interaction of a plane time-harmonic SH-wave with an elastic fiber of quasi-square or quasi-triangular cross section, when an interface crack is present between an infinite elastic matrix and the fiber, is considered. The modified null-field method taking into account the asymptotic behavior of the solution at crack tips is exploited for obtaining numerical results. The effects of fiber shape, fiber/matrix material combination, debonding (crack size), and direction of wave incidence on the scattering amplitude in the far zone are analyzed. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 44, No. 2, pp. 245–254, March–April, 2008.     

地下结构变形速度在抗爆结构分析中的作用

在地下抗爆结构动力计算中,结构变形速度的作用是非常突出的.考虑结构变形速度的作用,所建立的结构体系运动微分方程可以真实地描述结构振动的实际情况.文中通过一维平面波理论,导出了作用在结构周边上的荷载表达式,给出了地下曲杆结构动力分析的广义变分原理,同时还进行了数值计算结构的对比分析.     

Numerical grids used in a coastal ocean model with breaking wave effects

In coastal ocean modeling, traditional single-block rectangular (Cartesian) grids have been most commonly used for their simplicity. In many cases, these grids may be not well suited (even at very high resolutions) for regions with complicated physical fields, open boundaries, coastlines, and bottom bathymetry. The numerical curvilinear nearly orthogonal/orthogonal, single/multi-block coastline-following grids for the Mediterranean Sea, Monterey Bay and the South China Sea (SCS) are presented. These grids can be used in coastal ocean modeling to enhance model numerical solutions and save computer resources by giving better treatment of regions with high gradients such as areas of complicated coastlines and steep slopes of shelf breaks, complicated bottom topography, open boundaries, and multi-scale physical phenomenon. Grid generation techniques are used to designed these grids. This kind of grids can also easily increase horizontal resolutions in the subregion of the model domain, without increasing the computational expense, with a higher resolution over the entire domain.A three dimensional coastal ocean model with breaking wave effects is also presented and applied. The ocean system is a primitive equation modeling system with grid generation routines and a turbulent closure which is capable of taking surface breaking wave effects into account. The system also includes a grid package which allows model numerical grids to be coupled with the ocean model. The model code is written for multi-block grids, but a single-block grid is used for the South China Sea (SCS). The model with breaking wave effects and a grid of 121 × 121 grid points are used to simulate the winter circulation of the SCS as an example. The model output of the 60-day run shows the observed upwelling locations in the sea surface salinity field.     

Characteristic-based schemes for dispersive waves I. The method of characteristics for smooth solutions

In order to embark on the development of numerical schemes for stiff problems, we have studied a model of relaxing heat flow. To isolate those errors unavoidably associated with discretization, a method of characteristics is developed, containing three free parameters depending on the stiffness ratio. It is shown that such “decoupled” schemes do not take into account the interaction between the wave families and hence result in incorrect wave speeds. We also demonstrate that schemes can differ by up to two orders of magnitude in their rms errors even while maintaining second-order accuracy. We show that no method of characteristics solution can be better than second-order accurate. Next, we develop “coupled” schemes which account for the interactions, and here we obtain two additional free parameters. We demonstrate how coupling of the two wave families can be introduced in simple ways and how the results are greatly enhanced by this coupling. Finally, numerical results for several decoupled and coupled schemes are presented, and we observe that dispersion relationships can be a very useful qualitative tool for analysis of numerical algorithms for dispersive waves. © 1993 John Wiley & Sons, Inc.     

Scattering of an incident acoustic wave by an elastic sphere in a shallow water

The scattering of acoustic waves by an elastic sphere in a shallow ocean wave guide is investigated taking into account the shear waves which can exist in addition to compressional waves in scatterers of solid material. Expressions for the scattered waves are given. Numerical values for a quantity called the farfield form function for various depth are presented in graphical forms.      

Localized solutions for the finite difference semi-discretization of the wave equation

We study the propagation properties of the solutions of the finite difference space semi-discrete wave equation on a uniform grid of the whole Euclidean space. We provide a construction of high frequency wave packets that propagate along the corresponding bi-characteristic rays of Geometric Optics with a group velocity arbitrarily close to zero. Our analysis is motivated by control theoretical issues. In particular, the continuous wave equation has the so-called observability property: for a sufficiently large time, the total energy of its solutions can be estimated in terms of the energy concentrated in the exterior of a compact set. This fails to be true, uniformly on the mesh-size parameter, for the semi-discrete schemes and the observability constant blows-up at an arbitrarily large polynomial order. Our contribution consists in providing a rigorous derivation of those wave packets and in analyzing their behavior near that ray, by taking into account the subtle added dispersive effects that the numerical scheme introduces.     

Multivariate tail conditional expectation for elliptical distributions

In this paper we introduce a novel type of a multivariate tail conditional expectation (MTCE) risk measure and explore its properties. We derive an explicit closed-form expression for this risk measure for the elliptical family of distributions taking into account its variance–covariance dependency structure. As a special case we consider the normal, Student-t and Laplace distributions, important and popular in actuarial science and finance. The motivation behind taking the multivariate TCE for the elliptical family comes from the fact that unlike the traditional tail conditional expectation, the MTCE measure takes into account the covariation between dependent risks, which is the case when we are dealing with real data of losses. We illustrate our results using numerical examples in the case of normal and Student-t distributions.     

Numerical model for vibration damping resulting from the first-order phase transformations

A numerical model is constructed for modelling macroscale damping effects induced by the first-order martensite phase transformations in a shape memory alloy rod. The model is constructed on the basis of the modified Landau–Ginzburg theory that couples nonlinear mechanical and thermal fields. The free energy function for the model is constructed as a double well function at low temperature, such that the external energy can be absorbed during the phase transformation and converted into thermal form. The Chebyshev spectral methods are employed together with backward differentiation for the numerical analysis of the problem. Computational experiments performed for different vibration energies demonstrate the importance of taking into account damping effects induced by phase transformations.     

Analysis of damage and fracture mechanisms on the micro-scale

The paper deals with generalized damage and failure criteria taking into account different mechanisms on the micro-scale. Some micro-mechanical numerical results are discussed and compared with the macroscopic analysis based on the continuum damage model. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Wave Propagation in a beam with random material properties

Modern composite materials, e.g., carbon fibre reinforced plastics (CFRP), exhibit a complex micro structure due to their fabrication process. The latter, being usually omitted in mechanical models through the homogenization of elastic properties, has a strong influence on the propagation of ultrasonic guided waves [1, 2]. Though it is possible to model the wave phenomena deterministically, taking into account a realistic distribution of fibres and polymer matrix, it is desirable to develop an improved model for the finite element analysis (FEM), which consider the stochastic properties in a more general way. In the current work, an approach for the simulation of waves in a isotropic beam with random material properties is presented. For the numerical computations with the FEM the Young's modulus was discretized by the Karhunen-Loève Expansion (KLE). Numerical investigations on the excited and propagating guided waves are presented. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Augmented Galerkin schemes for the numerical solution of scattering by small obstacles

We are interested in the problem of a bidimensional acoustic wave propagation in a medium including a small obstacle with homogeneous Dirichlet boundary condition. We present and analyse a numerical scheme suitable for finite elements that does not suffer from numerical locking, and takes the presence of the small obstacle into account. It is based on the fictitious domain method combined with matched asymptotic expansions.     

A numerical method for calculating solitons of the nonlinear Schödinger equation in the axially symmetric case

A numerical method is proposed for determination of the eigenfunctions and eigenvalues of the nonlinear Schrödinger equation in the axially symmetric case. Optical solitons interpreted in the physical sense are found for various values of the nonlinearity coefficient by means of the developed method. As has previously been shown by other authors, such solitons are unstable under small perturbations of their shape. Since the considered problem finds numerous applications, methods providing for soliton stabilization are widely discussed in the literature. One of these methods involves strong modulation of the medium nonlinearity or even the reversal of the nonlinearity sign, which necessitates taking into account the wave reflected from irregularities and analyzing additionally the applicability of the mathematical model. We show that, theoretically, it is possible to stabilize a soliton via weak modulation of the cubic-nonlinearity coefficient. Such modulation ensures alternation of the length of nonlinear layers and enables one to increase the path length by a factor of 70 without a beam collapse.     

弹塑性杆在刚性块轴向撞击下的动力屈曲

基于能量原理,对弹塑性杆在刚性块轴向撞击下的动力屈曲问题进行了讨论．用特征线法分析了刚性块轴向撞击弹塑性直杆时应力波传播的过程．考虑了弹塑性应力波传播对屈曲的影响,建立了该问题横向扰动方程．用幂级数解法,理论上给出了该问题的级数解．分析解的性质,得到了发生屈曲时的临界条件．通过理论分析和数值计算,得到了临界速度与冲击质量、临界长度及线性强化模量间的关系．     

Natural and parametric vibrations of a cylindrical shell-filler system with tangential interaction taken into account

Conclusions 1. The spectrum of the eigenfrequencies and dynamical instability regions of a shell-filler system breaks up into two infinite spectra for each mode of wave formation; the first of them is determined by the shear modulus of the filler material, and the second — by the volume deformation modulus. The second spectrum is absent for incompressible fillers.2. It has been shown that taking tangential interactions into account has a strong effect on the arrangement and width of the dynamical instability regions belonging to the first spectrum and hardly changes the arrangement and widths of the regions of the second spectrum.3. As a result of the limiting transition as the frequency of the driving force tends to zero, expressions are obtained from the formulas of this paper for calculating the static stability of a shell with an incompressible filler. The numerical results obtained for this case correspond to results given in the literature.Institute of Polymer Mechanics, Academy of Sciences of the Latvian SSR, Riga. P. Stuchki Latvian State University, Riga. Translated from Mekhanika Polimerov, No. 3, pp. 503–509, May–June, 1977.     

Modeling and numerical simulations of rate and temperature-dependent damage and fracture in aluminum alloys

The paper deals with a continuum damage model and corresponding numerical simulations of the inelastic deformation and failure behavior of aluminum alloys. The focus is on generalized plasticity, damage and fracture approaches taking into account different tiaxialities as well as rate and temperature-dependence. Furthermore, the special intent of the present work are dynamic applications. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Méthodes d'optimisation pour l'approche de problèmes aux limites non linéaires mixtes elliptiques hyperboliques

Functional and numerical results are presented for a significative nonlinear boundary value problem of mixed type elliptic hyperbolic: a transonic flows model taking into account an entropy condition. Relation between continuous problem and its discretization in finite dimension is stated.     

On the application of the nonstationary form of the Tappert equation as an artificial boundary condition

In this paper, a nonstationary analog of the range refraction parabolic equation is derived. A new approach to the derivation of Tappert’s operator asymptotic formula with the use of noncommutative analysis is presented. The obtained nonstationary equation is proposed as an artificial boundary condition for the wave equation in underwater acoustics. This form of artificial boundary condition has low computational cost and systematically takes into account variations of sound speed. This is confirmed by various numerical experiments, including propagation of normal modes and wave fields produced by point source.