A high-efficiency second-order numerical scheme for time-fractional phase field models by using extended SAV method

Nonlinear Dynamics - In this paper, a second-order numerical scheme for the time-fractional phase field models is proposed. In this scheme, the fractional backward difference formula is used to...     

A numerical implementation using mid-node collocation for the hypersingular boundary contour method

A variant of the boundary element method, called the boundary contour method (BCM), offers a further reduction in dimensionality. Consequently, boundary contour analysis of two-dimensional problems does not require any numerical integration at all. In another development, a boundary contour implementation of a regularized hypersingular boundary integral equation (HBIE) using quadratic elements and end-node collocation was proposed and the technique is termed the hypersingular boundary contour method (HBCM). As reported in that work, the approach requires highly refined meshes in order to numerically enforce the stress continuity across boundary contour elements. This continuity requirement is very crucial since the regularized HBIE is only valid at collocation points where the stress tensor is continuous, while the computed stress at the endpoints of a boundary contour element, which is a non-conforming element, is generally not. This paper presents a new implementation of the HBCM for which the regularized HBIE is collocated at the mid-node of a boundary contour element. As the computed stress tensor is continuous at these mid-nodes, there is no need for unusually refined meshes. Some numerical tests herein show that, for the same mesh density, the HBCM using mid-node collocation has a comparable accuracy as the BCM.     

The numerical method for analysing the transient temperature field and transient phase transformation of metal materials in heat treating

In this paper,the Kirchhoffs transformation is popularized to the nonlinear heat conduction problem which the heat conductivity can be expressd as a multinomial of temperature firstly,the boundary condition of heat conduction problem is determined by analytics.Secondly,the incubation peroid superposition and the linear combination law is employed to simulate the transient phasses transformation in the process of heat treatment of materials.That the begin time of phase transformation,the type of phase transformation and the amount of phase constitution is determined simply.Finally,the three-dimension Dual Reciprocity Boundary Element Method is usedto analysis the total process of various heat treatment of component,the results of numerical calculation of examples show that the method provided in this paper is effectivce.     

Novel numerical implementation of asymptotic homogenization method for periodic plate structures

The present paper develops and implements finite element formulation for the asymptotic homogenization theory for periodic composite plate and shell structures, earlier developed in  and , and thus adopts this analytical method for the analysis of periodic inhomogeneous plates and shells with more complicated periodic microstructures. It provides a benchmark test platform for evaluating various methods such as representative volume approaches to calculate effective properties. Furthermore, the new numerical implementation (Cheng et al., 2013) of asymptotic homogenization method of 2D and 3D materials with periodic microstructure is shown to be directly applicable to predict effective properties of periodic plates without any complicated mathematical derivation. The new numerical implementation is based on the rigorous mathematical foundation of the asymptotic homogenization method, and also simplicity similar to the representative volume method. It can be applied easily using commercial software as a black box. Different kinds of elements and modeling techniques available in commercial software can be used to discretize the unit cell. Several numerical examples are given to demonstrate the validity of the proposed methods.     

On the extension of the integral length-scale approximation model to complex geometries

A new model for the unresolved stresses in large-eddy simulations was recently proposed by Piomelli et al. [J Fluid Mech 2015; 766:499–527] and Rouhi et al., [Phys Rev Fluids 2016; 1(4):0444011], in which the length scale is not related to the grid size, but determined based on turbulence properties. This model, the Integral Length-Scale Approximation (ILSA), has a single parameter, s_τ, which represents the contribution of the unresolved scales to the momentum transport, and is assigned by the user. We test ILSA in complex geometries using a low-dissipation finite-element method, and propose a rational method to determine s_τ on the basis of a grid-convergence study. The interaction of the model with the numerical method and grid topology is studied first; then, two cases are considered: the subcritical flow around a sphere, and the flow over the Ahmed body, a simplified car model. In each case calculations are performed using three grids and varying s_τ. With a consistent combination of grid size and s_τ the statistical results are in very good agreement with DNS data and experimental measurements. The eddy viscosity is insensitive to sudden variation of the mesh size, and the model adjusts to the different dissipation and diffusion characteristics associated with different grid topologies and numerical techniques.     

A numerical method for the evaluation of Fourier integrals

