Modelling the attenuation in the ATHENA finite elements code for the ultrasonic testing of austenitic stainless steel welds

Multipass welds made in austenitic stainless steel, in the primary circuit of nuclear power plants with pressurized water reactors, are characterized by an anisotropic and heterogeneous structure that disturbs the ultrasonic propagation and makes ultrasonic non-destructive testing difficult. The ATHENA 2D finite element simulation code was developed to help understand the various physical phenomena at play. In this paper, we shall describe the attenuation model implemented in this code to give an account of wave scattering phenomenon through polycrystalline materials. This model is in particular based on the optimization of two tensors that characterize this material on the basis of experimental values of ultrasonic velocities attenuation coefficients. Three experimental configurations, two of which are representative of the industrial welds assessment case, are studied in view of validating the model through comparison with the simulation results. We shall thus provide a quantitative proof that taking into account the attenuation in the ATHENA code dramatically improves the results in terms of the amplitude of the echoes. The association of the code and detailed characterization of a weld’s structure constitutes a remarkable breakthrough in the interpretation of the ultrasonic testing on this type of component.     

Measurement of ultrasonic scattering attenuation in austenitic stainless steel welds: Realistic input data for NDT numerical modeling

Multipass welds made of 316L stainless steel are specific welds of the primary circuit of pressurized water reactors in nuclear power plants. Because of their strong heterogeneous and anisotropic nature due to grain growth during solidification, ultrasonic waves may be greatly deviated, split and attenuated. Thus, ultrasonic assessment of the structural integrity of such welds is quite complicated. Numerical codes exist that simulate ultrasonic propagation through such structures, but they require precise and realistic input data, as attenuation coefficients. This paper presents rigorous measurements of attenuation in austenitic weld as a function of grain orientation. In fact attenuation is here mainly caused by grain scattering. Measurements are based on the decomposition of experimental beams into plane-wave angular spectra and on the modeling of the ultrasonic propagation through the material. For this, the transmission coefficients are calculated for any incident plane wave on an anisotropic plate. Two different hypotheses on the welded material are tested: first it is considered as monoclinic, and then as triclinic. Results are analyzed, and validated through comparison to theoretical predictions of related literature. They underline the great importance of well-describing the anisotropic structure of austenitic welds for UT modeling issues.     

Probabilistic approaches to compute uncertainty intervals and sensitivity factors of ultrasonic simulations of a weld inspection

For comprehension purpose, numerical computations are more and more used to simulate the propagation phenomena observed during experimental inspections. However, the good agreement between experimental and simulated data necessitates the use of accurate input data and thus a good characterization of the inspected material. Generally the input data are provided by experimental measurements and are consequently tainted with uncertainties. Thus, it becomes necessary to evaluate the impact of these uncertainties on the outputs of the numerical model. The aim of this study is to perform a probabilistic analysis of an ultrasonic inspection of an austenitic weld containing a manufactured defect based on advanced techniques such as polynomial chaos expansions and computation of sensitivity factors (Sobol, DGSM). The simulation of this configuration with the finite element code ATHENA2D was performed 6000 times with variations of the input parameters (the columnar grain orientation and the elastic constants of the material). The 6000 sets of input parameters were obtained from adapted statistical laws. The output parameters (the amplitude and the position of the defect echo) distributions were then analyzed and the 95% confidence intervals were determined.     

各向异性焊缝缺陷超声阵列全聚焦成像方法

针对奥氏体不锈钢焊接构件安全评价需要,对各向异性结构中超声波传播特性进行了研究,利用射线追踪法确定了各向异性介质中声波传播路径。在传统全聚焦成像基础上,发展了一种用于各向异性焊缝中缺陷检测的超声阵列全聚焦成像方法。通过数值仿真和实验,研究了介质材质对超声阵列成像的影响,结果表明,发展的全聚焦成像方法可以很好实现各向异性焊缝中缺陷检测,缺陷定位更准确。项目研究工作为奥氏体不锈钢焊缝检测提供了可行的技术方案.     

Ultrasonic wave propagation in heterogeneous solid media: theoretical analysis and experimental validation

This research deals with the ultrasonic characterization of thermal damage in concrete. This damage leads to the appearance of microcracks which then evolve in terms of volume rate and size in the material. The scattering of ultrasonic waves from the inclusions is present in this type of medium. The propagation of the longitudinal wave in the heterogeneous media is studied via a homogenization model that integrates the multiple scattering of waves. The model allows us to determine the phase velocity and the attenuation according to the elements which make the medium. Simulations adapted to the concrete are developed in order to test the responses of the model. These behaviors are validated by an experimental study: the measurements of phase velocity and attenuation are performed in immersion, with a comparison method, on a frequency domain which ranges from 160 kHz to 1.3 MHz. The analysis of different theoretical and experimental results obtained on cement-based media leads to the model validation, on the phase velocity behavior, in the case of a damage simulated by expanded polystyrene spheres in granular media. The application to the case of a thermally damaged concrete shows a good qualitative agreement for the changes in velocity and attenuation.     

