The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

Assessment of the residual fatigue strength in RC beams

In conventional analysis and design procedures of reinforced concrete structures, the ability of concrete to resist tension is neglected. Under cyclic loading, the tension-softening behavior of concrete influences its residual strength and subsequent crack propagation. The stability and the residual strength of a cracked reinforced concrete member under fatigue loading, depends on a number of factors such as, reinforcement ratio, specimen size, grade of concrete, and the fracture properties, and also on the tension-softening behavior of concrete. In the present work, a method is proposed to assess the residual strength of a reinforced concrete member subjected to cyclic loading. The crack extension resistance based approach is used for determining the condition for unstable crack propagation. Three different idealization of tension softening models are considered to study the effect of post-peak response of concrete. The effect of reinforcement is modeled as a closing force counteracting the effect of crack opening produced by the external moment. The effect of reinforcement percentage and specimen size on the failure of reinforced beams is studied. Finally, the residual strength of the beams are computed by including the softening behavior of concrete.

Fatigue crack propagation in rocker and roller–rocker bearings of railway steel bridges

In this work, a method is proposed for rolling contact fatigue crack propagation analysis using contact and fracture theories in conjunction with fatigue laws. The proposed method is used in the fatigue analysis of rocker and roller–rocker bearings of a railway open web girder bridge which is instrumented with strain gages. Using a contact algorithm based on the minimum energy principle for bodies in rolling contact with dry friction, the normal and tangential pressure distribution are computed. It is seen that the most critical location of a crack in bearings is at a point very close to the contact region, as expected.

Fatigue crack propagation in rocker and roller-rocker bearings of railway steel bridges

In this work, a method is proposed for rolling contact fatigue crack propagation analysis using contact and fracture theories in conjunction with fatigue laws. The proposed method is used in the fatigue analysis of rocker and roller-rocker bearings of a railway open web girder bridge which is instrumented with strain gages. Using a contact algorithm based on the minimum energy principle for bodies in rolling contact with dry friction, the normal and tangential pressure distribution are computed. It is seen that the most critical location of a crack in bearings is at a point very close to the contact region, as expected. (C) 2010 Elsevier Ltd. All rights reserved.

Pickup and delivery problem using metaheuristics techniques

Dial-a-ride problem (DARP) is an optimization problem which deals with the minimization of the cost of the provided service where the customers are provided a door-to-door service based on their requests. This optimization model presented in earlier studies, is considered in this study. Due to the non-linear nature of the objective function the traditional optimization methods are plagued with the problem of converging to a local minima. To overcome this pitfall we use metaheuristics namely Simulated Annealing (SA), Particle Swarm Optimization (PSO), Genetic Algorithm (GA) and Artificial Immune System (AIS). From the results obtained, we conclude that Artificial Immune System method effectively tackles this optimization problem by providing us with optimal solutions. Crown Copyright (C) 2011 Published by Elsevier Ltd. All rights reserved.

Studies on characteristics of lightning generated currents in an interconnected lightning protection system

An in-depth knowledge about the characteristics of lightning generated currents will facilitate evaluation of the interception efficacy of lightning protection systems. In addition, it would aid in extraction of valuable statistics (from measured current data) on local lightning parameters. Incidentally, present day knowledge on characteristics of lightning induced current in typical lightning protection systems is rather limited. This is particularly true with closely interconnected protection systems, like the one employed in Indian Satellite Launch Pad-II. This system is taken as a specific example in the present study. Various aspects suggest that theoretical modelling would be the best possible approach for the intended work. From the survey of pertinent literature, it is concluded that electromagnetic modelling of lightning return-stroke with current source at the channel base is best suited for this study. Numerical electromagnetic code was used for the required electromagnetic field solution and Fourier transform techniques were employed for computing time-domain results. A validation for the numerical modelling is provided by laborator experiments on a reduced scale model of the system. Apart from ascertaining the influence of various parameters, salient characteristics of tower base currents for different kinds of events are deduced. This knowledge can be used in identifying the type of event, as well as its approximate location. A method for estimation of injected stroke current has also been proposed.

Regularized numerical integration of multibody dynamics with the generalized alpha method*

This paper discusses the consistent regularization property of the generalized α method when applied as an integrator to an initial value high index and singular differential-algebraic equation model of a multibody system. The regularization comes from within the discretization itself and the discretization remains consistent over the range of values the regularization parameter may take. The regularization involves increase of the smallest singular values of the ill-conditioned Jacobian of the discretization and is different from Baumgarte and similar techniques which tend to be inconsistent for poor choice of regularization parameter. This regularization also helps where pre-conditioning the Jacobian by scaling is of limited effect, for example, when the scleronomic constraints contain multiple closed loops or singular configuration or when high index path constraints are present. The feed-forward control in Kane's equation models is additionally considered in the numerical examples to illustrate the effect of regularization. The discretization presented in this work is adopted to the first order DAE system (unlike the original method which is intended for second order systems) for its A-stability and same order of accuracy for positions and velocities.

