Jtam.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 4, 39, 2001 DETERMINATION OF PERMEABILITY AND TORTUOSITY OF PERMEABLE MEDIA BY ULTRASONIC METHOD. STUDIES FOR SINTERED BRONZE Jan Kochański Mariusz Kaczmarek Józef Kubik Department of Environmental Mechanics, Bydgoszcz University e-mail: mkk@rose.man.poznan.pl The paper presents a combined graphical and analytical method of de- termination of two structural parameters of a permeable medium: per- meability and tortuosity. Themethod is based on a two-phasemodel of permeable material and uses experimental data from ultrasonic studies of wave parameters in alcohol-saturated sintered bronze. Key words: saturated porousmaterials, ultrasonic waves, attenuation 1. Introduction The fundamental parameters characterizing geometrical structure of po- rous permeable materials are: volume porosity fv, permeability k, and tor- tuosity α. Such parameters play an important role in the motion equations of the two-phase model of permeable media and particularly, they determine interaction forces Rs and Rf between the pore fluid and the porous skeleton (Biot, 1956; Kaczmarek and Kubik, 1993) R s =−Rf = b(vf −vs)+ρf(α−1) ∂ ∂t (vf −vs) (1.1) where b= ηf2v/k, v f and vs denote the velocities of fluid and solid skeleton, ρf is the partial density of fluid, and η is the dynamic viscosity of fluid. The values of the above pore structure parameters of permeablematerials can be evaluated by the analysis of solutions of wave equations for the case of harmonic excitation, obtained within the two-phase model of a saturated 924 J.Kochański et al. porous material (Kaczmarek and Kubik, 1993; Kochański, 1998). Then, the following general relationships must be satisfied {Vj,Aj}= {Vj,Aj}{A,Q,R,N,η,ρ s,ρf,fv,α,k,ω} (1.2) where Vj and Aj denote wave velocities and coefficients of attenuation for the three modes propagating in porousmedia (fast, slow, and shear wave, re- spectively). A,Q,R,N are the coefficients which determine elastic properties of a saturatedmedium, ρf, ρs are partial densities of fluid and solid, and ω is the angular frequency.Using the relation (1.2) and the appropriate experimen- tal data for wave parameters (phase velocity and attenuation coefficient), one can identify the structural parameters of a porous material if an appropriate searching procedure is applied. This paper presents a procedure of determination of the two structural parameters: permeability and tortuosity for porous sintered bronze. 2. Experimental results Experimental results for themodel of porousmaterialsmade from sintered bronze particles of the average grain diameter equal to 70µm and saturated with ethylene alcohol are performed. The experimental technique is based on the broad-band ultrasonic transmission method (Kaczmarek et al., 1998; Kochański, 1998). A cylindrical sample of sintered bronze is placed between transmitting and receiving transducers immersed in liquid (alcohol). The sam- ple can be rotated to change the angle of wave incidencewhich serves as away to separate different wave modes (fast, slow and shear wave) (Plona, 1980; Kochański, 1998). From the three registrated wave modes the data obtained for the slow wave have been used to identify the structural parameters. The phase velocity and coefficient of attenuation of the wave are shown in Fig.1. 3. Determination of structural parameters In order to determine the permeability k and tortuosity α, a combined graphical and analytical method based on searching for the common point (point of intersection) of zero error functions of the wave velocity ∆V and attenuation ∆A is used. In Fig.2 the error functions for fixed frequency Determination of permeability and... 925 Fig. 1. Attenuation coefficient (a) and phase velocity (b) of a slow wave for alcohol-saturated bronze with the average grain size 70µm Fig. 2. Illustration of the steps of the proposedmethod of determination of structural parameters. Determination of the common point for ∆V =0 and ∆A=0 926 J.Kochański et al. (0.4Mhz) are plotted versus permeability k and tortuosity α. The error functions are defined as follow ∆V = [Vpom−Vmod(α,k)](ω=const) ∆A= [Apom−Amod(α,k)](ω=const) where Vpom, Vmod(α,k), Apom, Amod(α,k) are the velocities and attenuation coefficients determined from experiments and calculated from the two-phase model as the functions of α and k (see Kaczmarek and Kubik (1993) for the review of dynamical methods of determination of material parameters of the two-phase model). The proposed method applied with experimental data obtained for the sintered bronze and frequency equal to 0.4MHzallowed us to find the following values of permeability and tortuosity: k=2.04 ·10−12m2, α=1.85. 4. Analysis of the obtained results and discussion Fig. 3. Influence of the frequency on the values of structural parameters (the structural parameters should be identified with coordinates of the points) In order to test the sensitivity of the proposed method to the frequency, for which the structural parameters are identified, the same procedure was Determination of permeability and... 927 applied for various values of frequency from 0.05 to 1.2MHz. The results of calculations are presented in Fig.3. It can be seen that for relatively low frequencies (below 0.2MHz), the results for permeability and tortuosity are very close to each other and remain within the range from 2.45 · 10−12 to 2.52 · 10−12m2 for permeability, and from 1.8 to 1.84 for tortuosity. For the frequencies above 0.2MHz, the tortuosity changes are still insignificant (from 1.85 to 1.88) but permeability reaches the values from 0.33·10−12 to 2.45·10−12m2. Thus, it can be assumed that the values of permeability and tortuosity derived for frequencies below 0.2MHz are close to the real values of the material studied. One can expect that the observed dispersion of the parameters for higher frequencies are the result of scattering of the waves on grains of the material (Kaczmarek et al., 1998). The latter effect is not included into the classical model of the two- phasematerial and thus itmay essentially influence the errors of the identified structural parameters. In order to incorporate the influence of scattering on the wave parameters, we propose (at the present stage of themodelling it is a postulate) a modification of the model by the following frequency correction of the parameter b b= η f2v k [1+χ(f)] χ(f)= lf2+nf4 (4.1) The best fit to the experimental data is obtained when l = 10MHz−2 and n = 2.5MHz−4. Application of the function (4.1)1 to the two-phase model of porous medium allows one to minimize the dispersion of the structural parameters derived for various frequencies. The values of permeability are between 2.46 ·10−12 and 2.54 ·10−12m2 while the values of tortuosity range from 1.81 to 1.87, see Fig.4. Fig. 4. Distribution of structural parameters for corrected form of coefficient b 928 J.Kochański et al. The phase velocity and attenuation coefficient of the slow wave predicted by the model with uncorrected and corrected coefficient b, for k = 2.5 · 10−12m2, and α=1.85, were comparedwith the experimental data giving an excellent agreement of the model predictions with the experimental results. References 1. BiotM.A., 1956,TheoryofPropagationofElasticWaves in aFluid Saturated Porous Solid. I. LowFrequency Range, J. Acoust. Soc. Am., 28, 2, 168-178 2. Kaczmarek M., Kubik J., 1993, Dynamiczne metody wyznaczania charak- terystykmateriałowych nasyconych ośrodkówPorowatych,Prace IPPT, 16 3. Kaczmarek M., Kochański J., Kubik J., 1998, UltrasonicWaves in Satu- rated Porous Materials. Discussion of Modeling and Experimental Results, J. Theor. Appl. Mechanics, 36, 597-618 4. Kochański J., 1998, Wyznaczanie parametrów struktury ośrodka przepusz- czalnego na podstawie pomiarów fal ultradźxwiękowych, RozprawaDoktorska, Politechnika Poznańska 5. Plona T.J., 1980, Observation of a SecondCompressionalWaves in a Porous Medium at Ultrasonic Frequencies,Appl. Phys. Lett., 36, 259-261 Wyznaczanie przepuszczalności i krętości ośrodka przepuszczalnego przy wykorzystaniu wolnej fali ultradźwiękowej. Badania dla porowatego spiekanego brązu Streszczenie W pracy przedstawiono analityczno-graficzną metodę wyznaczania dwóch para- metrówstrukturyośrodkaprzepuszczalnego:przepuszczalności i krętości.Proponowa- na metoda bazuje na dwufazowym modelu ośrodka przepuszczalnego i na wynikach pomiaru parametrów propagacji fali ultradźwiękowej – wolnej w nasyconym ośrodku porowatym. Pomiary wykonano dla porowatego brązu nasyconego alkoholem. Manuscript received April 6, 2000; accepted for print February 5, 2001