Lamb waves propagation in elastic plane layers with a joint strip

The Lamb waves are used for the ultrasonic characterization of welds because of their relative long-range propagation. In this paper, a simplified model of a weld-strip between two identical semi-infinite elastic layers is investigated. The reflected and transmitted ultrasonic fields are expressed by modal series whose coefficients are obtained by application of orthogonality relation. Comparisons with solutions obtained by finite elements wave propagation simulations are made. The energy balance between the incident and the scattered waves is also used to verify the accuracy of the obtained modal amplitudes.     

奥氏体不锈钢焊缝中的声学仿真模型与全聚焦成像检测

声波在奥氏体不锈钢焊缝中传播时声束弯曲,为超声成像带来了困难。基于Ogilvy焊缝模型,建立了奥氏体不锈钢焊缝非均匀各向异性声场仿真模型,采用Dijkstra路径搜索算法对各向异性条件下的声传播路径和声传播时间进行了数值模拟和分析。在此基础上,采用Verasonics超声相控阵成像系统,进行了奥氏体不锈钢焊缝的全聚焦成像实验,采集全矩阵回波数据,并结合理论模型计算的声传播路径和相应的传播时间,进行了成像结果修正。结果表明,与均匀介质模型的全聚焦成像结果相比,基于该文模型的焊缝全聚焦成像检测结果具有更高的缺陷定位精度和分辨率,验证了该方法的可行性,为奥氏体不锈钢焊缝成像检测提供了新的思路。     

Numerical predictions and experiments on the free-plate edge mode

In this paper, te edge mode variation is studied with three different methods: the reciprocal work method, already used by Torvik [J. Acoust. Soc. Am. 41 (1967) 346] to model this phenomenon, the S-parameter method and a finite element model that are applied for the first time to the study of the edge resonance. Moreover, laser probe measurements of the edge mode have also been performed and compared to the numerical predictions. The good agreement between the numerical predictions and the experimental data allows full understanding of the resonant phenomenon. The edge resonance is linked to the strong increase in amplitude of two complex Lamb waves, and the edge mode is proved to radiate into the plate as the first symmetrical Lamb mode S(0). Displacements at the edge and away from the edge have been computed and measured to evaluate the spatial and temporal behaviour of the edge mode. The dependence of the edge resonance frequency and amplitude on the Poisson coefficient has also been studied.     

各向异性焊缝中超声相控阵瞬态声场仿真及缺陷定位方法

针对多层各向异性奥氏体不锈钢焊缝中超声相控阵瞬态声场的仿真问题,提出应用高斯声束等效点源模型计算宽带离散化的多个单频稳态声场,通过傅里叶变换将其拓展为瞬态声场,并分析了声场转换过程的主要影响参数。该方法可快速计算焊缝内部超声相控阵聚焦声场的瞬态能量分布和任意一点的时域波形信号。在此基础上针对多层奥氏体不锈钢焊缝内部缺陷的超声相控阵成像检测问题,提出利用上述时域高斯声束法对多通道缺陷散射信号进行时间反转计算,并根据时域声场焦点确定缺陷位置。最后通过实验,验证使用此方法检测实际奥氏体不锈钢焊缝试块内部缺陷的效果。结果表明,提出的方法能够确定缺陷位置,且计算速度快、运算量小,适合作为多层介质内部缺陷实时成像的声场仿真模型。      

Scattering of thermal phonon pulses by isotopes in silicon

The pure ballistic propagation of acoustic phonons in crystals at low temperatures can be described within anisotropic continuum acoustics. One needs only the elastic constants and mass densities to calculate the time-dependent spectral phonon irradition of the bolometer for a given radiator pulse power and detector/radiator geometry. We extend this treatment by including single isotope-scattering events for the phonons in a (111)-cut silicon disk on their flight from the radiator to the detector. Using the earlier experimentally determined polarization and frequency dependent phonon absorption in the bolometer metal, the instantaneous temperature of the bolometer can be calculated from this irradiation. This allows a direct comparison with measured bolometer temperatures using exactly the same transmission or reflection arrangement as in calculation. A very satisfying agreement is observed in the expected range of single phonon scattering.     

Directing ultrasound at the cemento-enamel junction (CEJ) of human teeth: I. Asymmetry of ultrasonic path lengths

