On upgrading the numerics in combustion chemistry codes

A method of updating and reusing legacy FORTRAN codes for combustion simulations is presented using the DAEPACK software package. The procedure is demonstrated on two codes that come with the CHEMKIN-II package, CONP and SENKIN, for the constant-pressure batch reactor simulation. Using DAEPACK generated code, analytical derivative calculations, sparsity pattern information, and hidden discontinuity information can be obtained for the models of interest. This information can be easily integrated with different solvers giving the modeler great flexibility in selecting the best solution procedure. Using the generated code, the CONP code was connected to three different solvers, and the SENKIN code was connected to two different solvers. The effect of model formulation, analytical derivatives, sparsity, and sensitivity equation solution method were analyzed for three large kinetic mechanisms for methane, acetylene, and n-heptane. For the n-heptane model, with 544 species and 2446 reactions, a factor of 10-speed improvement over the original solution procedure was found using analytical derivatives and sparse linear algebra. For sensitivity calculations, for a small number of parameters, a factor of 55 improvement over the original solution procedure was found for the n-heptane problem. Upon closer examination of results, no one method is found to always be superior to other methods, and selection of the appropriate solution procedure requires an examination of the specific kinetic mechanism, which is easily conducted using DAEPACK generated code.     

Network topology and the theory of rubber elasticity

The earliest investigations on rubber elasticity, commencing in the 19th century, were necessarily limited to phenomenological interpretations. The realisation that polymers consist of very long molecular chains. commencing c. 1930, gave impetus to the molecular theory of rubber elasticity (1932-). according to which the high deformability of an elastomer, and the elastic force generated by deformation, stem from the configurations accessible to long molecular chains. Theories of rubber elasticity put forward from 1934-1946 relied on the assumption that the junctions of the rubber network undergo displacements that are affine in macroscopic strain. The theory of James and Guth (1947) dispensed with this premise, and demonstrated instead that the mean positions of the junctions of a ‘phantom’ network consisting of Gaussian chains devoid of material properties are affine in the strain. The vital significance of the distinction between the actual distribution of chain vectors in a network and their distribution if the junctions would be fixed at their mean positions went unnoticed for nearly 30 years. Experimental investigations, commencing with the incisive work of Gee in 1946. revealed large departures from the relationship of stress to strain predicted by the theories cited. This discrepancy prompted extensive studies, theoretical and experimental, during succeeding years. Inquiry into the fundamentals of polymer networks, formed for example by interlinking very long polymer molecules, exposed the need to take account of network imperfections, typically consisting of chains attached at only one end to a network junction. Various means were advocated to make corrections for these imperfections. The cycle rank ζ of the network has been shown (1976) to be the fundamental measure of its connectivity, regardless of the junction functionality and pattern of imperfections. Often overlooked is the copious interpenetration of the chains comprising typical elastomeric networks. Theories that attempt to represent such networks on a lattice are incompatible with this universal feature. Moreover, the dense interpenetration of chains may limit the ability of junctions in real networks to accommodate the fluctuations envisaged in the theory of phantom networks. It was suggested in 1975 that departures from the form predicted for the elastic equation of state are due to constraints on the fluctuations of junctions whose effect diminishes with deformation and with dilation. Formulation of a self-consistent theory based on this suggestion required recognition of the non-affine connection between the chain vector distribution function and the macroscopic strain in a real network, which may partake of characteristics of a phantom network in some degree. Implementation of the idea was achieved through postulation of domains of constraint affecting the equilibrium distribution of fluctuations of network junctions from their mean positions. This led in due course to a theory that accounts for the relationship of stress to strain virtually throughout the ranges of strain accessible to measurement. The theory establishes connections between structure and elastic properties. This is achieved with utmost frugality in arbitrary parameters.     

Cooperation between memory systems: Acetylcholine release in the amygdala correlates positively with performance on a hippocampus-dependent task.

The present experiment examined the relationship between release of acetylcholine (ACh) in the amygdala and performance on a hippocampus-dependent spatial working memory task. Using in vivo microdialysis, the authors measured ACh release in rats during testing on a spontaneous alternation task. Amygdala ACh release was positively correlated with performance on the hippocampus-dependent task. These findings suggest that activation of the amygdala promotes processing in other neural systems important for learning and memory. (PsycINFO Database Record (c) 2010 APA, all rights reserved)     

An adaptive descent method for non-linear viscoplasticity

A forming model based on a viscoplastic flow formulation is derived which includes the effects of small elastic strains. A significant feature of the formulation is its reliance on the dominant inelastic material characteristics in the formation of the stiffness matrix for large strain problems. The resultant non-linear system of equations is solved by an adaptive descent method which combines the rapid convergence of Newton's method near the solution with the robustness of a method of successive approximations. The use of the adaptive descent method effectively extends the viscoplastic flow formulations into the nearly rate-insensitive range of behaviours exhibited, for example, by metals at low temperature, where slow convergence of the non-linear solution algorithm has traditionally hampered their use.     

The hot blow forming of AZ31 Mg sheet: Formability assessment and application development

The hot blow forming of magnesium sheet offers significant opportunity for forming complex, lightweight parts for automotive applications. This paper characterizes the elevated-temperature formability of AZ31 magnesium sheet materials and the effect of processing conditions on the performance of these materials. In addition, magnesium sheet application development at General Motors Corporation is reviewed.     

Analysis of polar compound classes in SRC-II liquids: Comparison of non-aqueous titrametric,i.r. spectrometric and h.p.l.c. methods

Results for various types of polar compounds in SRC-11 coal-derived liquids and other fuels were obtained by one or more of the three analytical methods: non-aqueous titration, i.r. spectroscopy and h.p.l.c. Practical aspects as well as precision, accuracy and assumptions necessary for effective application of each of the methods are discussed. H.p.l.c. is applicable to the widest variety of compound types, is the most rapid, most sensitive and shows the best promise for increased development and improvement. Application and development of these techniques is a logical step toward improving process monitoring, catalyst development, toxicological screening and general fuel analysis.     

Modeling and optimization of fringe capacitance of nanoscale DGMOS devices

We analyze the impact of gate electrode thickness and gate underlap on the fringe capacitance of nanoscale double-gate MOS (DGMOS) transistors. We propose an analytical fringe capacitance model considering gate underlap and finite source/drain length. A comparison with the simulation results show that the model can accurately estimate the fringe capacitance of the device. We show that an optimum gate underlap can significantly reduce the fringe capacitance resulting in higher performance and lower power consumption. Also, the effects of process variation in gate underlap devices are discussed. Simulation results on a three-stage ring oscillator show that with optimum gate underlap 32% improvement in delay can be achieved.     

Identification of microstructural zones and thermal cycles in a weldment of modified 9Cr-1Mo steel

This paper presents the results of a study on the microstructural and microchemical variations in a multipass Gas Tungsten Arc weld (GTAW) of modified 9Cr-1Mo steel. The changes brought about in the steel due to the heating and cooling cycles during welding and the subsequent effects due to reheating effects during multipass welding are described. Detailed analytical transmission electron microscopy has been carried to study the type and composition of the primary and secondary phases in this steel. The systematic changes in microstructural parameters such as Prior Austenite Grain Size, martensite lath size, number density, size and microchemistry of carbides, have been understood based on the different transformations that the steel undergoes during the heating and cooling process. Based on the observed microstructure, an attempt has been made to identify distinct microstructural zones and possible thermal cycles experienced by different regions of the weldment.