An Application of 3-D Kinematical Conservation Laws: Propagation of a 3-D Wavefront

Three-dimensional (3-D) kinematical conservation laws (KCL) are equations of evolution of a propagating surface Omega(t) in three space dimensions. We start with a brief review of the 3-D KCL system and mention some of its properties relevant to this paper. The 3-D KCL, a system of six conservation laws, is an underdetermined system to which we add an energy transport equation for a small amplitude 3-D nonlinear wavefront propagating in a polytropic gas in a uniform state and at rest. We call the enlarged system of 3-D KCL with the energy transport equation equations of weakly nonlinear ray theory (WNLRT). We highlight some interesting properties of the eigenstructure of the equations of WNLRT, but the main aim of this paper is to test the numerical efficacy of this system of seven conservation laws. We take several initial shapes for a nonlinear wavefront with a suitable amplitude distribution on it and let it evolve according to the 3-D WNLRT. The 3-D WNLRT is a weakly hyperbolic 7 x 7 system that is highly nonlinear. Here we use the staggered Lax-Friedrichs and Nessyahu-Tadmor central schemes and have obtained some very interesting shapes of the wavefronts. We find the 3-D KCL to be suitable for solving many complex problems for which there presently seems to be no other method capable of giving such physically realistic features.

Development of twist in an emerging magnetic flux tube by poloidal field accretion

Aims. Following an earlier proposal for the origin of twist in the magnetic fields of solar active regions, we model the penetration of a wrapped up background poloidal field into a toroidal magnetic flux tube rising through the solar convective zone.Methods. The rise of the straight, cylindrical flux tube is followed by numerically solving the induction equation in a comoving Lagrangian frame, while an external poloidal magnetic field is assumed to be radially advected onto the tube with a speed corresponding to the rise velocity.Results. One prediction of our model is the existence of a ring of reverse current helicity on the periphery of active regions. On the other hand, the amplitude of the resulting twist depends sensitively on the assumed structure ( diffuse vs. concentrated/intermittent) of the active region magnetic field right before its emergence, and on the assumed vertical profile of the poloidal field. Nevertheless, in the model with the most plausible choice of assumptions a mean twist comparable to the observations results.Conclusions. Our results indicate that the contribution of this mechanism to the twist can be quite significant, and under favourable circumstances it can potentially account for most of the current helicity observed in active regions.

Generalized Burgers equations and Euler–Painlevé transcendents. III

It was proposed earlier [P. L. Sachdev, K. R. C. Nair, and V. G. Tikekar, J. Math. Phys. 27, 1506 (1986); P. L. Sachdev and K. R. C. Nair, ibid. 28, 977 (1987)] that the Euler–Painlevé equations  y(d2y/dη2)+a(dy/dη)2 +f(η)y(dy/dη)+g(η)y2+b(dy/dη) +c=0 represent generalized Burgers equations (GBE’s) in the same way as Painlevé equations represent the Korteweg–de Vries type of equations. The earlier studies were carried out in the context of GBE’s with damping and those with spherical and cylindrical symmetry. In the present paper, GBE’s with variable coefficients of viscosity and those with inhomogeneous terms are considered for their possible connection to Euler–Painlevé equations. It is found that the Euler–Painlevé equation, which represents the GBE ut+uβux=(δ/2)g(t)uxx, g(t)=(1+t)n, β>0, has solutions, which either decay or oscillate at η=±∞, only when −1

Shaping of adherends in bonded joints

A method of analysis of adhesive bonded joints with arbitrary adherend shapes is developed and applied to their design for efficient load transfer. The formulation of the simplified physical model is essentially one dimensional. An efficient numerical procedure is developed and its applicability is demonstrated. Guide lines are indicated for choosing suitable shape functions.

Complexes of lanthanide nitrates with O,O′,N-triisopropyl phosphoramidate

A substituted phosphoramidate has been used as a ligand to lanthanides for the first time. New complexes of lanthanide nitrates with O,O′,N-triisopropyl phosphoramidate (TIP) of the general formula Ln(TIP)3(NO3)3 where Ln=La-Yb and Y have been synthesised and characterised by chemical analysis, infrared and visible electronic spectra and electrical conductance.Infrared spectra indicate the coordination of the ligand to the metal ions through the oxygen of the P=O group. IR and conductance show that the nitrate groups are all coordinated. Electronic spectral shapes have been interpreted in terms of an eight coordinate geometry around the metal ions.

Approximate multipole coefficients of RF ion traps as functions of aperture size