This paper suggests the use of spline function interpolation in the evaluation of Fourier integrals. At the same time, the numerical results of some common functions by various interpolation methods and a simplified method of construction of spline function for various boundary conditions are also presented.     

A phase field model of electromechanical fracture

Structural reliability analyses of piezoelectric solids need the modeling of failure under coupled electromechanical actions. However, the numerical simulation of failure due to fracture based on sharp crack discontinuities may suffer in situations with complex crack topologies. This can be overcome by a diffusive crack modeling based on the introduction of a crack phase field. In this work, we develop a framework of diffusive fracture in piezoelectric solids. We start our investigation with the definition of a crack surface functional of the phase field that Γ-converges for vanishing length-scale parameter to a sharp crack topology. This functional provides the basis for the definition of suitable dissipation functions which govern the evolution of the crack phase field. Based on experimental results available in the literature, we suggest a non-associative dissipative framework where the fracture phase field is driven by the mechanical part of the coupled electromechanical driving force. This accounts for a hierarchical view that considers (i) the decrease of stiffness due to mechanical rupture as the primary action that is followed by (ii) the decrease of electric permittivity due to the generated free space. The proposed definition of mechanical and electrical parts of the fracture driving force follows in a natural format from a kinematic assumption, that decomposes the total strains and the total electric field into energy-enthalpy-producing parts and fracture parts, respectively. Such an approach allows the insertion of well-known anisotropic piezoelectric storage functions without change. We end up with a three-field-problem that couples the displacement with the electric potential and the fracture phase field. The latter is governed by a micro-balance equation, which appears in a very transparent form in terms of a history field containing a maximum fracture source obtained in the time history of the electromechanical process. This representation allows the construction of a very robust algorithmic treatment based on a operator split scheme, which successively updates in a typical time step the history field, the crack phase field and finally the two piezoelectric fields. The proposed model is considered to be the canonically simple scheme for the simulation of diffusive electromechanical crack propagation in solids. We demonstrate its modeling capacity by means of representative numerical examples.     

A numerical method for the calculation of the stokes constants

A numerical method is introduced for the computation of the Stokes constants of certain linear differential equations. The method is tested on Whittaker's equation whose Stokes constants are known analytically. The factors which affect the accuracy of the results are shown. Some numerical results for an equation similar to Whittaker's equation are presented.     

A crystal plasticity analysis of length-scale dependent internal stresses with image effects

In this work, we present a stress functions approach to include image effects in continuum crystal plasticity arising from the long-range elastic interactions (LRI) between the GND density and free surfaces. The resulting length-scale dependent internal stresses augment those produced by the GND density variation. The formulation is applied to the case of a long, thin specimen subjected to uniform curvature. The analysis shows that under nominally uniform GND density distribution, internal stresses arise from two sources: (1) GND–GND LRI arising from the finite spatial extent of the uniform GND density field and (2) the LRI between the GND density and free surfaces appearing as image fields. A comparison with experimental results suggests that the length-scale for internal stresses, described as a correlation length-scale, should increase with decreasing specimen thickness. This observation is rationalized by associating the internal length-scale with the average slip-plane spacing, which may increase with decreasing specimen size due to paucity of dislocation sources. Finally, we also discuss the length-scale dependent image stress in terms of the Peach–Koehler force density proposed by Gurtin (2002).     

Nonlocal damage model using the phase field method: Theory and applications

A new nonlocal, gradient based damage model is proposed for isotropic elastic damage using the phase field method in order to show the evolution of damage in brittle materials. The general framework of the phase field model (PFM) is discussed and the order parameter is related to the damage variable in continuum damage mechanics (CDM). The time dependent Ginzburg–Landau equation which is also termed the Allen–Cahn equation is used to describe the damage evolution process. Specific length scale which addresses the interface region in which the process of changing undamaged solid to fully damaged material (microcracks) occurs is defined in order to capture the effect of the damaged localization zone. A new implicit damage variable is proposed through the phase field theory. Details of the different aspects and regularization capabilities are illustrated by means of numerical examples and the validity and usefulness of the phase field modeling approach is demonstrated.     

Two-dimensional numerical calculation of the flow field about a scoop inlet by the piecewise linear method

