基于格子 Bhatnagar-Gross-Krook模型的地震压力波模拟

给出一种新的用于模拟地震压力波的格子Boltzmann模型. 通过使 用Chapman-Enskog展开和多重尺度技术，得到了一系列的格子Boltzmann方程和时间 尺度$t_0$上守恒律，给出了满足地震压力波方程所要求的高阶矩以及简单的平衡态分布 函数表达式. 数值结果表明这种方法可以用来模拟地震压力波.     

Numerical study of the properties of the central moment lattice Boltzmann method

Central moment lattice Boltzmann method (LBM) is one of the more recent developments among the lattice kinetic schemes for computational fluid dynamics. A key element in this approach is the use of central moments to specify the collision process and forcing, and thereby naturally maintaining Galilean invariance, an important characteristic of fluid flows. When the different central moments are relaxed at different rates like in a standard multiple relaxation time (MRT) formulation based on raw moments, it is endowed with a number of desirable physical and numerical features. Because the collision operator exhibits a cascaded structure, this approach is also known as the cascaded LBM. While the cascaded LBM has been developed sometime ago, a systematic study of its numerical properties, such as the accuracy, grid convergence, and stability for well‐defined canonical problems is lacking, and the present work is intended to fulfill this need. We perform a quantitative study of the performance of the cascaded LBM for a set of benchmark problems of differing complexity, viz., Poiseuille flow, decaying Taylor–Green vortex flow, and lid‐driven cavity flow. We first establish its grid convergence and demonstrate second‐order accuracy under diffusive scaling for both the velocity field and its derivatives, that is, the components of the strain rate tensor, as well. The method is shown to quantitatively reproduce steady/unsteady analytical solutions or other numerical results with excellent accuracy. The cascaded MRT LBM based on the central moments is found to be of similar accuracy when compared with the standard MRT LBM based on the raw moments, when a detailed comparison of the flow fields are made, with both reproducing even the small scale vortical features well. Numerical experiments further demonstrate that the central moment MRT LBM results in significant stability improvements when compared with certain existing collision models at moderate additional computational cost. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of a vortex lattice numerical model in the calculation of inviscid incompressible flow around delta wings

The flow over a flat plate delta wing at incidence and in sideslip is studied using vortex lattice models based on streamwise penelling. For the attached flow problem the effect of sideslip is simulated by modifying the standard vortex lattice model for zero sideslip by aligning the trailing vortices aft of the wing along the resultant flow direction. For the separated flow problem a non-linear vortex lattice model is developed for both zero and non-zero sideslip angles in which the shape and position of the leading edge separation vortices are calculated by an iterative procedure starting from an assumed initial shape. The theoretical values are compared with available theoretical and experimental results.     

Incorporation of twinning into a crystal plasticity finite element model: Evolution of lattice strains and texture in Zircaloy-2

A crystal plasticity finite element code is developed to model lattice strains and texture evolution of HCP crystals. The code is implemented to model elastic and plastic deformation considering slip and twinning based plastic deformation. The model accounts for twinning reorientation and growth. Twinning, as well as slip, is considered to follow a rate dependent formulation. The results of the simulations are compared to previously published in situ neutron diffraction data. Experimental results of the evolution of the texture and lattice strains under uniaxial tension/compression loading along the rolling, transverse, and normal direction of a piece of rolled Zircaloy-2 are compared with model predictions. The rate dependent formulation introduced is capable of correctly capturing the influence of slip and twinning deformation on lattice strains as well as texture evolution.     

Parallel computation of turbulent flows using moment base lattice Boltzmann method

This paper describes parallel computing approach for simulating turbulent flows using a moment base lattice Boltzmann method. The distribution functions of the lattice Boltzmann method are expressed by corresponding moments. Choosing proper relaxation times for higher order moments, a minimum numerical dissipation is implicitly added to stabilise the method at high Reynolds numbers. Validation of the method is made by computing free decaying periodic turbulent flows and fully developed turbulent channel flows on a GPU platform. Though the present method requires additional work to calculate the higher order moments, it is shown that additional computational cost is negligible in the GPU computing. The numerical results stably obtained for the turbulent flows are in good agreement with those of a pseudo-spectral method and corresponding DNS database.     

A lattice deformation based model of metallic lattice sandwich plates subjected to impulsive loading

