A generalized chemistry version of SPARK

An extension of the reacting H2-air computer code SPARK is presented, which enables the code to be used on any reacting flow problem. Routines are developed calculating in a general fashion, the reaction rates, and chemical Jacobians of any reacting system. In addition, an equilibrium routine is added so that the code will have frozen, finite rate, and equilibrium capabilities. The reaction rate for the species is determined from the law of mass action using Arrhenius expressions for the rate constants. The Jacobian routines are determined by numerically or analytically differentiating the law of mass action for each species. The equilibrium routine is based on a Gibbs free energy minimization routine. The routines are written in FORTRAN 77, with special consideration given to vectorization. Run times for the generalized routines are generally 20 percent slower than reaction specific routines. The numerical efficiency of the generalized analytical Jacobian, however, is nearly 300 percent better than the reaction specific numerical Jacobian used in SPARK

The Effects of Mandatory Seatbelt Laws on Seatbelt Use, Motor Vehicle Fatalities, and Crash-Related Injuries among Youths

We provide the first comprehensive assessment of the effects of mandatory seatbelt laws on self-reported seatbelt use, highway fatalities, and crash-related injuries among high school age youths using data from the Centers for Disease Control's (CDC) national, state, and local Youth Risk Behavior Surveys (YRBS) and the Fatality Analysis Reporting System (FARS) from 1991 to 2005, a period spanning over 20 changes in state seatbelt laws. Our quasi-experimental approaches isolate the independent effects of seatbelt laws net of demographic characteristics, area and year fixed effects, and smooth area-specific trends. Across all data sources, we find consistent evidence that state mandatory seatbelt laws -- particularly those permitting primary enforcement -- significantly increased seatbelt use among high school age youths by 45-80 percent, primarily at the extensive margin. Unlike previous research for adults, however, we find evidence against the selective recruitment hypothesis: seatbelt laws had consistently larger effects on those most likely to be involved in traffic accidents (drinkers, alcohol-involved drivers). We also find that mandatory seatbelt laws significantly reduced traffic fatalities and serious injuries resulting from fatal crashes by 8 and 9 percent, respectively. Our results suggest that if all states had primary enforcement seatbelt laws then regular youth seatbelt use would be nearly universal and youth fatalities would fall by about 120 per year.

Decoding messages from the yield curve

Bank supervision ; Federal Reserve District, 5th

High-order cyclo-difference techniques: An alternative to finite differences

The summation-by-parts energy norm is used to establish a new class of high-order finite-difference techniques referred to here as 'cyclo-difference' techniques. These techniques are constructed cyclically from stable subelements, and require no numerical boundary conditions; when coupled with the simultaneous approximation term (SAT) boundary treatment, they are time asymptotically stable for an arbitrary hyperbolic system. These techniques are similar to spectral element techniques and are ideally suited for parallel implementation, but do not require special collocation points or orthogonal basis functions. The principal focus is on methods of sixth-order formal accuracy or less; however, these methods could be extended in principle to any arbitrary order of accuracy

The stability of numerical boundary treatments for compact high-order finite-difference schemes

The stability characteristics of various compact fourth and sixth order spatial operators are assessed using the theory of Gustafsson, Kreiss and Sundstrom (G-K-S) for the semi-discrete Initial Boundary Value Problem (IBVP). These results are then generalized to the fully discrete case using a recently developed theory of Kreiss. In all cases, favorable comparisons are obtained between the G-K-S theory, eigenvalue determination, and numerical simulation. The conventional definition of stability is then sharpened to include only those spatial discretizations that are asymptotically stable. It is shown that many of the higher order schemes which are G-K-S stable are not asymptotically stable. A series of compact fourth and sixth order schemes, which are both asymptotically and G-K-S stable for the scalar case, are then developed

Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes

We present a systematic method for constructing boundary conditions (numerical and physical) of the required accuracy, for compact (Pade-like) high-order finite-difference schemes for hyperbolic systems. First, a roper summation-by-parts formula is found for the approximate derivative. A 'simultaneous approximation term' (SAT) is then introduced to treat the boundary conditions. This procedure leads to time-stable schemes even in the system case. An explicit construction of the fourth-order compact case is given. Numerical studies are presented to verify the efficacy of the approach

Threshold Determination for ARTMAP-FD Familiarity Discrimination

The ARTMAP-FD neural network performs both identification (placing test patterns in classes encountered during training) and familiarity discrimination (judging whether a test pattern belongs to any of the classes encountered during training). ARTMAP-FD quantifies the familiarity of a test pattern by computing a measure of the degree to which the pattern's components lie within the ranges of values of training patterns grouped in the same cluster. This familiarity measure is compared to a threshold which can be varied to generate a receiver operating characteristic (ROC) curve. Methods for selecting optimal values for the threshold are evaluated. The performance of validation-set methods is compared with that of methods which track the development of the network's discrimination capability during training. The techniques are applied to databases of simulated radar range profiles.Advanced Research Projects Agency; Office of Naval Research (N00011-95-1-0657, N00011-95-0109, NOOOB-96-0659); National Science Foundation (IRI-94-01659