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 Dynamic Comprehensive Method for Landslide Control

A slope failure is developed due to progressive external loads and deteriorations of slope geomaterials, thus forming a progressive and dynamic development and occurrence of landslides. Site geological properties and other active factors such as hydrodynamic load and human activities are complex and usually unknown, thus this dynamic development and occurrence of landslides can only be understood through the progressive accumulation of knowledge on the landslides. For such a progressive process, this paper proposes a dynamic comprehensive control method for landslide control. This control method takes full advantage of updated monitoring data and site investigations of landslides, and emphasizes the implementation of possible measures for landslide control at reasonable stages and in different groups. These measures are to prevent the occurrence of a landslide disaster. As a case study, a landslide project at the Panluo open-pit iron mine is analyzed to illustrate this dynamic comprehensive control method.

A generalized self-consistent method accounting for fiber section shape

A three-phase confocal elliptical cylinder model is proposed for fiber-reinforced composites, in terms of which a generalized self-consistent method is developed for fiber-reinforced composites accounting for variations in fiber section shapes and randomness in fiber section orientation. The reasonableness of the fiber distribution function in the present model is shown. The dilute, self-consistent, differential and Mori-Tanaka methods are also extended to consider randomness in fiber section orientation in a statistical sense. A full comparison is made between various micromechanics methods and with the Hashin and Shtrikman's bounds. The present method provides convergent and reasonable results for a full range of variations in fiber section shapes (from circular fibers to ribbons), for a complete spectrum of the fiber volume fraction (from 0 to 1, and the latter limit shows the correct asymptotic behavior in the fully packed case) and for extreme types of the inclusion phases (from voids to rigid inclusions). A very different dependence of the five effective moduli on fiber section shapes is theoretically predicted, and it provides a reasonable explanation on the poor correlation between previous theory and experiment in the case of longitudinal shear modulus.

An Approximate Method For Calculation Of Fluid Force And Response Of A Circular Cylinder At Lock-In

In this paper, equations calculating lift force of a rigid circular cyclinder at lock-in uniform flow are deduced in detail. Besides, equations calculating the lift force on a long flexible circular cyclinder at lock-in are deduced based on mode analysis of a multi-degree freedom system. The simplified forms of these equations are also given. Furthermore, an approximate method to predict the forces and response of rigid circular cyclinders and long flexible circular cyclinders at lock-in is introduced in the case of low mass-damping ratio. A method to eliminate one deficiency of these equations is introduced. Comparison with experimental results show the effectiveness of this approximate method.

An Improved Energy Method For Determining Young'S Modulus By Instrumented Indentation Using A Berkovich Tip

We previously proposed a method for estimating Young's modulus from instrumented nanoindentation data based on a model assuming that the indenter had a spherical-capped Berkovich geometry to take account of the bluntness effect. The method is now further improved by releasing the constraint on the tip shape, allowing it to have a much broader arbitrariness to range from a conical-tipped shape to a flat-ended shape, whereas the spherical-capped shape is just a special case in between. This method requires two parameters to specify a tip geometry, namely, a volume bluntness ratio V-r and a height bluntness ratio h(r). A set of functional relationships correlating nominal hardness/reduced elastic modulus ratio (H-n/E-r) and elastic work/total work ratio (W-e/W) were established based on dimensional analysis and finite element simulations, with each relationship specified by a set of V-r and h(r). Young's modulus of an indented material can be estimated from these relationships. The method was shown to be valid when applied to S45C carbon steel and 6061 aluminum alloy.

Comparison Of The Quasi-Static Method And The Dynamic Method For Simulating Fracture Processes In Concrete

Concrete is heterogeneous and usually described as a three-phase material, where matrix, aggregate and interface are distinguished. To take this heterogeneity into consideration, the Generalized Beam (GB) lattice model is adopted. The GB lattice model is much more computationally efficient than the beam lattice model. Numerical procedures of both quasi-static method and dynamic method are developed to simulate fracture processes in uniaxial tensile tests conducted on a concrete panel. Cases of different loading rates are compared with the quasi-static case. It is found that the inertia effect due to load increasing becomes less important and can be ignored with the loading rate decreasing, but the inertia effect due to unstable crack propagation remains considerable no matter how low the loading rate is. Therefore, an unrealistic result will be obtained if a fracture process including unstable cracking is simulated by the quasi-static procedure.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

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.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Construction of habitat suitability models (HSMs) for benthic macroinvertebrate and their applications to instream environmental flows: A case study in Xiangxi River of Three Gorges Reservior region, China

