An Efficient Front-Tracking Method For Fully Nonlinear Interfacial Waves

A fully nonlinear and dispersive model within the framework of potential theory is developed for interfacial (2-layer) waves. To circumvent the difficulties arisen from the moving boundary problem a viable technique based on the mixed Eulerian and Lagrangian concept is proposed: the computing area is partitioned by a moving mesh system which adjusts its location vertically to conform to the shape of the moving boundaries but keeps frozen in the horizontal direction. Accordingly, a modified dynamic condition is required to properly compute the boundary potentials. To demonstrate the effectiveness of the current method, two important problems for the interfacial wave dynamics, the generation and evolution processes, are investigated. Firstly, analytical solutions for the interfacial wave generations by the interaction between the barotropic tide and topography are derived and compared favorably with the numerical results. Furthermore simulations are performed for the nonlinear interfacial wave evolutions at various water depth ratios and satisfactory agreement is achieved with the existing asymptotical theories. (c) 2008 Elsevier Inc. All rights reserved.

Wave-front analysis method of circular aperture sampling for collimation testing

A new type of wave-front analysis method for the collimation testing of laser beams is proposed. A concept of wave-front height is defined, and, on this basis, the wave-front analysis method of circular aperture sampling is introduced. The wave-front height of the tested noncollimated wave can be estimated from the distance between two identical fiducial diffraction planes of the sampled wave, and then the divergence is determined. The design is detailed, and the experiment is demonstrated. The principle and experiment results of the method are presented. Owing to the simplicity of the method and its low cost, it is a promising method for checking the collimation of a laser beam with a large divergence. © 2005 Optical Society of America.

Thermocapillary Migration Of Nondeformable Drops

In this paper, we present a numerical study on the thermocapillary migration of drops. The Navier-Stokes equations coupled with the energy conservation equation are solved by the finite-difference front-tracking scheme. The axisymmetric model is adopted in Our simulations, and the drops are assumed to be perfectly spherical and nondeformable. The benchmark simulation starts from the classical initial condition with a uniform temperature gradient. The detailed discussions and physical explanations of migration phenomena are presented for the different values of (1) the Marangoni numbers and Reynolds numbers of continuous phases and drops and (2) the ratios of drop densities and specific heats to those of continuous phases. It is found that fairly large Marangoni numbers may lead to fluctuations in drop velocities at the beginning part of simulations. Finally, we also discuss the influence of initial conditions on the thermocapillary migrations. (C) 2008 American Institute of Physics.

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 simulation of drop Marangoni migration under microgravity

Thermocapillary motion of a drop in a uniform temperature gradient is investigated numerically. The three-dimensional incompressible Navier-Stokes and energy equations are solved by the finite-element method. The front tracking technique is employed to describe the drop interface. To simplify the calculation, the drop shape is assumed to be a sphere. It has been verified that the assumption is reasonable under the microgravity environment. Some calculations have been performed to deal with the thermocapillary motion for the drops of different sizes. It has been verified that the calculated results are in good agreement with available experimental and numerical results. (C) 2003 Elsevier Ltd. All rights reserved.

Highly sensitive wave-front sensor with a non-zero-order phase plate

We propose a novel highly sensitive wave front detection method for a quick check of a flat wave front by taking advantage of a non-zero-order pi phase plate that yields a non-zero-order diffraction pattern. When a light beam with a flat wave front illuminates a phase plate, the zero-order intensity is zero. When there is a slight distortion of the wave front, the zero-order intensity increases. The ratio of first-order intensity to that of zero-order intensity is used as the criterion with which to judge whether the wave front under test is flat, eliminating the influence of background light. Experimental results demonstrate that this method is efficient, robust, and cost-effective and should be highly interesting for a quick check of a flat wave front of a large-aperture laser beam and adaptive optical systems. (c) 2005 Optical Society of America.

An N-shape thermal front in the western South Yellow Sea in winter

An N-shape thermal front in the western South Yellow Sea (YS) in winter was detected using Advanced Very High Resolution Radiation (AVHRR) Sea Surface Temperature data and in-situ observations with a merged front-detecting method. The front, which exists from late October through early March, consists of western and eastern wings extending roughly along the northeast-southwest isobaths with a southeastward middle segment across the 20-50 m isobaths. There are north and south inflexions connecting the middle segment with the western and eastern wings, respectively. The middle segment gradually moves southwestward from November through February with its length increasing from 62 km to 107 km and the southern inflexion moving from 36.2A degrees N to 35.3A degrees N. A cold tongue is found to coexist with the N-shape front, and is carried by the coastal jet penetrating southward from the tip of the Shandong Peninsula into the western South YS as revealed by a numerical simulation. After departing from the coast, the jet flows as an anti-cyclonic recirculation below 10 m depth, trapping warmer water originally carried by the compensating Yellow Sea Warm Current (YSWC). A northwestward flowing branch of the YSWC is also found on the lowest level south of the front. The N-shape front initially forms between the cold tongue and warm water involved in the subsurface anti-cyclonical recirculation and extends upwards to the surface through vertical advection and mixing. Correlation analyses reveal that northerly and easterly winds tend to be favorable to the formation and extension of the N-shape front probably through strengthening of the coastal jet and shifting the YSWC pathway eastward, respectively.

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.

