A multistep transversal linearization (MTL) method in non-linear dynamics through a Magnus characterization

A new form of a multi-step transversal linearization (MTL) method is developed and numerically explored in this study for a numeric-analytical integration of non-linear dynamical systems under deterministic excitations. As with other transversal linearization methods, the present version also requires that the linearized solution manifold transversally intersects the non-linear solution manifold at a chosen set of points or cross-section in the state space. However, a major point of departure of the present method is that it has the flexibility of treating non-linear damping and stiffness terms of the original system as damping and stiffness terms in the transversally linearized system, even though these linearized terms become explicit functions of time. From this perspective, the present development is closely related to the popular practice of tangent-space linearization adopted in finite element (FE) based solutions of non-linear problems in structural dynamics. The only difference is that the MTL method would require construction of transversal system matrices in lieu of the tangent system matrices needed within an FE framework. The resulting time-varying linearized system matrix is then treated as a Lie element using Magnus’ characterization [W. Magnus, On the exponential solution of differential equations for a linear operator, Commun. Pure Appl. Math., VII (1954) 649–673] and the associated fundamental solution matrix (FSM) is obtained through repeated Lie-bracket operations (or nested commutators). An advantage of this approach is that the underlying exponential transformation could preserve certain intrinsic structural properties of the solution of the non-linear problem. Yet another advantage of the transversal linearization lies in the non-unique representation of the linearized vector field – an aspect that has been specifically exploited in this study to enhance the spectral stability of the proposed family of methods and thus contain the temporal propagation of local errors. A simple analysis of the formal orders of accuracy is provided within a finite dimensional framework. Only a limited numerical exploration of the method is presently provided for a couple of popularly known non-linear oscillators, viz. a hardening Duffing oscillator, which has a non-linear stiffness term, and the van der Pol oscillator, which is self-excited and has a non-linear damping term.

Beta-structure of polypeptides in non-aqueous solutions. I. Spectral characteristics of the polypeptide backbone

The optical rotatory features of the beta-structure of the polypeptides in non-aqueous solutions and films cast from these solutions have been investigated. The beta-structure of poly-S-benzyl-L-cysteine, poly-S-carbobenzoxy-L-cysteine and poly-S-benzyl-L-cysteine, poly-S-carbobenzoxy-L-cysteine and poly-O-carbo-bands of their films. The optical rotatory dispersion (ORD) and circular dichroism (CD) spectra of these polypeptides are found to be very similar in both film and solution. In solvents promoting the beta-structure, the polypeptides are characterized by CD troughs in the n-pi* transition region of the peptide chromophore. The ORD spectra are found to be positive in sign throughout the visible and accessible ultraviolet regions and are interpreted in terms of the possible existence of a relatively much larger positive pi-pi* CD bands as compared with the negative n-pi* band. The rotatory data obtained in the non-aqueous solution are compared with those obtained for other poly peptides in aqueous solutions, with respect to the type and extent of beta-structure present.

A detailed study of Li-DNA fibres at various salt concentrations reveals a non-helical B-DNA and a possible similarity of solution and solid state structures

A thorough investigation of salt concentration dependence of lithium DNA fibres is made using X-ray diffraction. While for low salt the C-form pattern is obtained, crystalline B-type diffraction patterns result on increasing the salt concentration. The salt content in the gel (from which fibres are drawn) is estimated by equilibrium dialysis using the Donnan equilibrium principle. The salt range giving the best crystalline B pattern is determined. It is found that in this range meridional reflections occur on the fourth and sixth layer lines. In addition, the tenth layer meridian is absent at a particular salt concentration. These results strongly suggest the presence of non-helical features in the DNA molecule. Preliminary analysis of the diffraction patterns indicates a structural variability within the B-form itself. Further, the possibility of the structural parameters of DNA being similar in solid state and in solution is discussed.

On the non-existence of hydrodynamic solution of surface plasmons at the bimetallic interface

It is shown that for an abrupt bimetallic interface a hydrodynamic solution for interface plasmons does not exist. It appears that this result is valid irrespective of the choice of of the additional boundary condition, thereby suggesting a careful look at the use of usual hydrodynamic equations for a bimetallic interface.

Application of Laplace transform technique to the solution of certain third-order non-linear systems

A number of papers have appeared on the application of operational methods and in particular the Laplace transform to problems concerning non-linear systems of one kind or other. This, however, has met with only partial success in solving a class of non-linear problems as each approach has some limitations and drawbacks. In this study the approach of Baycura has been extended to certain third-order non-linear systems subjected to non-periodic excitations, as this approximate method combines the advantages of engineering accuracy with ease of application to such problems. Under non-periodic excitations the method provides a procedure for estimating quickly the maximum response amplitude, which is important from the point of view of a designer. Limitations of such a procedure are brought out and the method is illustrated by an example taken from a physical situation.

A solution framework for scheduling a BPM with non-identical job dimensions

