Computing with functions in spherical and polar geometries I. The sphere

A collection of algorithms is described for numerically computing with smooth functions defined on the unit sphere. Functions are approximated to essentially machine precision by using a structure-preserving iterative variant of Gaussian elimination together with the double Fourier sphere method. We show that this procedure allows for stable differentiation, reduces the oversampling of functions near the poles, and converges for certain analytic functions. Operations such as function evaluation, differentiation, and integration are particularly efficient and can be computed by essentially one-dimensional algorithms. A highlight is an optimal complexity direct solver for Poisson's equation on the sphere using a spectral method. Without parallelization, we solve Poisson's equation with $100$ million degrees of freedom in one minute on a standard laptop. Numerical results are presented throughout. In a companion paper (part II) we extend the ideas presented here to computing with functions on the disk.Comment: 23 page

Cost of Forward Contracting Hard Red Winter Wheat

Marketing,

COST OF FORWARD CONTRACTING HARD RED WINTER WHEAT

Two methods were used to estimate the cost of forward contracting hard red winter wheat. One hundred days before delivery, the estimated cost of forward contracting ranged from six cents/bu. To eight cents/bu. Thus, further evidence is provided that the cost of forward contracting grain is not zero.forward contracting, nonparametric regression, wheat, Marketing,

Dynamic measurements of gear tooth friction and load

As part of a program to study fundamental mechanisms of gear noise, static and dynamic gear tooth strain measurements were made on the NASA gear-noise rig. Tooth-fillet strains from low-contact ratio-spur gears were recorded for 28 operating conditions. A method is introduced whereby strain gage measurements taken from both the tension and compression sides of a gear tooth can be transformed into the normal and frictional loads on the tooth. This technique was applied to both the static and dynamic strain data. The static case results showed close agreement with expected results. For the dynamic case, the normal-force computation produced very good results, but the friction results, although promising, were not as accurate. Tooth sliding friction strongly affected the signal from the strain gage on the tensionside of the tooth. The compression gage was affected by friction to a much lesser degree. The potential of the method to measure friction force was demonstrated, but further refinement will be required before this technique can be used to measure friction forces dynamically with an acceptable degree of accuracy

Nucleon and deuteron scattering cross sections from 25 MV/Nucleon to 22.5 GeV/Nucleon

Within the context of a double-folding optical potential approximation to the exact nucleus-nucleus multiple-scattering series, eikonal scattering theory is used to generate tables of nucleon and deuteron total and absorption cross sections at kinetic energies between 25 MeV/nucleon and 22.5 GeV/nucleon for use in cosmic-ray transport and shielding studies. Comparisons of predictions for nucleon-nucleus and deuteron-nucleus absorption and total cross sections with experimental data are also made

A comparison between theoretical prediction and experimental measurement of the dynamic behavior of spur gears

A comparison was made between computer model predictions of gear dynamics behavior and experimental results. The experimental data were derived from the NASA gear noise rig, which was used to record dynamic tooth loads and vibration. The experimental results were compared with predictions from the DSTO Aeronautical Research Laboratory's gear dynamics code for a matrix of 28 load speed points. At high torque the peak dynamic load predictions agree with the experimental results with an average error of 5 percent in the speed range 800 to 6000 rpm. Tooth separation (or bounce), which was observed in the experimental data for light torque, high speed conditions, was simulated by the computer model. The model was also successful in simulating the degree of load sharing between gear teeth in the multiple tooth contact region

Heavy-ion total and absorption cross sections above 25 MeV/nucleon

Within the context of a double-folding optical potential approximation to the exact nucleus-nucleus multiple-scattering series, eikonal scattering theory is used to generate tables of heavy ion total and absorption cross sections at incident kinetic energies above 25 MeV/nucleon for use in cosmic ray high-energy heavy ion transport and shielding studies. Comparisons of predictions with nucleus-nucleus experimental data show excellent agreement except at the lowest energies, where the eikonal approximation may not be completely valid. Even at the lowest energies, however, agreement is typically within 20 percent

Phenomenological optical potential analysis of proton-carbon elastic scattering at 200 MeV

Differential cross sections for 200 MeV protons elastically scattered from C-12 were analyzed utilizing a local, complex, spin-dependent optical potential with a harmonic well radial dependence. Analyses were performed using the WKB and eikonal approximations. For the latter, first-order corrections to he phase shifts were incorporated to account for the spin-orbit contribution. Large disagreement between theory and experiment was observed when the usual Thomas form for the spin-orbit potential was utilized. Substantial improvement was obtained by allowing the parameters in the central and spin-orbit potential terms to vary independently

Analytic determinations of single-folding optical potentials

A simple analytic method for calculating nucleon-nucleus optical potentials using a single folding of a Gaussian two body interaction with an arbitrary nuclear distribution is presented. When applied to proton-lead elastic scattering, the predicted real part of the Woods-Saxon potential is in substantial agreement with the experimentally determined phenomenological potential, although there are no adjustable parameters. In addition, the volume integrals of both real potentials are nearly identical

Membrane solitons in eight-dimensional hyper-Kaehler backgrounds

We derive the BPS equations satisfied by lump solitons in $(2+1)$-dimensional sigma models with toric 8-dimensional hyper-K\"ahler (${HK}_8$) target spaces and check they preserve 1/2 of the supersymmetry. We show how these solitons are realised in M theory as M2-branes wrapping holomorphic 2-cycles in the \bE^{1,2}\times {HK}_8 background. Using the $\kappa$-symmetry of a probe M2-brane in this background we determine the supersymmetry they preserve, and note that there is a discrepancy in the fraction of supersymmetry preserved by these solitons as viewed from the low energy effective sigma model description of the M2-brane dynamics or the full M theory. Toric ${HK}_8$ manifolds are dual to a Hanany-Witten setup of D3-branes suspended between 5-branes. In this picture the lumps correspond to vortices of the three dimensional ${\mathcal N}=3$ or ${\mathcal N}=4$ theory.Comment: 12+1 pages. LaTex. v2: Typos corrected and references adde