Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.     

Regular wave integral approach to numerical simulation of radiation and diffraction of surface waves

A regular wave integral method is developed in the discretisation of a linear hydrodynamic problem on radiation and diffraction of surface waves by a floating or submerged body. The velocity potential of the problem is expressed as a solution of a body boundary integral equation involving the pulsating free surface Green function or pulsating free surface sources distributed on the body surface. With the use of a discretisation on the regular wave integral rather than discretisations on the singular wave integral of the Green function as in earlier investigations, the singular wave integral is approximated as an expansion of regular (or nonirregular) wave potentials. Influence coefficients between pulsating free surface source points are computed by the approximate expansion together with Hess–Smith panel integral formulas. Thus the velocity potential solution is evaluated by a boundary element algorithm. The numerical results produced from the proposed method agree well with semi-analytic solution results.     

Mode decomposition of nonlinear eigenvalue problems and application in flow stability

Direct numerical simulations are carried out with different disturbance forms introduced into the inlet of a flat plate boundary layer with the Mach number 4.5. According to the biorthogonal eigenfunction system of the linearized Navier-Stokes equations and the adjoint equations, the decomposition of the direct numerical simulation results into the discrete normal mode is easily realized. The decomposition coefficients can be solved by doing the inner product between the numerical results and the eigenfunctions of the adjoint equations. For the quadratic polynomial eigenvalue problem, the inner product operator is given in a simple form, and it is extended to an Nth-degree polynomial eigenvalue problem. The examples illustrate that the simplified mode decomposition is available to analyze direct numerical simulation results.     

Calculation of Transport‐Equation Eigenvalues

The eigenvalue problem for the transport equation with variable coefficients in an arbitrary domain with a smooth boundary is considered. A saturationfree numerical algorithm is constructed. Examples of numerical calculations are given, which prove the effectiveness of the proposed procedure.     

The asymptotics for the eigenelements of the Laplacian in a cylinder with frequently alternating boundary conditions

The singular perturbed eigenvalue problem for the Laplace operator in a cylindrical body with frequently alternating boundary conditions is considered. The asymptotics for the eigenelements are constructed.     

V形切口应力强度因子的一种边界元分析方法

将V形切口结构分成围绕切口尖端的小扇形和剩余结构两部分. 尖端处扇形域应力场表示成关于尖端距离$\rho$的渐近级数展开式,从线弹性理论方程推导出了一组分析平面V形切口奇异性的常微分方程特征值问题,通过求解特征方程,得到前若干个奇性指数和相应的特征向量. 再将切口尖端的位移和应力表示为有限个奇性阶和特征向量的组合. 然后用边界元法分析挖去小扇形后的剩余结构. 将位移和应力的线性组合与边界积分方程联立,求解获得切口根部区域的应力场、应力幅值系数和整体结构的位移和应力. 从而准确计算出平面V形切口的奇异应力场和应力强度因子.      

Sliding friction across the scales: Thermomechanical interactions and dissipation partitioning

A homogenization framework is developed for determining the complete macroscopic thermomechanical sliding contact response of soft interfaces with microscopic roughness. To this end, a micro–macro mechanical dissipation equality is first established which enables defining a macroscopic frictional traction. The derivation allows both contacting bodies to be deformable, thereby extending the commonly adopted setting where one of the bodies is rigid. Moreover, it forms a basis for the second step, where a novel micro–macro thermal dissipation equality is established which enables defining partitioning coefficients that are associated with the frictional dissipation as it is perceived on the macroscale. Finally, a comparison of the temperature fields from the original heterogeneous thermomechanical contact problem and an idealized homogeneous one reveals an identification of the macroscopic temperature jump. The computational implementation of the framework is carried out within an incrementally two-phase micromechanical test which delivers a well-defined macroscopic response that is not influenced by purely algorithmic choices such as the duration of sliding. Two-dimensional numerical investigations on periodic and random samples from thermo-viscoelastic boundary layers with unilateral and bilateral roughness demonstrate the temperature–velocity–pressure dependence of the macroscopic contact response.     