A Batch Processing Machine (BPM) is one which processes a number of jobs simultaneously as a batch with common beginning and ending times. Also, a BPM, once started cannot be interrupted in between (Pre-emption not allowed). This research is motivated by a BPM in steel casting industry. There are three main stages in any steel casting industry viz., pre-casting stage, casting stage and post-casting stage. A quick overview of the entire process, is shown in Figure 1. There are two BPMs : (1) Melting furnace in the pre-casting stage and (2) Heat Treatment Furnace (HTF) in the post casting stage of steel casting manufacturing process. This study focuses on scheduling the latter, namely HTF. Heat-treatment operation is one of the most important stages of steel casting industries. It determines the final properties that enable components to perform under demanding service conditions such as large mechanical load, high temperature and anti-corrosive processing. In general, different types of castings have to undergo more than one type of heat-treatment operations, where the total heat-treatment processing times change. To have a better control, castings are primarily classified into a number of job-families based on the alloy type such as low-alloy castings and high alloy castings. For technical reasons such as type of alloy, temperature level and the expected combination of heat-treatment operations, the castings from different families can not be processed together in the same batch.

Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

Exact Solution of a One-Dimensional Multicomponent Lattice Gas with Hyperbolic Interaction

We present the exact solution to a one-dimensional multicomponent quantum lattice model interacting by an exchange operator which falls off as the inverse sinh square of the distance. This interaction contains a variable range as a parameter and can thus interpolate between the known solutions for the nearest-neighbor chain and the inverse-square chain. The energy, susceptibility, charge stiffness, and the dispersion relations for low-lying excitations are explicitly calculated for the absolute ground state, as a function of both the range of the interaction and the number of species of fermions.

Non-Marcus energy gap dependence of electron transfer rate in contact ion pairs. Novel interplay between relaxation and reaction in solution

Electron transfer reactions between donor-acceptor pairs in solution and in organized media exhibit diverse behaviour. Recent experiments have indicated an interesting breakdown of the Marcus parabolic energy gap dependence in the normal regime for back electron transfer from contact ion pairs. A novel explanation of this breakdown has recently been proposed (M. Tachiya and S. Murata, J. Am. Chem. Sec., 116(1994) 2434) which attributes the breakdown to the interplay between the relaxation in the reactant well and the reaction. A particularly interesting aspect of the model is that it envisages the electron transfer in the normal regime to take place from a completely non-equilibrium condition. In this article a time dependent solution of the model is presented for the first time, after generalizing it to include a realistic initial population distribution. The decay of the contact ion pair population is completely non-exponential. This can be used to check the validity of the Tachiya-Murata model. The dynamics of electron transfer from the solvent separated ion pair, which seem to obey the Marcus relation, is exponential.

Solution of a class of surface water wave scattering problems involving non-reflecting vertical barriers

A large class of scattering problems of surface water waves by vertical barriers lead to mixed boundary value problems for Laplace equation. Specific attentions are paid, in the present article, to highlight an analytical method to handle this class of problems of surface water wave scattering, when the barriers in question are non-reflecting in nature. A new set of boundary conditions is proposed for such non-reflecting barriers and tile resulting boundary value problems are handled in the linearized theory of water waves. Three basic poblems of scattering by vertical barriers are solved. The present new theory of non-reflecting vertical barriers predict new transmission coefficients and tile solutions of tile mathematical problems turn out to be extremely simple and straight forward as compared to the solution for other types of barriers handled previously.

Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.

The nature of non-equilibrium solidification of undercooled Agsingle bondCu alloys in contact with primary copper dendrites

Liquids of silver-copper alloys with near eutectic compositions embedded in a copper matrix were undercooled. The structural and microstructural investigations of these alloys solidified from undercooled temperature indicated the absence of both the eutectic reaction and diffusionless transformation below the equal free energy curve (T0). Instead the liquid maintained local equilibrium with the copper dendrites continuously until it intersected the extended solidus of the silver rich solid solution.

Non-linear resonance in strings under narrow band random excitation, part I: Planar response and stability

Non-linear planar response of a string to planar narrow band random excitation is investigated in this paper. A response equation for the mean square deflection σ2 is obtained under a single mode approximation by using the equivalent linearization technique. It is shown that the response is triple valued, as in the case of harmonic excitation, if the centre frequency of excitation Ω lies in a certain specified range. The triple valued response occurs only if the excitation bandwidth β is smaller than a critical value βcrit which is a monotonically increasing function of the intensity of excitation. An approximate method of investigating the almost sure asymptotic stability of the solution is presented and regions of instability in the Ω-σ2 plane have been charted. It is shown that planar response can become unstable either due to an unbounded growth of the in-plane component of motion or due to a spontaneous appearance of an out-of-plane component.

Sol → gel transition with rings in a finite polycondensing system in solution

We report a theoretical formulation for the mean cluster size distribution in a finite polycondensing system. Expressions for the mean number of n-mers with j bonds ( nj) are developed. Numerical calculations show that while the non-cyclic molecules make the dominant contribution to the small clusters, the large clusters are dominated by cyclic structures. The number of particles in ringless chains, n n,n-1, decays monotonically with n at all extents of reaction, but n n becomes bimodal near the gel point. We also find that the solvent plays an important role in the cluster size distribution.