The piecewise linear method (PLM) based on time operator splitting is used to solve the unsteady compressible Euler equations describing the two-dimensional flow around and through a straight wall inlet placed stationary in a rapidly rotating supersonic flow. The PLM scheme is formulated as a Lagrangian step followed by an Eulerian remap. The inhomogeneous terms in the Euler equations written in cylindrical coordinates are first removed by Sod's method and the resulting set of equations is further reduced to two sets of one-dimensional Lagrangian equations, using time operator splitting. The numerically generated flow fields are presented for different values of the back pressure imposed at the downstream exit of the inlet nozzle. An oblique shock wave is formed in front of the almost whole portion of the inlet entrance, the incoming streamlines being deflected towards the higher pressure side after passing through the oblique shock wave and then bending down to the lower pressure side. A reverse flow appears inside the inlet nozzle owing to the recovery pressure of the incoming streams being lower than the back pressure of the inlet nozzle.     

陡峭山坡风场数值模拟方法研究

地形起伏是造成山地风场复杂多变的主要原因,其风场特性与基于均匀粗糙平面的风场有很大的区别.为准确模拟山坡地形风场,以坡角45°的简化陡峭山坡为研究对象,采用CFD数值模拟方法进行了流场分析,通过与风洞试验和各国规范对比,详细分析了网格分辨率、湍流模型和坡顶局部光滑处理等因素对数值模拟风场精度的影响.结果 表明,采用基本网格尺度布置以及Realizable k-ε湍流模型,在坡顶位置的风压系数值与风洞试验存在一定偏差;在坡顶分离点处采用具有二阶连续性的曲线进行局部光滑,可使得数值模拟所得风速比和风压系数与试验结果更好地吻合,且光滑曲线的过渡段水平距离越短,模拟效果相对越好,坡顶位置的地形加速效应模拟结果与文献试验、中国以及澳大利亚规范具有更好的一致性.     

A decoupling numerical method for fluid flow

A first-order non-conforming numerical methodology, Separation method, for fluid flow problems with a 3-point exponential interpolation scheme has been developed. The flow problem is decoupled into multiple one-dimensional subproblems and assembled to form the solutions. A fully staggered grid and a conservational domain centred at the node of interest make the decoupling scheme first-order-accurate. The discretization of each one-dimensional subproblem is based on a 3-point interpolation function and a conservational domain centred at the node of interest. The proposed scheme gives a guaranteed first-order accuracy. It is shown that the traditional upwind (or exponentially weighted upstream) scheme is less than first-order-accurate. The pressure is decoupled from the velocity field using the pressure correction method of SIMPLE. Thomas algorithm (tri-diagonal solver) is used to solve the algebraic equations iteratively. The numerical advantage of the proposed scheme is tested for laminar fluid flows in a torus and in a square-driven cavity. The convergence rates are compared with the traditional schemes for the square-driven cavity problem. Good behaviour of the proposed scheme is ascertained.     

3D two phase flows numerical simulations by SL‐VOF method

This paper describes the development of a semi‐Lagrangian computational method for simulating complex 3D two phase flows. The Navier–Stokes equations are solved separately in both fluids using a robust pseudo‐compressibility method able to deal with high density ratio. The interface tracking is achieved by the segment Lagrangian volume of fluid (SL‐VOF) method. The 2D SL‐VOF method using the concepts of VOF, piecewise linear interface calculation (PLIC) and Lagrangian advection of the interface is herein extended to 3D flows. Three different test cases of SL‐VOF 3D are presented for validation and comparison either with 2D flows or with other numerical methods. A good agreement is observed in each case. Copyright © 2004 John Wiley & Sons, Ltd.     

非网格重剖分模拟宏观裂纹体的扩展有限单元法(2:数值实现)

探讨了扩展有限单元法的具体实现过程,包括裂纹体几何结构的拓扑分析、广义节点的选取及详细的单元数值计算。并针对前文提出的扩展有限单元平衡方程的统一矩阵实现模式,提出了采用虚拟层合元的思想来处理被裂纹横贯单元的子域积分问题,自然地解决了原方法中由于特殊的位移插值场在裂纹两侧不连续造成的单元刚度阵求解困难。同时依托比较成熟的虚拟层合单元法,可以方便地考虑域内及裂纹面上分布载荷影响。此外,一、二维算例较高精度的数值结果验证了本文算法的有效性和精度。