The diagnosis of degenerative changes in human teeth is of general interest because early detections can avoid greater health problems and further weakening effects. Since the wear of teeth determines their stability and lifetime in relation to the physiological load, an ultrasonic survey of any dimensional changes of the enamel layer and especially of the dentin wall thickness may be very helpful. However, an ultrasonographic diagnosis requires first to determine the anisotropic human tooth properties at clinically relevant locations and to simulate wave propagation phenomena in inhomogeneous tooth models with proper dimensions. The first article of a series that provides modular data of mineralized tissues in human teeth at the cemento-enamel junction (CEJ) deals with an ultrasonic method for measuring the asymmetry of dimensional characteristics of extracted human teeth and their ultrasonic path lengths (UPL). Heavily attenuating tooth halves were investigated with respect to the symmetry of normal and inclined oppositely directed radial ultrasonic paths. The measured UPLs ranged from 1.2 mm to 4.4 mm. The relative difference in inclined UPLs between the left and the right tooth halves reaches almost 30%. This reveals a large asymmetry. The mean difference of angles that represent fastest path lengths was 2.2+/-8.1 degrees, which indicates large asymmetry and anisotropy. Several aspects, which are required for a proper integration of asymmetric data into models designed for medical element engineering and simulation (MEES), are discussed.     

Backscattering grain noise modelling in ultrasonic non-destructive testing

Abstract

Many materials present an internal grain microstructure. When these materials are subjected to ultrasonic non-destructive testing, the grains behave like scattering centres producing unwanted backscattered noise that can make the detection of true defects difficult. This paper is devoted to the modelling of the probability density and the spacetime correlation functions of the grain noise complex envelope. Assuming statistical independence between any pair of grains, the authors derive analytical expressions for the above functions. Specifically, the envelope comes to be K-distributed, the parameters of the distribution may be related, under reasonable simplifying assumptions, to the material characteristics (grain density, grain size distribution, propagation velocity). The spacetime correlation function is a separable function. It may be expressed as the product of a spatial factor due to the spatial correlation introduced by the non-zero beamwidth, and a time factor due to the time correlation introduced by the non-zero pulse duration. The analytical expressions are verified by means of real data measured in austenitic stainless steel specimens.     

Correction of ultrasonic array images to improve reflector sizing and location in inhomogeneous materials using a ray-tracing model

The use of ultrasonic arrays has increased dramatically within recent years due to their ability to perform multiple types of inspection and to produce images of the structure through post-processing of received signals. Phased arrays offer many advantages over conventional transducers in the inspection of materials that are inhomogeneous with spatially varying anisotropic properties. In this paper, the arrays are focused on austenitic steel welds as a representative inhomogeneous material. The method of ray-tracing through a previously developed model of an inhomogeneous weld is shown, with particular emphasis on the difficulties presented by material inhomogeneity. The delay laws for the structure are computed and are used to perform synthetic focusing at the post-processing stage of signal data acquired by the array. It is demonstrated for a simulated austenitic weld that by taking material inhomogeneity and anisotropy into account, superior reflector location (and hence, superior sizing) results when compared to cases where these are ignored. The image is thus said to have been corrected. Typical images are produced from both analytical data in the frequency domain and data from finite element simulations in the time domain in a variety of wave modes, including cases with mode conversion and reflections.     

Comparison between experimental and 2-D numerical studies of multiple scattering in Inconel600® by means of array probes

Ultrasonic non-destructive testing of polycrystalline structures can be disturbed by scattering at grain boundaries. Understanding and modeling this so-called “structural noise” is crucial for characterization as well as detection purposes. Structural noise can be considered as a fingerprint of the material under investigation, since it contains information about its microstructure. The interpretation of experimental data necessitates an accurate comprehension of complex phenomena that occur in multiple scattering media and thus robust scattering models. In particular, numerical models can offer the opportunity to realize parametrical studies on controlled microstructures. However, the ability of the model to simulate wave propagation in complex media must be validated. In that perspective, the main objective of the present work is to evaluate the ability of the finite-element code ATHENA 2D to reproduce typical features of multiple wave scattering in the context of ultrasonic non-destructive evaluation, with an array of sources and receivers. Experiments were carried out with a 64-element array, around 2 MHz. The sample was a mock-up of Inconel600® exhibiting a coarse grain structure with a known grain size distribution. The numerical model of this microstructure is based on Voronoi diagrams. Two physical parameters were used to compare numerical and experimental data: the coherent backscattering peak, and the singular value distribution of the array response matrix. Though the simulations are 2-D, a good agreement was found between simulated and experimental data.     