In this study we present approximate analytical expressions for estimating the variation in multipole expansion coefficients as a function of the size of the apertures in the electrodes in axially symmetric (3D) and two-dimensional (2D) ion trap ion traps. Following the approach adopted in our earlier studies which focused on the role of apertures to fields within the traps, here too, the analytical expression we develop is a sum of two terms, A(n,noAperiure), the multipole expansion coefficient for a trap with no apertures and A(n,dueToAperture), the multipole expansion coefficient contributed by the aperture. A(n,noAperture) has been obtained numerically and A(n,dueToAperture) is obtained from the n th derivative of the potential within the trap. The expressions derived have been tested on two 3D geometries and two 2D geometries. These include the quadrupole ion trap (QIT) and the cylindrical ion trap (CIT) for 3D geometries and the linear ion trap (LIT) and the rectilinear ion trap (RIT) for the 2D geometries. Multipole expansion coefficients A(2) to A(12), estimated by our analytical expressions, were compared with the values obtained numerically (using the boundary element method) for aperture sizes varying up to 50% of the trap dimension. In all the plots presented, it is observed that our analytical expression for the variation of multipole expansion coefficients versus aperture size closely follows the trend of the numerical evaluations for the range of aperture sizes considered. The maximum relative percentage errors, which provide an estimate of the deviation of our values from those obtained numerically for each multipole expansion coefficient, are seen to be largely in the range of 10-15%. The leading multipole expansion coefficient, A(2), however, is seen to be estimated very well by our expressions, with most values being within 1% of the numerically determined values, with larger deviations seen for the QIT and the LIT for large aperture sizes. (C) 2010 Elsevier B.V. All rights reserved.

Design of Constant Accuracy Linear Proportional Weir

Five different shaped weirs were designed and pertinent data for their use are given. One of these weir shapes had the least “sharp edge” at the junction of the base weir and “complementary weir.” Two other types of weirs had equal slopes at the junction of the base weir and complementary weir. Another shape, for which neither the indication accuracy was constant nor the slope was equal at the junction of the base weir and complementary weir, was also tested. The results of the four weir shapes hydraulically tested give consistent values for the coefficient of discharge varying between 0.625 to 0.631. The indication accuracies of all the previously designed linear proportional weirs (includig Sutro weir) are neither constant nor unity, as is believed.

Correlation of helicopter rotor aeroelastic response with HART-II wind tunnel test data

Purpose - This paper aims to validate a comprehensive aeroelastic analysis for a helicopter rotor with the higher harmonic control aeroacoustic rotor test (HART-II) wind tunnel test data. Design/methodology/approach - Aeroelastic analysis of helicopter rotor with elastic blades based on finite element method in space and time and capable of considering higher harmonic control inputs is carried out. Moderate deflection and coriolis nonlinearities are included in the analysis. The rotor aerodynamics are represented using free wake and unsteady aerodynamic models. Findings - Good correlation between analysis and HART-II wind tunnel test data is obtained for blade natural frequencies across a range of rotating speeds. The basic physics of the blade mode shapes are also well captured. In particular, the fundamental flap, lag and torsion modes compare very well. The blade response compares well with HART-II result and other high-fidelity aeroelastic code predictions for flap and torsion mode. For the lead-lag response, the present analysis prediction is somewhat better than other aeroelastic analyses. Research limitations/implications - Predicted blade response trend with higher harmonic pitch control agreed well with the wind tunnel test data, but usually contained a constant offset in the mean values of lead-lag and elastic torsion response. Improvements in the modeling of the aerodynamic environment around the rotor can help reduce this gap between the experimental and numerical results. Practical implications - Correlation of predicted aeroelastic response with wind tunnel test data is a vital step towards validating any helicopter aeroelastic analysis. Such efforts lend confidence in using the numerical analysis to understand the actual physical behavior of the helicopter system. Also, validated numerical analyses can take the place of time-consuming and expensive wind tunnel tests during the initial stage of the design process. Originality/value - While the basic physics appears to be well captured by the aeroelastic analysis, there is need for improvement in the aerodynamic modeling which appears to be the source of the gap between numerical predictions and HART-II wind tunnel experiments.

Analytical prediction of break-out noise from a reactive rectangular plenum with four flexible walls

This paper describes an analytical calculation of break-out noise from a rectangular plenum with four flexible walls by incorporating three-dimensional effects along with the acoustical and structural wave coupling phenomena. The breakout noise from rectangular plenums is important and the coupling between acoustic waves within the plenum and structural waves in the flexible plenum walls plays a critical role in prediction of the transverse transmission loss. The first step in breakout noise prediction is to calculate the inside plenum pressure field and the normal flexible plenum wall vibration by using an impedance-mobility approach, which results in a compact matrix formulation. In the impedance-mobility compact matrix (IMCM) approach, it is presumed that the coupled response can be described in terms of finite sets of the uncoupled acoustic subsystem and the structural subsystem. The flexible walls of the plenum are modeled as an unfolded plate to calculate natural frequencies and mode shapes of the uncoupled structural subsystem. The second step is to calculate the radiated sound power from the flexible walls using Kirchhoff-Helmholtz (KH) integral formulation. Analytical results are validated with finite element and boundary element (FEM-BEM) numerical models. (C) 2010 Acoustical Society of America. DOI: 10.1121/1.3463801]