Parametric solutions of the eigenvalue equation for a convectively cooled slab

The convective cooling of a slab by an ambient fluid under the most general linear and homogeneous boundary conditions is considered. For the roots of the corresponding transcendental eigenvalue equation an explicit formula is written down in a parametric form. The practical consequences of this representation, among them certain “singular solutions” which cannot be obtained by a direct numerical treatment of the original eigenvalue equation, are discussed in detail.     

纤维端部的界面裂纹分析

基于弹性力学空间轴对称问题的通解,研究了短纤维增强复合材料中纤维端部的轴对称币形和柱形界面裂纹尖端的应力奇异性,得到了裂纹尖端附近的奇异应力场．研究结果表明,这两种轴对称界面裂纹尖端的应力奇异性相同,并且与平面应变状态下相应模型的应力奇异性一致,材料性能对裂纹尖端附近奇异应力场的影响可用三个组合参数描述     

Efficient computation of the eigenspectrum of viscoelastic flows using submatrix-based transformation of the linearized equations

Linear stability analysis of incompressible viscoelastic flows based on normal mode expansions of the eigenfunctions requires the numerical solution of a generalized eigenvalue problem (GEVP). The complex boundary layer structure of the leading eigenfunctions and the singular character of the continuous set of eigenvalues, necessitate the use of fine mesh sizes, leading to large algebraic GEVPs. In this paper, we present a submatrix-based transformation of the linearized equations (SubTLE) that converts the GEVP into a simple eigenvalue problem (EVP) of half the original dimension for the purely elastic isothermal and non-isothermal flows of an Oldroyd-B liquid. This leads to significant (up to an order of magnitude) reduction in the CPU time and memory required for the solution of the EVP. This is illustrated in the context of isothermal and non-isothermal shear flows.     

Karhunen–Loeve representations of turbulent channel flows using the method of snapshots

For three‐dimensional flows with one inhomogeneous spatial coordinate and two periodic directions, the Karhunen–Loeve procedure is typically formulated as a spatial eigenvalue problem. This is normally referred to as the direct method (DM). Here we derive an equivalent formulation in which the eigenvalue problem is formulated in the temporal coordinate. It is shown that this so‐called method of snapshots (MOS) has some numerical advantages when compared to the DM. In particular, the MOS can be formulated purely as a matrix composed of scalars, thus avoiding the need to construct a matrix of matrices as in the DM. In addition, the MOS avoids the need for so‐called weight functions, which emerge in the DM as a result of the non‐uniform grid typically employed in the inhomogeneous direction. The avoidance of such weight functions, which may exhibit singular behaviour, guarantees satisfaction of the boundary conditions. The MOS is applied to data sets recently obtained from the direct simulation of turbulence in a channel in which viscoelasticity is imparted to the fluid using a Giesekus model. The analysis reveals a steep drop in the dimensionality of the turbulence as viscoelasticity is increased. This is consistent with the results that have been obtained with other viscoelastic models, thus revealing an essential generic feature of polymer‐induced drag reduced turbulent flows. Published in 2006 by John Wiley & Sons, Ltd.     

比例边界等几何分析方法Ⅰ:波导本征问题

提出比例边界等几何方法 (scaled boundary isogeometric analysis, SBIGA), 并用以求解波导本征值问题. 在比例边界等几何坐标变换的基础上, 利用加权余量法将控制偏微分方程进行离散处理, 半弱化为关于边界控制点变量的二阶常微分方程, 即 TE 波或 TM 波波导的比例边界等几何分析的频域方程以及波导动刚度方程, 同时利用连分式求解波导动刚度矩阵. 通过引入辅助变量进一步得出波导本征方程. 该方法只需在求解域的边界上进行等几何离散, 使问题降低一维, 计算工作量大为节约, 并且由于边界的等几何离散, 使得解的精度更高, 进一步节省求解自由度. 以矩形和 L 形波导的本征问题分析为例, 通过与解析解和其他数值方法比较, 结果表明该方法具有精度高、计算工作量小的优点.     

