A DNS study of jet control with microjets using an immersed boundary method

In this work, a microjet arrangement to control a turbulent jet is studied by means of direct numerical simulation. A customised numerical strategy was developed to investigate the interactions between the microjets and the turbulent jet. This approach is based on an improved immersed boundary method in order to reproduce realistically the control device while being compatible with the accuracy and the parallel strategy of the in-house code Incompact3d. The 16 converging microjets, so-called fluidevrons, lead to an increase of the turbulent kinetic energy in the near-nozzle region through an excitation at small scale caused by the interaction between the fluidevrons and the main jet. As a consequence, very intense unstable ejections are produced from the centre of the jet toward its surrounding. Further downstream, the turbulent kinetic energy levels are lower with a lengthening of the potential core compared to a natural jet, in agreement with experimental results.     

Thermal modulation and breakup of liquid jets

In this paper, we study the breakup behavior of Newtonian liquid and non‐Newtonian liquid jets with an arbitrary variation surface tension imposed along its length. The effect of duty cycle, fluid properties, and the various profiles of the surface tension is investigated. It is shown that the breakup behavior of a jet can be constructed by using the Fourier expansion of the surface tension profile. When the dimensionless wavenumber k is larger than 0.5, the jet breakup behavior is determined by the lowest frequency of the Fourier series expansion of the surface tension profile. As k decreases, higher frequency Fourier modes come to play. In general, for k between, 1∕(n+ 1) and 1∕n,n Fourier modes are needed to determine the jet breakup behavior. The current nonlinear model differs from the existing linear slender jet model in the literature in several ways. While the principle of superposition is valid for the linear model, it is not generally valid for the current nonlinear model. For the linear model, the jet will never break up when the wavenumber is larger than 1. The current model, however, shows clearly that the jet can indeed break up when the wavenumber is larger than 1. Furthermore, the current nonlinear model predicts a breakup time substantially higher than that from the linear model.Copyright © 2011 John Wiley & Sons, Ltd.     

Direct numerical simulation of thermocapillary flow based on the Volume of Fluid method

A numerical method for direct simulation of thermal Marangoni effects at dynamically deformable interface of two-phase incompressible fluids is developed. The approach is based on the Volume of Fluid (VOF) method with special focus on the numerical treatment of the temperature surface gradient because of its decisive role as the driving force of the flow. The surface gradient calculation is split into computing its length and direction in order to satisfy the correct thermal boundary condition at the interface without losing mobility of the interface. The method is applied to three different types of thermocapillary flow, namely thermocapillary migration of a droplet in an ambient fluid with linear temperature gradient, thermocapillary convection in a liquid layer under linear temperature gradient along the interface, and Marangoni convection due to Bénard–Marangoni instability. In the first case, different aspects of the dynamics of the migration are considered for validation such as the terminal migration velocity, the initial acceleration and quantification of the wall effects. The simulation also considers high convective heat transfer and covers a wide range of Marangoni numbers up to 5000, where good agreement with both theoretical and experimental results is achieved. In the second case, the convection velocity in the liquid layer is compared with an analytical result. In the final application, pattern formation due to the Bénard–Marangoni instability in a liquid layer in square geometry of small aspect ratio is investigated for realistic Biot number and dynamically deformable fluid interface. The results show good agreement with experiments from literature, where our numerical simulation also predicts cell pattern for a particular aspect ratio which is outside the limitation of the above cited experimental work.     

Application of the TVD scheme to the nonlinear instability analysis of a capillary jet

In the past, when either the perturbation‐type method or direct‐simulation approach was used to analyse capillary jets, the governing equations, which are parabolic in time and elliptic in space, were simplified or linearized. In the present study, the convective derivative term and a full, nonlinear form of the capillary pressure term are retained in the governing equations to investigate nonlinear effects on the break‐up of capillary jets. In this work, the TVD (i.e. total variation diminishing) scheme with flux‐vector splitting is applied to obtain the solutions of the system of nonlinear equations in a matrix form. Numerical results show that the present nonlinear model predicts longer jet break‐up lengths and slower growth rates for capillary jets than the previous linear model does. Comparing with other measurements from past literatures, the nonlinear results are consistent with the experimental data and appear more accurate than the linear analysis. In the past, the classic perturbation‐type analyses assumed constant growth rates for the fundamental and all harmonic components. By contrast, the present model is able to capture the local features of growth rates, which are not spatially and temporally constant. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical analysis of moving interfaces using a level set method coupled with adaptive mesh refinement

