Stability analysis of spectral methods for hyperbolic initial-boundary value systems

A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations

Convergence of spectral methods for hyperbolic initial-boundary value systems

A convergence proof for spectral approximations is presented for hyperbolic systems with initial and boundary conditions. The Chebyshev collocation is treated in detail, but the final result is readily applicable to other spectral methods, such as Legendre collocation or tau-methods

Quadrature imposition of compatibility conditions in Chebyshev methods

Often, in solving an elliptic equation with Neumann boundary conditions, a compatibility condition has to be imposed for well-posedness. This condition involves integrals of the forcing function. When pseudospectral Chebyshev methods are used to discretize the partial differential equation, these integrals have to be approximated by an appropriate quadrature formula. The Gauss-Chebyshev (or any variant of it, like the Gauss-Lobatto) formula can not be used here since the integrals under consideration do not include the weight function. A natural candidate to be used in approximating the integrals is the Clenshaw-Curtis formula, however it is shown that this is the wrong choice and it may lead to divergence if time dependent methods are used to march the solution to steady state. The correct quadrature formula is developed for these problems. This formula takes into account the degree of the polynomials involved. It is shown that this formula leads to a well conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied

Staggered Fermion Thermodynamics using Anisotropic Lattices

Numerical simulations of full QCD on anisotropic lattices provide a convenient way to study QCD thermodynamics with fixed physics scales and reduced lattice spacing errors. We report results from calculations with 2-flavors of dynamical fermions where all bare parameters and hence the physics scales are kept constant while the temperature is changed in small steps by varying only the number of the time slices. The results from a series of zero-temperature scale setting simulations are used to determine the Karsch coefficients and the equation of state at finite temperatures.Comment: Lattice2002(nonzerot), 3 pages, 2 figure

Urban Resurgence and the Consumer City

Cities make it easier for humans to interact, and one of the main advantages of dense, urban areas is that they facilitate social interactions. This paper provides evidence suggesting that the resurgence of big cities in the 1990s is due, in part, to the increased demand for these interactions and due to the reduction in big city crime, which had made it difficult for urban residents to enjoy these social amenities. However, while density is correlated with consumer amenities, we show that it is not correlated with social capital and that there is no evidence that sprawl has hurt civic engagement.

A quasi-linear control theory analysis of timesharing skills

The compliance of the human ankle joint is measured by applying 0 to 50 Hz band-limited gaussian random torques to the foot of a seated human subject. These torques rotate the foot in a plantar-dorsal direction about a horizontal axis at a medial moleolus of the ankle. The applied torques and the resulting angular rotation of the foot are measured, digitized and recorded for off-line processing. Using such a best-fit, second-order model, the effective moment of inertia of the ankle joint, the angular viscosity and the stiffness are calculated. The ankle joint stiffness is shown to be a linear function of the level of tonic muscle contraction, increasing at a rate of 20 to 40 Nm/rad/Kg.m. of active torque. In terms of the muscle physiology, the more muscle fibers that are active, the greater the muscle stiffness. Joint viscosity also increases with activation. Joint stiffness is also a linear function of the joint angle, increasing at a rate of about 0.7 to 1.1 Nm/rad/deg from plantar flexion to dorsiflexion rotation

QCD Thermodynamics at Nt=8N_t=8 and 12

We present results from studies of high temperature QCD with two flavors of Kogut-Susskind quarks on $16^3\times 8$ lattices at a quark mass of $am_q=0.00625$ and on $24^3\times 12$ lattices at quark masses $am_q=0.008$ and 0.016. The value of the crossover temperature is consistent with that obtained on coarser lattices and/or at larger quark masses. Results are presented for the chiral order parameter and for the baryon number susceptibility.Comment: 3-pages, uuencoded compressed postscript file, contribution to Lattice'94 conferenc

Further observations on the relationship of EMG and muscle force

Human skeletal muscle may be regarded as an electro-mechanical transducer. Its physiological input is a neural signal originating at the alpha motoneurons in the spinal cord and its output is force and muscle contraction, these both being dependent on the external load. Some experimental data taken during voluntary efforts around the ankle joint and by direct electrical stimulation of the nerve are described. Some of these experiments are simulated by an analog model, the input of which is recorded physiological soleus muscle EMG. The output is simulated foot torque. Limitations of a linear model and effect of some nonlinearities are discussed

Anisotropic Lattices and Dynamical Fermions

We report results from full QCD calculations with two flavors of dynamical staggered fermions on anisotropic lattices. The physical anisotropy as determined from spatial and temporal masses, their corresponding dispersion relations, and spatial and temporal Wilson loops is studied as a function of the bare gauge anisotropy and the bare velocity of light appearing in the Dirac operator. The anisotropy dependence of staggered fermion flavor symmetry breaking is also examined. These results will then be applied to the study of 2-flavor QCD thermodynamics.Comment: Lattice2001(spectrum

Strong Stability Preserving Two-Step Runge-Kutta Methods

We investigate the strong stability preserving (SSP) property of two-step Runge– Kutta (TSRK) methods. We prove that all SSP TSRK methods belong to a particularly simple\ud subclass of TSRK methods, in which stages from the previous step are not used. We derive simple order conditions for this subclass. Whereas explicit SSP Runge–Kutta methods have order at most four, we prove that explicit SSP TSRK methods have order at most eight. We present TSRK methods of up to eighth order that were found by numerical search. These methods have larger SSP coefficients than any known methods of the same order of accuracy, and may be implemented in a form with relatively modest storage requirements. The usefulness of the TSRK methods is demonstrated through numerical examples, including integration of very high order WENO discretizations