Sharif Digital Repository

Hydrodynamic instabilities at the interface of stratified shear layers could occur in various modes and have an important role in the mixing process. In this work, the linear stability analysis in the temporal framework is used to study the stability characteristics of a particle-laden stratified two-layer flow for two different background density profiles: smooth (hyperbolic tangent) and piecewise linear. The effect of parameters, such as bed slope, viscosity, and particle size, on the stability is also considered. The pseudospectral collocation method employing Chebyshev polynomials is used to solve two coupled eigenvalue equations. Based on the results, there are some differences in the... 

Various researchers have contributed to the identification of the mass and stiffness matrices of two dimensional (2-D) shear building structural models for a given set of vibratory frequencies. The suggested methods are based on the specific characteristics of the Jacobi matrices, i.e., symmetric, tri-diagonal and semi-positive definite matrices. However, in case of three dimensional (3-D) structural models, those methods are no longer applicable, since their stiffness matrices are not tri-diagonal. In this paper the inverse problem for a special class of vibratory structural systems, i.e., 3-D shear building models, is investigated. A practical algorithm is proposed for solving the inverse... 

We use the Fokker-Planck equation and its moment equations to study the collective behavior of interacting particles in unsteady one-dimensional flows. Particles interact according to a longrange attractive and a short-range repulsive potential field known as Morse potential. We assume Stokesian drag force between particles and their carrier fluid and find analytic single-peaked traveling solutions for the spatial density of particles in the catastrophic phase. In steady flow conditions the streaming velocity of particles is identical to their carrier fluid, but we show that particle streaming is asynchronous with an unsteady carrier fluid. Using linear perturbation analysis, the stability... 

Elastic stability of shell structures under certain loading conditions is characterized by a dramatically unstable postbuckling behavior. The presence of simultaneous 'competing' buckling modes (corresponding to the same critical buckling load) is understood to be largely responsible for such behavior. In this paper, within the framework of linear bifurcation eigenvalue analysis and Donnell shallow shell theory, the presence of simultaneous buckling modes in axially compressed isotropic cones is determined using the semi-analytical method of Galerkin. The results are presented in the plane of the dimensionless reciprocal meridional and circumferential buckling half wavelengths, and are... 

During the year, due to weather conditions, the temperature fluctuations at surface level cause problems in underground pipes as a result of freezing water. One of the best prevention strategies is the use of polyurethane floor insulation for keeping the temperature of clay above zero degrees Celsius. In this study to calculate the minimum thickness of polyurethane insulation layer, the differential equation of energy is solved based on principle of separation of variables using imaginary eigenvalues for consistency with the temperature distribution in multi-layer consist of asphalt, gravel and polyurethane with finite thickness and clay as a semiinfinite medium with periodic thermal... 

In this article, OPSEM (Orthonormal Polynomial Series Expansion Method) is developed as a new computational approach for the evaluation of thin beams of variable thickness transverse vibration. Capability of the OPSEM in assessing the free vibration frequencies and mode shapes of an Euler-Bernoulli beam with varying thickness is discussed. Multispan continuous beams with various classical boundary conditions are included. Contribution of BOPs (Basic Orthonormal Polynomials) in capturing the beam vibrations is also illustrated in numerical examples to give a quantitative measure of convergence rate. Furthermore, OPSEM is adopted for the forced vibration of a thin beam caused by a moving mass.... 

Recent researches on automatic skill acquisition in reinforcement learning have focused on subgoal discovery methods. Among them, algorithms based on graph partitioning have achieved higher performance. In this paper, we propose a new automatic skill acquisition framework based on graph partitioning approach. The main steps of this framework are identifying subgoals and discovering useful skills. We propose two subgoal discovery algorithms, which use spectral analysis on the transition graph of the learning agent. The first proposed algorithm, incorporates k′-means algorithm with spectral clustering. In the second algorithm, eigenvector centrality measure is utilized and options are... 

Stability analysis of a horizontal cantilevered pipe conveying fluid with an inclined terminal nozzle is considered in this paper. The pipe is modelled as a cantilevered Euler-Bernoulli beam, and the flow-induced inertia, Coriolis and centrifugal forces along the pipe as well as the follower force induced by the jet-flow are taken into account. The governing equations of the coupled bending-torsional vibrations of the pipe are obtained using extended Hamilton's principle and are then discretized via the Galerkin method. The resulting eigenvalue problem is then solved, and several cases are examined to determine the effect of nozzle inclination angle, nozzle aspect ratio, mass ratio and... 

The interaction of the dynamics of doubly fed wind generators with the electromechanical mode of nearby synchronous generators (SGs) can affect the small signal stability of power systems with high penetration levels of wind power. In this study, a novel approach is developed to investigate these interactions and their impact on the damping of power system oscillations. In this approach it is not necessary to model the dynamic behaviour of system SGs and only the frequencies of system oscillations are important. This approach is based on the sensitivity of SGs electromechanical eigenvalue with respect to variations in the Jacobian matrix of power system. By applying this approach to a test... 

Introduction of large amounts of new wind generation can affect the small signal stability of power systems with three mechanisms: displacing synchronous generators (SGs); reducing SGs power generation; and the dynamics of wind farms (WFs) interacting with the electromechanical mode of SGs. In this paper a novel approach is developed to investigate the impact of the latter mechanism on existing power systems oscillations. In this approach, the dynamic behavior of grid connected WFs is studied independent of the dynamic behavior of system SGs. This approach helps to identify the conditions in which the dynamics of WFs may interact with the electromechanical mode of SGs. Also it helps to... 