A lattice structure deformation mechanism based theoretical model is developed to predict the dynamic response of square lattice sandwich plates under impulsive loading. The analytical model is established on the basis of the three-stage framework proposed by Fleck and Deshpande (2004). In the first stage, the impulse transmitted from air shock loading to the sandwich plates by fluid-structure interaction is analytically calculated. The lattice core suffers non-uniform compression in the second stage due to the clamped boundary conditions. The structure deformation mechanism is introduced in the lattice core compression and the analytical nominal stress–strain curve of core compression accords well with previous experimental results. In the final stage, the sandwich plate is analyzed as a continuum plate with non-uniform thickness deduced by inconsistent deformation of the front and back sheets.The experiment results of square metallic sandwich plates with tetrahedral lattice core are presented and compared with analytical prediction to validate the theoretical model. Good agreements are found between the predicted and testing results for both the impulse transmitted to the sandwich plates and the maximum deflection of the back face sheet.     

Tuning of band gaps for a two-dimensional piezoelectric phononic crystal with a rectangular lattice

In this paper, the elastic wave propagation in a two-dimensional piezoelectric phononic crystal is studied by considering the mechanic-electric coupling. The generalized eigenvalue equation is obtained by the relation of the mechanic and electric fields as well as the Bloch-Floquet theorem. The band structures of both the in-plane and anti-plane modes are calculated for a rectangular lattice by the planewave expansion method. The effects of the lattice constant ratio and the piezoelectricity with different filling fractions are analyzed. The results show that the largest gap width is not always obtained for a square lattice. In some situations, a rectangular lattice may generate larger gaps. The band gap characteristics are influenced obviously by the piezoelectricity with the larger lattice constant ratios and the filling fractions.     

Hydroelastic interaction of a wedge-shaped construction entering a liquid

The aim of the paper is to create a mathematical model of the hydroelastic interaction of a wedge-shaped construction with a liquid that it enters. The various relevant factors — velocity and acceleration of elastic vibrations, influence of catinary forces, initial bending, and mobility of brackets — are to be taken into account as fully as possible.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 118–125, May–June, 1993.     

Continuum slip viscoplasticity with the Haasen constitutive model: Application to CdTe single crystal inelasticity

Small strain constitutive equations are developed for the thermomechanical behavior of semiconductor single crystals, including dislocation density as an evolving parameter. The model of Haasen, Alexander and coworkers is modified (extended) to include evolution of coefficients in the definition of internal stress. These account for an evolving dislocation substructure. The resulting model is applied in a continuum slip framework to allow multiple slip orientations. Slip system interaction rules are adapted to include slip system interaction for multiple slip conditions and to suppress secondary slip and dislocation density generation for single slip orientations. The approach is discussed relative to other models for viscoplasticity of single crystals and is examined in the context of thermodynamics with internal state variables. The framework is used to correlate experimental data from compression tests of single crystals of the compound semiconductor CdTe from room temperature to near the melting point. Sensitivity of the model to uncertainties such as initial dislocation density is explored.     

Crucial properties of the moment closure model FENE-QE

We investigate four crucial properties for testing and evaluating a moment closure approximation of the FENE dumbbell model for dilute polymer solutions: non-negative configuration distribution function, energy dissipation, accuracy of approximation and computational expense. Through mathematical analysis, numerical experiments and comparisons with closure model FENE-P and FENE-YDL, we prove that the FENE-QE approximation has non-negative configuration distribution function, approximates the energy dissipation behavior of original kinetic theory and provides good accuracy. To improve the efficiency of this closure approximation, we introduce a piecewise linear approximation technique that greatly reduces the computational cost. This extension of FENE-QE, FENE-QE-PLA, is the closure model we recommend for simulating dilute polymer solutions.     

Intrinsic instability of the lattice BGK model

Based on the stability analysis with no linearization and expansion, it is argued that instability in the lattice BGK model is originated from the linear relaxation hypothesis of collision in the model. The hypothesis stands up only when the deviation from the local equilibrium is weak. In this case the computation is absolutely stable for real fluids. But for flows of high Reynolds number, this hypothesis is violated and then instability takes place physically. By performing a transformation a quantified stability criteria is put forward without those approximation. From the criteria a sufficient condition for stability can be obtained and serve as an estimation of the limited Reynolds number as high as possible.     

Analytical model for a rigid, linear-hardening plastic tube subjected to impulsive torsional moment

