Shocked Flows Induced by Supersonic Projectiles Moving in Tubes

A numerical study on shocked flows induced by a supersonic projectile moving in tubes is described in this paper. The dispersion-controlled scheme was adopted to solve the Euler equations implemented with moving boundary conditions. Four test cases were carried out in the present study: the first two cases are for validation of numerical algorithms and verification of moving boundary conditions, and the last two cases are for investigation into wave dynamic processes induced by the projectile moving at Mach numbers of M-p = 2.0 and 2.4, respectively, in a short time duration after the projectile was released from a shock tube into a big chamber. It was found that complex shock phenomena exist in the shocked flow, resulting from shock-wave/projectile interaction, shock-wave focusing, shock-wave reflection and shock-wave/contact-surface interactions, from which turbulence and vortices may be generated. This is a fundamental study on complex shock phenomena, and is also a useful investigation for understanding on shocked flows in the ram accelerator that may provide a highly efficient facility for launching hypersonic projectiles.

Pattern Recognition With Weighted Complex Networks

In this paper we introduce a weighted complex networks model to investigate and recognize structures of patterns. The regular treating in pattern recognition models is to describe each pattern as a high-dimensional vector which however is insufficient to express the structural information. Thus, a number of methods are developed to extract the structural information, such as different feature extraction algorithms used in pre-processing steps, or the local receptive fields in convolutional networks. In our model, each pattern is attributed to a weighted complex network, whose topology represents the structure of that pattern. Based upon the training samples, we get several prototypal complex networks which could stand for the general structural characteristics of patterns in different categories. We use these prototypal networks to recognize the unknown patterns. It is an attempt to use complex networks in pattern recognition, and our result shows the potential for real-world pattern recognition. A spatial parameter is introduced to get the optimal recognition accuracy, and it remains constant insensitive to the amount of training samples. We have discussed the interesting properties of the prototypal networks. An approximate linear relation is found between the strength and color of vertexes, in which we could compare the structural difference between each category. We have visualized these prototypal networks to show that their topology indeed represents the common characteristics of patterns. We have also shown that the asymmetric strength distribution in these prototypal networks brings high robustness for recognition. Our study may cast a light on understanding the mechanism of the biologic neuronal systems in object recognition as well.

Thermocapillary Motion Of Droplets At Large Marangoni Numbers

In this paper, the thermocapillary motion problem of drops is investigated using the axisymmetric model. The front-tracking method is employed to capture the drop interface. We find that the migration velocity of the drop is greatly influenced by the temperature field in the drop when Ma is fairly large (>100), which leads to an increase-decrease migration velocity at the beginning of our simulations. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Numerical investigation of thermocapillary migration of the drop for large Marangoni numbers

An axisymmetric model is adopted to simulate the problem of unsteady drop thermocapillary motion for large Marangoni numbers. Front tracking methods are used in the investigation. It is found that the non-dimensional drop migration velocity will decrease with increasing Marangoni number. This agrees well with the experimental results obtained from the 4th Shen-Zhou space ship. In the meanwhile, this is also the first time for numerical simulations to verify the experimental phenomenon under large Marangoni numbers.

A high order accurate difference scheme for complex flow fields

A high order accurate finite difference method for direct numerical simulation of coherent structure in the mixing layers is presented. The reason for oscillation production in numerical solutions is analyzed, It is caused by a nonuniform group velocity of wavepackets. A method of group velocity control for the improvement of the shock resolution is presented. In numerical simulation the fifth-order accurate upwind compact difference relation is used to approximate the derivatives in the convection terms of the compressible N-S equations, a sixth-order accurate symmetric compact difference relation is used to approximate the viscous terms, and a three-stage R-K method is used to advance in time. In order to improve the shock resolution the scheme is reconstructed with the method of diffusion analogy which is used to control the group velocity of wavepackets. (C) 1997 Academic Press.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Complex stress functions and plastic zone sizes for notch cracks subjected to various loading conditions

In this paper, the conformal mapping method is used to solve the plane problem of an infinite plate containing a central lip-shaped notch subjected to biaxial loading at a remote boundary or a surface uniform pressure on the notch. The stress intensity factors KI and KII are obtained by the derived complex stress functions. The simple analytical expressions can be applied to the situation of cracks originating from a circular or an elliptical notch. The plastic zone sizes for such notch cracks are subsequently evaluated in light of the Dugdale strip yield concept. The results are consistent with available numerical data.