Because of the complicated buckling characteristics of Y-shaped bracings, calculation of their strength is beyond routine engineering procedures. A method based on slope-deflection equations incorporating stability functions has been used for computation of buckling eigenvalues and effective length factors. Lateral strength is computed based on the least buckling strength of bracing members. This study is limited to bracings with similar sections and fixed end connections resisting out-of-plane rotation. Lateral strengths predicted by the analytical method for two cases are compared with experimental results on full-scale specimens. It is shown that the proposed analytical method predicts... 

This paper presents a scheme to enumerate the incident waves impinging on a high dimensional uniform linear array using relatively few samples. The approach is based on Minimum Description Length (MDL) criteria and statistical properties of eigenvalues of the Sample Covariance Matrix (SCM). We assume that several models, with each model representing a certain number of sources, will compete and MDL criterion will select the best model with the minimum model complexity and maximum model decision. Statistics of noise eigenvalue of SCM can be approximated by the distributional properties of the eigenvalues given by Marcenko-Pastur distribution in the signal-free SCM. In this paper we use random... 

This letter presents new sufficient conditions for robust Hurwitz stability of interval matrices. The proposed conditions are based on two approaches: (i) finding a common Lyapunov matrix for the interval family and (ii) converting the robust stability problem into a robust non-singularity problem using Kronecker operations. The main contribution of the letter is to derive accurate and computationally simple optimal estimates of the robustness margin and spectral bound of general interval matrices. The evaluation of the condition relies on the solutions of linear matrix inequalities (LMIs) and eigenvalue problems, both of which are solved very efficiently. The improvements gained by using... 

Thermally developing mixed electroosmotically and pressure-driven flow in a parallel plate microchannel with a step change in wall temperature is considered in the framework of an extended Graetz problem. Both Joule heating and viscous dissipation effects are taken into consideration. Expressions for the dimensionless temperature and Nusselt number in the form of infinite series are presented. In general, the associated eigenvalue problem is solved numerically. Nevertheless, an analytical solution is also presented for axial locations close to the entrance. Comparisons are made between the present results and those obtained by approximating the electroosmotic velocity with the... 

Shape analysis plays an important role in many medical imaging studies. One of the recent shape analysis methods uses the Laplace Beltrami eigenvalues which is also used in this paper for global shape comparison of hippocampus of normal subjects and epileptic patients. Popularity of the Laplace Beltrami operator in this field is due to its isometry-invariance which avoids pre-processing steps like mapping, registration, and alignment. In addition, it is capable of revealing fine details in shapes that makes this method a good choice for deformation detecting purposes like epilepsy diagnosis. To examine capability of the proposed method, statistical analysis and two ways of classification,... 

To solve the controversy, regarding the existence of an analytic solution to the 1-D Ising model with nearest-neighbor (NN) and next-nearest-neighbor (NNN) interactions in the presence of a magnetic field, we apply the transfer matrix method to solve the 1-D Ising model in the presence of a magnetic field, taking both NN and NNN interactions into account. We show that it is possible to write a transfer matrix only if the number of sites is even. Even in such a case, it is impossible to diagonalize the transfer matrix in an analytic form. Therefore, we employ a numerical method to obtain the eigenvalues of the transfer matrix. Moreover, the heat capacity, magnetization, and magnetic... 

Controlling networked structures has many applications in science and engineering. In this paper, we consider the problem of pinning control (pinning the dynamics into the reference state), and optimally placing the driver nodes, i.e., the nodes to which the control signal is fed. Considering the local controllability concept, a metric based on the eigenvalues of the Laplacian matrix is taken into account as a measure of controllability. We show that the proposed optimal placement strategy considerably outperforms heuristic methods including choosing hub nodes with high degree or betweenness centrality as drivers. We also study properties of optimal drivers in terms of various centrality... 

In this paper, a robotic manipulator modelled as a cantilever rotating Euler-Bernoulli beam is considered. Control objective is achieving a desired angular rotation of the manipulator tip while its lateral vibration is suppressed. An external driving torque is the control input of the system. Two dynamic transfer functions are derived to describe beam tip motion and angular rotation in terms of the desired angular rotation. After state-space representation of the problem, an observer is designed to estimate state variables of the system. Then, a feedback control is designed for both regulation and tracking objectives. Eigenvalues are chosen such that an appropriate response is achieved while... 

It is has been long known that the curved space in the presence of gravitation can be described as a non-homogeneous anisotropic medium in flat geometry with different constitutive equations. In this article, we show that the eigenpolarizations of such medium can be exactly solved, leading to a pseudo-isotropic description of curved vacuum with two refractive index eigenvalues having opposite signs, which correspond to forward and backward travel in time. We conclude that for a rotating universe, time-reversal symmetry is broken. We also demonstrate the applicability of this method to Schwarzschild metric and derive exact forms of refractive index. We derive the subtle optical anisotropy of... 

In this paper, we consider the problem of spectrum sensing by using multiple antenna in cognitive radios when the noise and the primary user signal are assumed as independent complex zero-mean Gaussian random signals. The optimal multiple antenna spectrum sensing detector needs to know the channel gains, noise variance, and primary user signal variance. In practice some or all of these parameters may be unknown, so we derive the Generalized Likelihood Ratio (GLR) detectors under these circumstances. The proposed GLR detector, in which all the parameters are unknown, is a blind and invariant detector with a low computational complexity. We also analytically compute the missed detection and...