A novel numerical scheme is developed by coupling the level set method with the adaptive mesh refinement in order to analyse moving interfaces economically and accurately. The finite element method (FEM) is used to discretize the governing equations with the generalized simplified marker and cell (GSMAC) scheme, and the cubic interpolated pseudo‐particle (CIP) method is applied to the reinitialization of the level set function. The present adaptive mesh refinement is implemented in the quadrangular grid systems and easily embedded in the FEM‐based algorithm. For the judgement on renewal of mesh, the level set function is adopted as an indicator, and the threshold is set at the boundary of the smoothing band. With this criterion, the variation of physical properties and the jump quantity on the free surface can be calculated accurately enough, while the computation cost is largely reduced as a whole. In order to prove the validity of the present scheme, two‐dimensional numerical simulation is carried out in collapse of a water column, oscillation and movement of a drop under zero gravity. As a result, its effectiveness and usefulness are clearly shown qualitatively and quantitatively. Among them, the movement of a drop due to the Marangoni effect is first simulated efficiently with the present scheme. Copyright © 2004 John Wiley & Sons, Ltd.     

多跨连续长索在支座处存在滑移的非线性静力分析

提出了一种考虑悬索结构在中支座存在滑移时的非线性计算方法，推导出关于索元无应力长度发生变化时的滑移刚度公式，并进行了实例计算。计算结果表明，本文的方法正确，计算精度满足实际工程要求，并能使计算工作且大大减少，可供悬索结构设计、施工时参考。     

Numerical modelling and validation of Marangoni and surface tension phenomena using the finite volume method

Surface tension induced flow is implemented into a numerical modelling framework and validated for a number of test cases. Finite volume unstructured mesh techniques are used to discretize the mass, momentum and energy conservation equations in three dimensions. An explicit approach is used to include the effect of surface tension forces on the flow profile and final shape of a liquid domain. Validation of this approach is made against both analytical and experimental data. Finally, the method is used to model the wetting balance test for solder alloy material, where model predictions are used to gain a greater insight into this process. Copyright © 2000 John Wiley & Sons, Ltd.     

Formulations of the three-dimensional discontinuous deformation analysis method

This paper extends the original 2D discontinuous deformation analysis (DDA) method proposed by Shi to 3D cases, and presents the formulations of the 3D DDA. The formulations maintain the characteristics of the original 2D DDA approach. Contacts between the blocks are detected by using Common-Plane (C-P) approach and the non-smooth contact, such as of vertex-to-vertex, vertex-to-edge and edge-to-edge types, can be handled easily based on the C-P method. The matrices of equilibrium equations have been given in detail for programming purposes. TheC program codes for the 3D DDA are developed. The ability and accuracy of the formulations and the program are verified by the analytical solutions of several dynamic examples. The robustness and versatility of the algorithms presented in this paper are demonstrated with the aid of an example of scattering of densely packed cubes. Finally, implications and future extensions are discussed. The project supported by the National Natural Science Foundation of China (50139010)     

A finite amplitude analysis of tunnel deformation by the finite element method with highly viscous fluid

A finite element method for highly viscous fluid is used to calculate the velocity and stress fields in the surrounding soft rock of a tunnel. In order to fit the calculated values with the measured displacement of tunnel wall, we inverted the boundary forces and the mechanical parameters of the surrounding rocks.     

在中支座可滑动的两跨连续索的静力分析

基于非线性弹性理论,从单索的平衡出发,在索元的无应力长度已知情况下,提出了两跨连续索在中间节点滑移的初始平衡位置确定及均布荷载作用下的静力响应分析方法,建立了其非缄性计算的基本理论、公式和求解方法,并进行了实例计算。     

Nonlinear coupling modeling and dynamics analysis of rotating flexible beams with stretching deformation effect

Dynamic coupling modeling and analysis of rotating beams based on the nonlinear Green-Lagrangian strain are introduced in this work. With the reservation of the axial nonlinear strain, there are more coupling terms for axial and transverse deformations. The discretized dynamic governing equations are obtained by using the finite element method and Lagrange’s equations of the second kind. Time responses are conducted to compare the proposed model with other previous models. The stretching deforma...     

基于梁平截面假设的复杂细长结构动力模型简化方法