Quantitative evaluation of ultrasonic C-scan image in acoustically homogeneous and layered anisotropic materials using three dimensional ray tracing method

Quantitative evaluation of ultrasonic C-scan images in homogeneous and layered anisotropic austenitic materials is of general importance for understanding the influence of anisotropy on wave fields during ultrasonic non-destructive testing and evaluation of these materials. In this contribution, a three dimensional ray tracing method is presented for evaluating ultrasonic C-scan images quantitatively in general homogeneous and layered anisotropic austenitic materials. The directivity of the ultrasonic ray source in general homogeneous columnar grained anisotropic austenitic steel material (including layback orientation) is obtained in three dimensions based on Lamb’s reciprocity theorem. As a prerequisite for ray tracing model, the problem of ultrasonic ray energy reflection and transmission coefficients at an interface between (a) isotropic base material and anisotropic austenitic weld material (including layback orientation), (b) two adjacent anisotropic weld metals and (c) anisotropic weld metal and isotropic base material is solved in three dimensions. The influence of columnar grain orientation and layback orientation on ultrasonic C-scan image is quantitatively analyzed in the context of ultrasonic testing of homogeneous and layered austenitic steel materials. The presented quantitative results provide valuable information during ultrasonic characterization of homogeneous and layered anisotropic austenitic steel materials.     

Analysis of laser-generated ultrasonic force source at specimen surface and display of bulk wave in transversely isotropic plate by numerical method

Taking into account the effects of thermal diffusion and optical penetration, as well as the finite width and duration of the laser source, the laser-generated ultrasonic force source at surface vicinity is presented. The full acoustic fields of laser-generated ultrasonic bulk wave are obtained and displayed in transversely isotropic plate. The features of laser-generated ultrasound bulk waves are analyzed. The features of laser-generated ultrasonic bulk wave are in good agreement with the theoretical results (the phase velocity surfaces), demonstrating the validity of this simulation. The numerical results indicate that the features of laser-generated ultrasound waveforms in anisotropic specimen, different from the case in isotropic materials, have a close relation with the propagating plane and propagation direction. This method can provide insight to the generation and propagation of laser-generated ultrasonic bulk wave in transversely isotropic material.     

A semi-analytical model for predicting multiple propagating axially symmetric modes in cylindrical waveguides

A semi-analytical model for multiple mode axially symmetric wave propagation in finite solid cylindrical waveguides is presented. The model is designed as a tool for predicting and interpreting experimental signals. The model is based on a common experimental configuration and considers the excitation, propagation and reception of the ultrasonic signal in the waveguide. The Pochhammer-Chree solution for an infinite cylinder is the basis for the model. Extensions are made to enable comparison to experimental results. Comparisons with experiment are performed in the time, frequency and joint-time frequency domain for both narrow band and broad band excitation of the piezo-electric transducer.     

Numerical simulation of transfer and attenuation characteristics of soft-tissue conducted sound originating from vocal tract

A non-audible murmur (NAM), a very weak speech sound produced without vocal cord vibration, can be detected by a special NAM microphone attached to the neck, thereby providing a new speech communication tool for functional speech disorders as well as human-to-machine and human-to-human interfaces with inaudible voice input for use with unimpaired. The NAM microphone is a condenser microphone covered with soft-silicone impression material that provides good impedance matching with the soft tissues of the neck. Because higher-frequency components are suppressed severely, however, the NAM detected with this device can be insufficiently clear. To improve NAM clarity, the mechanism of NAM production as well as the transfer characteristics of the NAM in soft neck tissues must be clarified. We have investigated sound propagation from the vocal tract to the neck surface, using a finite difference time domain method and a head model based on magnetic resonance imaging scans. Numerical results show that, compared to air-conducted sound detected in front of a mouth, soft-tissue-conducted sound attenuates 50 dB at 1 kHz, which consists of 30 dB full-range attenuation due to air-to-soft-tissues transmission loss and −10 dB/octave spectral decay due to a propagation loss in soft tissues. The decay agrees well with the spectral characteristics of the measured NAM.     

Experimental study of light diffraction by standing ultrasonic wave with cylindrical symmetry

We report a new acousto-optic arrangement based on ultrasonic wave with cylindrical symmetry. The theory of light interaction with standing cylindrical ultrasonic wave is experimentally verified in the Fraunhofer region. A very good agreement of experimental results with numerical calculations based on the proposed theory is found. The diffraction pattern consists of ring-shaped diffraction orders which posses a fine structure. The time average light intensity of the whole zeroth diffraction order as a function of the Raman-Nath parameter is investigated. The modulation properties of presented system are examined by means of single photon counting technique. Finally, some potentially useful applications in the laser and fibre technology are suggested.