Scattering of SHPlease define abbreviation SH waves by arc-shaped interface cracks between a cylindrical magneto-electro-elastic inclusion and matrix: near fields

Summary In this paper, the scattering of SH waves by a magneto-electro-elastic cylindrical inclusion partially debonded from its surrounding magneto-electro-elastic material is investigated by using the wavefunction expansion method and a singular integral equation technique. The debonding regions are modeled as multiple arc-shaped interface cracks with non-contacting faces. The magneto-electric impermeable boundary conditions are adopted. By expressing the scattered fields as wavefunction expansions with unknown coefficients, the mixed boundary-value problem is firstly reduced to a set of simultaneous dual-series equations. Then, dislocation density functions are introduced as unknowns to transform these dual-series equations to Cauchy singular integral equations of the first type,which can be numerically solved easily. The solution is valid for arbitrary number and size of the arc-shaped interface cracks. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond. The effects of incident direction, crack configuration and various material parameters on the dynamic stress intensity factors are discussed. The solution of this problem is expected to have applications in the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.The work was supported by the National Natural Science Fund of China (Project No. 19772029) and the Research Fund for Doctors of Hebei Province, China (Project No. B2001213).     

基于新型裂尖杂交元的压电材料断裂力学研究

提出了一种裂尖邻域杂交元模型，将其与标准杂交应力元结合来求解压电材料裂纹尖 端的奇性电弹场和断裂参数的数值解．裂纹尖端杂交元的建立步骤为：1) 利用高次内插有限元特征法求解特征问题，得到反映裂尖奇异性电弹场状况的特 征值和特征角分布函数；2) 利用广义Hellinger-Reissner变分泛函以及特征问题的解来建立裂尖邻域杂交元模型．该 方法求解电弹场时，摒弃了传统有限元方法中裂尖奇异性场需要借助解析解的做法，也避免 了单纯有限元方法中需要在裂尖端部进行高密度单元划分．采用PZT5板中心裂纹问题 作为考核例，数值结果显示了良好的精确性．作为进一步应用，求解了含中心界面裂纹 的PZT4-PZT5两相压电材料的应力强度因子和电位移强度因子．所有的算例都考虑 了3种裂纹面电边界条件．     

Nonlinear forced vibration analysis of postbuckled beams

A numerical solution methodology is proposed herein to investigate the nonlinear forced vibrations of Euler–Bernoulli beams with different boundary conditions around the buckled configurations. By introducing a set of differential and integral matrix operators, the nonlinear integro-differential equation that governs the buckling of beams is discretized and then solved using the pseudo-arc-length method. The discretized governing equation of free vibration around the buckled configurations is also solved as an eigenvalue problem after imposing the boundary conditions and some complicated matrix manipulations. To study forced and nonlinear vibrations that take place around a buckled configuration, a Galerkin-based numerical method is applied to reduce the partial integro-differential equation into a time-varying ordinary differential equation of Duffing type. The Duffing equation is then discretized using time differential matrix operators, which are defined based on the derivatives of a periodic base function. Finally, for any given magnitude of axial load, the pseudo -arc-length method is used to obtain the nonlinear frequencies of buckled beams. The effects of axial load on the free vibration, nonlinear, and forced vibrations of beams in both prebuckling and postbuckling domains for the lowest three vibration modes are analyzed. This study shows that the nonlinear response of beams subjected to periodic excitation is complex in the postbuckling domain. For example, the type of boundary conditions significantly affects the nonlinear response of the postbuckled beams.     

The effect of an additional boundary condition on the plane creeping flow of a second-order fluid

A similarity solution is obtained for the two-dimensional creeping flow of a second-order fluid with non-parallel porous walls. The resulting ordinary differential equation is fifth order. Thus, an additional velocity boundary condition is needed, the other four being due to the usual no-slip conditions. Having chosen to prescribe the rate of shear at the wall, the problem is solved by a standard numerical routine. A singular perturbation analysis is developed for small values of the Deborah number. A type of boundary layer forms for which the viscous Newtonian case is the outer solution.     

