BOOTSTRAPPING YOUR FISH OR FISHING FOR BOOTSTRAPS?: PRECISION OF WELFARE LOSS ESTIMATES FROM A GLOBALLY CONCAVE INVERSE DEMAND MODEL OF COMMERCIAL FISH LANDINGS IN THE U.S. GREAT LAKES

Replaced with revised version of paper 06/30/04.Demand and Price Analysis, Resource /Energy Economics and Policy,

Systematic Inclusion of High-Order Multi-Spin Correlations for the Spin-121\over2 XXZXXZ Models

We apply the microscopic coupled-cluster method (CCM) to the spin-$1\over2$ $XXZ$ models on both the one-dimensional chain and the two-dimensional square lattice. Based on a systematic approximation scheme of the CCM developed by us previously, we carry out high-order {\it ab initio} calculations using computer-algebraic techniques. The ground-state properties of the models are obtained with high accuracy as functions of the anisotropy parameter. Furthermore, our CCM analysis enables us to study their quantum critical behavior in a systematic and unbiased manner.Comment: (to appear in PRL). 4 pages, ReVTeX, two figures available upon request. UMIST Preprint MA-000-000

Integrity static analysis of COTS/SOUP

This paper describes the integrity static analysis approach developed to support the justification of commercial off-the-shelf software (COTS) used in a safety-related system. The static analysis was part of an overall software qualification programme, which also included the work reported in our paper presented at Safecomp 2002. Integrity static analysis focuses on unsafe language constructs and “covert” flows, where one thread can affect the data or control flow of another thread. The analysis addressed two main aspects: the internal integrity of the code (especially for the more critical functions), and the intra-component integrity, checking for covert channels. The analysis process was supported by an aggregation of tools, combined and engineered to support the checks done and to scale as necessary. Integrity static analysis is feasible for industrial scale software, did not require unreasonable resources and we provide data that illustrates its contribution to the software qualification programme

Linearized solutions of the Einstein equations within a Bondi-Sachs framework, and implications for boundary conditions in numerical simulations

We linearize the Einstein equations when the metric is Bondi-Sachs, when the background is Schwarzschild or Minkowski, and when there is a matter source in the form of a thin shell whose density varies with time and angular position. By performing an eigenfunction decomposition, we reduce the problem to a system of linear ordinary differential equations which we are able to solve. The solutions are relevant to the characteristic formulation of numerical relativity: (a) as exact solutions against which computations of gravitational radiation can be compared; and (b) in formulating boundary conditions on the $r=2M$ Schwarzschild horizon.Comment: Revised following referee comment

Lightweight ducts fabricated from reinforced plastics and elastomers

Method has been developed for fabrication of lightweight ducts that are three times stronger than aluminum ducts. Method can be used to produce either flexible or rigid ducts

Numerical relativity with characteristic evolution, using six angular patches

The characteristic approach to numerical relativity is a useful tool in evolving gravitational systems. In the past this has been implemented using two patches of stereographic angular coordinates. In other applications, a six-patch angular coordinate system has proved effective. Here we investigate the use of a six-patch system in characteristic numerical relativity, by comparing an existing two-patch implementation (using second-order finite differencing throughout) with a new six-patch implementation (using either second- or fourth-order finite differencing for the angular derivatives). We compare these different codes by monitoring the Einstein constraint equations, numerically evaluated independently from the evolution. We find that, compared to the (second-order) two-patch code at equivalent resolutions, the errors of the second-order six-patch code are smaller by a factor of about 2, and the errors of the fourth-order six-patch code are smaller by a factor of nearly 50.Comment: 12 pages, 5 figures, submitted to CQG (special NFNR issue

High-powered Gravitational News

We describe the computation of the Bondi news for gravitational radiation. We have implemented a computer code for this problem. We discuss the theory behind it as well as the results of validation tests. Our approach uses the compactified null cone formalism, with the computational domain extending to future null infinity and with a worldtube as inner boundary. We calculate the appropriate full Einstein equations in computational eth form in (a) the interior of the computational domain and (b) on the inner boundary. At future null infinity, we transform the computed data into standard Bondi coordinates and so are able to express the news in terms of its standard $N_{+}$ and $N_{\times}$ polarization components. The resulting code is stable and second-order convergent. It runs successfully even in the highly nonlinear case, and has been tested with the news as high as 400, which represents a gravitational radiation power of about $10^{13}M_{\odot}/sec$.Comment: 24 pages, 4 figures. To appear in Phys. Rev.

Research study for determination of liquid surface profile in a cryogenic tank during gas injection Quarterly progress report no. 9, Jun. 18 - Sep. 17, 1966

Determining liquid surface profiles in cryogenic tank during gas injectio

Direct calculation of the spin stiffness on square, triangular and cubic lattices using the coupled cluster method

We present a method for the direct calculation of the spin stiffness by means of the coupled cluster method. For the spin-half Heisenberg antiferromagnet on the square, the triangular and the cubic lattices we calculate the stiffness in high orders of approximation. For the square and the cubic lattices our results are in very good agreement with the best results available in the literature. For the triangular lattice our result is more precise than any other result obtained so far by other approximate method.Comment: 5 pages, 2 figure