Reaction Mechanism Reduction for Ozone-Enhanced CH4/Air Combustion by a Combination of Directed Relation Graph with Error Propagation, Sensitivity Analysis and Quasi-Steady State Assumption

In this study, an 18-steps, 22-species reduced global mechanism for ozone-enhanced CH4/air combustion processes was derived by coupling GRI-Mech 3.0 and a sub-mechanism for ozone decomposition. Three methods, namely, direct relation graphics with error propagation, (DRGRP), sensitivity analysis (SA), and quasi-steady-state assumption (QSSA), were used to downsize the detailed mechanism to the global mechanism. The verification of the accuracy of the skeletal mechanism in predicting the laminar flame speeds and distribution of the critical components showed that that the major species and the laminar flame speeds are well predicted by the skeletal mechanism. However, the pollutant NO was predicated inaccurately due to the precursors for generating NO were removed as redundant components. The laminar flame speeds calculated by the global mechanism fit the experimental data well. The comparisons of simulated results between the detailed mechanism and global mechanism were investigated and showed that the global mechanism could accurately predict the major and intermediate species and significantly reduced the time cost by 72%Peer reviewe

Simulations of Incompressible MHD Turbulence

We simulate incompressible MHD turbulence in the presence of a strong background magnetic field. Our major conclusions are: 1) MHD turbulence is most conveniently described in terms of counter propagating shear Alfven and slow waves. Shear Alfven waves control the cascade dynamics. Slow waves play a passive role and adopt the spectrum set by the shear Alfven waves, as does a passive scalar. 2) MHD turbulence is anisotropic with energy cascading more rapidly along k_perp than along k_parallel, where k_perp and k_parallel refer to wavevector components perpendicular and parallel to the local magnetic field. Anisotropy increases with increasing k_perp. 3) MHD turbulence is generically strong in the sense that the waves which comprise it suffer order unity distortions on timescales comparable to their periods. Nevertheless, turbulent fluctuations are small deep inside the inertial range compared to the background field. 4) Decaying MHD turbulence is unstable to an increase of the imbalance between the flux of waves propagating in opposite directions along the magnetic field. 5) Items 1-4 lend support to the model of strong MHD turbulence by Goldreich & Sridhar (GS). Results from our simulations are also consistent with the GS prediction gamma=2/3. The sole notable discrepancy is that 1D power law spectra, E(k_perp) ~ k_perp^{-alpha}, determined from our simulations exhibit alpha ~ 3/2, whereas the GS model predicts alpha = 5/3.Comment: 56 pages, 30 figures, submitted to ApJ 59 pages, 31 figures, accepted to Ap

A semi-explicit multi-step method for solving incompressible navier-stokes equations

The fractional step method is a technique that results in a computationally-efficient implementation of Navier–Stokes solvers. In the finite element-based models, it is often applied in conjunction with implicit time integration schemes. On the other hand, in the framework of finite difference and finite volume methods, the fractional step method had been successfully applied to obtain predictor-corrector semi-explicit methods. In the present work, we derive a scheme based on using the fractional step technique in conjunction with explicit multi-step time integration within the framework of Galerkin-type stabilized finite element methods. We show that under certain assumptions, a Runge–Kutta scheme equipped with the fractional step leads to an efficient semi-explicit method, where the pressure Poisson equation is solved only once per time step. Thus, the computational cost of the implicit step of the scheme is minimized. The numerical example solved validates the resulting scheme and provides the insights regarding its accuracy and computational efficiency.Peer ReviewedPostprint (published version

Universally Rigid Framework Attachments

A framework is a graph and a map from its vertices to R^d. A framework is called universally rigid if there is no other framework with the same graph and edge lengths in R^d' for any d'. A framework attachment is a framework constructed by joining two frameworks on a subset of vertices. We consider an attachment of two universally rigid frameworks that are in general position in R^d. We show that the number of vertices in the overlap between the two frameworks must be sufficiently large in order for the attachment to remain universally rigid. Furthermore, it is shown that universal rigidity of such frameworks is preserved even after removing certain edges. Given positive semidefinite stress matrices for each of the two initial frameworks, we analytically derive the PSD stress matrices for the combined and edge-reduced frameworks. One of the benefits of the results is that they provide a general method for generating new universally rigid frameworks.Comment: 16 pages, 4 figure

Instabilities of the Hubbard chain in a magnetic field

We find and characterize the instabilities of the repulsive Hubbard chain in a magnetic field by studing all response functions at low frequency \omega and arbitrary momentum. The instabilities occur at momenta which are simple combinations of the (U=0) \sigma =\uparrow ,\downarrow Fermi points, \pm k_{F\sigma}. For finite values of the on-site repulsion U the instabilities occur for single \sigma electron adding or removing at momenta \pm k_{F\sigma}, for transverse spin-density wave (SDW) at momenta \pm 2k_F (where 2k_F=k_{F\uparrow}+k_{F\downarrow}), and for charge-density wave (CDW) and SDW at momenta \pm 2k_{F\uparrow} and \pm 2k_{F\downarrow}. While at zero magnetic field removing or adding single electrons is dominant, the presence of that field brings about a dominance for the transverse \pm 2k_F SDW over all the remaining instabilities for a large domain of $U$ and density n values. We go beyond conformal-field theory and study divergences which occur at finite frequency in the one-electron Green function at half filling and in the transverse-spin response function in the fully-polarized ferromagnetic phase.Comment: LaTeX file, 15 pages plus 9 figures. Accepted for publication in Phys. Rev. B. The figures can be obtained upon request from Pedro Sacramento at [email protected]