Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.

Unsteady Model of Drop Marangoni Migration in Microgravity

The experiments of drop Marangoni migration have been performed by the drop shift facility of short period of 4.5 s, and the drop accelerates gradually to an asymptotic velocity during the free fall. The unsteady and axisymmetric model is developed to study the drop migration for the case of moderate Reynolds number Re = O(1), and the results are compared with the experimental ones in the present paper. Both numerical and experimental results show that the migration velocity for moderate Reynolds number is several times smaller than that given by the linear YGB theory.

Non-Linear Wave-Induced Transient Response of Soil Around a Trenched Pipeline

Submarine pipelines are always trenched within a seabed for reducing wave loads and thereby enhancing their stability. Based on Biot’s poroelastic theory, a two-dimensional finite element model is developed to investigate non-linear wave-induced responses of soil around a trenched pipeline, which is verified with the flume test results by Sudhan et al. [Sudhan, C.M., Sundar, V., Rao, S.N., 2002. Wave induced forces around buried pipeline. Ocean Engineering, 29, 533–544] and Turcotte et al. [Turcotte, B.R., Liu, P.L.F., Kulhawy, F.H., 1984. Laboratory evaluation of wave tank parameters for wave-sediment interaction. Joseph H. Defree Hydraulic Laboratory Report 84-1, School of Civil and Environmental Engineering, Cornell University]. Non-linear wave-induced transient pore pressure around pipeline at various phases of wave loading is examined firstly. Unlike most previous investigations, in which only a single sediment layer and linear wave loading were concerned, in this study, the influences of the non-linearity of wave loading, the physical properties of backfill materials and the geometry profile of trenches on the excess pore pressures within the soil around pipeline, respectively, were explored, taking into account the in situ conditions of buried pipeline in the shallow ocean zones. Based on the parametric study, it is concluded that the shear modulus and permeability of backfill soils significantly affect the wave-induced excess pore pressures around trenched pipeline, and that the effect of wave non-linearity becomes more pronounced and comparable with that of trench depth, especially at high wave steepness in shallow water.

Linear Stability Analysis of Convection in Two-Layer System with an Evaporating Vapor-Liquid Interface

Classical theories have successfully provided an explanation for convection in a liquid layer heated from below without evaporation. However, these theories are inadequate to account for the convective instabilities in an evaporating liquid layer, especially in the case when it is cooled from below. In the present paper, we study the onset of Marangoni convection in a liquid layer being overlain by a vapor layer.A new two-sided model is put forward instead of the one-sided model in previous studies. Marangoni-Bénard instabilities in evaporating liquid thin layers are investigated with a linear instability analysis. We define a new evaporation Biot number, which is different from that in previous studies and discuss the influences of reference evaporating velocity and evaporation Biot number on the vapor-liquid system. At the end, we explain why the instability occurs even when an evaporating liquid layer is cooled from below.

Bayesian model selection for time series using Markov chain Monte Carlo

We present a stochastic simulation technique for subset selection in time series models, based on the use of indicator variables with the Gibbs sampler within a hierarchical Bayesian framework. As an example, the method is applied to the selection of subset linear AR models, in which only significant lags are included. Joint sampling of the indicators and parameters is found to speed convergence. We discuss the possibility of model mixing where the model is not well determined by the data, and the extension of the approach to include non-linear model terms.

Model for solidification cracking in low alloy steel weld metals

Data on the occurrence of solidification cracking in low alloy steel welds have been analysed using a classification neural network based on a Bayesian framework. It has thereby been possible to express quantitatively the effect of variables such as the chemical composition, welding conditions, and weld geometry, on the tendency for solidification cracking during solidification. The ability of the network to express the relationship in a suitably non-linear form is shown to be vital in reproducing known experimental phenomena. © 1996 The Institute of Materials.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

A consistent model of isolated force-free magnetic arches

The magnetic flux tube concentrating strong magnetic field is the basic configuration of magneticfield in the solar atmosphere. In the present paper, the equilibrium of isolated magnetic flux tube inthe solar atmosphere is discussed. In the viewpoint of mathematics, the boundary condition is nonlinearand the position of boundary needs to be determined by the physical condition although the equation ofmagnetic potential is linear for the linear force-free field. Analytical solutions to the arches of bothuniform circular cross-section and non-uniform cross section have been obtained. The results show thatthe nonlinear problem may have or not have any solution according to different azimuthal components of the magnetic field; the number of solutions to the nonlinear problem is four at most, and two in some cases. In the present paper, the analytical solutions to the approximations of both fat and slender arches are given in detail, and the general features of magnetic arch structure are shown.

A Supply Function Competition Model for the Spanish Wholesale Electricity Market

