Learning-Based Fast Iterative Convergence of 3-D MoM via Eigen-AGMRES Method

Three-dimensional (3-D) full-wave electromagnetic simulation using method of moments (MoM) under the framework of fast solver algorithms like fast multipole method (FMM) is often bottlenecked by the speed of convergence of the Krylov-subspace-based iterative process. This is primarily because the electric field integral equation (EFIE) matrix, even with cutting-edge preconditioning techniques, often exhibits bad spectral properties arising from frequency or geometry-based ill-conditioning, which render iterative solvers slow to converge or stagnate occasionally. In this communication, a novel technique to expedite the convergence of MoMmatrix solution at a specific frequency is proposed, by extracting and applying Eigen-vectors from a previously solved neighboring frequency in an augmented generalized minimum residual (AGMRES) iterative framework. This technique can be applied in unison with any preconditioner. Numerical results demonstrate up to 40% speed-up in convergence using the proposed Eigen-AGMRES method.

Modeling of dynamic material behavior in hot deformation: Forging of Ti-6242

A new method of modeling material behavior which accounts for the dynamic metallurgical processes occurring during hot deformation is presented. The approach in this method is to consider the workpiece as a dissipator of power in the total processing system and to evaluate the dissipated power co-contentJ = ∫o σ ε ⋅dσ from the constitutive equation relating the strain rate (ε) to the flow stress (σ). The optimum processing conditions of temperature and strain rate are those corresponding to the maximum or peak inJ. It is shown thatJ is related to the strain-rate sensitivity (m) of the material and reaches a maximum value(J max) whenm = 1. The efficiency of the power dissipation(J/J max) through metallurgical processes is shown to be an index of the dynamic behavior of the material and is useful in obtaining a unique combination of temperature and strain rate for processing and also in delineating the regions of internal fracture. In this method of modeling, noa priori knowledge or evaluation of the atomistic mechanisms is required, and the method is effective even when more than one dissipation process occurs, which is particularly advantageous in the hot processing of commercial alloys having complex microstructures. This method has been applied to modeling of the behavior of Ti-6242 during hot forging. The behavior of α+ β andβ preform microstructures has been exam-ined, and the results show that the optimum condition for hot forging of these preforms is obtained at 927 °C (1200 K) and a strain rate of 1CT•3 s•1. Variations in the efficiency of dissipation with temperature and strain rate are correlated with the dynamic microstructural changes occurring in the material.

On an analytical-numerical procedure for the analysis of cylindrical shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

Bi-criteria single machine scheduling problem with a learning effect: Aneja-Nair method to obtain the set of optimal sequences

In this paper, we consider the bi-criteria single machine scheduling problem of n jobs with a learning effect. The two objectives considered are the total completion time (TC) and total absolute differences in completion times (TADC). The objective is to find a sequence that performs well with respect to both the objectives: the total completion time and the total absolute differences in completion times. In an earlier study, a method of solving bi-criteria transportation problem is presented. In this paper, we use the methodology of solvin bi-criteria transportation problem, to our bi-criteria single machine scheduling problem with a learning effect, and obtain the set of optimal sequences,. Numerical examples are presented for illustrating the applicability and ease of understanding.

CSSI-PRO: a method for secondary structure type editing, assignment and estimation in proteins using linear combination of backbone chemical shifts

Estimation of secondary structure in polypeptides is important for studying their structure, folding and dynamics. In NMR spectroscopy, such information is generally obtained after sequence specific resonance assignments are completed. We present here a new methodology for assignment of secondary structure type to spin systems in proteins directly from NMR spectra, without prior knowledge of resonance assignments. The methodology, named Combination of Shifts for Secondary Structure Identification in Proteins (CSSI-PRO), involves detection of specific linear combination of backbone H-1(alpha) and C-13' chemical shifts in a two-dimensional (2D) NMR experiment based on G-matrix Fourier transform (GFT) NMR spectroscopy. Such linear combinations of shifts facilitate editing of residues belonging to alpha-helical/beta-strand regions into distinct spectral regions nearly independent of the amino acid type, thereby allowing the estimation of overall secondary structure content of the protein. Comparison of the predicted secondary structure content with those estimated based on their respective 3D structures and/or the method of Chemical Shift Index for 237 proteins gives a correlation of more than 90% and an overall rmsd of 7.0%, which is comparable to other biophysical techniques used for structural characterization of proteins. Taken together, this methodology has a wide range of applications in NMR spectroscopy such as rapid protein structure determination, monitoring conformational changes in protein-folding/ligand-binding studies and automated resonance assignment.

On an analytical-numerical procedure for the analysis of cylindrical-shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

The transient response of certain third-ordernon-linear systems

In this paper a method of solving certain third-order non-linear systems by using themethod of ultraspherical polynomial approximation is proposed. By using the method of variation of parameters the third-order equation is reduced to three partial differential equations. Instead of being averaged over a cycle, the non-linear functions are expanded in ultraspherical polynomials and with only the constant term retained, the equations are solved. The results of the procedure are compared with the numerical solutions obtained on a digital computer. A degenerate third-order system is also considered and results obtained for the above system are compared with numerical results obtained on the digital computer. There is good agreement between the results obtained by the proposed method and the numerical solution obtained on digital computer.

A unified approach to the solution of plane problems of Magneto-elasticity with special reference to a hole in a thin infinite conducting plate