基于栅格地图的多化学羽流源定位

本文研究了移动机器人在湍流环境中定位多个化学羽流源的问题。利用粒子随机行走方法建立动态羽流模型。在此环境中机器人进行梳状搜索并采集羽流数据,使用一种基于后验概率独立假设的改进贝叶斯算法融合这些传感器数据建立一张羽流源位置的概率地图,地图中具有高概率值的栅格指出了羽流源可能的位置。仿真结果说明了该方法的有效性,通过与标准贝叶斯算法的比较说明了该方法优点。

On a novel method of impact by a front-end-coated bullet to evaluate the interface adhesion between film and substrate

A preliminary experiment was carried out to validate the feasibility of the method of impact by a front-end-coated bullet to evaluate the interface adhesion between film and substrate. The theoretical description of the initiation, propagation and evolution of the stress pulse during impact was generalized and formulized. The effects of the crucial parameters on the interface stress were further investigated with FEM. The results found the promising prospect of the application of such a method and provided useful guidance for experimental design.

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Stress state in front of a crack perpendicular to bimaterial interface

Using a dislocation simulation approach, the basic equation for a crack perpendicular to a bimaterial interface is formulated in this paper. A novel expansion method is proposed for solving the problem. The complete solution for the problem, including the T stress ahead of the crack tip and the stress intensity factors are presented. The stress field characteristics are analyzed in detail. It is found that ahead of the crack tip and near the interface the normal stress, perpendicular to the crack plane, sigma(x), is characterized by the K fields and the normal stress sigma(y) is dominated by the K field plus T stress in the region of 0 < r/b < 0.4 for b/a(0) less than or equal to 0.1, where b is the distance from the crack tip to the interface.

Three-dimensional crack problem analysis using boundary element method with finite-part integrals

A set of hypersingular integral equations of a three-dimensional finite elastic solid with an embedded planar crack subjected to arbitrary loads is derived. Then a new numerical method for these equations is proposed by using the boundary element method combined with the finite-part integral method. According to the analytical theory of the hypersingular integral equations of planar crack problems, the square root models of the displacement discontinuities in elements near the crack front are applied, and thus the stress intensity factors can be directly calculated from these. Finally, the stress intensity factor solutions to several typical planar crack problems in a finite body are evaluated.

Pipeline Architecture and Parallel Computation-Based Real-Time Stereovision Tracking System for Surgical Navigation

This paper studies the development of a real-time stereovision system to track multiple infrared markers attached to a surgical instrument. Multiple stages of pipeline in field-programmable gate array (FPGA) are developed to recognize the targets in both left and right image planes and to give each target a unique label. The pipeline architecture includes a smoothing filter, an adaptive threshold module, a connected component labeling operation, and a centroid extraction process. A parallel distortion correction method is proposed and implemented in a dual-core DSP. A suitable kinematic model is established for the moving targets, and a novel set of parallel and interactive computation mechanisms is proposed to position and track the targets, which are carried out by a cross-computation method in a dual-core DSP. The proposed tracking system can track the 3-D coordinate, velocity, and acceleration of four infrared markers with a delay of 9.18 ms. Furthermore, it is capable of tracking a maximum of 110 infrared markers without frame dropping at a frame rate of 60 f/s. The accuracy of the proposed system can reach the scale of 0.37 mm RMS along the x- and y-directions and 0.45 mm RMS along the depth direction (the depth is from 0.8 to 0.45 m). The performance of the proposed system can meet the requirements of applications such as surgical navigation, which needs high real time and accuracy capability.

Pipeline Architecture and Parallel Computation-Based Real-Time Stereovision Tracking System for Surgical Navigation

This paper studies the development of a real-time stereovision system to track multiple infrared markers attached to a surgical instrument. Multiple stages of pipeline in field-programmable gate array (FPGA) are developed to recognize the targets in both left and right image planes and to give each target a unique label. The pipeline architecture includes a smoothing filter, an adaptive threshold module, a connected component labeling operation, and a centroid extraction process. A parallel distortion correction method is proposed and implemented in a dual-core DSP. A suitable kinematic model is established for the moving targets, and a novel set of parallel and interactive computation mechanisms is proposed to position and track the targets, which are carried out by a cross-computation method in a dual-core DSP. The proposed tracking system can track the 3-D coordinate, velocity, and acceleration of four infrared markers with a delay of 9.18 ms. Furthermore, it is capable of tracking a maximum of 110 infrared markers without frame dropping at a frame rate of 60 f/s. The accuracy of the proposed system can reach the scale of 0.37 mm RMS along the x- and y-directions and 0.45 mm RMS along the depth direction (the depth is from 0.8 to 0.45 m). The performance of the proposed system can meet the requirements of applications such as surgical navigation, which needs high real time and accuracy capability.