Numerical analysis of the initial conditions in fractional systems

Fractional dynamics is a growing topic in theoretical and experimental scientific research. A classical problem is the initialization required by fractional operators. While the problem is clear from the mathematical point of view, it constitutes a challenge in applied sciences. This paper addresses the problem of initialization and its effect upon dynamical system simulation when adopting numerical approximations. The results are compatible with system dynamics and clarify the formulation of adequate values for the initial conditions in numerical simulations.     

Numerical study on combustion of H2—O2 in a supersonic reactive expansion in a nozzle

A study on expansion flow inside a nozzle considering full mechanism chemistry of hydrogen and oxygen was carried out. In this study, a full implicit scheme for turbulent reactive flow was obtained by combining the second order TVD scheme of Yee and Harten (1987, Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates. American Institute of Aeronautics and Astronautics Journal, 25(2), 266–274) with the efficient implicit lower-upper scheme of Shuen and Yoon (1989, Numerical study of chemically reacting flows using a lower-upper symmetric successive overrelaxation scheme. American Institute of Aeronautics and Astronautics Journal, 27(12), 1752–1756). The species equations, Navier–Stokes equations and turbulence model were implemented in the numerical scheme and solved in conjunction with full detailed finite rate chemistry. The numerical scheme is verified by comparison with experimental results of a converging–diverging nozzle. Effects of inlet pressure, temperature and fuel-oxidant mass ratio on nozzle flow field were studied. Variation of chemical species under different conditions was investigated by considering a chemical mechanism. Results show that increasing inlet pressure increases the rate of reactions due to increasing the concentration of reactants. For lower inlet pressure the radical H increases slightly in the diverging part of the nozzle, while for higher pressures it decreases along the nozzle. Inlet fuel–oxidant mass ratio affects the variation of all species with a greater effect for a near stoichiometric ratio. It was also shown that a higher inlet temperature provides a more enhanced reaction zone in the diverging part of the nozzle.     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程，结合广义Betti-Rayleigh互易等式与时域下的边界积分方程，推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解，将经过主部分析的积分核函数分解为静态和动态部分，其中动态积分核不具有奇异性。在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换，采用配置点法计算超奇异积分，获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序，得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律，数值结果稳定且收敛速度快。     

Modeling and numerical simulation of the resistance spot welding of zinc coated steel sheets using rounded tip electrode: Analysis of required conditions

This paper investigates the essential conditions to improve the accuracy of a resistance spot welding computational study of advanced zinc coated steel sheets using rounded tip electrode. An experimental analysis is performed to highlight the required considerations for a suitable simulation. A sequential Electrical-Thermal-Metallurgical and Mechanical (ETMM) finite element analysis with appropriate precautions of the contact conditions enables to accurately simulate the nugget development during the welding. A critical smooth evolution of the contact radius is required. A fine meshing with an interfacial mesh size of at least 0.05 × 10⁻³ m combined with a coupling time step of 0.0025 s between the electrical-thermal-metallurgical and the mechanical analysis allows a regular incrementation of the contact radius, without burdening the time computing. Accurate values of the contact resistance depending on the interfacial pressure and temperature are essential for a good simulation of the nugget size. The ETMM calculation is successfully extended to the simulation of the welding of a typical two sheets assembly.     

A study on Brinkman number variation on water based nanofluid heat transfer in partially heated tubes

Even if the variation of Nusselt number with Reynolds number has been observed before, little to no studies are concerning nanofluids. The present numerical data in the laminar regime were also found to correlate well with the Brinkman number for the individual sets and much better globally, by combining all the data sets for nanofluids and water.Thus, this work proposes a new correlation between heat transfer rate and Brinkman number which is a nondimensional number of viscosity, flow velocity and temperature. The results showed a good empirical equation that Nu/(Re^0.62 Pr^0.33) is dependent on a power law function to the Brinkman number in laminar flow regime. The proposed correlation can be applied both to water and nanofluid flow, with respect to the nanoparticle loading. It is expected that the equation suggested by this work can be useful to design heating/cooling devices.     

非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

On the use of edge cracked short bend beam specimen for PMMA fracture toughness testing under mixed-mode I/II

In this paper an inclined edge cracked short beam specimen subjected to symmetric three-point bend loading was designed and examined for conducting mixed-mode I/II fracture toughness experiments. The aspect ratio (i.e. length to width ratio) and the loading span distance are considered much lower than the other conventional cracked bend beam samples. Crack tip parameters such as stress intensity factors and T-stress were computed numerically for this specimen by several finite element analyses and it was demonstrated that the specimen is able to produce full combinations of mode I and II including pure mode II. The practical capability of the short bend beam specimen was studied experimentally by conducting a set of mixed-mode fracture tests on PolymethylMethacrylate (PMMA) as a well-known model brittle material. The critical stress intensity factors, the direction of fracture kinking and the path of fracture trajectory were investigated both experimentally and theoretically using two stress and strain-based fracture criteria. The fracture toughness of tested PMMA was decreased by moving towards mode II case due to the effect of T-stress on the fracture mechanism of the short bend beam specimen.