Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

A time domain computation method for dynamic behavior of mooring system

A quasi-steady time domain method is developed for the prediction of dynamic behavior of a mooring system under the environmental disturbances, such as regular or irregular waves, winds and currents. The mooring forces are obtained in a static sense at each instant. The dynamic feature of the mooring cables can be obtained by incorporating the extended 3-D lumped-mass method with the known ship motion history. Some nonlinear effects, such as the influence of the instantaneous change of the wetted hull surface on the hydrostatic restoring forces and Froude-Krylov forces, are included. The computational results show a satisfactory agreement with the experimental ones.

Statistical Simulation of Rarefied Gas Flows in Micro-Channels

Rarefied gas flows through micro-channels are simulated using particle approaches, named as the information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method. In simulating the low speed flows in long micro-channels the DSMC method encounters the problem of large sample size demand and the difficulty of regulating boundary conditions at the inlet and outlet. Some important computational issues in the calculation of long micro-channel flows by using the IP method, such as the use the conservative form of the mass conservation equation to guarantee the adjustment of the inlet and outlet boundary conditions and the super-relaxation scheme to accelerate the convergence process, are addressed. Stream-wise pressure distributions and mass fluxes through micro-channels given by the IP method agree well with experimental data measured in long micro-channels by Pong et al. (with a height to length ratio of 1.2:3000), Shih et al. (l.2:4800), Arkilic et al. and Arkilic (l.3:7500), respectively. The famous Knudsen minimum of normalized mass flux is observed in IP and DSMC calculations of a short micro-channel over the entire flow regime from continuum to free molecular, whereas the slip Navier-Stokes solution fails to predict it.

Statistical and stochastic description of microvoid evolution observed in dynamic failure of ductile materials

A brief review is presented of statistical approaches on microdamage evolution. An experimental study of statistical microdamage evolution in two ductile materials under dynamic loading is carried out. The observation indicates that there are large differences in size and distribution of microvoids between these two materials. With this phenomenon in mind, kinetic equations governing the nucleation and growth of microvoids in nonlinear rate-dependent materials are combined with the balance law of void number to establish statistical differential equations that describe the evolution of microvoids' number density. The theoretical solution provides a reasonable explanation of the experimentally observed phenomenon. The effects of stochastic fluctuation which is influenced by the inhomogeneous microscopic structure of materials are subsequently examined (i.e. stochastic growth model). Based on the stochastic differential equation, a Fokker-Planck equation which governs the evolution of the transition probability is derived. The analytical solution for the transition probability is then obtained and the effects of stochastic fluctuation is discussed. The statistical and stochastic analyses may provide effective approaches to reveal the physics of damage evolution and dynamic failure process in ductile materials.

Molecular/Cluster Statistical Thermodynamics Methods To Simulate Quasi-Static Deformations At Finite Temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their "seamlessness" or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation. (C) 2008 Elsevier Ltd. All rights reserved.

Statistical Interpolation Method of Turbulent Phase Screen

A relative displacement between the grid points of optical fields and those of phase screens may occur in the simulation of light propagation through the turbulent atmosphere. A statistical interpolator is proposed to solve this problem in this paper. It is evaluated by the phase structure function and numerical experiments of light propagation through atmospheric turbulence with/without adaptive optics (AO) and it is also compared with the well-known linear interpolator under the same condition. Results of the phase structure function show that the statistical interpolator is more accurate in comparison with the linear one, especially in the high frequency region. More importantly, the long-exposure results of light propagation through the turbulent atmosphere with/without AO also show that the statistical interpolator is more accurate and reliable than the linear one. (C) 2009 Optical Society of America.

Avalanche behavior and statistical properties in a microcrack coalescence process

Fracture owing to the coalescence of numerous microcracks can be described by a simple statistical model, where a coalescence event stochastically occurs as the number density of nucleated microcracks increases. Both numerical simulation and statistical analysis reveal that a microcrack coalescence process may display avalanche behavior and that the final failure is catastrophic. The cumulative distribution of coalescence events in the vicinity of critical fracture follows a power law and the fracture profile has self-affine fractal characteristic. Some macromechanical quantities may be traced back and extracted from the mesoscopic process based on the statistical analysis of coalescence events.

Statistical analysis of synthetic earthquake catalogs generated by models with various levels of fault zone disorder

The stress release model, a stochastic version of the elastic rebound theory, is applied to the large events from four synthetic earthquake catalogs generated by models with various levels of disorder in distribution of fault zone strength (Ben-Zion, 1996) They include models with uniform properties (U), a Parkfield-type asperity (A), fractal brittle properties (F), and multi-size-scale heterogeneities (M). The results show that the degree of regularity or predictability in the assumed fault properties, based on both the Akaike information criterion and simulations, follows the order U, F, A, and M, which is in good agreement with that obtained by pattern recognition techniques applied to the full set of synthetic data. Data simulated from the best fitting stress release models reproduce, both visually and in distributional terms, the main features of the original catalogs. The differences in character and the quality of prediction between the four cases are shown to be dependent on two main aspects: the parameter controlling the sensitivity to departures from the mean stress level and the frequency-magnitude distribution, which differs substantially between the four cases. In particular, it is shown that the predictability of the data is strongly affected by the form of frequency-magnitude distribution, being greatly reduced if a pure Gutenburg-Richter form is assumed to hold out to high magnitudes.

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, K-I and K-II, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.

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.

The degree of predictability of earthquakes in several regions of China: statistical analysis of historical data

The stress release model, a stochastic version of the elastic-rebound theory, is applied to the historical earthquake data from three strong earthquake-prone regions of China, including North China, Southwest China, and the Taiwan seismic regions. The results show that the seismicity along a plate boundary (Taiwan) is more active than in intraplate regions (North and Southwest China). The degree of predictability or regularity of seismic events in these seismic regions, based on both the Akaike information criterion (AIC) and fitted sensitivity parameters, follows the order Taiwan, Southwest China, and North China, which is further identified by numerical simulations. (c) 2004 Elsevier Ltd. All rights reserved.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

用Symbolic Computation计算毛细重力波六阶解

本文首先运用Symbolic Computation在半物理平面(x,)上计算了毛细重力波的六阶解,得到了波形与色散关系,低阶解与 Hogan 结果一致。

Analysis of spallation based upon statistical microdamage mechanics

A relationship between the cumulative length of microcracks and the amplitude and duration of tensile impulse in spallation was established based on the application of statistical microdamage mechanics, which included a statistical formulation and dynamic laws of microdamage under loading. Since the degrees of spallation, called incipient, intermediate and complete spallation, can be characterized by the cumulative length of microcracks, a physical interpretation of an empirical criterion to spallation was presented.

Statistical formulation and experimental determination of growth rate of micrometre cracks under impact loading

A new statistical formulation and a relevant experimental approach to determine the growth rate of microcracks were proposed. The method consists of experimental measurements and a statistical analysis' on the basis of the conservation law of number density of microcracks in phase space. As a practical example of the method, the growth rate of microcracks appearing in an aluminium alloy subjected to planar impact loading was determined to be ca. 10 mu m/mu s under a tensile stress of 1470 MPa and load duration between 0.26 mu s and 0.80 mu s.