曲线裂纹和反平面圆形夹杂相交问题

建立了和反平面圆夹杂界面相交的曲线裂纹的弱奇异积分方程,利用Cauchy型奇异积分方程主部分析方法研究了穿过反平面圆夹杂界面的曲线裂纹在交点处的奇性应力指数以及交点处角形域内的奇性应力,并根据奇性应力定义了交点处的应力强度因子。通过对弱奇异积分方程的数值求解,可得裂纹端点和交点处的应力强度因子。     

Elastic stability of all edges clamped stepped and stiffened rectangular plate under uni-axial, bi-axial and shearing forces

The stability of clamped stepped and stiffened rectangular plate subjected to in-plane forces is examined. The plate is divided into 900 rectangular meshes and the partial derivatives are approximated using central difference formula. Altogether 841 equations of equilibrium and 248 equations representing boundary conditions are formed, finally leading to the solution of eigenvalue problem. The buckling coefficients are calculated for various types of stepped plates and the results are presented in tables for ready use by designers. The results are compared with the published results and they are in close agreement.     

Transient 1D transport equation simulated by a mixed Green element formulation

New discrete element equations or coefficients are derived for the transient 1D diffusion–advection or transport equation based on the Green element replication of the differential equation using linear elements. The Green element method (GEM), which solves the singular boundary integral theory (a Fredholm integral equation of the second kind) on a typical element, gives rise to a banded global coefficient matrix which is amenable to efficient matrix solvers. It is herein derived for the transient 1D transport equation with uniform and non-uniform ambient flow conditions and in which first-order decay of the containment is allowed to take place. Because the GEM implements the singular boundary integral theory within each element at a time, the integrations are carried out in exact fashion, thereby making the application of the boundary integral theory more utilitarian. This system of discrete equations, presented herein for the first time, using linear interpolating functions in the spatial dimensions shows promising stable characteristics for advection-dominant transport. Three numerical examples are used to demonstrate the capabilities of the method. The second-order-correct Crank–Nicolson scheme and the modified fully implicit scheme with a difference weighting value of two give superior solutions in all simulated examples. © 1997 John Wiley & Sons, Ltd.     

Scattering of Tollmien-Schlichting waves by localized roughness in transonic boundary layers

The laminar-turbulent transition in boundary-layer flows is often affected by wall imperfections, because the latter may interact with either the freestream perturbations or the oncoming boundary-layer instability modes, leading to a modification of the accumulation of the normal modes. The present paper particularly focuses on the latter mechanism in a transonic boundary layer, namely, the effect of a two-dimensional(2 D) roughness element on the oncoming Tollmien-Schlichting(T-S) modes when they propagate through the region of the rapid mean-flow distortion induced by the roughness. The wave scattering is analyzed by adapting the local scattering theory developed for subsonic boundary layers(WU, X. S. and DONG, M. A local scattering theory for the effects of isolated roughness on boundary-layer instability and transition: transmission coefficient as an eigenvalue. Journal of Fluid Mechanics, 794, 68–108(2006)) to the transonic regime, and a transmission coefficient is introduced to characterize the effect of the roughness. In the sub-transonic regime, in which the Mach number is close to, but less than, 1, the scattering system reduces to an eigenvalue problem with the transmission coefficient being the eigenvalue; while in the super-transonic regime, in which the Mach number is slightly greater than 1, the scattering system becomes a high-dimensional group of linear equations with the transmission coefficient being solved afterward. In the largeReynolds-number asymptotic theory, the K′arm′an-Guderley parameter is introduced to quantify the effect of the Mach number. A systematical parametric study is carried out,and the dependence of the transmission coefficient on the roughness shape, the frequency of the oncoming mode, and the K′arm′an-Guderley parameter is provided.