A CSP and tabulation-based adaptive chemistry model

We demonstrate the feasibility of a new strategy for the construction of an adaptive chemistry model that is based on an explicit integrator stabilized by an approximation of the Computational Singular Perturbation (CSP)-slow-manifold projector. We examine the effectiveness and accuracy of this technique first using a model problem with variable stiffness. We assess the effect of using an approximation of the CSP-slow-manifold by either reusing the CSP vectors calculated in previous steps or from a pre-built tabulation. We find that while accuracy is preserved, the associated CPU cost was reduced substantially by this method. We used two ignition simulations – hydrogen–air and heptane–air mixtures – to demonstrate the feasibility of using the new method to handle realistic kinetic mechanisms. We test the effect of utilizing an approximation of the CSP-slow-manifold and find that its use preserves the order of the explicit integrator, produces no degradation in accuracy, and results in a scheme that is competitive with traditional implicit integration. Further analysis on the performance data demonstrates that the tabulation of the CSP-slow-manifold provides an increasing level of efficiency as the size of the mechanism increases. From the software engineering perspective, all the machinery developed is Common Component Architecture compliant, giving the software a distinct advantage in the ease of maintainability and flexibility in its utilization. Extension of this algorithm is underway to implement an automated tabulation of the CSP-slow-manifold for a detailed chemical kinetic system either off-line, or on-line with a reactive flow simulation code.     

Chemically shortened multi-walled carbon nanotubes used as anode materials for lithium-ion batteries

An easy chemically cutting process, modified Hummers' method, was proposed to treat multi-walled carbon nanotubes, successfully cutting pristine long, entangled carbon nanotubes into hydrosoluble pieces, mostly less than 200 nm. This short, chemically oxidized carbon nanotube was then applied as an anode material for lithium-ion batteries. The as-prepared material possessed higher reversible capacity and coulombic efficiency. The intrinsic factors were explored by X-ray photoelectron spectroscopy and cyclic voltammetry.     

厄尔尼诺和南方涛动海气耦合模型中参数估计的变分方法

提出一种基于变分原理估计厄尔尼诺和南方涛动海气耦合模型中未知参数的方法. 首先将所研究的非线性海气耦合动力方程引入到目标泛函中; 接着利用变分方法导出伴随方程和待辨识参数泛函梯度的公式; 然后设计了估计未知参数的算法.数值试验结果表明变分方法是一 种能有效估计海气耦合非线性系统未知参数的方法.     

Some MnIII-porphyrins with de-polymerization activity toward humic acid

Four Mn^III-porphyrin complexes, chloro(tetraphenylporphinato)Mn^III(1,6-diaminohexane), bromo(tetraphenylporphinato)Mn^III(1,6-diaminohexane), azido(tetraphenylporphinato)Mn^III(1,6-diaminohexane), and thiocyanato(tetraphenylporphinato)Mn^III(1,6-diaminohexane), have been synthesized. These complexes have been characterized using UV-Vis, IR, ESI-mass spectra, elemental analyses, magnetic susceptibility measurements, and conductivity measurement. The molar conductance values of these complexes in ethanol indicate non-electrolytes. The utility of these complexes in de-polymerization of coal using humic acid as the coal model has been tested by the optical density method.     

An improved progressive preconditioning method for steady non‐cavitating and sheet‐cavitating flows

An improved progressive preconditioning method for analyzing steady inviscid and laminar flows around fully wetted and sheet‐cavitating hydrofoils is presented. The preconditioning matrix is adapted automatically from the pressure and/or velocity flow‐field by a power‐law relation. The cavitating calculations are based on a single fluid approach. In this approach, the liquid/vapour mixture is treated as a homogeneous fluid whose density is controlled by a barotropic state law. This physical model is integrated with a numerical resolution derived from the cell‐centered Jameson's finite volume algorithm. The stabilization is achieved via the second‐and fourth‐order artificial dissipation scheme. Explicit four‐step Runge–Kutta time integration is applied to achieve the steady‐state condition. Results presented in the paper focus on the pressure distribution on hydrofoils wall, velocity profiles, lift and drag forces, length of sheet cavitation, and effect of the power‐law preconditioning method on convergence speed. The results show satisfactory agreement with numerical and experimental works of others. The scheme has a progressive effect on the convergence speed. The results indicate that using the power‐law preconditioner improves the convergence rate, significantly. Copyright © 2010 John Wiley & Sons, Ltd.     

On the usage of NURBS as interface representation in free‐surface flows