Based on a long-term ecological monitoring, the present study chose the most dominant benthic macroinvertebrate (Baetis spp.) as target organisms in Xiangxi River, built the habitat suitability models (HSMs) for water depth, current velocity and substrate, respectively, which is the first aquatic organisms model for habitat suitability in the Chinese Mainland with a long-term consecutive in situ measurement. In order to protect the biointegrity and function of the river ecosystem, the theory system of instream environmental flow should be categorized into three hierarchies, namely minimum required instream flow (hydrological level), minimum instream environmental flow (biospecies level), and optimum instream environmental flow (ecosystem level). These three hierarchies of instream environmental flow models were then constructed with the hydrological and weighted usable area (WUA) method. The results show that the minimum required instream flow of Xiangxi River calculated by the Tennant method (10% of the mean annual flow) was 0.615 m(3) s(-1); the minimum instream environmental flow accounted for 19.22% of the mean annual flow (namely 1.182 m(3) s(-1)), which was the damaged river channel. ow in the dry season; and 42.91% of the mean annual flow (namely 2.639 m(3) s(-1)) should be viewed as the optimum instream environmental flow in order to protect the health of the river ecosystem, maintain the instream biodiversity, and reduce the impact of small hydropower stations nearby the Xiangxi River. We recommend that the hydrological and biological methods can help establish better instream environmental. ow models and design best management practices for use in the small hydropower station project. (C) 2008 National Natural Science Foundation of China and Chinese Academy of Sciences. Published by Elsevier Limited and Science in China Press. All rights reserved.

Development of EST-SSRs by an efficient FIASCO-based strategy: A case study in rare minnow (Gobiocyrpis rarus)

Type I markers are useful for comparative mapping and other genetic analyses, but relatively difficult to develop. In the present study a microsatellite (SSR)-enriched cDNA library was constructed for the first time using the fast isolation by AFLP of sequences containing repeats (FIASCO) method in a small fish, Chinese rare minnow (Gobiocypris rarus). A total of 97.4% of the expressed sequence tags (ESTs) contained targeted CA-repeats, in which 29 unique EST-SSRs were identified. Ten out of the 28 loci for which primer pairs were designed were polymorphic with alleles ranging from three to seven (mean 4.50). Some of these EST-SSRs can be amplified in other species. These results proved that cDNA-FIASCO is an efficient way to isolate novel EST-SSRs in a fish.

A method for analyzing the coupling of 3-D tapered rib waveguides with fibers

The fields in 3-dimensional tapered waveguides are unstable compared with the fields in the straight waveguides. In the case of waveguide-to-fiber coupling and fiber-to-waveguide coupling, a sequence of short straight waveguides has been modeled to approximate the 3-dimensional tapered waveguide; and the unstable incident and reflected fields, as well as their derivatives, were determined by the beam propagation method(BPM). Then free space radiation mode(FSRM) was employed to calculate the reflected and transmitted powers. Analysis results of the coupling of fiber with silicon-on-insulator(SOI) tapered rib waveguides showed the feasibility of the method.

encoding-based method for pattern matching

A new algorithm for compiling pattern matching is presented. Different from the traditional traversal-based approaches, it can represent a sequence of patterns as an integer by an encoding method and translate equations into case-expressions. The algorithm is simple to implement, and efficient for a kind of patterns, i.e. simple and dense patterns. This method can be complementary to traditional approaches.

automated test case generation for the stress testing of multimedia systems

With the advancement in network bandwidth and computing power, multimedia systems have become a popular means for information delivery. However, general principles of system testing cannot be directly applied to testing of multimedia systems on account of their stringent temporal and synchronization requirements. In particular, few studies have been made on the stress testing of multimedia systems with respect to their temporal requirements under resource saturation. Stress testing is important because erroneous behavior is most likely to occur under resource saturation. This paper presents an automatable method of test case generation for the stress testing of multimedia systems. It adapts constraint solving techniques to generate test cases that lead to potential resource saturation in a multimedia system. Coverage of the test cases is defined upon the reachability graph of a multimedia system. The proposed stress testing technique is supported by tools and has been successfully applied to a real-life commercial multimedia system. Although our technique focuses on the stress testing of multimedia systems, the underlying issues and concepts are applicable to other types of real-time systems.