A nonlinear model for marble sulphation including surface rugosity and mechanical damage

Here we propose and analyze a mathematical model that aims to describe the marble sulphation process occurring in a given material. The model accounts for rugosity as well as for damaging effects. This model is characterized by some technical difficulties that seem hard to overcome from a theoretical viewpoint. Therefore, we introduce some physically reasonable modifications in order to establish the existence of a suitable notion of solution on a given time interval. Numerical simulations are presented and discussed, also in view of further research.     

Time reversal for obstacle location in elastodynamics from acoustic recording

The Note is concerned with a feasibility study of time reversal in a non-homogeneous elastic medium, from data recorded in an acoustic medium. Our aim here is to determine the presence and some physical properties of elastic “inclusions” (unknown, not observable solid objects, characterized by their elastic properties) from partial observations of acoustic waves scattered by these inclusions. A finite element numerical method, based on a variational acousto-elastodynamics formulation, is derived and used to solve the forward, and then, the time-reversed problem. A criterion, derived from the reverse time migration framework, is introduced, to help construct images of the inclusions to be determined. Numerical illustrations on configurations that mimic the breast cancer configuration are proposed, and show that one can differentiate between two inclusions, even with different properties.     

Lamb波频域逆散射全聚焦成像方法研究

Lamb波因其传播距离远、衰减小常被用于板状结构的无损检测中,在基于Lamb波损伤检测的诸多成像技术中,全聚焦方法(Total Focus Method,TFM)方法因其成像分辨率高、信噪比高而受青睐。然而Lamb波的频散效应导致时域延时量不能被准确计算,进而影响传统TFM方法对损伤定位及成像的精度;此外,既有的TFM方法仅以回波幅值作为成像指标,忽略了Lamb波与损伤的相互作用,故而不能通过损伤表面的物理参数增强成像质量。针对这两个问题,本文首先在时域TFM基础上发展了频域TFM,在计算中纳入频散关系以规避频散的影响;其次以包含明确的损伤特征参数--反射率为成像指标,结合频域TFM方法建立损伤逆散射模型,以实现对损伤的准确成像。仿真和实验结果表明:频域逆散射TFM成像方法能够实现对铝板结构中的损伤检测,在工程实践中具有应用价值。     

基于NICE-OHMS技术进行大气压气样直接检测的理论分析

噪声免疫腔增强光外差分子光谱技术(NICE-OHMS)是目前世界上最灵敏的激光吸收光谱技术,其在低压环境中具有极高的探测灵敏度。然而当测量样品处于大气压时,NICE-OHMS系统的探测灵敏度会大幅下降。主要原因之一是大气压下获取最大NICE-OHMS信号幅度的条件与低气压下不同。通过对大气压NICE-OHMS理论进行分析,分析了影响信号幅度的参数,并通过数值模拟来寻找最佳的实验条件。本文着重讨论影响信号的主要参数包括光学腔腔长L,调制系数β,探测相位θ。其中,由于在NICE-OHMS中使用DeVoe-Brewer技术将调制频率ν_m锁定到Fabry-Parot(FP)腔的自由光谱区(FSR)。因此FP腔的腔长决定了ν_m,同时还作用于信号幅度S■。模拟结果显示,当腔长增大时,由于ν_m随之减小,载波和边带的光谱成分相互重叠部分增大,因此线型函数的幅度逐渐减小。而吸收信号幅度随着腔长的增加而逐渐增加,色散信号幅度先增大后减小,并且在腔长等于8 cm时达到最大值。调制系数β会影响频率调制后激光载波和边带的幅度大小,并且影响信号线型。随着腔长的增加,最大信号幅度对应的β值也随之增加。在相同腔长下,色散信号的最佳β值小于吸收信号,更容易使用电光调制器实现。最后分析了参数的可实现性,分析了不同种类激光器的频率调谐能力,压电陶瓷的扫描宽度等。以乙炔气体为例,大气压下NICE-OHMS的谱线半宽达到~3 GHz,而光谱覆盖范围大于10 GHz。分布反馈式半导体激光器(DFB)与外腔二极管激光器(ECDL)的频率调谐范围可以达到30 GHz以上,但是由于激光线宽宽,得到的PDH锁定性能欠佳。回音壁模式激光器(WGM)和掺饵光纤激光器(EDFL)线宽为百Hz量级,是目前高灵敏NICE-OHMS系统中常用的光源。但是WGM目前可以实现了5 GHz的激光频率调谐范围,而EDFL的外部电压可控制的调谐范围仅为3 GHz。使用精细度为55000的腔进行模拟,调制系数β=1,腔长大于8 cm时,可使用WGM激光器实现,腔长大于25 cm时,可以使用EDFL激光器实现。而对于在设计光学腔中常用的伸缩长度为25μm的PZT,随着腔长的增加,对应的腔模频移范围逐渐减小,在腔长为典型的40 cm时,扫描范围大于12 GHz。     