When simulating free‐surface flows using the finite element method, there are many cases where the governing equations require information which must be derived from the available discretized geometry. Examples are curvature or normal vectors. The accurate computation of this information directly from the finite element mesh often requires a high degree of refinement—which is not necessarily required to obtain an accurate flow solution. As a remedy and an option to be able to use coarser meshes, the representation of the free surface using non‐uniform rational B‐splines (NURBS) curves or surfaces is investigated in this work. The advantages of a NURBS parameterization in comparison with the standard approach are discussed. In addition, it is explored how the pressure jump resulting from surface tension effects can be handled using doubled interface nodes. Numerical examples include the computation of surface tension in a two‐phase flow as well as the computation of normal vectors as a basis for mesh deformation methods. For these examples, the improvement of the numerical solution compared with the standard approaches on identical meshes is shown. Copyright © 2011 John Wiley & Sons, Ltd.     

A least‐square weighted residual method for level set formulation

In this paper, a least‐square weighted residual method (LSWRM) for level set (LS) formulation is introduced to achieve interface capturing in two‐dimensional (2D) and three‐dimensional (3D) problems. An LSWRM was adopted for two semi‐discretized advection and reinitialization equations of the LS formulation. The present LSWRM provided good mathematical properties such as natural numerical diffusion and the symmetry of the resulting algebraic systems for the advection and reinitialization equations. The proposed method was validated by solving some 2D and 3D benchmark problems such as those involving a rotating slotted disk, the rotation of a slotted sphere, and a time‐reversed single‐vortex flow and a deformation problem of a spherical fluid. The numerical results were compared with those obtained from essentially non‐oscillatory type formulations and particle LS methods. Further, the proposed LSWRM for the LS formulation was coupled with a splitting finite element method code to solve the incompressible Navier–Stokes equations, and then, the collapse of a 3D broken dam flow was well simulated; in the simulation, the entrapping of air and the splashing of the surge front of water were reproduced. The mass conservation of the present method was found to be satisfactory during the entire simulation. Copyright © 2011 John Wiley & Sons, Ltd.     

Three‐dimensional computation for flow‐induced vibrations in in‐line and cross‐flow directions of a circular cylinder

We present numerical results for in‐line and cross‐flow vibrations of a circular cylinder, which is immersed in a uniform flow and is elastically supported by damper‐spring systems to compute vibrations of a rigid cylinder. In the case of a circular cylinder with a low Scruton number, it is well‐known that two types of self‐excited vibrations appear in the in‐line direction in the range of low reduced velocities. On the other hand, a cross‐flow vibration of the circular cylinder can be excited in the range of high reduced velocities. Therefore, we compute the flow‐induced vibrations of the circular cylinder in the wide range of the reduced velocities at low and high Scruton numbers and discuss about excitation mechanisms in the in‐line and cross‐flow directions. Copyright © 2011 John Wiley & Sons, Ltd.     

A vectorial algorithm with finite element method for prediction of powder segregation in metal injection molding

The design of the mold and the choice of the injection parameters for metal injection molding (MIM) is required to maintain homogeneity of the filled mixture. However, powder segregation is unavoidable in MIM because of the significant difference in densities of the metallic powder and the polymer binder. To achieve an effective prediction of segregation effect, a biphasic model based on mixture theory is employed. The viscous behaviors of each phase and the interaction coefficient between the flows of the two phases should be determined. The solution of two coupled Navier–Stokes equations results in a tremendous computation effort. The previous development of an explicit algorithm makes the biphasic simulation much faster than that of the classic methods. However, it is strongly desired to reduce or even eliminate the numerous global solutions for pressure fields at each time step. Hence, a new vectorial algorithm is proposed and developed to perform the simulation only by vectorial operations. It provides the anticipated efficiency in the simulation of biphasic modeling, and the advantage to use the classic elements of equal‐order interpolations. Some results produced by the two algorithms are compared with the experimental values to validate the new vectorial algorithm. Copyright © 2012 John Wiley & Sons, Ltd.     

The spherical-model limit in a random field

The spherical-model limitn of then-vector model in a random field, with either a statistically independent distribution or with long-range correlated random fields, is studied to demonstrate the correctness of the replica method in which then and replica limits limits are interchanged, provided the replica and thermodynamic limits are taken in the right order, in the case of long-range correlated random fields. A scaling form for the two-point correlation function relevant to the first-order phase transition below the lower critical dimensionality of the random system is also obtained.