Resolution estimates for simple one-dimensional flames

The accurate resolution of flame structure is critically important to the direct numerical simulation of reacting flows. While grid-resolution estimates are readily available for cold flows, premixed flames appear to have received relatively little attention. In this paper, a premixed flame characterised by single-step chemistry at moderate-to-high Zel'dovitch numbers (β) is analysed, and its structure is used to provide estimates for sufficient grid resolution. It is found that the critical region of the flame exhibits a weak inverse Zel'dovitch number dependence, and that heuristic methods of resolution estimation based on flame thickness grow relatively less meaningful as β → ∞. Resolution estimates for second- and fourth-order finite-difference schemes are provided.     

Issue Information

Three‐dimensional direct numerical simulation results of flow past a circular cylinder are influenced by numerical aspects, for example the spanwise domain length and the lateral boundary condition adopted for the simulation. It is found that inappropriate numerical set‐up, which restricts the development of intrinsic wake structure, leads to an over‐prediction of the onset point of the secondary wake instability (Re_cr). A best practice of the numerical set‐up is presented for the accurate prediction of Re_cr by direct numerical simulation while minimizing the computational cost. The cylinder span length should be chosen on the basis of the intrinsic wavelength of the wake structure to be simulated, whereas a long span length is not necessary. For the wake transitions above Re_cr, because the wake structures no longer follow particular wavelengths but become disordered and chaotic, a span length of more than 10 cylinder diameters (approximately three times the intrinsic wavelength) is recommended for the simulations to obtain wake structures and hydrodynamic forces that are not strongly restricted by the numerical set‐up. The performances of the periodic and symmetry lateral boundary conditions are compared and discussed. The symmetry boundary condition is recommended for predicting Re_cr, while the periodic boundary condition is recommended for simulating the wake structures above Re_cr. The general conclusions drawn through a circular cylinder are expected to be applicable to other bluff body configurations. Copyright © 2017 John Wiley & Sons, Ltd.     

Simulation of three‐dimensional flows over moving objects by an improved immersed boundary–lattice Boltzmann method

An improved immersed boundary–lattice Boltzmann method (IB–LBM) developed recently [28] was applied in this work to simulate three‐dimensional (3D) flows over moving objects. By enforcing the non‐slip boundary condition, the method could avoid any flow penetration to the wall. In the developed IB–LBM solver, the flow field is obtained on the non‐uniform mesh by the efficient LBM that is based on the second‐order one‐dimensional interpolation. As a consequence, its coefficients could be computed simply. By simulating flows over a stationary sphere and torus [28] accurately and efficiently, the proposed IB–LBM showed its ability to handle 3D flow problems with curved boundaries. In this paper, we further applied this method to simulate 3D flows around moving boundaries. As a first example, the flow over a rotating sphere was simulated. The obtained results agreed very well with the previous data in the literature. Then, simulation of flow over a rotating torus was conducted. The capability of the improved IB–LBM for solving 3D flows over moving objects with complex geometries was demonstrated via the simulations of fish swimming and dragonfly flight. The numerical results displayed quantitative and qualitative agreement with the date in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Impacts of air pressure on the evolution of nanosecond pulse discharge products

Based on the nonequilibrium plasma dynamics of air discharge, a dynamic model of zero-dimensional plasma is established by combining the component density equation, the Boltzmann equation, and the energy transfer equation. The evolution properties of nanosecond pulse discharge (NPD) plasma under different air pressures are calculated. The results show that the air pressure has significant impacts on the NPD products and the peak values of particle number density for particles such as O atoms, O3 molecules, N2(A3) molecules in excited states, and NO molecules. It increases at first and then decreases with the increase of air pressure. On the other hand, the peak values of particle number density for N2(B3) and N2(C3) molecules in excited states are only slightly affected by the air pressure.     

Hybrid Particle–Continuum Simulations of Polymer Brushes in Shear Flow

A hybrid particle–continuum method is used to study the shear flow confined between two opposing walls, one of which is coated with polymer chains. Molecular dynamics (MD) is used in the particle region near the brush and Navier–Stokes (NS) equations are applied in the remaining region where the continuum assumption holds. The information exchange from the continuum region to the particle region is implemented using the constrained particle dynamics. Both Couette shear flow and oscillatory flow are considered in the present work. The effect of the shear flow on the conformational characteristics of polymer brushes is analyzed. In the overlap region, the velocities obtained from MD simulations are smoothly connected with those from NS equations. Our investigations demonstrate that the hybrid particle–continuum model is valid in exploring the shear behavior of polymer brushes.     

Long-time simulation of a highly oscillatory Vlasov equation with an exponential integrator

We change a previous time-stepping algorithm for solving a multi-scale Vlasov–Poisson system within a Particle-In-Cell method, in order to perform accurate long-time simulations. As an exponential integrator, the new scheme allows us to use large time steps compared to the size of the oscillations in the solution.     

Highly Selective Zeolite Topologies for Flue Gas Separation

The separation of carbon dioxide from flue gas is essential for the reduction of greenhouse gas emissions. In adsorptive methods, the challenge lies in the choice of suitable porous materials. Among all zeolite topologies, a number of adsorbents with pore dimensions in the range of the guest molecules were identified to allow an excellent separation by diffusion, and MRE and AFO zeolite topologies appear to be the best candidates based on equilibrium adsorption. Also, it was found that the behavior of this gas mixture in DFT and APD zeolites differed from the normal behavior.     

Computational analysis of glycolytic reaction in open spatial reactor

We consider the computational analysis of processes within the spatially-distributed model simulating the glycolytic reaction taking place in the one-side fed open chemical reactor. The main point of the simulation is the decomposition of the reaction–diffusion system into unidirectional reaction in a bulk supplied by feedback terms stated as boundary conditions on the lower boundary of the reactor, i.e. the unique plane where an exchange with an outer medium is possible within the real experimental situation. Analysis of the curvature of the reagents distribution curves proves kinematic character of the observed lateral waves corresponding to the picture of experimentally observed glycolytic traveling waves. At the same time, their origin relates to diffusion of the reagents in a vertical cross-section of the reactor. Study of the solutions for the concerned reaction–diffusion model in the case of stochastically different diffusion coefficients reveals the Turing structures.     

Modeling and numerical simulation of the resistance spot welding of zinc coated steel sheets using rounded tip electrode: Analysis of required conditions

This paper investigates the essential conditions to improve the accuracy of a resistance spot welding computational study of advanced zinc coated steel sheets using rounded tip electrode. An experimental analysis is performed to highlight the required considerations for a suitable simulation. A sequential Electrical-Thermal-Metallurgical and Mechanical (ETMM) finite element analysis with appropriate precautions of the contact conditions enables to accurately simulate the nugget development during the welding. A critical smooth evolution of the contact radius is required. A fine meshing with an interfacial mesh size of at least 0.05 × 10⁻³ m combined with a coupling time step of 0.0025 s between the electrical-thermal-metallurgical and the mechanical analysis allows a regular incrementation of the contact radius, without burdening the time computing. Accurate values of the contact resistance depending on the interfacial pressure and temperature are essential for a good simulation of the nugget size. The ETMM calculation is successfully extended to the simulation of the welding of a typical two sheets assembly.