Vibration and Stability of Fluid Conveying Pipes with Stochastic Parameters

Flexible cantilever pipes conveying fluids with high velocity are analysed for their dynamic response and stability behaviour. The Young's modulus and mass per unit length of the pipe material have a stochastic distribution. The stochastic fields, that model the fluctuations of Young's modulus and mass density are characterized through their respective means, variances and autocorrelation functions or their equivalent power spectral density functions. The stochastic non self-adjoint partial differential equation is solved for the moments of characteristic values, by treating the point fluctuations to be stochastic perturbations. The second-order statistics of vibration frequencies and mode shapes are obtained. The critical flow velocity is-first evaluated using the averaged eigenvalue equation. Through the eigenvalue equation, the statistics of vibration frequencies are transformed to yield critical flow velocity statistics. Expressions for the bounds of eigenvalues are obtained, which in turn yield the corresponding bounds for critical flow velocities.

Reactivation kinetics of boron acceptors in hydrogenated silicon during zero bias anneal

The reactivation kinetics of passivated boron accepters in hydrogenated silicon during zero bias annealing in the temperature range of 65-130 degrees C are reported, For large annealing times and high annealing temperatures, the reactivation process follows second-order kinetics and is rate limited by a thermally activated <(H)over tilde (2)> complex formation process, For short annealing times and low annealing temperatures, the reactivation rate is found to be larger than that due to <(H)over tilde (2)> complex formation alone. We conclude that the faster reactivation is caused by the diffusion of the liberated hydrogen atoms into the bulk as well as <(H)over tilde (2)> complex formation. The effective diffusion coefficient of hydrogen is measured and found to obey the Arrhenius relation with an activation energy (1.41 +/- 0.1) eV. (C) 1997 American Institute of Physics.

Magnetization properties of some quantum spin ladders

The experimental realization of various spin ladder systems has prompted their detailed theoretical investigations. Hen we study the evolution of ground-state magnetization with an external magnetic field for two different antiferromagnetic systems: a three-legged spin-1/2 ladder, and a two-legged spin-1/2 ladder with an additional diagonal interaction. The finite system density-matrix renormalization-group method is employed for numerical studies of the three-chain system, and an effective low-energy Hamiltonian is used in the limit of strong interchain coupling to study the two- and three-chain systems. The three-chain system has a magnetization plateau at one-third of the saturation magnetization. The two-chain system has a plateau at zero magnetization due to a gap above the singlet ground state. It also has a plateau at half of the saturation magnetization for a certain range of values of the couplings. We study the regions of transitions between plateaus numerically and analytically, and find that they are described, at first order in a strong-coupling expansion, by an XXZ spin-1/2 chain in a magnetic field; the second-order terms give corrections to the XXZ model, We also study numerically some low-temperature properties of the three-chain system, such as the magnetization, magnetic susceptibility and specific heat. [S0163-1829(99)303001-5].

A differential geometric method for kinematic analysis of two- and three-degree-of-freedom rigid body motions

In this paper, we present a novel differential geometric characterization of two- and three-degree-of-freedom rigid body kinematics, using a metric defined on dual vectors. The instantaneous angular and linear velocities of a rigid body are expressed as a dual velocity vector, and dual inner product is defined on this dual vector, resulting in a positive semi-definite and symmetric dual matrix. We show that the maximum and minimum magnitude of the dual velocity vector, for a unit speed motion, can be obtained as eigenvalues of this dual matrix. Furthermore, we show that the tip of the dual velocity vector lies on a dual ellipse for a two-degree-of-freedom motion and on a dual ellipsoid for a three-degree-of-freedom motion. In this manner, the velocity distribution of a rigid body can be studied algebraically in terms of the eigenvalues of a dual matrix or geometrically with the dual ellipse and ellipsoid. The second-order properties of the two- and three-degree-of-freedom motions of a rigid body are also obtained from the derivatives of the elements of the dual matrix. This results in a definition of the geodesic motion of a rigid body. The theoretical results are illustrated with the help of a spatial 2R and a parallel three-degree-of-freedom manipulator.

Assessment of conditional moment closure models of turbulent autoignition using DNS data

