A Nonstationary Negative Binomial Time Series with Time-Dependent Covariates: Enterococcus Counts in Boston Harbor

Boston Harbor has had a history of poor water quality, including contamination by enteric pathogens. We conduct a statistical analysis of data collected by the Massachusetts Water Resources Authority (MWRA) between 1996 and 2002 to evaluate the effects of court-mandated improvements in sewage treatment. Motivated by the ineffectiveness of standard Poisson mixture models and their zero-inflated counterparts, we propose a new negative binomial model for time series of Enterococcus counts in Boston Harbor, where nonstationarity and autocorrelation are modeled using a nonparametric smooth function of time in the predictor. Without further restrictions, this function is not identifiable in the presence of time-dependent covariates; consequently we use a basis orthogonal to the space spanned by the covariates and use penalized quasi-likelihood (PQL) for estimation. We conclude that Enterococcus counts were greatly reduced near the Nut Island Treatment Plant (NITP) outfalls following the transfer of wastewaters from NITP to the Deer Island Treatment Plant (DITP) and that the transfer of wastewaters from Boston Harbor to the offshore diffusers in Massachusetts Bay reduced the Enterococcus counts near the DITP outfalls.

Linear Approach to the Orbiting Spacecraft Thermal Problem

A linear method is developed for solving the nonlinear differential equations of a lumped-parameter thermal model of a spacecraft moving in a closed orbit. This method, based on perturbation theory, is compared with heuristic linearizations of the same equations. The essential feature of the linear approach is that it provides a decomposition in thermal modes, like the decomposition of mechanical vibrations in normal modes. The stationary periodic solution of the linear equations can be alternately expressed as an explicit integral or as a Fourier series. This method is applied to a minimal thermal model of a satellite with ten isothermal parts (nodes), and the method is compared with direct numerical integration of the nonlinear equations. The computational complexity of this method is briefly studied for general thermal models of orbiting spacecraft, and it is concluded that it is certainly useful for reduced models and conceptual design but it can also be more efficient than the direct integration of the equations for large models. The results of the Fourier series computations for the ten-node satellite model show that the periodic solution at the second perturbative order is sufficiently accurate.

Fourier analysis of the aerodynamic behavior of cup anemometers

The calibration results (the transfer function) of an anemometer equipped with several cup rotors were analyzed and correlated with the aerodynamic forces measured on the isolated cups in a wind tunnel. The correlation was based on a Fourier analysis of the normal-to-the-cup aerodynamic force. Three different cup shapes were studied: typical conical cups, elliptical cups and porous cups (conical-truncated shape). Results indicated a good correlation between the anemometer factor, K, and the ratio between the first two coefficients in the Fourier series decomposition of the normal-to-the-cup aerodynamic force

Determination of the distance between two spin labels attached to a macromolecule.

An EPR "spectroscopic ruler" was developed using a series of alpha-helical polypeptides, each modified with two nitroxide spin labels. The EPR line broadening due to electron-electron dipolar interactions in the frozen state was determined using the Fourier deconvolution method. These dipolar spectra were then used to estimate the distances between the two nitroxides separated by 8-25 A. Results agreed well with a simple alpha-helical model. The standard deviation from the model system was 0.9 A in the range of 8-25 A. This technique is applicable to complex systems such as membrane receptors and channels, which are difficult to access with high-resolution NMR or x-ray crystallography, and is expected to be particularly useful for systems for which optical methods are hampered by the presence of light-interfering membranes or chromophores.

Serie di Fourier e altre rappresentazioni analitiche delle funzioni di una variabile reale.

Available on demand as hard copy or computer file from Cornell University Library.

A. L. Cauchy's Vorlesungen über die Differenzialrechnung : mit Fourier's Auflösungsmethode der bestimmten Gleichungen verbunden /

Mode of access: Internet.

An Analytical Solution for the Elastic Lateral-Distortional Buckling of I-section Beams

Lateral-distortional buckling may occur in I-section beams with slender webs and stocky flanges. A computationally efficient method is presented in this paper to study this phenomenon. Previous studies on distortional buckling have been on the use of 3(rd) and 5(th) order polynomials to model the displacements. The present study provides an alternative way, using Fourier Series, to model the behaviour. Beams of different cross-sectional dimensions, load cases and restraint conditions are examined and compared. The accuracy and versatility of the method are verified by calibrating against the results of other published studies. The present method is believed to be a simple and efficient way of determining the buckling load and mode shapes of I-section beams that are susceptible to lateral-distortional buckling modes.