针对纵向尺度显著大于横向的复杂细长结构,提出了一个基于梁平截面假设的模型简化方法,建立了具有物理意义的基向量。同时考虑这类结构截面变形的翘曲模式,增加了翘曲基向量,得到精度较高的简化动力模型。这样的简化模型可以用于原结构在不同边界条件的动力分析。结构频率计算的两个具体算例表明了简化方法的有效性。     

Nonlocal continuum model for large deformation analysis of SLGSs using the kp-Ritz element-free method

A geometrically nonlinear large deformation analysis of SLGSs is presented using the element-free kp-Ritz method. Classical plate theory (CLP) is applied to describe the geometrically nonlinear behavior of SLGSs. Nonlocal elasticity theory is incorporated into CLP to take the small-scale effect into consideration. The system nonlinear equations are derived from the Ritz procedure based on the total Lagrangian formulation. The modified Newton–Raphson method and arc-length continuation are employed to solve the nonlinear equations. The efficiency of the element-free kp-Ritz method is verified through comparison with results reported in previous research. Numerical cases are studied to examine the influence of boundary conditions, aspect ratio, side length and nonlocal parameters on the nonlinear large deformation behavior of SLGSs. An interesting phenomenon is observed in that the nonlocal parameter effect is related to the mathematical expression of the transverse load.     

应用有限差分法的轴向流作用下叠层板的稳定性和模态分析

分析了轴向流作用下两端简支和固支叠层板的稳定性。基于势流理论建立轴向流作用下叠层板的流固耦合系统连续型运动方程,基于有限差分法建立了流场网格和结构网格统一的离散化格式,流场势函数用板的横向振动位移变量来表示,得到关于叠层板的横向振动位移变量的控制方程。求解控制方程的广义特征值,计算分析结果表明,两端简支和两端固支模型发生屈曲失稳,且得到了屈曲失稳临界速度与叠层板的层数和无量纲板间距的关系。此外,轴向流作用下叠层板的一阶模态并不是叠层板的同相弯曲模态。     

Vibration analysis of cracked plates using the extended finite element method

In the present paper, the extended finite element method (X-FEM) is adopted to analyze vibrations of cracked plates. Mindlin’s plate theory taking into account the effects of shear deformation and rotatory inertia is included in the development of the model. First, conventional FEM without any discontinuity is carried out, then the enrichment proposed by Moës et al. (Int J Numer Methods Eng 46, 131–150, 1999) of nodal elements containing cracks is added to the FEM formulation. Numerical implementation of enriched elements by discontinuous functions is performed, and thus dynamic equations (stiffness and mass matrices) are established. A FORTRAN computer code based on the X-FEM formulation is hence developed. Rectangular and square plates containing through-edge and central cracks with different boundary conditions are considered. The subspace iteration method is used to solve the eigenvalue problem. Natural frequencies as well as the corresponding eigenfunctions are consequently calculated as a function of the crack length. The obtained results show that the X-FEM is an efficient method in the dynamic analysis of plates containing discontinuities.     

Cracking analysis of fracture mechanics by the finite element method of lines （FEMOL）

The Finite Element Method of Lines (FEMOL) is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor (SIF) and crack opening displacement (COD) are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems. The English text was polished by Yunming Chen.     

Sensitivity analysis of low Reynolds number channel flow using the finite volume method

An unsteady finite volume‐based fractional step algorithm solved on a staggered grid has been developed for computing design sensitivity parameters in two‐dimensional flows. Verification of the numerical code is performed for the case of low Reynolds number, pressure‐driven flow through a straight channel, which has an exact steady‐state solution to the Navier–Stokes equations. Sensitivity of the flow to the channel height, fluid viscosity, and imposed pressure gradient is considered. Three different numerical techniques for computing the design sensitivity parameters: finite difference, complex‐step differentiation, and sensitivity equation method (SEM), are compared in terms of numerical error (relative to the exact solution), computational expense, and ease of implementation. Results indicate that, of all the three methods, complex step is the most accurate and requires the least computational time. In addition, treatment of the boundary conditions in SEM is addressed, within the framework of the present finite volume approach, with special attention given to parameter dependence in the boundary conditions. Error estimation based on the Grid Convergence Index provides a good indication of the exact error in the SEM solutions. An example of application of the use of sensitivity parameters to estimate the propagation of input uncertainty through the numerical simulation is also provided. Copyright © 2007 John Wiley & Sons, Ltd.