A method for the experimental evaluation of the acoustic characteristics of an engine exhaust system in the presence of mean flow

The paper deals with an exact analysis of standing waves in an impedance tube with mean flow. A method is offered for the experimental evaluation of the various wave parameters. Navier–Stokes equations have been solved for evaluating the volume velocity taking into account mean flow, viscosity, etc. The engine exhaust system has been characterized as an acoustic source with an acoustic pressure and internal impedance. A method is suggested for the evaluation of these hypothetical parameters using the exhaust pipe as an impedance tube.Subject Classification: [43]85.20; [43]20.40.

Numerical study of heat transfer in laminar film boiling by the finite-difference method

The finite-difference form of the basic conservation equations in laminar film boiling have been solved by the false-transient method. By a judicious choice of the coordinate system the vapour-liquid interface is fitted to the grid system. Central differencing is used for diffusion terms, upwind differencing for convection terms, and explicit differencing for transient terms. Since an explicit method is used the time step used in the false-transient method is constrained by numerical instability. In the present problem the limits on the time step are imposed by conditions in the vapour region. On the other hand the rate of convergence of finite-difference equations is dependent on the conditions in the liquid region. The rate of convergence was accelerated by using the over-relaxation technique in the liquid region. The results obtained compare well with previous work and experimental data available in the literature.

Landau-Zener problem with waiting at the minimum gap and related quench dynamics of a many-body system

We discuss a technique for solving the Landau-Zener (LZ) problem of finding the probability of excitation in a two-level system. The idea of time reversal for the Schrodinger equation is employed to obtain the state reached at the final time and hence the excitation probability. Using this method, which can reproduce the well-known expression for the LZ transition probability, we solve a variant of the LZ problem, which involves waiting at the minimum gap for a time t(w); we find an exact expression for the excitation probability as a function of t(w). We provide numerical results to support our analytical expressions. We then discuss the problem of waiting at the quantum critical point of a many-body system and calculate the residual energy generated by the time-dependent Hamiltonian. Finally, we discuss possible experimental realizations of this work.

Effect of static deflection on natural frequency of non-linear spring mass system by direct linearization method

An exact expression for the frequency of a non-linear cubic spring mass system is obtained considering the effect of static deflection. An alternative expression for the approximate frequency is also obtained by the direct linearization procedure; it is shown that this is very accurate as compared with the exact method. This approximate frequency equation is used to explain a “dual behaviour” of the frequency amplitude curves.

On an asymptotic method of Krylov-Bogoliubov for overdamped nonlinear systems

A general asymptotic method based on the work of Krylov-Bogoliubov is developed to obtain the response of nonlinear over damped systems. A second-order system with both roots real is treated first and the method is then extended to higher-order systems. Two illustrative examples show good agreement with results obtained by numerical integration.

The influence of Zr/Ti content on the morphotropic phase boundary in the PZT-PZN system

A quantitative structural investigation was carried out on (1-y)PbZrxTi1-xO3-yPbZn(1/3)Nb(2/3)O(3) where y=0.1 and 0.2 ((1-y)PZT-yPZN). High resolution XRD data have been used for quantitative phase analysis. The nominal compositions were prepared by a two-step low temperature calcining solid-state method. The sintered samples show an average grain size of 1-2 mu m. It is demonstrated that the increase in the concentration of PZN leads to the shift of the morphotropic phase boundary (MPB) of PZT towards the PbZrO3 end member. In the present work, an effort has been made to quantitatively determine the MPB phase contents and to regain the coexistence of tetragonal and monoclinic phases by varying the value of x(i.e. Zr/Ti ratio). The width of the MPB becomes considerably larger for y=0.10 and 0.20 as compared to pure PZT. This is attributed to the considerably lower grain size of our samples resulting from the adopted preparation method. (C) 2010 Elsevier B.V. All rights reserved.

Time-domain-finite-wave analysis of the engine exhaust system by means of the stationary-frame method of characteristics. Part I. Theory

Time-domain-finite-wave analysis of the engine exhaust system is usually done using the method of characteristics. This makes use of either the moving frame method, or the stationary frame method. The stationary frame method is more convenient than its counterpart inasmuch as it avoids the tedium of graphical computations. In this paper (part I), the stationary-frame computational scheme along with the boundary conditions has been implemented. The analysis of a uniform tube, cavity-pipe junction including the engine and the radiation ends, and also the simple area discontinuities has been presented. The analysis has been done accounting for wall friction and heat-transfer for a one-dimensional unsteady flow. In the process, a few inconsistencies in the formulations reported in the literature have been pointed out and corrected. In the accompanying paper (part II) results obtained from the simulation are shown to be in good agreement with the experimental observations.