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.     

Structural transformation from shish‐kebab crystals to micro‐fibrils through hot stretching process of gel‐spun ultra‐high molecular weight polyethylene fibers with high concentration solution

Structural evolution of gel‐spun ultra‐high molecular weight polyethylene fibers with high concentration solution via hot stretching process was investigated by in situ small‐angle X‐ray scattering, in situ wide‐angle X‐ray diffraction measurements, scanning electron microscopy, and differential scanning calorimetry. With the increase of stretching strain, the long period continuously increases at relative lower stretching temperature, while it first increases and then decreases rapidly at relative higher stretching temperature. The kebab thickness almost keeps constant during the whole hot‐stretching process and the kebab diameter continually decreases for all stretching temperatures. Moreover, the length of shish decreases slightly and the shish quantity increases although there is almost no change in the diameter of shish crystals during the hot stretching process. The degree of crystal orientation at different temperatures is as high as above 0.9 during the whole stretching process. These results indicate that the shish‐kebab crystals in ultra‐high molecular weight polyethylene fibers can transform continuously into the micro‐fibril structure composed mostly of shish crystals through the hot stretching process. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 225–238     

Mixed convection unsteady stagnation point flow over a stretching sheet with heat transfer in the presence of variable free stream

An analysis is performed for the unsteady mixed convection flow of an incompressible viscous fluid about a stagnation point on a stretching sheet in the presence of a variable free stream. The equations of motion and energy are transformed into the ordinary differential equations by using similarity transformations. Homotopy analysis method is used for the solution of the governing problem. The results have been discussed by plots. The present values of the function are shown very close to the previous limiting solutions. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Study of Polymer-Cement Composite Containing Portland Cement and Aqueous Poly (methyl methacrylate) Latex Polymer by Fourier-Transform Infrared (FT-IR) Spectroscopy

The polymer-cement composite is an immiscible band having modified morphology. Polymer phase improves the toughness of brittle cement and has a reinforcing effect, and interactions between cement components and polymers provide stability. In this study, polymer latex [poly(methyl methacrylate), PMMA] – cement paste was prepared with different weight proportions of polymer content. The main scientific objective of this report has been the analysis and the prediction of high strength with the aim of understanding the interactions of PMMA polymer with Portland cement after hydration reaction of cement pastes and curing the specimens for 28 days. The samples were studied by Fourier transform infrared (FT-IR) spectroscopy. The results are critically examined. It is shown that bands support the interaction of PMMA with cement in the composite.     

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.     

高β相聚偏氟乙烯基复合体系的制备及压电性能

为了提高聚偏氟乙烯(PVDF)的压电性能,需要寻找有效的途径来提高PVDF的电活性相(β相)含量。 通过水热法成功合成了Ag、ZnO以及二者复合(Ag-ZnO)的3种类型纳米粒子,并与PVDF共混形成PVDF复合薄膜。 通过表征PVDF复合材料的形貌,结晶性能和压电性能,可以发现Ag-ZnO复合纳米粒子的协同作用可以有效提高PVDF的结晶性能和压电性能。 此外,通过单轴拉伸可以使得所有PVDF膜的β相含量得到进一步提高,其中拉伸后的PVDF/Ag-ZnO纳米颗粒(P-C)的β相物质的量分数最高,达到70.0％,最佳的压电系数(d₃₃)达到了31.0 pC/N。