Application of a line implicit method to fully coupled system of equations for turbulent flow problems

For unstructured finite volume methods, we present a line implicit Runge–Kutta method applied as smoother in an agglomerated multigrid algorithm to significantly improve the reliability and convergence rate to approximate steady-state solutions of the Reynolds-averaged Navier–Stokes equations. To describe turbulence, we consider a one-equation Spalart–Allmaras turbulence model. The line implicit Runge–Kutta method extends a basic explicit Runge–Kutta method by a preconditioner given by an approximate derivative of the residual function. The approximate derivative is only constructed along predetermined lines which resolve anisotropies in the given grid. Therefore, the method is a canonical generalisation of point implicit methods. Numerical examples demonstrate the improvements of the line implicit Runge–Kutta when compared with explicit Runge–Kutta methods accelerated with local time stepping.     

Toward low‐noise synthetic turbulent inflow conditions for aeroacoustic calculations

The accuracy of boundary conditions for computational aeroacoustics is a well‐known challenge, due in part to the necessity of truncating the flow domain and replacing the analytical boundary conditions at infinity with numerical boundary conditions. In particular, the inflow boundary condition involving turbulent velocity or scalar fields is likely to introduce spurious waves into the domain, therefore degrading the flow behavior and deteriorating the physical acoustic waves. In this work, a method to generate low‐noise, divergence‐free, synthetic turbulence for inflow boundary conditions is proposed. It relies on the classical view of turbulence as a superposition of random eddies convected with the mean flow. Within the proposed model, the vector potential and the requirement that the individual eddies must satisfy the linearized momentum equations about the mean flow are used. The model is tested using isolated eddies convected through the inflow boundary and an experimental benchmark data for spatially decaying isotropic turbulence. Copyright © 2013 John Wiley & Sons, Ltd.     

Residual‐based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods

In this paper, we present a finite element method with a residual‐based artificial viscosity for simulation of turbulent compressible flow, with adaptive mesh refinement based on a posteriori error estimation with sensitivity information from an associated dual problem. The artificial viscosity acts as a numerical stabilization, as shock capturing, and as turbulence capturing for large eddy simulation of turbulent flow. The adaptive method resolves parts of the flow indicated by the a posteriori error estimates but leaves shocks and turbulence under‐resolved in a large eddy simulation. The method is tested for examples in 2D and 3D and is validated against experimental data. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Modelling effects of subgrid-scale mixture fraction variance in LES of a piloted diffusion flame

A posteriori analysis of the statistics of two large-eddy simulation (LES) solutions describing a piloted methane–air (Sandia D) flame is performed on a series of grids with progressively increased resolution reaching about 10.5 million cells. Chemical compositions, density and temperature fields are modelled with a steady flamelet approach and parametrised by the mixture fraction. The difference between the LES solutions arises from a different numerical treatment of the subgrid scale (SGS) mixture fraction variance – an important quantity of interest in non-premixed combustion modelling. In the first case (model I), the variance transport equation is solved directly, while in the second (model II), an equation for the square of the mixture fraction is solved, and the variance is computed from its definition. The comparison of the LES solutions is based on the convergence properties of their statistics with respect to the turbulence resolution length scale. The dependence of the LES statistics is analysed for velocity and the mixture fraction fields, and tested for convergence. For the most part, the statistics converge for the finest grids, but the variance of the mixture fraction shows some residual grid dependence in the high-gradient regions of the jet near field. The SGS variance given by model I exhibits realisability everywhere, whereas in regions of the flame model II is non-realisable, predicting negative variances. Furthermore, the LES statistics of model I exhibit superior convergence behaviour.     

Experimental investigation on the effects of the standoff distance and the initial radius on the dynamics of a single bubble near a rigid wall in an ultrasonic field

Bubble behaviors near a boundary in an ultrasonic field are the fundamental forms of acoustic cavitation and of substantial importance in various applications, such as industry cleaning, chemical engineering and food processing. The effects of two important factors that strongly affect the dynamics of a single acoustic cavitation bubble, namely, the initial bubble radius and the standoff distance, were investigated in this work. The temporal evolution of the bubble was recorded using high speed microphotography. Meanwhile, the time of bubble collapse and the characteristics of the liquid jets were analyzed. The results demonstrate that the intensity of the acoustic cavitation, which is characterized by the time of bubble collapse and the liquid jet speed, reaches the optimum level under suitable values of the initial bubble radius and the normalized standoff distance. As the initial bubble radius and the normalized standoff distance increase or decrease from the optimal values, the time of the bubble collapse increases, and the first liquid jet’s speed decreases substantially, whereas the speeds of the second and third liquid jets exhibit no substantial changes. These results on bubble dynamics in an ultrasonic field are important for identifying or correcting the mechanisms of acoustic cavitation and for facilitating its optimization and application.     

Numerical modeling of chlorine concentration in water storage tanks

In this paper, we describe a numerical model to simulate the evolution in time of the hydrodynamics of water storage tanks, with particular emphasis on the time evolution of chlorine concentration. The mathematical model contains several ingredients particularly designed for this problem, namely, a boundary condition to model falling jets on free surfaces, an arbitrary Lagrangian–Eulerian formulation to account for the motion of the free surface because of demand and supply of water, and a coupling of the hydrodynamics with a convection–diffusion–reaction equation modeling the time evolution of chlorine. From the numerical point of view, the equations resulting from the mathematical model are approximated using a finite element formulation, with linear continuous interpolations on tetrahedra for all the unknowns. To make it possible, and also to be able to deal with convection‐dominated flows, a stabilized formulation is used. In order to capture the sharp gradients present in the chlorine concentration, particularly near the injection zone, a discontinuity capturing technique is employed. Copyright © 2015 John Wiley & Sons, Ltd.     

A puzzle of the radio-quiet quasar PG 1407+265: Are optical and X-ray emissions produced from a relativistic jet?

PG 1407+265 is a radio quiet quasar that has a relativistic jet.In this report,we show some peculiar properties of its optical and X-ray emissions,which indicate possible non-thermal origins produced from the jet.We use a simple synchrotron + synchrotron self Compton (SSC) model to fit the emissions with different ratios of energy densities between the magnetic field and electrons (η≡U B/U e),which predicts a different γ-ray luminosity.The First LAT AGN Catalog (1LAC) did not include PG 1407+265,which indic...