_{1}

^{*}

The examination of wave motions is traditionally based on the differential equation of D’Alambert , the solution of which describes the motion along a single dimension, while its bidimensional extension takes on the concept of plane waves. Considering these elements and/or limits, the research is divided into two parts: in the first are written the dif ferential equations relating on the conditions two/three-dimensional for which the exact solutions are found; in the sec ond the concepts are extended to the analysis of the propagation of wave motions in porous media both artificial and natural. In the end the work is completed by a series of tests, which show the high reliability of the ph y s i cal-mat h ema t i cal models proposed.

The mathematical theory, that analysed the propagation of elastic impulses, has its origin from the mechanics of the wave motions the latter are defined as the integral of a force generated naturally (example: earthquake) or artificially (example: geophysical prospecting; vehicular traffic), time-averaged:

Also, the wave nature assumed for the analysis of the phenomenon implies the propagations time of the impulse in accordance with the relationships between the dynamic elastic modules—longitudinal (E_{d}), tangential (G_{d}) and volumetric (K_{d})—through the dynamic Poisson’s ratio (ν_{d}) in turn dependent on the speed of compressions waves (v_{S}):

Known these elements (

the vibratory phenomenon is attributable to the study of elastic behaviour of the individual waves, that exploit the properties of sine and cosine functions and are repeated periodically in the initial characteristics of the motion, expressed in terms of amplitude A, frequency f and phase costant δ (

the angular pulsation (or angular frequency):

And the number of the wave:

In turn the Equations (4) and (5) can be combined to derive the expression of velocity:

Given that the real motions (example: earthquake) can be simulated by the addition of a suitable number of monochromatic waves (

Under these conditions, a mathematical model of general validity and its extension to the porous media will be described in the continuation of research.

The analysis of vibratory motions originates from the

D’Alembert’s differential equation [4,5] which analyze the propagation of monodimensional waves in function of the vectorial component u of the movement (variable between A = 0 and A = 1—_{x}:

The solution of this problem, provided by the author, has been perfected by Euler [

In fact, calculated the second derivatives of the Equation (8) as a function of space and time:

you can replace the Equations (9b), (9d) and (6) in the Equation (7) to obtain the identity:

Alternatively, can be introduced only Equations (9b) and (9d) to obtain the Equation (7).

The extension of D’Alembert equation to the two-dimensional case is an unsolved problem that has required the introduction of isotropic homogeneous media in which the motions are propagated as circular waves (

Each monochromatic components of a vibratory motion which propagates in a homogeneous medium travels at the same speed, in turn, dependent on the angular fre-

quency and the wave number through Equation (6); in this case, indicated with ρ the density of the medium, it proves [

Besides, it is also the relationship v_{P}_{ }> v_{S} for the Equation (2a). In inhomogeneous media, for which the speed varies along the path, each frequency component of the vibratory motion is instead equipped with its own speed, such as to arrive to a generic receiver at different times. The described phenomenon, known as the dispersion of the phases, implies that in the means of non-dispersive (homogeneous) the vibratory motion arrives at the receivers with the same initial shape; on the contrary, in dispersive media is possible to identify a group velocity (if you choose to characterize the group through the analysis of the maximum amplitude) and a phase velocity of each component, so that a signal emitted by the source arrives distorted at the receiver. Obviously, means in the two non-dispersive group velocity and phase are.

Denote by s(x,y,z) the displacement vector in space, defined with respect to a reference system coordinated, having components u, v e w; consequently, if the Equation (7) describes the motion of a disturbance which is propagated only in the direction x, the generalized form of the equation of motion will be of the type:

(12a)

The Equation (12), expressed in matrix and extended form, in the case of the means non-dispersive, reduces to the form:

in which v_{j} identifies the column vector of the speed.

Then, the solution can still be expressed using the Equation (8) on condition to take account of the threedimensionality of the phenomenon:

By repeating the procedure previously seen with Equation (7) of the motion 1D, you get in sequence:

In Equation (16), replace the components of the velocity, expressed according to Equation (6), finally to prove the validity of the proposed solution:

The latter also demonstrates that the fundamental period of the wave (and consequently the frequency fundamental) corresponds to the period (and frequency) of the individual components.

In conclusion, for orthotropic media and from the structure of Equation (14), which is easily seen that the waves propagate according ellipsoids scalene (_{z}_{ }> v_{y} = v_{x} we are seeing waves with the form of ellipsoids of revolution (_{z} = v_{y} = v_{x} leads to the development of spherical waves (

obtaining:

Place also:

you get:

Finally, putting too:

we come to the equation of an ellipsoid:

On the other hand, the Laplacian content in Equation general (13a)—Ñ^{2}s—is the quadratic form of the divergence of the vector functions that, for positives values in a generic point P denotes the existence of an outgoing flow in the neighborhood of P; so, to the condition imposed v_{z}_{ }> v_{y}_{ }> v_{x} corresponds w > v > u which describes an ellipsoid. The discussed model is independent of temperature of the system assumed adiabatic, that respects the law of conservation of energy.

In case of the plane x z the Equation (13b) is reduced on the form:

whose solution is a special case of the Equation (14):

Again, the procedure is applied now known we arrive at the identity relation:

which proves the accuracy of the solution according to the propagation of the waves on a plane with elliptical or circular forms.

The problem of waves propagation in porous media has been addressed previously in [

In the Equation (28) appear porosity (n), the degree of saturation (S) and the density of the solid phase (ρ_{s}), fluida (ρ_{w}) and gaseous (ρ_{a}) capable of describing a structure consisting of a skeleton solid with interstitial interconnected pores between them (

Since the mass is connected with the speed, through the Equations (11), (2a) and (2b) shows a dependence of the type:

In other words, must be considered as a porous solid consists of three continuous elastic means that occupy the same region of space and which interact between them carving up the propagation of the same elastic pulse.

The problem outlined can be mathematically simplified if one analyzes the propagation of three waves decoupled from the same impulse, through the solid skeleton and phases fluid and gases contained in the pores (

Analyzing the Equation (30) analyzing the equation it turns out that the first bracket identifies the anisotropic component, relative to the propagation through the solid skeleton, linked with only the elastic modulus being able to neglect the effect of the transverse contraction [_{w} e K_{a}).

The symbol E^{*}_{d,i} identifies a law of variation of the longitudinal elastic modulus dynamic, determined from [

In the Equation (31) appear the parameterawhich represents a material constant and must be determinated experimentally, meanwhile the value of ρ, presents in the Equations (30) changes as a function of the porosity according to the law (28) of the mass balance.

Once defined the general structure of the equations, the problem can be further simplified if one considers that the speed of propagation of the impulses in water and air are approximately v_{w} ≈ 1500 m/s and v_{a} ≈ 340 m/s which would alter the Equations (30):

In conclusion, the same way as seen with the Equations (13) ÷ (19) related to the generalized theory of waves, the Equations (30) and (32) describe scalene ellipsoidal waves, ellipsoid of revolution or spherical with respect to the relationships between the elastic modulus of the soil skeleton that affect the dynamic response of the anisotropic component.

In the determination of the speed of shear waves is necessary to consider: 1) the physical impossibility of waves to propagate in fluids; 2) their polarization on mutually orthogonal planes (waves S_{V} e S_{H}); 3) their development in all directions with respect to the source; 4) to Equation (11b); 5) a law of variation of Poisson’s ratio:

Finally, the Equations (2a), (31) and (34) can be combined with each other:

Using the

Equations (35), compared with the Equations (32), demonstrate the existence of time differences of arrival of the shear wave—compared to compression—which must necessarily increase with the increase of the distance from the source.

The graph of

function of porosity. In detail, as reported in [

At the same graph have been superimposed the theoretical laws of variation of P-waves velocity, and their tabulated data, calculated using Equations (28), (31) and (32c); as can be noted from the same theoretical model fails to approximate, for a = 2.126 introduced in the calculation of, the speed difference between the dry and saturated rocks with the increase of the porosity.

Final interpretation of the results with the theoretical model must be considered that the theoretical behavior of lava rocks (lavas of Etna and Campi Flegrei) and pyroclastics (tuffs and ignimbrites) was standardized in terms of ρ and of, although the different mineralogical composition (Basalts and Trachytes alkaline in the first case, the variable in the second case [

The second test was performed using the experimental data contained in [_{P}-n and v_{S}-n measured in samples of alumina ceramic (Al_{2}O_{3}) known for the acid resistance and low thermal conductivity so as to be used as a catalyst in the chemical industry and as a graft material in the biomedical; therefore, with reference to

The

in the test: 1) the P-waves velocity for saturated samples is slightly lower than that responsible for the dry samples n < 60%; 2) for n > 60% the P-waves velocity of saturated samples becomes greater than that of the dry samples while the model manifests a marked non-linearity; 3) for n = 100% the speed of the saturated samples is reduced to that of water (point A) while that of the samples dried in the air (point B). Finally, as expected, n = 100% for the S-waves velocity vanishes.

The trends of the P-waves velocity associated to the fields 0 < n < 0.6 and 0.6 < n < 1depend on the laws of variation adopted for ρ and E^{*}_{d} that, described by Equations (28) and (31), respectively are linear and non-linear; which means that, while the density varies linearly for 0 < n < 1 the elastic modulus assumes a-heating an approximately linear for n < 0.6 and not linear for n > 0.6. Consequently, the increase of speed setting attributable to fluid phases fails to compensate—for n < 0.6— the effects adducts by variations in density and elastic modulus.

The last test was conducted using the experimental results contained in [

The analysis of wave motions takes origin by the differential equation of D’Alembert [4,5], the solution of which, based on the properties of the trigonometric sine and cosine functions, simulates the propagation of an elastic wave in a continuous medium 1-dimensional, and its extension to the floor is brought back to the simplified concept of circular waves reduced to plane waves.

After a brief review of the main elements that make up the rolling waves of D’Alembert, was written the differential equation governing the wave motion in space (generalized theory) to follow, it was found that the exact solution shows that the waves take the form of ellipsoids scalene in orthotropic means which, in turn, reduce to the ellipsoids of revolution in the means transversely isotropic and spherical waves in isotropic media; in the same way, which step consequent, the 3D equation has been reduced to the 2D field whose solution leads to waves having elliptical shapes that are reduced to circular waves. But be aware that these equations are valid for continuous media, which can be described by assigning them appropriate scalar and vector fields defined by means of functions regular and continuous over the entire domain configuration.

Which next step, given the nature of the particle actually real media (with particular reference to geomaterial), the equations have been applied to porous media, the description of which is based on concepts known in Geotechnical that predict the existence of three continuous media (solid skeleton with interparticle voids paths from water and air) that interact between them carving up the propagation of the same impulse elastic divided into the components of compression and shear. Finally, the search has been completed with some tests based on known experimental data relating to the propagation of waves in different types of rocks and in samples of alumina ceramics.

Ultimately, the tests showed that the proposed models fail to accurately simulate the behavior of a particular material or group of materials that share the same origin while further applications may be later developed: 1) in

the field seismological, considered that the Equations (30) reduce to Equations (11a) for means infinitely porous, as in the case in which the seismic compressional waves are reduced to acoustic waves when intercept the earth’s surface and propagate into the atmosphere; 2) still in the field seismological or in the analysis of vibrational motions artificial, through the introduction of damping function of the signal; 3) for the determination of the porosity of geomaterials coarse loose that can not be sampled and measured in the laboratory.

The author wishes to thank Luca Lussari, Department of Mathematics and Physics of “University of Sacro Cuore of Brescia”—Italy, for his suggestions and observations.