“Deborah Numbers”, Coupling Multiple Space and Time Scales and Governing Damage Evolution to Failure

Two different spatial levels are involved concerning damage accumulation to eventual failure. nucleation and growth rates of microdamage nN* and V*. It is found that the trans-scale length ratio c*/L does not directly affect the process. Instead, two independent dimensionless numbers: the trans-scale one * * ( V*)including the * **5 * N c V including mesoscopic parameters only, play the key role in the process of damage accumulation to failure. The above implies that there are three time scales involved in the process: the macroscopic imposed time scale tim = /a and two meso-scopic time scales, nucleation and growth of damage, (* *4) N N t =1 n c and tV=c*/V*. Clearly, the dimensionless number De*=tV/tim refers to the ratio of microdamage growth time scale over the macroscopically imposed time scale. So, analogous to the definition of Deborah number as the ratio of relaxation time over external one in rheology. Let De be the imposed Deborah number while De represents the competition and coupling between the microdamage growth and the macroscopically imposed wave loading. In stress-wave induced tensile failure (spallation) De* < 1, this means that microdamage has enough time to grow during the macroscopic wave loading. Thus, the microdamage growth appears to be the predominate mechanism governing the failure. Moreover, the dimensionless number D* = tV/tN characterizes the ratio of two intrinsic mesoscopic time scales: growth over nucleation. Similarly let D be the “intrinsic Deborah number”. Both time scales are relevant to intrinsic relaxation rather than imposed one. Furthermore, the intrinsic Deborah number D* implies a certain characteristic damage. In particular, it is derived that D* is a proper indicator of macroscopic critical damage to damage localization, like D* ∼ (10–3~10–2) in spallation. More importantly, we found that this small intrinsic Deborah number D* indicates the energy partition of microdamage dissipation over bulk plastic work. This explains why spallation can not be formulated by macroscopic energy criterion and must be treated by multi-scale analysis.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

MODIFIED DIRECT TWOS-COMPLEMENT PARALLEL ARRAY MULTIPLICATION ALGORITHM FOR COMPLEX MATRIX OPERATION

A direct twos-complement parallel array multiplication algorithm is introduced and modified for digital optical numerical computation. The modified version overcomes the problems encountered in the conventional optical twos-complement algorithm. In the array, all the summands are generated in parallel, and the relevant summands having the same weights are added simultaneously without carries, resulting in the product expressed in a mixed twos-complement system. In a two-stage array, complex multiplication is possible with using four real subarrays. Furthermore, with a three-stage array architecture, complex matrix operation is straightforwardly accomplished. In the experiment, parallel two-stage array complex multiplication with liquid-crystal panels is demonstrated.

Fuzzy inferring mathematical morphology and optical implementation

Fuzzification is introduced into gray-scale mathematical morphology by using two-input one-output fuzzy rule-based inference systems. The fuzzy inferring dilation or erosion is defined from the approximate reasoning of the two consequences of a dilation or an erosion and an extended rank-order operation. The fuzzy inference systems with numbers of rules and fuzzy membership functions are further reduced to a simple fuzzy system formulated by only an exponential two-input one-output function. Such a one-function fuzzy inference system is able to approach complex fuzzy inference systems by using two specified parameters within it-a proportion to characterize the fuzzy degree and an exponent to depict the nonlinearity in the inferring. The proposed fuzzy inferring morphological operators tend to keep the object details comparable to the structuring element and to smooth the conventional morphological operations. Based on digital area coding of a gray-scale image, incoherently optical correlation for neighboring connection, and optical thresholding for rank-order operations, a fuzzy inference system can be realized optically in parallel. (C) 1996 Society of Photo-Optical Instrumentation Engineers.

Three-dimensional superresolution by three-zone complex pupil filters

Complex pupil filters are introduced to improve the three-dimensional resolving power of an optical imaging system. Through the design of the essential parameters of such filters, the transmittance and radius of the first zone, three-dimensional superresolution is realized. The Strehl ratio and the transverse and axial gains of such filters are analyzed in detail. A series of simulation examples of such filters are also presented that prove that three-dimensional superresolution can be realized. The advantage of such filters is that it is easy to realize three-dimensional superresolution, and the disadvantage is that the sidelobes of the axial intensity distribution are too high. But this can be overcome by the application of a confocal system. (C) 2005 Optical Society of America.