Fourth order accurate compact scheme with group velocity control (GVC)

For solving complex flow field with multi-scale structure higher order accurate schemes are preferred. Among high order schemes the compact schemes have higher resolving efficiency. When the compact and upwind compact schemes are used to solve aerodynamic problems there are numerical oscillations near the shocks. The reason of oscillation production is because of non-uniform group velocity of wave packets in numerical solutions. For improvement of resolution of the shock a parameter function is introduced in compact scheme to control the group velocity. The newly developed method is simple. It has higher accuracy and less stencil of grid points.

Scheme Construction with Numerical Flux Residual Correction (NFRC) and Group Velocity Control (GVC)

For simulating multi-scale complex flow fields like turbulent flows, the high order accurate schemes are preferred. In this paper, a scheme construction with numerical flux residual correction (NFRC) is presented. Any order accurate difference approximation can be obtained with the NFRC. To improve the resolution of the shock, the constructed schemes are modified with group velocity control (GVC) and weighted group velocity control (WGVC). The method of scheme construction is simple, and it is used to solve practical problems.

Hybrid Finite Compact-WENO Schemes for Shock Calculation

Hybrid finite compact (FC)-WENO schemes are proposed for shock calculations. The two sub-schemes (finite compact difference scheme and WENO scheme) are hybridized by means of the similar treatment as in ENO schemes. The hybrid schemes have the advantages of FC and WENO schemes. One is that they possess the merit of the finite compact difference scheme, which requires only bi-diagonal matrix inversion and can apply the known high-resolution flux to obtain high-performance numerical flux function; another is that they have the high-resolution property of WENO scheme for shock capturing. The numerical results show that FC-WENO schemes have better resolution properties than both FC-ENO schemes and WENO schemes. In addition, some comparisons of FC-ENO and artificial compression method (ACM) filter scheme of Yee et al. are also given.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

Iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact

We present in this paper an iterative method using consistent mass matrix in axisymmetrical finite element analysis of hypervelocity impact. To retain the advantage of integration on an element-by-element basis which is at the heart of modern hydrocodes, we suggest that the first step should be to solve for accelerations at an advanced time step by using the lumped mass approach, then iterate using a consistent mass matrix to improve the estimate. Examples are given to show the improved resolution with the new method.

An elasto-plastic constitutive relation of whisker-reinforced metal-matrix composite

A material model for whisker-reinforced metal-matrix composites is constructed that consists of three kinds of essential elements: elastic medium, equivalent slip system, and fiber-bundle. The heterogeneity of material constituents in position is averaged, while the orientation distribution of whiskers and slip systems is considered in the structure of the material model. Crystal and interface sliding criteria are addressed. Based on the stress-strain response of the model material, an elasto-plastic constitutive relation is derived to discuss the initial and deformation induced anisotropy as well as other fundamental features. Predictions of the present theory for unidirectional-fiber-reinforced aluminum matrix composites are favorably compared with FEM results.

Assessment of Microscale Test Methods of Peeling and Splitting along Surface of Thin-Film/Substrates

Peel test methods are assessed through being applied to a peeling analysis of the ductile film/ceramic substrate system. Through computing the fracture work of the system using the either beam bend model (BB model) or the general plane analysis model (GPA model), surprisingly, a big difference between both model results is found. Although the BB model can capture the plastic dissipation phenomenon for the ductile film case as the GPA model can, it is much sensitive to the choice of the peeling criterion parameters, and it overestimates the plastic bending effect unable to capture crack tip constraint plasticity. In view of the difficulty of measuring interfacial toughness using peel test method when film is the ductile material, a new test method, split test, is recommended and analyzed using the GPA model. The prediction is applied to a wedge-loaded experiment for Al-alloy double-cantilever beam in literature.

Optical Factors Determined by the T-matrix Method in Turbidity Measurement of Absolute Coagulation Rate Constants