真空条件下LiCl溶液蒸发量随浓度变化规律研究

除湿溶液再生是除湿空调系统的核心过程,为了提高除湿溶液再生效率,以蒸发量作为评价标准,在不同的真空度下对不同浓度LiCl溶液进行再生。对不同再生条件下,溶液的再生情况进行数学模拟及实验研究,并将模拟结果与实验结果进行对比分析,结果表明,实验结果同模拟结果具有相同的变化趋势。在6 kPa和8 kPa的压强下,LiCl溶液的蒸发量随浓度的增大先上升后下降,而在1 kPa的压强下,LiCl溶液的蒸发量随浓度的增大呈现出先下降再升高接着又下降的趋势。     

Assembling hydrodynamic cloaks to conceal complex objects from drag

Transformation hydrodynamics and the corresponding metamaterials have been proposed as a means to exclude the drag force acting on an object. Here, we report a strategy to deploy the hydrodynamic cloaks in a more practical manner by assembling different-shaped cloaking parts. Our strategy is to first model a square-shaped cloak and a carpet cloak and then combine them to conceal a more complex-shaped space in the three-dimensional hydrodynamic flow. With the derivation of transformation hydrodynamics, the coordinate transformations for each hydrodynamic cloaking are demonstrated with the calculated viscosity tensors. The pressure and velocity fields of the square, triangular (carpet), and exemplary three-dimensional house-shaped cloaks are numerically simulated, thus showing a cloaking effect and reduced drag. This study suggests an efficient way of cloaking complex architectures from fluid-dynamic forces.     

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines. Friction factors play an important role in the accurate calculation of pressure drop. Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall, gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions. In this paper by modification of a friction model available in the literature, an improved semiempirical model is proposed. The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved. Pressure gradient and velocity profiles are validated against experimental data. Using the improved model, the pressure gradient deviation from experiments diminishes by about 3%; the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.     

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.     

Simple Equations Method (SEsM): Algorithm,Connection with Hirota Method,Inverse Scattering Transform Method,and Several Other Methods

The goal of this article is to discuss the Simple Equations Method (SEsM) for obtaining exact solutions of nonlinear partial differential equations and to show that several well-known methods for obtaining exact solutions of such equations are connected to SEsM. In more detail, we show that the Hirota method is connected to a particular case of SEsM for a specific form of the function from Step 2 of SEsM and for simple equations of the kinds of differential equations for exponential functions. We illustrate this particular case of SEsM by obtaining the three- soliton solution of the Korteweg-de Vries equation, two-soliton solution of the nonlinear Schrödinger equation, and the soliton solution of the Ishimori equation for the spin dynamics of ferromagnetic materials. Then we show that a particular case of SEsM can be used in order to reproduce the methodology of the inverse scattering transform method for the case of the Burgers equation and Korteweg-de Vries equation. This particular case is connected to use of a specific case of Step 2 of SEsM. This step is connected to: (i) representation of the solution of the solved nonlinear partial differential equation as expansion as power series containing powers of a “small” parameter ϵ; (ii) solving the differential equations arising from this representation by means of Fourier series, and (iii) transition from the obtained solution for small values of ϵ to solution for arbitrary finite values of ϵ. Finally, we show that the much-used homogeneous balance method, extended homogeneous balance method, auxiliary equation method, Jacobi elliptic function expansion method, F-expansion method, modified simple equation method, trial function method and first integral method are connected to particular cases of SEsM.     

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.