Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

Experimental Study on Vortex-Induced Vibrations of a Pipeline Near a Plane Boundary in Steady Flow

对单向水流作用下近壁管道横向涡激振动进行了实验模拟,重点探讨了管道与壁面间隙比(e/D)对管道涡激振动幅值和涡激振动频率响应特性的影响规律.实验结果表明,管道与壁面间隙宽度对管道涡激振动特性有较明显影响.在较大间隙比(e/D＞0.66)下,管道振幅随着Vr数的增大先快速增长到最大值,然后平缓下降;在振动初期(即Vr数较小时),管道振动频率变化基本符合Strouhal规律;在振动中后期(即Vr数较大时),管道振动频率变化不符合Strouhal规律,而在管道固有频率附近缓慢增长.在较小间隙比(e/D＜0.30)下,管道振幅随Vr数的增大先平缓上升到最大值,随后较快速下降;在振动初期,管道振动频率变化不遵循Strouhal规律;在整个振动范围内,与较大间隙比情况相比,随着Vr数增加,管道振动频率增长幅度明显较大.

Influence of grain boundary on melting

The temperature behaviour of an Al bicrystal with surfaces consisting of (110) and (111) crystals is simulated using molecular dynamics. The result shows that the (110) crystal losses its crystalline order at 820K, whereas the disorder does not propagate through the (111) crystal at this temperature. Instead, some disordered atoms are recrystallized into the (111) crystal and the initial grain boundary changes into a stable order-disorder interface. Thus, it was discovered that at a temperature near its melting point, the (111) crystal grew and obstructed the propagation of disorder. Such an obstruction is helpful for understanding melting.

Flow Characteristics and Its Friction Coefficient of Laminar Near-Wall Boundary Layer in Dusty Atmosphere

采用双向耦合的双流体模型,研究了大气悬浮沙尘的存在对大气边界层中层流底层流动特性及摩阻系数的影响,计算并讨论了不同沙尘含量下含尘大气相对于无尘大气摩阻系数的变化。结果表明：摩阻系数的变化取决于悬浮沙尘的初始运动状态和质量载荷率。

Coherent Structures in Transition of a Flat-Plate Boundary Layer at Ma=0.7

Direct numerical simulation (DNS) of a spatially evolving flat-plate boundary layer transition process at free stream Mach number 0.7 is performed. Tollmien-Schlichting (T-S) waves are added on the inlet boundary as the disturbances before transition. Typical coherent structures in the transition process are investigated based on the second invariant of velocity gradient tensor. The instantaneous shear stress and the mean velocity profile in the transition region are studied. In our view, the fact that the peak value of shear stress in the stress concentration area increases and exceeds a threshold value during the later stage of the transition process plays an important role in the laminar breakdown process.

Experimental Study on Vortex-Induced Vibrations of Submarine Pipeline Near Seabed Boundary in Ocean Currents

Unlike most previous studies on vortex- induced vibrations of a cylinder far from a boundary, this paper focuses On the influences of close proximity of a submarine pipeline to a rigid seabed boundary upon the dynamic responses of the pipeline in ocean currents. The effects of gap-to-diameter ratio and those of the stability parameter on the amplitude and frequency responses of a pipeline are investigated experimentally with a novel hydro-elastic facility. A comparison is made between the present experimental results Of the amplitude and frequency responses for the pipes with seabed boundary effects and those for wall-free cylinders given by Govardhan and Williamson (2000) and Anand ( 1985). The comparison shows that the close proximity of a pipeline to seabed has much influence on the vortex- induced vibrations of the pipeline. Both the width of the lock-in ranges in terms of V, and the dimensionless amplitude ratio A(max)/D become larger with the decrease of the gap-to-diameter ratio e/D. Moreover, the vibration of the pipeline becomes easier to occur and its amplitude response becomes more intensive with the decrease of the stability parameter, while the pipeline frequency responses are affected slightly by the stability parameter.

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes, equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tolimien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.

DNS of Compressible Turbulent Boundary Layer Over a Blunt Wedge

Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suction is located downstream of the sound-speed point. Statistics are studied and compared with the results in incompressible flat-plate boundary layer. The mean pressure gradient effects on the vortex structure are studied.

Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0.similar to 4.7) x 10(7)/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30 degrees for a sweepback angle of 67.6 degrees. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10 degrees and secondary separation were detected at deflection of theta greater than or equal to 20 degrees. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements at M greater than or equal to 6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers.

Finite Boundary Effects in Problem of a Crack Perpendicular to and Terminating at a Bimaterial Interface

In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731-2755) are discussed also.

Flow-Induced Vibrations Of A Cylinder With Two Degrees Of Freedom Near Rigid Plane Boundary

Based on similarity analyses, the flow-induced vibrations of a near-wall cylinder with 2 degrees of freedom are investigated experimentally by employing a hydroelastic apparatus in conjunction with a flume. The cylinder's vibration amplitude, vibration frequency and vortex shedding frequency were measured and analyzed. The effects of gap-to-diameter ratio (e,ID) upon the vibration responses are further investigated. The experimental results indicate that, when the reduced velocity (Vr) is small (e.g. Vr = 1.2 similar to 2.6), only streamwise vibration occurs, and its frequency is quite close to its natural frequency in still water. When increasing Vr (e.g. Vr > 3.4), both streamwise and transverse vibrations of the near-wall cylinder may occur. In the examined range of gap-to-diameter ratio (0.42 < e(0)/D < 2.68), 2 vibration stages (in terms of Vr) of streamwise vibrations usually exist: First Streamwise Vibration (FSV) and Second Streamwise Vibration (SSV). In the SSV stage, the vortex shedding frequency may either undergo a jump to that of the streamwise vibration, or stay consistent with that of the transverse vibration. The amplitudes of transverse vibration are usually much larger than those of streamwise vibration for the same value of e(0)/D. The maximum amplitudes of both streamwise and transverse vibration get larger with the increase of e(0)/D (0.42 < e(0)/D < 2.68).

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Grain Boundary Effects On Plastic Deformation And Fracture Mechanisms In Cu Nanowires: Molecular Dynamics Simulations

Metallic nanowires have many attractive properties such as ultra-high yield strength and large tensile elongation. However, recent experiments show that metallic nanowires often contain grain boundaries, which are expected to significantly affect mechanical properties. By using molecular dynamics simulations, here, we demonstrate that polycrystalline Cu nanowires exhibit tensile deformation behavior distinctly different from their single-crystal counterparts. A significantly lowered yield strength was observed as a result of dislocation emission from grain boundaries rather than from free surfaces, despite of the very high surface to volume ratio. Necking starts from the grain boundary followed by fracture, resulting in reduced tensile ductility. The high stresses found in the grain boundary region clearly play a dominant role in controlling both inelastic deformation and fracture processes in nanoscale objects. These findings have implications for designing stronger and more ductile structures and devices on nanoscale.