Effects of Subgrid-Scale Modeling on Time Correlations in Large Eddy Simulation

The effects of the unresolved subgrid-scale (SGS) motions on the energy balance of the resolved scales in large eddy simulation (LES) have been investigated actively because modeling the energy transfer between the resolved and unresolved scales is crucial to constructing accurate SGS models. But the subgrid scales not only modify the energy balance, they also contribute to temporal decorrelation of the resolved scales. The importance of this effect in applications including the predictability problem and the evaluation of sound radiation by turbulent flows motivates the present study of the effect of SGS modeling on turbulent time correlations. This paper compares the two-point, two-time Eulerian velocity correlation in isotropic homogeneous turbulence evaluated by direct numerical simulation (DNS) with the correlations evaluated by LES using a standard spectral eddy viscosity. It proves convenient to express the two-point correlations in terms of spatial Fourier decomposition of the velocity field. The LES fields are more coherent than the DNS fields: their time correlations decay more slowly at all resolved scales of motion and both their integral scales and microscales are larger than those of the DNS field. Filtering alone is not responsible for this effect: in the Fourier representation, the time correlations of the filtered DNS field are identical to those of the DNS field itself. The possibility of modeling the decorrelating effects of the unresolved scales of motion by including a random force in the model is briefly discussed. The results could have applications to the problem of computing sound sources in isotropic homogeneous turbulence by LES

Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Modeling and Simulation of Liquid Pulsed Particulate Fluidized Beds

Pulsed fluidization is of considerable interest in process engineering for improving fluidization quality. Quantitative understanding of the pulsed two-phase flow behaviors is very important for proper design and optimum operation of such contactors. The

Theoretical simulation studies of pulsed field magnetisation of (RE)BCO bulk superconductors

The magnetisation of bulk high temperature superconductors (HTS), such as RE-Ba-Cu-O [(RE)BCO, where RE is a rare earth element or Y], by a practical technique is essential for their application in high field, permanent magnet-like devices. Research to-date into the pulsed field magnetisation (PFM) of these materials, however, has been limited generally to experimental techniques, with relatively little progress in the development of theoretical models. This is because not only is a multi-physics approach needed to take account of the heating of the samples but also the high electric fields generated are well above the regime in which there are reliable experimental results. This paper describes a framework of theoretical simulation using the finite element method (FEM) that is applicable to both single- and multi-pulse magnetisation processes of (RE)BCO bulk superconductors. The model incorporates the heat equation and provides a convenient way of determining the distribution of trapped field, current density and temperature change within a bulk superconductor at each stage of the magnetisation process. An example of the single-pulse magnetisation of a (RE)BCO bulk is described. Potentially, the model may serve as a cost-effective tool for the optimisation of the bulk geometry and the magnetisation profile in multi-pulse magnetisation processes. © 2010 IOP Publishing Ltd.

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Parameter Estimation for Hidden Markov Models with Intractable Likelihoods

Approximate Bayesian computation (ABC) is a popular technique for analysing data for complex models where the likelihood function is intractable. It involves using simulation from the model to approximate the likelihood, with this approximate likelihood then being used to construct an approximate posterior. In this paper, we consider methods that estimate the parameters by maximizing the approximate likelihood used in ABC. We give a theoretical analysis of the asymptotic properties of the resulting estimator. In particular, we derive results analogous to those of consistency and asymptotic normality for standard maximum likelihood estimation. We also discuss how sequential Monte Carlo methods provide a natural method for implementing our likelihood-based ABC procedures.

Micro-scale Mechanics of the Surface-Nanocrystalline Al-Alloy Material

Based on the microscopic observations and measurements, the mechanical behavior of the surface-nanocrystallized Al-alloy material at microscale is investigated experimentally and theoretically. In the experimental research, the compressive stress-strain curves and the hardness depth curves are measured. In the theoretical simulation, based on the material microstructure characteristics and the experimental features of the compression and indentation, the microstructure cell models are developed and the strain gradient plasticity theory is adopted. The material compressive stress-strain curves and the hardness depth curves-are predicted and simulated. Through comparison of the experimental results with the simulation results, the material and model parameters are determined.

Effects Of Subgrid-Scale Modeling On Lagrangian Statistics In Large-Eddy Simulation