Free vibration of a stochastic beam-column using stochastic fem

The free vibrational characteristics of a beam-column, which is having randomly varying Young's modulus and mass density and subjected to randomly distributed axial loading is analysed. The material property fluctuations and axial loadings are considered to constitute independent one-dimensional, uni-variate, homogeneous real, spatially distributed stochastic fields. Hamilton's principle is used to formulate the problem using stochastic FEM. Vibration frequencies and mode shapes are analysed for their statistical descriptions. A numerical example is shown.

Numerical estimation of hotspot temperatures in electrodes subjected to pulsed electron beam heating in vacuum

A numerical solution for the transient temperature distribution in a cylindrical disc heated on its top surface by a circular source is presented. A finite difference form of the governing equations is solved by the Alternating Direction Implicit (ADI) time marching scheme. This solution has direct applications in analyzing transient electron beam heating of target materials as encountered in the prebreakdown current enhancement and consequent breakdown in high voltage vacuum gaps stressed by alternating and pulsed voltages. The solution provides an estimate of the temperature for pulsed electron beam heating and the size of thermally activated microparticles originating from anode hot spots. The calculated results for a typical 45kV (a.c.) electron beam of radius 2.5 micron indicate that the temperature of such spots can reach melting point and could give rise to microparticles which could initiate breakdown.

The Single Langmuir Probe under Orbital-Limited Conditions in a Low-Density Collisional Plasma

The cylindrical Langmuir probe under orbital-limited conditions was used to determine the charge density in a low-density collisional plasma. The Langmuir's theory was applied to both electron and ion saturation currents in their respective accelerating regions. Present study indicates that the length of the probe significantly affects the probe characteristics. A probe of suitable length under orbital-limited conditions may be useful under the experimental conditions where the radius of the probe is much smaller than the Debye lengt.

Epoch extraction from linear prediction residual for identification of closed glottis interval

In voiced speech analysis epochal information is useful in accurate estimation of pitch periods and the frequency response of the vocal tract system. Ideally, linear prediction (LP) residual should give impulses at epochs. However, there are often ambiguities in the direct use of LP residual since samples of either polarity occur around epochs. Further, since the digital inverse filter does not compensate the phase response of the vocal tract system exactly, there is an uncertainty in the estimated epoch position. In this paper we present an interpretation of LP residual by considering the effect of the following factors: 1) the shape of glottal pulses, 2) inaccurate estimation of formants and bandwidths, 3) phase angles of formants at the instants of excitation, and 4) zeros in the vocal tract system. A method for the unambiguous identification of epochs from LP residual is then presented. The accuracy of the method is tested by comparing the results with the epochs obtained from the estimated glottal pulse shapes for several vowel segments. The method is used to identify the closed glottis interval for the estimation of the true frequency response of the vocal tract system.

A geodesic-based triangulation technique for damage location in metallic and composite plates

Lamb-wave-based damage detection methods using the triangulation technique are not suitable for handling structures with complex shapes and discontinuities as the parametric/analytical representation of these structures is very difficult. The geodesic concept is used along with the triangulation technique to overcome the above problem. The present work is based on the fundamental fact that a wave takes the minimum energy path to travel between two points on any multiply connected surface and this reduces to the shortest distance path or geodesic. The geodesics are computed on the meshed surface of the structure using the fast marching method. The wave response matrix of the given sensor configuration for the healthy and the damaged structure is obtained experimentally. The healthy and damage response matrices are compared and their difference gives the time information about the reflection of waves from the damage. A wavelet transform is used to extract the arrival time information of the wave scattered by the damage from the acquired Lamb wave signals. The computed geodesics and time information are used in the ellipse algorithm of triangulation formulation to locate the loci of possible damage location points for each actuator-sensor pair. The results obtained for all actuator-sensor pairs are combined and the intersection of multiple loci gives the damage location result. Experiments were conducted in aluminum and composite plate specimens to validate this method.

A geodesic constant method for computing high-frequency mutual coupling between antennas on general quadric cylinders

A Geodesic Constant Method (GCM) is outlined which provides a common approach to ray tracing on quadric cylinders in general, and yields all the surface ray-geometric parameters required in the UTD mutual coupling analysis of conformal antenna arrays in the closed form. The approach permits the incorporation of a shaping parameter which permits the modeling of quadric cylindrical surfaces of desired sharpness/flatness with a common set of equations. The mutual admittance between the slots on a general parabolic cylinder is obtained as an illustration of the applicability of the GCM.