Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

A ghost-cell immersed boundary method for large-eddy simulations of compressible turbulent flows

A methodology to perform a ghost-cell-based immersed boundary method (GCIBM) is presented for simulating compressible turbulent flows around complex geometries. In this method, the boundary condition on the immersed boundary is enforced through the use of ‘ghost cells’ that are located inside the solid body. The computations of variables on these ghost cells are achieved using linear interpolation schemes. The validity and applicability of the proposed method is verified using a three-dimensional (3D) flow over a circular cylinder, and a large-eddy simulation of fully developed 3D turbulent flow in a channel with a wavy surface. The results agree well with the previous numerical and experimental results, given that the grid resolution is reasonably fine. To demonstrate the capability of the method for higher Mach numbers, supersonic turbulent flow over a circular cylinder is presented. While more work still needs to be done to demonstrate higher robustness and accuracy, the present work provides interesting insights using the GCIBM for the compressible flows.     

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.     

Investigation of the equality constraint effect on the reduction of the rotational ambiguity in three-component system using a novel grid search method

The obtained results by soft modeling multivariate curve resolution methods often are not unique and are questionable because of rotational ambiguity. It means a range of feasible solutions equally fit experimental data and fulfill the constraints. Regarding to chemometric literature, a survey of useful constraints for the reduction of the rotational ambiguity is a big challenge for chemometrician. It is worth to study the effects of applying constraints on the reduction of rotational ambiguity, since it can help us to choose the useful constraints in order to impose in multivariate curve resolution methods for analyzing data sets. In this work, we have investigated the effect of equality constraint on decreasing of the rotational ambiguity. For calculation of all feasible solutions corresponding with known spectrum, a novel systematic grid search method based on Species-based Particle Swarm Optimization is proposed in a three-component system.     

Wave body interaction in 2D using smoothed particle hydrodynamics (SPH) with variable particle mass

Wave interaction with bodies is an important practical application for smoothed particle hydrodynamics (SPH) which in principle applies to steep and breaking waves without special treatment. However, few detailed tests have been undertaken even with small amplitude waves. In order to reduce computer time a variable particle mass distribution is tested here with fine resolution near the body and coarse resolution further away, while maintaining a uniform kernel size. We consider two well‐defined test cases, in two dimensions, of waves generated by a heaving semi‐immersed cylinder and progressive waves interacting with a fixed cylinder. But first, still water with hydrostatic pressure is tested. The open‐source code SPHysics ( http://www.sphysics.org )^{§Update made here after initial online publication.} is used with a Riemann solver in an Arbitrary Lagrangian–Eulerian formulation. For the heaving cylinder, SPH results for far field wave amplitude and cylinder force show good agreement with the data of Yu and Ursell (J. Fluid Mech. 1961; 11 :529–551). For wave loading on a half‐submerged cylinder the agreement with the experimental data of Dixon et al. (J. Waterway Port Coastal Ocean Div. 1979; 105 :421–438) for the root mean square force is within 2%. For more submerged cases, the results show some discrepancy, but this was also found with other modelling approaches. The sensitivity of results to the value of the slope limiter used in the MUSCL‐based Riemann solver is demonstrated. The variable mass distribution leads to a computer run speedup of nearly 200% in these cases. Copyright © 2011 John Wiley & Sons, Ltd.     

Free convection of a dusty‐gas flow along a semi‐infinite vertical cylinder

An analysis is performed to study the free convection of a dusty‐gas flow along a semi‐infinite isothermal vertical cylinder. The governing equations of the flow problem are transformed into non‐dimensional form and the resulting nonlinear, coupled parabolic partial differential equations have been solved numerically using an implicit finite difference scheme of Crank–Nicholson type. The flow variables such as gas–velocity, dust‐particle velocity and temperature, shearing stress and heat transfer coefficients are calculated numerically for various parameters occurring in the problem. It is observed that due to the presence of dust particles, the gas velocity is found to decrease. Copyright © 2009 John Wiley & Sons, Ltd.     

A NEW SYSTEM EXHIBITING A HIGH BARRIER TO ROTATION ABOUT THE P—N BOND. SYNTHESES AND SPECTRAL STUDIES OF 1,3,2-DIAZAPHOSPHETIDIN-4-ONE DERIVATIVES

Abstract

1,3,2-Diazaphosphetidin-4-one derivatives were synthesized. Their ¹H NMR signals at ambient temperature displaed a broad doublet which eventually coalesced before becoming sharp upon warming, and the ¹³C NMR spectra also showed that the α-and β-carbons of the dialkylamino group were non-equivalent. The phenomena were ascribed to restricted rotation about the exocyclic P—N bond.     

The Allen–Cahn equation with dynamic boundary conditions and mass constraints

The Allen–Cahn equation, coupled with dynamic boundary conditions, has recently received a good deal of attention. The new issue of this paper is the setting of a rather general mass constraint, which may involve either the solution inside the domain or its trace on the boundary. The system of nonlinear partial differential equations can be formulated as a variational inequality. The presence of the constraint in the evolution process leads to additional terms in the equation and the boundary condition containing a suitable Lagrange multiplier. A well‐posedness result is proved for the related initial value problem. Copyright © 2014 John Wiley & Sons, Ltd.     

Start‐up of electrophoresis of an arbitrarily oriented dielectric cylinder

An analytical study is presented for the transient electrophoretic response of a circular cylindrical particle to the step application of an electric field. The electric double layer adjacent to the particle surface is thin but finite compared with the radius of the particle. The time‐evolving electroosmotic velocity at the outer boundary of the double layer is utilized as a slip condition so that the transient momentum conservation equation for the bulk fluid flow is solved. Explicit formulas for the unsteady electrophoretic velocity of the particle are obtained for both axially and transversely applied electric fields, and can be linearly superimposed for an arbitrarily‐oriented applied field. If the cylindrical particle is neutrally buoyant in the suspending fluid, the transient electrophoretic velocity is independent of the orientation of the particle relative to the applied electric field and will be in the direction of the applied field. If the particle is different in density from the fluid, then the direction of electrophoresis will not coincide with that of the applied field until the steady state is attained. The growth of the electrophoretic mobility with the elapsed time for a cylindrical particle is substantially slower than for a spherical particle.