Turbidity measurement for the absolute coagulation rate constant of suspensions has been extensively adopted because of its simplicity and easy implementation. A key factor to derive the rate constant from experimental data is how to theoretically evaluate the so-called optical factor involved in calculating the extinction cross section of doublets formed in the aggregation. In a previous paper, we have shown that compared with other theoretical approaches, the T-matrix method provides a robust solution to this problem and is effective in extending the applicability range of the turbidity methodology as well as increasing measurement accuracy. This paper will provide a more comprehensive discussion about the physical insight of using the T-matrix method in turbidity measurement and associated technical details. In particular, the importance of ensuring the correct value for the refractive indices for colloidal particles and the surrounding medium used in the calculation is addressed because the indices generally vary with the wavelength of the incident light. The comparison of calculated results with experiments shows that the T-matrix method can correctly calculate optical factors even for large particles, whereas other existing theories cannot. In addition, the calculated data of the optical factor by the T-matrix method for a range of particle radii and incident light wavelengths are listed.

Onset of Adiabatic Shear Instability in Strain Gradient Dependent Metal Matrix Composites

In this paper, effect of strain gradient on adiabatic shear instability in particle reinforced metal matrix composites is investigated by making use of the strain gradient dependent constitutive equation developed by Dai et al. [9] and the linear perturbation analysis presented by Bai [10]. The results have shown that the onset of adiabatic shear instability in metal matrix composites reinforced with small particles is more prone to occur than in the composites reinforced with large particles. This means that the strain gradient provides a strong deriving force for onset of adiabatic shear instability in metal matrix composites.

Particulate Size Effects in the Particle-Reinforced Metal-Matrix Composites

The influences of I,article size on the mechanical properties of the particulate metal matrix composite;are obviously displayed in the experimental observations. However, the phenomenon can not be predicted directly using the conventional elastic-plastic theory. It is because that no length scale parameters are involved in the conventional theory. In the present research, using the strain gradient plasticity theory, a systematic research of the particle size effect in the particulate metal matrix composite is carried out. The roles of many composite factors, such as: the particle size, the Young's modulus of the particle, the particle aspect ratio and volume fraction, as well as the plastic strain hardening exponent of the matrix material, are studied in detail. In order to obtain a general understanding for the composite behavior, two kinds of particle shapes, ellipsoid and cylinder, are considered to check the strength dependence of the smooth or non-smooth particle surface. Finally, the prediction results will be applied to the several experiments about the ceramic particle-reinforced metal-matrix composites. The material length scale parameter is predicted.

Failure of SiC particulate-reinforced metal matrix composites induced by laser thermal shock

Thermal failure of SiC particulate-reinforced 6061 aluminum alloy composites induced by both laser thermal shock and mechanical load has been investigated. The specimens with a single-edge notch were mechanically polished to 0.25 mm in thickness. The notched-tip region of the specimen is subjected to laser beam rapid heating. In the test, a pulsed Nd:glass laser beam is used with duration 1.0 ms or 250 mu s, intensity 15 or 70 kW/cm(2), and spot size 5.0 mm in diameter. Threshold intensity was tested and fracture behavior was studied. The crack-tip process zone development and the microcrack formation were macroscopically and microscopically observed. It was found that in these materials, the initial crack occurred in the notched-tip region, wherein the initial crack was induced by either void nucleation, growth, and subsequent coalescence of the matrix materials or separation of the SiC particulate-matrix interface. It was further found that the process of the crack propagation occurred by the fracture of the SiC particulates.

Dynamic Effective Properties of Particle-Reinforced Composites with Viscoelastic Matrix

The frequency-dependent dynamic effective properties of the particle-reinforced composites with the viscoelastic matrix are studied. Several equations to predict the effective wavenumber of the coherent plane waves propagating through particle-reinforced

Debonding and cracking energy release rate of the fiber/matrix interface

An investigation of fiber/matrix interfacial fracture energy is presented in this paper. Several existing theoretical expressions for the fracture energy of interfacial debonding are reviewed. For the single-fiber/matrix debonding and pull-out experimental model, a study is carried out on the effect of interfacial residual compressive stress and friction on interface cracking energy release rate.