A systematic assessment of the submodels of conditional moment closure (CMC) formalism for the autoignition problem is carried out using direct numerical simulation (DNS) data. An initially non-premixed, n-heptane/air system, subjected to a three-dimensional, homogeneous, isotropic, and decaying turbulence, is considered. Two kinetic schemes, (1) a one-step and (2) a reduced four-step reaction mechanism, are considered for chemistry An alternative formulation is developed for closure of the mean chemical source term , based on the condition that the instantaneous fluctuation of excess temperature is small. With this model, it is shown that the CMC equations describe the autoignition process all the way up to near the equilibrium limit. The effect of second-order terms (namely, conditional variance of temperature excess sigma(2) and conditional correlations of species q(ij)) in modeling is examined. Comparison with DNS data shows that sigma(2) has little effect on the predicted conditional mean temperature evolution, if the average conditional scalar dissipation rate is properly modeled. Using DNS data, a correction factor is introduced in the modeling of nonlinear terms to include the effect of species fluctuations. Computations including such a correction factor show that the species conditional correlations q(ij) have little effect on model predictions with a one-step reaction, but those q(ij) involving intermediate species are found to be crucial when four-step reduced kinetics is considered. The "most reactive mixture fraction" is found to vary with time when a four-step kinetics is considered. First-order CMC results are found to be qualitatively wrong if the conditional mean scalar dissipation rate is not modeled properly. The autoignition delay time predicted by the CMC model compares excellently with DNS results and shows a trend similar to experimental data over a range of initial temperatures.

Non-Fermi-liquid theory for disordered metals near two dimensions

We consider the Finkelstein action describing a system of spin-polarized or spinless electrons in 2+2epsilon dimensions, in the presence of disorder as well as the Coulomb interactions. We extend the renormalization-group analysis of our previous work and evaluate the metal-insulator transition of the electron gas to second order in an epsilon expansion. We obtain the complete scaling behavior of physical observables like the conductivity and the specific heat with varying frequency, temperature, and/or electron density. We extend the results for the interacting electron gas in 2+2epsilon dimensions to include the quantum critical behavior of the plateau transitions in the quantum Hall regime. Although these transitions have a very different microscopic origin and are controlled by a topological term in the action (theta term), the quantum critical behavior is in many ways the same in both cases. We show that the two independent critical exponents of the quantum Hall plateau transitions, previously denoted as nu and p, control not only the scaling behavior of the conductances sigma(xx) and sigma(xy) at finite temperatures T, but also the non-Fermi-liquid behavior of the specific heat (c(v)proportional toT(p)). To extract the numerical values of nu and p it is necessary to extend the experiments on transport to include the specific heat of the electron gas.

Optimum design of a helicopter rotor for low vibration using aeroelastic analysis and response surface methods

An aeroelastic analysis based on finite elements in space and time is used to model the helicopter rotor in forward flight. The rotor blade is represented as an elastic cantilever beam undergoing flap and lag bending, elastic torsion and axial deformations. The objective of the improved design is to reduce vibratory loads at the rotor hub that are the main source of helicopter vibration. Constraints are imposed on aeroelastic stability, and move limits are imposed on the blade elastic stiffness design variables. Using the aeroelastic analysis, response surface approximations are constructed for the objective function (vibratory hub loads). It is found that second order polynomial response surfaces constructed using the central composite design of the theory of design of experiments adequately represents the aeroelastic model in the vicinity of the baseline design. Optimization results show a reduction in the objective function of about 30 per cent. A key accomplishment of this paper is the decoupling of the analysis problem and the optimization problems using response surface methods, which should encourage the use of optimization methods by the helicopter industry. (C) 2002 Elsevier Science Ltd. All rights reserved.

A differential-geometric analysis of singularities of point trajectories of serial and parallel manipulators