We model the Spanish wholesale market as a multiplant linear supply function competition model. According to the theory, the larger generators should have supply curves for each plant which are to the left of the supply curves of plants owned by smaller generators. We test this prediction for fuel plants using data from the Spanish Market Operator (OMEL) from May 2001 to December 2003. Our results indicate that the prediction of the model holds.

Cyclical Features of Uzawa-Lucas Endogenous Growth Model

This paper analyzes the cyclical properties of a generalized version of Uzawa-Lucas endogenous growth model. We study the dynamic features of different cyclical components of this model characterized by a variety of decomposition methods. The decomposition methods considered can be classified in two groups. On the one hand, we consider three statistical filters: the Hodrick-Prescott filter, the Baxter-King filter and Gonzalo-Granger decomposition. On the other hand, we use four model-based decomposition methods. The latter decomposition procedures share the property that the cyclical components obtained by these methods preserve the log-linear approximation of the Euler-equation restrictions imposed by the agent’s intertemporal optimization problem. The paper shows that both model dynamics and model performance substantially vary across decomposition methods. A parallel exercise is carried out with a standard real business cycle model. The results should help researchers to better understand the performance of Uzawa-Lucas model in relation to standard business cycle models under alternative definitions of the business cycle.

A rational fraction polynomials model to study vertical dynamic wheel-rail interaction

This paper presents a model designed to study vertical interactions between wheel and rail when the wheel moves over a rail welding. The model focuses on the spatial domain, and is drawn up in a simple fashion from track receptances. The paper obtains the receptances from a full track model in the frequency domain already developed by the authors, which includes deformation of the rail section and propagation of bending, elongation and torsional waves along an infinite track. Transformation between domains was secured by applying a modified rational fraction polynomials method. This obtains a track model with very few degrees of freedom, and thus with minimum time consumption for integration, with a good match to the original model over a sufficiently broad range of frequencies. Wheel-rail interaction is modelled on a non-linear Hertzian spring, and consideration is given to parametric excitation caused by the wheel moving over a sleeper, since this is a moving wheel model and not a moving irregularity model. The model is used to study the dynamic loads and displacements emerging at the wheel-rail contact passing over a welding defect at different speeds.

Asymptotic Hyperstability of a Class of Linear Systems under Impulsive Controls Subject to an Integral Popovian Constraint

This paper is focused on the study of the important property of the asymptotic hyperstability of a class of continuous-time dynamic systems. The presence of a parallel connection of a strictly stable subsystem to an asymptotically hyperstable one in the feed-forward loop is allowed while it has also admitted the generation of a finite or infinite number of impulsive control actions which can be combined with a general form of nonimpulsive controls. The asymptotic hyperstability property is guaranteed under a set of sufficiency-type conditions for the impulsive controls.

Downscaling of surface moisture flux and precipitation in the Ebro Valley (Spain) using analogues and analogues followed by random forests and multiple linear regression

In this paper, reanalysis fields from the ECMWF have been statistically downscaled to predict from large-scale atmospheric fields, surface moisture flux and daily precipitation at two observatories (Zaragoza and Tortosa, Ebro Valley, Spain) during the 1961-2001 period. Three types of downscaling models have been built: (i) analogues, (ii) analogues followed by random forests and (iii) analogues followed by multiple linear regression. The inputs consist of data (predictor fields) taken from the ERA-40 reanalysis. The predicted fields are precipitation and surface moisture flux as measured at the two observatories. With the aim to reduce the dimensionality of the problem, the ERA-40 fields have been decomposed using empirical orthogonal functions. Available daily data has been divided into two parts: a training period used to find a group of about 300 analogues to build the downscaling model (1961-1996) and a test period (19972001), where models' performance has been assessed using independent data. In the case of surface moisture flux, the models based on analogues followed by random forests do not clearly outperform those built on analogues plus multiple linear regression, while simple averages calculated from the nearest analogues found in the training period, yielded only slightly worse results. In the case of precipitation, the three types of model performed equally. These results suggest that most of the models' downscaling capabilities can be attributed to the analogues-calculation stage.

Linear Chromatic Adaptation Transform Based on Delaunay Triangulation

Computer vision algorithms that use color information require color constant images to operate correctly. Color constancy of the images is usually achieved in two steps: first the illuminant is detected and then image is transformed with the chromatic adaptation transform ( CAT). Existing CAT methods use a single transformation matrix for all the colors of the input image. The method proposed in this paper requires multiple corresponding color pairs between source and target illuminants given by patches of the Macbeth color checker. It uses Delaunay triangulation to divide the color gamut of the input image into small triangles. Each color of the input image is associated with the triangle containing the color point and transformed with a full linear model associated with the triangle. Full linear model is used because diagonal models are known to be inaccurate if channel color matching functions do not have narrow peaks. Objective evaluation showed that the proposed method outperforms existing CAT methods by more than 21%; that is, it performs statistically significantly better than other existing methods.