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.

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation-based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1-D as well as the 3-D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1-D problems.

New colorimetric method for the estimation of aldehyde derivatives of vitamin A

A new colorimetric method for the estimation of aldehyde derivatives of vitamin A is described. The assay depends upon the formation of colored derivatives of the vitamin A aldehydes by reaction with p-aminobenzoic acid or p-aminosalicylic acid in the presence of 2.0 N hydrochloric acid.

A colorimetric method for the estimation of phloroglucinol

1. 1. A simple method has been devised for the estimation of phloroglucinol based on the formation of an intense colored compound with a modified Ehrlich reagent in the presence of trichloroacetic acid. Some factors affecting the formation of color have been studied. 2. 2. Careful regulation of trichloroacetic acid content in the system permits its estimation in 1–15 μg in the micro range and in 10–50 μg in the macro range. Phloroglucinol lends itself to ready separation by paper chromatography and estimation after clution from paper.

A new method for the titrimetric determination of perchlorate

A new method for the titrimetric determination of perchlorate has been developed, based on its reduction to chloride by iron(II) in a strong sulphuric acid medium at high temperature. The effect of variables, such as the sulphuric acid concentration, the temperature and the period of heating, on the extent of reduction has been studied and the optimal conditions for analytical determination of perehlorate derived.

A colorimetric method for determination of pyrocatechol and related substances

A colorimetric assay for the quantitative determination of catecholic compounds was developed. The method was based on the observation that a red color was formed when nitrite was added to a solution containing pyrocatechol and sodium tungstate. Aromatic amines interfere with the reaction but this could be overcome by the addition of formaldehyde. When interfering substances are present along with pyrocatechol, it can be readily separated by paper chromatography and estimated after elution from the filter paper.

An acid method for the volumetric estimation of water-soluble dithiocarbamates

Dithiocarbamates have been estimated previously by reaction with a strong acid, the carbon disulfide evolved being converted into a xanthate and the latter estimated iodimetrically. In the present method, a water-soluble dithiocarbamate is reacted with a decinormal mineral acid and the excess acid is determined to compute the amount of dithiocarbamate present. This method is applicable for the determination of a dithiocarbamate in a mixture containing thiuram disulfide.

On Application Of The Method Of Kinetic Invariants To The Description Of Dissolution Accompanied By A Chemical Reaction

The dissolution, accompanied by chemical reaction, of monodisperse solid particles has been analysed. The resulting model, which accounts for the variation of mass transfer coefficient with the size of the dissolving particles, yields an approximate analytical form of a kinetic function. Rigorous numerical and approximate analytical solutions have been obtained for the governing system of nonlinear ordinary differential equations. The transient nature of the dissolution process as well as the accuracy of the analytical solution is brought out by the rigorous numerical solution. The analytical solution is fairly accurate for the major part of the range of operational times encountered in practice.

A simplified method for the purification of riboflavin-binding protein from hen's egg yolk

A simple and efficient procedure for the purification of the riboflavin-binding protein from hen's egg yolk is described. This method involves the removal by exclusion of lipoproteins and subsequent fractionation of soluble yolk proteins held on a DEAE-cellulose column by a salt gradient which is followed by purification by gel filtration on Sephadex G-100. The protein thus isolated is homogeneous by various physicoehemical, immunological, and functional criteria.

Plant aspartate transcarbamylase: An affinity chromatographic method for the purification of the enzyme from germinated seedlings

A simple and rapid affinity chromatographic method for the isolation of aspartate transcarbamylase from germinated seedlings of mung bean (Phaseolus aureus) was developed. A partially purified preparation of the enzyme was chromatographed on an affinity column containing aspartate linked to CNBr-activated Sepharose 4B. Aspartate transcarbamylase was specifically eluted from the column with 10 mImage aspartate or 0.5 Image KCl. The enzyme migrated as a single sharp band during disc electrophoresis at pH 8.6 on polyacrylamide gels. Electrophoresis of the sodium dodecyl sulfate-treated enzyme showed two distinct protein bands, suggesting that the mung bean aspartate transcarbamylase was made up of nonidentical subunits. Like the enzyme purified by conventional procedures, this enzyme preparation also exhibited positive homotropic interactions with carbamyl phosphate and negative heterotropic interactions with UMP. This method was extended to the purification of aspartate transcarbamylase from Lathyrus sativus, Eleucine coracona, and Trigonella foenum graecum.

Biorthogonal series method for Oseen type flows

A biorthogonal series method is developed to solve Oseen type flow problems. The theory leads to a new set of eigenfunctions for a specific class of linear non-selfadjoint operators containing the biharmonic one. These eigenfunctions differ from those given earlier in the literature for the biharmonic operator. The method is applied to the problem of thermocapillary flow in a cylindrical liquid bridge of finite length with axial through flow. Flow and temperature distributions are obtained at leading order of an expansion for small surface tension Reynolds number and Prandtl number. Another related problem considered is that of cylindrical cavity flow. Solutions for both cases are presented in terms of biorthogonal series. The effect of axial through flow on velocity and temperature fields is discussed by numerical evaluation of the truncated analytical series. The presence of axial through flow not only convectively shifts the vortices induced by surface forces in the direction of the through flow, but also moves their centers toward the outer cylindrical boundary. This process can lead to significantly asymmetric flow structures.

Transverse plane wave analysis of short elliptical chamber mufflers: An analytical approach

Short elliptical chamber mufflers are used often in the modern day automotive exhaust systems. The acoustic analysis of such short chamber mufflers is facilitated by considering a transverse plane wave propagation model along the major axis up to the low frequency limit. The one dimensional differential equation governing the transverse plane wave propagation in such short chambers is solved using the segmentation approaches which are inherently numerical schemes, wherein the transfer matrix relating the upstream state variables to the downstream variables is obtained. Analytical solution of the transverse plane wave model used to analyze such short chambers has not been reported in the literature so far. This present work is thus an attempt to fill up this lacuna, whereby Frobenius solution of the differential equation governing the transverse plane wave propagation is obtained. By taking a sufficient number of terms of the infinite series, an approximate analytical solution so obtained shows good convergence up to about 1300 Hz and also covers most of the range of muffler dimensions used in practice. The transmission loss (TL) performance of the muffler configurations computed by this analytical approach agrees excellently with that computed by the Matrizant approach used earlier by the authors, thereby offering a faster and more elegant alternate method to analyze short elliptical muffler configurations. (C) 2010 Elsevier Ltd. All rights reserved.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Simulation of the infrared spectra of thioacetamide by the extended molecular mechanics method

The infrared spectrum of the matrix-isolated species of thioacetamide has been simulated using the extended molecular mechanics method. The equilibrium structure, vibrational frequencies, dipole moment and infrared absorption intensities of thioacetamide have been calculated in good agreement with the experiment. The vibrational frequencies and infrared absorption intensities for the isotopic molecules (CH2CSNH2)-C-13, (CH3CSNH2)-N-15 and (CH2CSND2)-C-13 have also been calculated consistent with the experiment. The infrared spectra of the matrix isolated species of N- and C- deuterated isotopomers of thioacetamide, CH3CSND2 and CD3CSNH2 have also been simulated in satisfactory agreement with the experimental spectra.