SCI Publications

Multigrid methods for solving the steady-state incompressible Navier-Stokes equations require an appropriate smoother and coarse grid solution strategy to be effective. Classical pressure-correction methods, such as SIMPLE and SIMPLER, are widely used as solvers in engineering analysis codes, but can also be used as effective multigrid smoothers. An inexact Newton method preconditioned by a linear multigrid method with a pressurecorrection smoother can serve as a coarse grid solver. A hybrid nonlinear multigrid scheme based on combinations of these components is described. A standard benchmark problem is used to demonstrate the effectiveness of SIMPLER smoothing and the impact an inexact Newton coarse grid solver has on the resulting nonlinear multigrid scheme.

Neutron scattering and computer simulations are powerful tools for studying structural and dynamical properties of condensed matter systems in general and of polymer melts in particular. When neutron scattering studies and quantitative atomistic molecular dynamics simulations of the same material are combined, synergy between the methods can result in exciting new insights into polymer melts not obtainable from either method separately. We present here an overview of our recent efforts to combine neutron scattering and atomistic simulations in the study of melt dynamics of polyethylene and polybutadiene. Looking at polymer segmental motion on a picosecond time scale, we show how atomistic simulations can be used to identify molecular motions giving rise to relaxation processes observed in experimental dynamic susceptibility spectra. Examining larger length and longer time scale polymer dynamics involving chain self-diffusion and overall conformational relaxation, we show how simulation results can motivate experiment and how combined results of scattering and simulation can be used to critically test theories that attempt to describe melt dynamics of short polymer chains.

A concept and algorithm are presented for non-iterative robust estimation of piecewise smooth curves of maximal edge strength in small image windows-typically 8/spl times/8 to 32/spl times/32. This boundary-estimation algorithm has the nice properties that it uses all the data in the window and thus can find locally weak boundaries embedded in noise or texture and boundaries when there are more than two regions to be segmented in a window; it does not require step edges-but handles ramp edges well. The curve-estimates found are among the level sets of a dth degree polynomial fit to "suitable" weightings of the image gradient vector at each pixel in the image window. Since the polynomial fitting is linear least squares, the computation to this point is very fast. Level sets then chosen to be appropriate boundary curves are those having the highest differences in average gray level in regions to either side. This computation is also fast. The boundary curves and segmented regions found are suitable for all purposes but especially for indexing using algebraic curve invariants in this form.

We present a new method called Reaction Class Transition State Theory (RC-TST) for estimating thermal rate constants of a large number of reactions in a class. This method is based on the transition state theory framework within the reaction class approach. Thermal rate constants of a given reaction in a class relative to those of its principal reaction can be efficiently predicted from only its differential barrier height and reaction energy. Such requirements are much less than what is needed by the conventional TST method. Furthermore, we have shown that the differential energetic information can be calculated at a relatively low level of theory. No frequency calculation beyond those of the principal reaction is required for this theory. The new theory was applied to a number of hydrogen abstraction reactions. Excellent agreement with experimental data shows that the RC-TST method can be very useful in design of fundamental kinetic models of complex reactions.

We present a new practical computational methodology for predicting thermal rate constants of reactions involving large molecules or a large number of elementary reactions in the same class. This methodology combines the integrated molecular orbital+molecular orbital (IMOMO) approach with our recently proposed reaction class models for tunneling. With the new methodology, we show that it is possible to significantly reduce the computational cost by several orders of magnitude while compromising the accuracy in the predicted rate constants by less than 40% over a wide range of temperatures. Another important result is that the computational cost increases only slightly as the system size increases.

A statistical procedure is proposed for estimating realistic confidence intervals for the activation energy determined by using an advanced isoconversional method. Nine sets of five thermogravimetric measurements have been produced for the process of gassification of ammonium nitrate at five different heating rates. Independent estimates of the confidence intervals for the activation energy have been obtained from these data sets. Agreement with these independent estimates demonstrates that the proposed statistical procedure is capable of adequately estimating the actual uncertainty in the activation energy determined from a small number of measurements. The resulting averaged relative errors in the activation energy were found to be 26, 21, and 17% for three, four, and five heating rate estimates, respectively.