Sigmoidal cosine series on the interval

We construct a set of functions, say, psi([r])(n) composed of a cosine function and a sigmoidal transformation gamma(r) of order r > 0. The present functions are orthonormal with respect to a proper weight function on the interval [-1, 1]. It is proven that if a function f is continuous and piecewise smooth on [-1, 1] then its series expansion based on psi([r])(n) converges uniformly to f so long as the order of the sigmoidal transformation employed is 0 < r

A New Algorithm for Monte Carlo for American Options

2000 Mathematics Subject Classification: 91B28, 65C05.

Hybrid formulation solution for stress analysis of curved pipes with welded bending joints

A method is presented for evaluating the stress intensity factor of part-through cracks in a thin pipe elbow. A hybrid formulation solution is used to evaluate the stress field close to the crack area. The stress field values are then inputted into a previously developed method published in the literature to evaluate the stress intensity factor in cylindrical shells. Results from cylindrical shells with part-through cracks are extended to double-curvature pipe configurations that contain the same kind of flaw.

Complete spin and orbital evolution of close-in bodies using a Maxwell viscoelastic rheology

In this paper, we present a formalism designed to model tidal interaction with a viscoelastic body made of Maxwell material. Our approach remains regular for any spin rate and orientation, and for any orbital configuration including high eccentricities and close encounters. The method is to integrate simultaneously the rotation and the position of the planet as well as its deformation. We provide the equations of motion both in the body frame and in the inertial frame. With this study, we generalize preexisting models to the spatial case and to arbitrary multipole orders using a formalism taken from quantum theory. We also provide the vectorial expression of the secular tidal torque expanded in Fourier series. Applying this model to close-in exoplanets, we observe that if the relaxation time is longer than the revolution period, the phase space of the system is characterized by the presence of several spin-orbit resonances, even in the circular case. As the system evolves, the planet spin can visit different spin-orbit configurations. The obliquity is decreasing along most of these resonances, but we observe a case where the planet tilt is instead growing. These conclusions derived from the secular torque are successfully tested with numerical integrations of the instantaneous equations of motion on HD 80606 b. Our formalism is also well adapted to close-in super-Earths in multiplanet systems which are known to have non-zero mutual inclinations.

Voltage and current spectra for matrix power converters

Matrix power converters are used for transforming one alternating-current power supply to another, with different peak voltage and frequency. There are three input lines, with sinusoidally varying voltages which are 120◦ out of phase one from another, and the output is to be delivered as a similar three-phase supply. The matrix converter switches rapidly, to connect each output line in sequence to each of the input lines in an attempt to synthesize the prescribed output voltages. The switching is carried out at high frequency and it is of practical importance to know the frequency spectra of the output voltages and of the input and output currents. We determine in this paper these spectra using a new method, which has significant advantages over the prior default method (a multiple Fourier series technique), leading to a considerably more direct calculation. In particular, the determination of the input current spectrum is feasible here, whereas it would be a significantly more daunting procedure using the prior method instead.

Large-Eddy Simulation of physiological pulsatile non-newtonian flow in a channel with double constriction

Physiological pulsatile flow in a 3D model of arterial double stenosis, using the modified Power-law blood viscosity model, is investigated by applying Large Eddy Simulation (LES) technique. The computational domain has been chosen is a simple channel with biological type stenoses. The physiological pulsation is generated at the inlet of the model using the first four harmonics of the Fourier series of the physiological pressure pulse. In LES, a top-hat spatial grid-filter is applied to the Navier-Stokes equations of motion to separate the large scale flows from the subgrid scale (SGS). The large scale flows are then resolved fully while the unresolved SGS motions are modelled using the localized dynamic model. The flow Reynolds numbers which are typical of those found in human large artery are chosen in the present work. Transitions to turbulent of the pulsatile non-Newtonian along with Newtonian flow in the post stenosis are examined through the mean velocity, wall shear stress, mean streamlines as well as turbulent kinetic energy and explained physically along with the relevant medical concerns.

Bayesian models for the determination of resonant frequencies in a DI diesel engine

A time series method for the determination of combustion chamber resonant frequencies is outlined. This technique employs the use of Markov-chain Monte Carlo (MCMC) to infer parameters in a chosen model of the data. The development of the model is included and the resonant frequency is characterised as a function of time. Potential applications for cycle-by-cycle analysis are discussed and the bulk temperature of the gas and the trapped mass in the combustion chamber are evaluated as a function of time from resonant frequency information.