An analytical model of a rigid, linear-hardening plastic tube subjected to impulsive torsion at its free end is suggested. It is shown that an impulsive torsional moment at the free end of the tube gives rise to a torsional disturbance of strong discontinuity that propagate at a constant speed and the plastic deformation is concentrated at the front of the disturbances. Two examples are given to illustrate the characterization of the structure under the intense dynamic pure torsional loading condition.     

Implementation of a high order lattice spring model for elasticity

In a lattice spring model (LSM), the material is discretised into particles linked by springs. However, LSMs always adopt linear springs, which results in a stiff approximation of the corresponding elastic solution. In this work, a high order LSM is proposed to overcome this limitation by introducing additional degrees of freedoms (DOFs) to the particles. Based on the energy minimisation principle and the local strain technique, equations for the stiffness matrices of high order LSM are derived. Relationships between micro spring parameters and macro material constants are derived from the Cauchy-born rules and the hyperelastic theory. Numerical examples show that the high order LSM can provide a better solution than that of the linear LSM and that the LSM is more suitable for modelling singularity and fracture problems.     

Simulation of shockwave propagation with a thermal lattice Boltzmann model

A two‐dimensional 19‐velocity (D2Q19) lattice Boltzmann model which satisfies the conservation laws governing the macroscopic and microscopic mass, momentum and energy with local equilibrium distribution order O(u⁴) rather than the usual O(u³) has been developed. This model is applied to simulate the reflection of shockwaves on the surface of a triangular obstacle. Good qualitative agreement between the numerical predictions and experimental measurements is obtained. As the model contains the higher‐order terms in the local equilibrium distribution, it performs much better in terms of numerical accuracy and stability than the earlier 13‐velocity models with the local equilibrium distribution accurate only up to the second order in the velocity u. Copyright © 2003 John Wiley & Sons, Ltd.     

A dislocation density based crystal plasticity finite element model: Application to a two-phase polycrystalline HCP/BCC composites

We present a multiscale model for anisotropic, elasto-plastic, rate- and temperature-sensitive deformation of polycrystalline aggregates to large plastic strains. The model accounts for a dislocation-based hardening law for multiple slip modes and links a single-crystal to a polycrystalline response using a crystal plasticity finite element based homogenization. It is capable of predicting local stress and strain fields based on evolving microstructure including the explicit evolution of dislocation density and crystallographic grain reorientation. We apply the model to simulate monotonic mechanical response of a hexagonal close-packed metal, zirconium (Zr), and a body-centered cubic metal, niobium (Nb), and study the texture evolution and deformation mechanisms in a two-phase Zr/Nb layered composite under severe plastic deformation. The model predicts well the texture in both co-deforming phases to very large plastic strains. In addition, it offers insights into the active slip systems underlying texture evolution, indicating that the observed textures develop by a combination of prismatic, pyramidal, and anomalous basal slip in Zr and primarily {110}〈111〉 slip and secondly {112}〈111〉 slip in Nb.     

Numerical simulation of internal gravity waves using a lattice gas model

Internal waves are modelled in two different circumstances: in a continuously stratified fluid and at the interface between two immiscible fluids. This is done using the lattice gas approach. The standard single phase model and an immiscible two-phase model are both modified to incorporate gravitational interactions. Standing internal waves are set up in both models and are seen to oscillate under the action of the gravitational interaction. The results obtained suggest that the lattice gas approach can be a useful tool in the modelling of such phenomena. © 1998 John Wiley & Sons, Ltd.     

Application of a contact model due to L. A. galin to problems of contact interaction between stringers and massive deformable bodies

L. A. Galin’s contact model for a narrow beam bending on an elastic half-space and Melan’s contact model for a stringer are used to consider two problems of contact interaction between one or two identical symmetrically loaded stringers with small rectangular cross-sections and an elastic half-space. The basic characteristics of these problems are expressed by explicit formulas, and the results of their numerical analysis are given as well.     

Quasi-one-dimensional model of the rod-target interaction

Numerical techniques are proposed for determining the integral characteristics of penetration of a rod into an target. An algorithm for solving two-dimensional elastoplastic problems is employed. To construct the solution, a one-dimensional finite-element column is used (a two-dimensional domain is replaced by a one-dimensional domain). Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 5, pp. 205–210, September–October, 2000.