Under certain specific assumption it has been observed that the basic equations of magneto-elasticity in the case of plane deformation lead to a biharmonic equation, as in the case of the classical plane theory of elasticity. The method of solving boundary value problems has been properly modified and a unified approach in solving such problems has been suggested with special reference to problems relating thin infinite plates with a hole. Closed form expressions have been obtained for the stresses due to a uniform magnetic field present in the plane of deformation of a thin infinite conducting plate with a circular hole, the plate being deformed by a tension acting parallel to the direction of the magnetic field.

Application of ultraspherical polynomials to forced oscillations of a third-order non-linear system

In this paper the method of ultraspherical polynomial approximation is applied to study the steady-state response in forced oscillations of a third-order non-linear system. The non-linear function is expanded in ultraspherical polynomials and the expansion is restricted to the linear term. The equation for the response curve is obtained by using the linearized equation and the results are presented graphically. The agreement between the approximate solution and the analog computer solution is satisfactory. The problem of stability is not dealt with in this paper.

Study of third-order non-linear systems- Variation of parameters approach

In this paper, the transient response of a third-order non-linear system is obtained by first reducing the given third-order equation to three first-order equations by applying the method of variation of parameters. On the assumption that the variations of amplitude and phase are small, the functions are expanded in ultraspherical polynomials. The expansion is restricted to the constant term. The resulting equations are solved to obtain the response of the given third-order system. A numerical example is considered to illustrate the method. The results show that the agreement between the approximate and digital solution is good thus vindicating the approximation.

Analysis of vibration of clamped square plates by the Rayleigh-Ritz method with asymptotic solutions from a modified Bolotin method

Estimates of flexural frequencies of clamped square plates are initially obtained by the modified Bolotin's method. The mode shapes in “each direction” are then determined and the product functions of these mode shapes are used as admissible functions in the Rayleigh-Ritz method. The data for the first twenty eigenvalues in each of the three (four) symmetric groups obtained by the (i) Bolotin, (ii) Rayleigh and (iii) Rayleigh-Ritz methods are reported here. The Rayleigh estimates are found to be much closer to the true eigenvalues than the Bolotin estimates. The present product functions are found to be much superior to the conventional beam eigenmodes as admissible functions in the Rayleigh-Ritz method of analysis.

On Multi-Dimensional Packets of Surface Waves

By using the method of operators of multiple scales, two coupled nonlinear equations are derived, which govern the slow amplitude modulation of surface gravity waves in two space dimensions. The equations of Davey and Stewartson, which also govern the two-dimensional modulation of the amplitude of gravity waves, are derived as a special case of our equations. For a fully dispersed wave, symmetric about a point which moves with the group velocity, the coupled equations reduce to a nonlinear Schrödinger equation with extra terms representing the effect of the curvature of the wavefront.

The non-linear response of a gyrostabilized platform: Effects of correction torque time delay and of random inputs

First, the non-linear response of a gyrostabilized platform to a small constant input torque is analyzed in respect to the effect of the time delay (inherent or deliberately introduced) in the correction torque supplied by the servomotor, which itself may be non-linear to a certain extent. The equation of motion of the platform system is a third order nonlinear non-homogeneous differential equation. An approximate analytical method of solution of this equation is utilized. The value of the delay at which the platform response becomes unstable has been calculated by using this approximate analytical method. The procedure is illustrated by means of a numerical example. Second, the non-linear response of the platform to a random input has been obtained. The effects of several types of non-linearity on reducing the level of the mean square response have been investigated, by applying the technique of equivalent linearization and solving the resulting integral equations by using laguerre or Gaussian integration techniques. The mean square responses to white noise and band limited white noise, for various values of the non-linear parameter and for different types of non-linearity function, have been obtained. For positive values of the non-linear parameter the levels of the non-linear mean square responses to both white noise and band-limited white noise are low as compared to the linear mean square response. For negative values of the non-linear parameter the level of the non-linear mean square response at first increases slowly with increasing values of the non-linear parameter and then suddenly jumps to a high level, at a certain value of the non-linearity parameter.

The pulse response of non-linear third order systems

The response of a third order non-linear system subjected to a pulse excitation is analysed. A transformation of the displacement variable is effected. The transformation function chosen is the solution of the linear problem subjected to the same pulse. With this transformation the equation of motion is brought into a form in which the method of variation of parameters is applicable for the solution of the problem. The method is applied to a single axis gyrostabilized platform subjected to an exponentially decaying pulse. The analytical results are compared with digital and analog computer solutions.

Entire Solutions of Hydrodynamical Equations with Exponential Dissipation

We consider a modification of the three-dimensional Navier-Stokes equations and other hydrodynamical evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is replaced by an operator whose Fourier symbol grows exponentially as e(vertical bar k vertical bar/kd) at high wavenumbers vertical bar k vertical bar. Using estimates in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately entire in the space variables and that the Fourier coefficients decay faster than e-(C(k/kd) ln(vertical bar k vertical bar/kd)) for any C < 1/(2 ln 2). The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der Hoeven, indicate that the leading-order asymptotics is precisely of the above form with C = C-* = 1/ ln 2. The same behavior with a universal constant C-* is conjectured for the Navier-Stokes equations with exponential dissipation in any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which is obtained for ordinary Navier-Stokes turbulence. Possible applications to improved spectral simulations are briefly discussed.