Sidewall inclined slot in a rectangular waveguide: Theory and experiment

A novel analysis to compute the admittance characteristics of the slots cut in the narrow wall of a rectangular waveguide, which includes the corner diffraction effects and the finite waveguide wall thickness, is presented. A coupled magnetic field integral equation is formulated at the slot aperture which is solved by the Galerkin approach of the method of moments using entire domain sinusoidal basis functions. The externally scattered fields are computed using the finite difference method (FDM) coupled with the measured equation of invariance (MEI). The guide wall thickness forms a closed cavity and the fields inside it are evaluated using the standard FDM. The fields scattered inside the waveguide are formulated in the spectral domain for faster convergence compared to the traditional spatial domain expansions. The computed results have been compared with the experimental results and also with the measured data published in previous literature. Good agreement between the theoretical and experimental results is obtained to demonstrate the validity of the present analysis.

Some limit theorems for regenerative queues

We consider a single server queue with the interarrival times and the service times forming a regenerative sequence. This traffic class includes the standard models: lid, periodic, Markov modulated (e.g., BMAP model of Lucantoni [18]) and their superpositions. This class also includes the recently proposed traffic models in high speed networks, exhibiting long range dependence. Under minimal conditions we obtain the rates of convergence to stationary distributions, finiteness of stationary moments, various functional limit theorems and the continuity of stationary distributions and moments. We use the continuity results to obtain approximations for stationary distributions and moments of an MMPP/GI/1 queue where the modulating chain has a countable state space. We extend all our results to feedforward networks where the external arrivals to each queue can be regenerative. In the end we show that the output process of a leaky bucket is regenerative if the input process is and hence our results extend to a queue with arrivals controlled by a leaky bucket.

A Quasi-Newton Adaptive Algorithm for Generalized Symmetric Eigenvalue Problem

In this paper, we first recast the generalized symmetric eigenvalue problem, where the underlying matrix pencil consists of symmetric positive definite matrices, into an unconstrained minimization problem by constructing an appropriate cost function, We then extend it to the case of multiple eigenvectors using an inflation technique, Based on this asymptotic formulation, we derive a quasi-Newton-based adaptive algorithm for estimating the required generalized eigenvectors in the data case. The resulting algorithm is modular and parallel, and it is globally convergent with probability one, We also analyze the effect of inexact inflation on the convergence of this algorithm and that of inexact knowledge of one of the matrices (in the pencil) on the resulting eigenstructure. Simulation results demonstrate that the performance of this algorithm is almost identical to that of the rank-one updating algorithm of Karasalo. Further, the performance of the proposed algorithm has been found to remain stable even over 1 million updates without suffering from any error accumulation problems.

Hybrid stiff-string-polynomial basis functions for vibration analysis of high speed rotating beams

A new finite element is developed for free vibration analysis of high speed rotating beams using basis functions which use a linear combination of the solution of the governing static differential equation of a stiff-string and a cubic polynomial. These new shape functions depend on rotation speed and element position along the beam and account for the centrifugal stiffening effect. The natural frequencies predicted by the proposed element are compared with an element with stiff-string, cubic polynomial and quintic polynomial shape functions. It is found that the new element exhibits superior convergence compared to the other basis functions.

Hingeless-rotor aeromechanical stability in axial and forward flight with wake dynamics

A finite-state wake model is used to investigate aeromechanical stability of hingeless-rotor helicopters in the ground-contact, hover and trimmed-night conditions. The investigation covers three items: (1) the convergence of the damping with increasing number of wake harmonics for the lag regressing, and body pitch and roll modes; (2) a parametric study of the damping over a range of thrust level, advance ratio and number of blades; and (3) correlations, primarily with the damping and frequency measurements of these lag and body modes. The convergence and parametric studies are conducted in the hover and trimmed-flight conditions; they include predictions from the widely used dynamic inflow model. The correlations are conducted in the ground-contact conditions and include predictions from the dynamic inflow and vortex models; recently, this vortex model is proposed for the axial-flight conditions and is used to investigate the coupled free vibrations of rotor flapping and body modes. The convergence and parametric studies show that a finite-state wake model that goes well beyond the dynamic inflow model is required for fairly converged damping, Moreover, the correlations from the finite-state wake, dynamic inflow and vortex models are generally satisfactory.

Convergent evolution of marine mammals is associated with distinct substitutions in common genes