The application of large-eddy simulation (LES) to turbulent transport processes requires accurate prediction of the Lagrangian statistics of flow fields. However, in most existing SGS models, no explicit consideration is given to Lagrangian statistics. In this paper, we focus on the effects of SGS modeling on Lagrangian statistics in LES ranging from statistics determining single-particle dispersion to those of pair dispersion and multiparticle dispersion. Lagrangian statistics in homogeneous isotropic turbulence are extracted from direct numerical simulation (DNS) and the LES with a spectral eddy-viscosity model. For the case of longtime single-particle dispersion, it is shown that, compared to DNS, LES overpredicts the time scale of the Lagrangian velocity correlation but underpredicts the Lagrangian velocity fluctuation. These two effects tend to cancel one another leading to an accurate prediction of the longtime turbulent dispersion coefficient. Unlike the single-particle dispersion, LES tends to underestimate significantly the rate of relative dispersion of particle pairs and multiple-particles, when initial separation distances are less than the minimum resolved scale due to the lack of subgrid fluctuations. The overprediction of LES on the time scale of the Lagrangian velocity correlation is further confirmed by a theoretical analysis using a turbulence closure theory.

On the computation of space-time correlations by large-eddy simulation

The effect of subgrid-scale (SGS) modeling on velocity (space-) time correlations is investigated in decaying isotropic turbulence. The performance of several SGS models is evaluated, which shows superiority of the dynamic Smagorinsky model used in conjunction with the multiscale large-eddy simulation (LES) procedure. Compared to the results of direct numerical simulation, LES is shown to underpredict the (un-normalized) correlation magnitude and slightly overpredict the decorrelation time scales. This can lead to inaccurate solutions in applications such as aeroacoustics. The underprediction of correlation functions is particularly severe for higher wavenumber modes which are swept by the most energetic modes. The classic sweeping hypothesis for stationary turbulence is generalized for decaying turbulence and used to analyze the observed discrepancies. Based on this analysis, the time correlations are determined by the wavenumber energy spectra and the sweeping velocity, which is the square root of the total energy. Hence, an accurate prediction of the instantaneous energy spectra is most critical to the accurate computation of time correlations. (C) 2004 American Institute of Physics.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Computer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap

Collisions of a particle pair induced by optical tweezers have been employed to study colloidal stability. In order to deepen insights regarding the collision-sticking dynamics of a particle pair in the optical trap that were observed in experimental approaches at the particle level, the authors carry out a Brownian dynamics simulation. In the simulation, various contributing factors, including the Derjaguin-Landau-Verwey-Overbeek interaction of particles, hydrodynamic interactions, optical trapping forces on the two particles, and the Brownian motion, were all taken into account. The simulation reproduces the tendencies of the accumulated sticking probability during the trapping duration for the trapped particle pair described in our previous study and provides an explanation for why the two entangled particles in the trap experience two different statuses. (c) 2007 American Institute of Physics.

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

Processing Simulation of Virtual Manufacture in Laser Surface Treatmenthang

Processing simulation is at the bottom of the coral technology of VM and is also difficult due to the complexity of mechanism and diversity of parameters. Previously much research has been mainly carried out on the geometrical simulation or physical simulation respectively. The aim of this paper is to study the processing simulation in laser surface treatment based on the mechanism, put forward the architecture of the whole processing simulation and give the models of the processing. As a result the data structure layers in the whole simulation is presented.

Comparisons of static, quasi-static and dynamic 3D porous media scale network models for two-phase immiscible flow in porous media

A dynamic 3D pore-scale network model is formulated for investigating the effect of interfacial tension and oil-water viscosity during chemical flooding. The model takes into account both viscous and capillary forces in analyzing the impact of chemical properties on flow behavior or displacement configuration, while the static model with conventional invasion percolation algorithm incorporates the capillary pressure only. From comparisons of simulation results from these models. it indicates that the static pore scale network model can be used successfully when the capillary number is low. With the capillary increases due to the enhancement of water viscosity or decrease of interfacial tension, only the quasi-static and dynamic model can give insight into the displacement mechanisms.

A complex systems modelling approach for decentralized simulation of electrical microgrids

ICECCS 2010