In this paper, we present a differential-geometric approach to analyze the singularities of task space point trajectories of two and three-degree-of-freedom serial and parallel manipulators. At non-singular configurations, the first-order, local properties are characterized by metric coefficients, and, geometrically, by the shape and size of a velocity ellipse or an ellipsoid. At singular configurations, the determinant of the matrix of metric coefficients is zero and the velocity ellipsoid degenerates to an ellipse, a line or a point, and the area or the volume of the velocity ellipse or ellipsoid becomes zero. The degeneracies of the velocity ellipsoid or ellipse gives a simple geometric picture of the possible task space velocities at a singular configuration. To study the second-order properties at a singularity, we use the derivatives of the metric coefficients and the rate of change of area or volume. The derivatives are shown to be related to the possible task space accelerations at a singular configuration. In the case of parallel manipulators, singularities may lead to either loss or gain of one or more degrees-of-freedom. For loss of one or more degrees-of-freedom, ther possible velocities and accelerations are again obtained from a modified metric and derivatives of the metric coefficients. In the case of a gain of one or more degrees-of-freedom, the possible task space velocities can be pictured as growth to lines, ellipses, and ellipsoids. The theoretical results are illustrated with the help of a general spatial 2R manipulator and a three-degree-of-freedom RPSSPR-SPR parallel manipulator.

System Dy-Fe-O: thermodynamic properties of ternary oxides using Calvet calorimetry and solid-state electrochemical cell

The enthalpy increments and the standard molar Gibbs energies of formation-of DyFeO3(s) and Dy3Fe5O12(s) have been measured using a Calvet micro-calorimeter and a solid oxide galvanic cell, respectively. A co-operative phase transition, related to anti-ferromagnetic to paramagnetic transformation, is apparent. from the heat capacity data for DyFeO3 at similar to 648 K. A similar type of phase transition has been observed for Dy3Fe5O12 at similar to 560 K which is related to ferrimagnetic to paramagnetic transformation. Enthalpy increment data for DyFeO3(s) and Dy3Fe5O12(s), except in the vicinity of the second-order transition, can be represented by the following polynomial expressions:{H(0)m(T) - H(0)m(298.15 K)) (Jmol(-1)) (+/-1.1%) = -52754 + 142.9 x (T (K)) + 2.48 x 10(-3) x (T (K))(2) + 2.951 x 10(6) x (T (K))(-1); (298.15 less than or equal to T (K) less than or equal to 1000) for DyFeO3(s), and {H(0)m(T) - H(0)m(298.15 K)} (Jmol(-1)) (+/-1.2%) = -191048 + 545.0 x (T - (K)) + 2.0 x 10(-5) x (T (K))(2) + 8.513 x 10(6) x (T (K))(-1); (208.15 less than or equal to T (K) less than or equal to 1000)for Dy3Fe5O12(s). The reversible emfs of the solid-state electrochemical cells: (-)Pt/{DyFeO3(s) + Dy2O3(s) + Fe(s)}/YDT/CSZ//{Fe(s) + Fe0.95O(s)}/Pt(+) and (-)Pt/{Fe(s) + Fe0.95O(s)}//CSZ//{DyFeO3(s) + Dy3Fe5O12(s) + Fe3O4(s)}/Pt(+), were measured in the temperature range from 1021 to 1250 K and 1035 to 1250 K, respectively. The standard Gibbs energies of formation of solid DyFeO3 and Dy3Fe5O12 calculated by the least squares regression analysis of the data obtained in the present study, and data for Fe0.95O and Dy2O3 from the literature, are given by Delta(f)G(0)m(DyFeO3,s)(kJmol(-1))(+/-3.2)= -1339.9 + 0.2473 x (T(K)); (1021 less than or equal to T (K) less than or equal to 1548)and D(f)G(0)m(Dy3Fe5O12,s) (kJmol(-1)) (+/-3.5) = -4850.4 + 0.9846 x (T (K)); (1035 less than or equal to T (K) less than or equal to 1250) The uncertainty estimates for Delta(f)G(0)m include the standard deviation in the emf and uncertainty in the data taken from the literature. Based on the thermodynamic information, oxygen potential diagram and chemical potential diagrams for the system Dy-Fe-O were developed at 1250 K. (C) 2002 Editions scientifiques et medicales Elsevier SAS. All rights reserved.

Direct numerical simulation of turbulent spots

A group of high-order finite-difference schemes for incompressible flow was implemented to simulate the evolution of turbulent spots in channel flows. The long-time accuracy of these schemes was tested by comparing the evolution of small disturbances to a plane channel flow against the growth rate predicted by linear theory. When the perturbation is the unstable eigenfunction at a Reynolds number of 7500, the solution grows only if there are a comparatively large number of (equispaced) grid points across the channel. Fifth-order upwind biasing of convection terms is found to be worse than second-order central differencing. But, for a decaying mode at a Reynolds number of 1000, about a fourth of the points suffice to obtain the correct decay rate. We show that this is due to the comparatively high gradients in the unstable eigenfunction near the walls. So, high-wave-number dissipation of the high-order upwind biasing degrades the solution especially. But for a well-resolved calculation, the weak dissipation does not degrade solutions even over the very long times (O(100)) computed in these tests. Some new solutions of spot evolution in Couette flows with pressure gradients are presented. The approach to self-similarity at long times can be seen readily in contour plots.