Phenotypic convergence is thought to be driven by parallel substitutions coupled with natural selection at the sequence level. Multiple independent evolutionary transitions of mammals to an aquatic environment offer an opportunity to test this thesis. Here, whole genome alignment of coding sequences identified widespread parallel amino acid substitutions in marine mammals; however, the majority of these changes were not unique to these animals. Conversely, we report that candidate aquatic adaptation genes, identified by signatures of likelihood convergence and/or elevated ratio of nonsynonymous to synonymous nucleotide substitution rate, are characterized by very few parallel substitutions and exhibit distinct sequence changes in each group. Moreover, no significant positive correlation was found between likelihood convergence and positive selection in all three marine lineages. These results suggest that convergence in protein coding genes associated with aquatic lifestyle is mainly characterized by independent substitutions and relaxed negative selection.

Reliable estimation via simulation

Let a and s denote the inter arrival times and service times in a GI/GI/1 queue. Let a (n), s (n) be the r.v.s, with distributions as the estimated distributions of a and s from iid samples of a and s of sizes n. Let w be a r.v. with the stationary distribution lr of the waiting times of the queue with input (a, s). We consider the problem of estimating E [w~], tx > 0 and 7r via simulations when (a (n), s (n)) are used as input. Conditions for the accuracy of the asymptotic estimate, continuity of the asymptotic variance and uniformity in the rate of convergence to the estimate are obtained. We also obtain rates of convergence for sample moments, the empirical process and the quantile process for the regenerative processes. Robust estimates are also obtained when an outlier contaminated sample of a and s is provided. In the process we obtain consistency, continuity and asymptotic normality of M-estimators for stationary sequences. Some robustness results for Markov processes are included.

Implicit and multigrid procedures for steady-state computations with upwind algorithms

A computational study for the convergence acceleration of Euler and Navier-Stokes computations with upwind schemes has been conducted in a unified framework. It involves the flux-vector splitting algorithms due to Steger-Warming and Van Leer, the flux-difference splitting algorithms due to Roe and Osher and the hybrid algorithms, AUSM (Advection Upstream Splitting Method) and HUS (Hybrid Upwind Splitting). Implicit time integration with line Gauss-Seidel relaxation and multigrid are among the procedures which have been systematically investigated on an individual as well as cumulative basis. The upwind schemes have been tested in various implicit-explicit operator combinations such that the optimal among them can be determined based on extensive computations for two-dimensional flows in subsonic, transonic, supersonic and hypersonic flow regimes. In this study, the performance of these implicit time-integration procedures has been systematically compared with those corresponding to a multigrid accelerated explicit Runge-Kutta method. It has been demonstrated that a multigrid method employed in conjunction with an implicit time-integration scheme yields distinctly superior convergence as compared to those associated with either of the acceleration procedures provided that effective smoothers, which have been identified in this investigation, are prescribed in the implicit operator.

A Jacobian-Based Algorithm for Planning Attitude Maneuvers Using Forward and Reverse Rotations

Algorithms for planning quasistatic attitude maneuvers based on the Jacobian of the forward kinematic mapping of fully-reversed (FR) sequences of rotations are proposed in this paper. An FR sequence of rotations is a series of finite rotations that consists of initial rotations about the axes of a body-fixed coordinate frame and subsequent rotations that undo these initial rotations. Unlike the Jacobian of conventional systems such as a robot manipulator, the Jacobian of the system manipulated through FR rotations is a null matrix at the identity, which leads to a total breakdown of the traditional Jacobian formulation. Therefore, the Jacobian algorithm is reformulated and implemented so as to synthesize an FR sequence for a desired rotational displacement. The Jacobian-based algorithm presented in this paper identifies particular six-rotation FR sequences that synthesize desired orientations. We developed the single-step and the multiple-step Jacobian methods to accomplish a given task using six-rotation FR sequences. The single-step Jacobian method identifies a specific FR sequence for a given desired orientation and the multiple-step Jacobian algorithm synthesizes physically feasible FR rotations on an optimal path. A comparison with existing algorithms verifies the fast convergence ability of the Jacobian-based algorithm. Unlike closed-form solutions to the inverse kinematics problem, the Jacobian-based algorithm determines the most efficient FR sequence that yields a desired rotational displacement through a simple and inexpensive numerical calculation. The procedure presented here is useful for those motion planning problems wherein the Jacobian is singular or null.

Temporal dependence of in vivo USPIO-enhanced MRI signal changes in human carotid atheromatous plaques

