A Lagrangian lattice Boltzmann method for Euler equations

A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.

A Lagrangian for von Karman Equations of Large Deflection Problem of Thin Circular Plate

By the semi-inverse method proposed by He, a Lagrangian is established for the large deflection problem of thin circular plate. Ritz method is used to obtain an approximate analytical solution of the problem. First order approximate solution is obtained, which is similar to those in open literature. By Mathematica a more accurate solution can be deduced.

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.

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.

Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence

A scale-similarity model for Lagrangian two-point, two-time velocity correlations LVCs in isotropic turbulence is developed from the Kolmogorov similarity hypothesis. It is a second approximation to the isocontours of LVCs, while the Smith-Hay model is only a first approximation. This model expresses the LVC by its space correlation and a dispersion velocity. We derive the analytical expression for the dispersion velocity from the Navier-Stokes equations using the quasinormality assumption. The dispersion velocity is dependent on enstrophy spectra and shown to be smaller than the sweeping velocity for the Eulerian velocity correlation. Therefore, the Lagrangian decorrelation process is slower than the Eulerian decorrelation process. The data from direct numerical simulation of isotropic turbulence support the scale-similarity model: the LVCs for different space separations collapse into a universal form when plotted against the separation axis defined by the model.

Subgrid-scale contributions to Lagrangian time correlations in isotropic turbulence

The application of large-eddy simulation (LES) to particle-laden turbulence raises such a fundamental question as whether the LES with a subgrid scale (SGS) model can correctly predict Lagrangian time correlations (LTCs). Most of the currently existing SGS models are constructed based on the energy budget equations. Therefore, they are able to correctly predict energy spectra, but they may not ensure the correct prediction on the LTCs. Previous researches investigated the effect of the SGS modeling on the Eulerian time correlations. This paper is devoted to study the LTCs in LES. A direct numerical simulation (DNS) and the LES with a spectral eddy viscosity model are performed for isotropic turbulence and the LTCs are calculated using the passive vector method. Both a priori and a posteriori tests are carried out. It is observed that the subgrid-scale contributions to the LTCs cannot be simply ignored and the LES overpredicts the LTCs than the DNS. It is concluded from the straining hypothesis that an accurate prediction of enstrophy spectra is most critical to the prediction of the LTCs.

Some new bounds on eigenvalues of the Hadamard product and the Fan product of matrices

Let A and B be nonsingular M-matrices. A lower bound on the minimum eigenvalue q(B circle A(-1)) for the Hadamard product of A(-1) and B, and a lower bound on the minimum eigenvalue q(A star B) for the Fan product of A and B are given. In addition, an upper bound on the spectral radius rho(A circle B) of nonnegative matrices A and B is also obtained. These bounds improve several existing results in some cases and the estimating formulas are easier to calculate for they are only depending on the entries of matrices A and B. (C) 2009 Elsevier Inc. All rights reserved.

Selection for variables in quantitative structure-activity/property relationship studies - Comparison between methods of orthogonal descriptors and leaps-and-bounds regression analysis

In this paper, the comparison of orthogonal descriptors and Leaps-and-Bounds regression analysis is performed. The results obtained by using orthogonal descriptors are better than that obtained by using Leaps-and-Bounds regression for the data set of nitrobenzenes used in this study. Leaps-and-Bounds regression can be used effectively for selection of variables in quantitative structure-activity/property relationship(QSAR/QSPR) studies. Consequently, orthogonalisation of descriptors is also a good method for variable selection for studies on QSAR/QSPR.

Studies on the method of leaps and bounds regression

In this paper, we introduce the method of leaps and bounds regression which can be used to select variables quickly and obtain the best regression models. These models contain one variable, two variables, three variables and so on. The results obtained by using leaps and bounds regression were compared with those achieved by using stepwise regression to lead to the conclusion that leaps and bounds regression is an effective method.

Flow properties of a dusty-gas point source in a supersonic free stream

By using Lagrangian method, the flow properties of a dusty-gas point source in a supersonic free stream were studied and the particle parameters in the near-symmetry-axis region were obtained. It is demonstrated that fairly inertial particles travel along oscillating and intersecting trajectories between the bow and termination shock waves. In this region,formation of "multi-layer structure" in particle distribution with alternating low- and highdensity layers is revealed. Moreover, sharp accumulation of particles occurs near the envelopes of particle trajectories.

裂纹在冲击载荷作用下起裂的临界载荷面

采用有限元方法研究裂纹在I型短脉冲载荷作用下应力强度因子随时间的变化，用应力强度 因子的初始上升时间T_r对时间坐标无量纲化，对应力强度因子初始上升段进行曲线拟合， 得到了上升段的曲线表达式。运用简单弹性梁理论和Lagrangian运动方程，获得载荷与时间对裂纹作用的关系式，结合有限元的结果，得到了上升时间T_r的计算表达式，并进一步推 出了裂纹在冲击载荷作用下起裂的临界载荷面。