Cavitation in holographic sQGP

We study the possibility of cavitation in the non-conformal N = 2* SU(N) theory which is a mass deformation of N = 4 SU(N) Yang-Mills theory. The second order transport coefficients are known from the numerical work using AdS/CFT by Buchel and collaborators. Using these and the approach of Rajagopal and Tripuraneni, we investigate the flow equations in a (1 + 1)-dimensional boost invariant set up. We find that the string theory model does not exhibit cavitation before phase transition is reached. We give a semi-analytic explanation of this finding. (C) 2011 Elsevier B.V. All rights reserved.

Dielectric Properties of Pulsed Excimer Laser Ablated BaBi2Nb2O9 Thin Films

The dielectric response of BaBi2Nb2O9 (BBN) thin films has been studied as a function of frequency over a wide range of temperatures. Both dielectric constant and loss tangent of BBN thin films showed a ‘power law’ dependence with frequency, which was analyzed using the Jonscher's universal dielectric response model. Theoretical fits were utilized to compare the experimental results and also to estimate the value of temperature dependence parameters such as n(T) and a(T) used in the Jonscher's model. The room temperature dielectric constant (ε') of the BBN thin films was 214 with a loss tangent (tanδ) of 0.04 at a frequency of 100 kHz. The films exhibited the second order dielectric phase transition from ferroelectric to paraelectric state at a temperature of 220 °C. The nature of phase transition was confirmed from the temperature dependence of dielectric constant and sponteneous polarization,respectively. The calculated Currie constant for BBN thin films was 4 × 105°C.

Turbulent flow computations using upwind RANS technique for a flat plate and an airfoil

Reynolds Averaged Navier Stokes (RANS) equations are solved using third order upwind biased Roe's scheme for the inviscid fluxes and second order central difference scheme for the viscous fluxes. The Baldwin & Lomax turbulence model is employed for Reynolds stresses. The governing equations are solved using finite-volume implicit scheme in body fitted curvilinear coordinate O-grid system. Computations axe reported for a flat plate apart from RAE 2822 and NACA 0012 airfoils. Results for the flat plate at M = 0.3, R-c = 4.0 x 10(6) compare favourably with the analytical solution. Results for the two airfoils are compared with experiment. There is a good agreement in C-p distribution between experiment and computation for both the airfoils. Comparison of C-f distribution with experiment for RAE 2822 airfoil is reasonable.

Acquisition time improvement of PLLs using some aiding functions

Analytical studies are carried out to minimize acquisition time in phase-lock loop (PLL) applications using aiding functions. A second order aided PLL is realized with the help of the quasi-stationary approach to verify the acquisition behavior in the absence of noise. Time acquisition is measured both from the study of the LPF output transient and by employing a lock detecting and indicating circuit to crosscheck experimental and analytical results. A closed form solution is obtained for the evaluation of the time acquisition using different aiding functions. The aiding signal is simple and economical and can be used with state of the art hardware.

Robust Aeroelastic Optimization Of Composite Helicopter Rotor

A robust aeroelastic optimization is performed to minimize helicopter vibration with uncertainties in the design variables. Polynomial response surfaces and space-¯lling experimental designs are used to generate the surrogate model of aeroelastic analysis code. Aeroelastic simulations are performed at the sample inputs generated by Latin hypercube sampling. The response values which does not satisfy the frequency constraints are eliminated from the data for model ¯tting. This step increased the accuracy of response surface models in the feasible design space. It is found that the response surface models are able to capture the robust optimal regions of design space. The optimal designs show a reduction of 10 percent in the objective function comprising six vibratory hub loads and 1.5 to 80 percent reduction for the individual vibratory forces and moments. This study demonstrates that the second-order response surface models with space ¯lling-designs can be a favorable choice for computationally intensive robust aeroelastic optimization.