Introduction: Ultrasmall superparamagnetic iron oxide (USPIO)-enhanced MRI has been shown to be a useful modality to image activated macrophages in vivo, which are principally responsible for plaque inflammation. This study determined the optimum imaging time-window to detect maximal signal change post-USPIO infusion using T1-weighted (T1w), T2*- weighted (T2*w) and quantitative T2*(qT 2*) imaging. Methods: Six patients with an asymptomatic carotid stenosis underwent high resolution T1w, T2*w and qT2*MR imaging of their carotid arteries at 1.5 T. Imaging was performed before and at 24, 36, 48, 72 and 96 h after USPIO (Sinerem™, Guerbet, France) infusion. Each slice showing atherosclerotic plaque was manually segmented into quadrants and signal changes in each quadrant were fitted to an exponential power function to model the optimum time for post-infusion imaging. Results: The power function determining the mean time to convergence for all patients was 46, 41 and 39 h for the T1w, T 2*w and qT2*sequences, respectively. When modelling each patient individually, 90% of the maximum signal intensity change was observed at 36 h for three, four and six patients on T1w, T 2*w and qT2*, respectively. The rates of signal change decrease after this period but signal change was still evident up to 96 h. Conclusion: This study showed that a suitable imaging window for T 1w, T2*w and qT2*signal changes post-USPIO infusion was between 36 and 48 h. Logistically, this would be convenient in bringing patients back for one post-contrast MRI, but validation is required in a larger cohort of patients.

A Doubly curved Quadrilateral finite-element for the analysis of laminated anisotropic thin shells of revoution

The details of development of the stiffness matrix for a doubly curved quadrilateral element suited for static and dynamic analysis of laminated anisotropic thin shells of revolution are reported. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first order Hermite polynomials, it is possible to ensure that the displacement state for the assembled set of such elements, is geometrically admissible. Monotonic convergence of total potential energy is therefore possible as the modelling is successively refined. Systematic evaluation of performance of the element is conducted, considering various examples for which analytical or other solutions are available.

Ocean atmosphere coupling over monsoon regions

In monsoon regions, the seasonal migration of the intertropical convergence zone (ITCZ) is manifested as a seasonal reversal of winds. Most of the summer monsoon rainfall over India occurs owing to synoptic and large-scale convection associated with the continental ITCZ (Fig. 1). We have investigated the interaction between these large-scale convective systems and the ocean over which they are generated1â3, concentrating on the relationship between organized convection over the Indian Ocean and sea surface temperature (SST). We report here that on a monthly basis the degree of cloudiness correlates well with SST for the relatively colder oceans, but when SST is maintained above 28 °C it ceases to be an important factor in determining the variability of cloudiness. Over the major regions of convection east of 70°E, which are warm year after year, the observed cloudiness cannot be correlated with variations in SST.

Standard errors and covariance matrices for smoothed rank estimators

A 'pseudo-Bayesian' interpretation of standard errors yields a natural induced smoothing of statistical estimating functions. When applied to rank estimation, the lack of smoothness which prevents standard error estimation is remedied. Efficiency and robustness are preserved, while the smoothed estimation has excellent computational properties. In particular, convergence of the iterative equation for standard error is fast, and standard error calculation becomes asymptotically a one-step procedure. This property also extends to covariance matrix calculation for rank estimates in multi-parameter problems. Examples, and some simple explanations, are given.

A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

Mesh-free approximations via the error reproducing kernel method and applications to nonlinear systems developing shocks

Error estimates for the error reproducing kernel method (ERKM) are provided. The ERKM is a mesh-free functional approximation scheme [A. Shaw, D. Roy, A NURBS-based error reproducing kernel method with applications in solid mechanics, Computational Mechanics (2006), to appear (available online)], wherein a targeted function and its derivatives are first approximated via non-uniform rational B-splines (NURBS) basis function. Errors in the NURBS approximation are then reproduced via a family of non-NURBS basis functions, constructed using a polynomial reproduction condition, and added to the NURBS approximation of the function obtained in the first step. In addition to the derivation of error estimates, convergence studies are undertaken for a couple of test boundary value problems with known exact solutions. The ERKM is next applied to a one-dimensional Burgers equation where, time evolution leads to a breakdown of the continuous solution and the appearance of a shock. Many available mesh-free schemes appear to be unable to capture this shock without numerical instability. However, given that any desired order of continuity is achievable through NURBS approximations, the ERKM can even accurately approximate functions with discontinuous derivatives. Moreover, due to the variation diminishing property of NURBS, it has advantages in representing sharp changes in gradients. This paper is focused on demonstrating this ability of ERKM via some numerical examples. Comparisons of some of the results with those via the standard form of the reproducing kernel particle method (RKPM) demonstrate the relative numerical advantages and accuracy of the ERKM.