Optical resonators such as silica micro-spheres and micro-toroids can support whispering gallery modes (WGMs), where light circulates near the resonator surface and is confined by the total internal reflection at the dielectric boundary. Such resonators can exhibit very high quality (Q) factors, since the resonator surface can maintain atomic level smoothness. The combination of high Q factors and small resonator volumes has led to a wide range of applications in sensing, optomechanics, nonlinear optics, and quantum optics. In this dissertation, we introduce a new type of whispering gallery resonators (WGRs) based on micro-droplets in an immiscible liquid-liquid system. Within such an all-liquid platform, it is possible to achieve highly nonlinear coupling between light and liquid that can potentially lead to single-photon level optical nonlinearity. Specifically, we experimentally characterize a droplet (D~500um) of index matching fluid submerged in the water as a high-Q optical resonator, where we use an optical fiber taper to couple light into the droplet through non-contact evanescent coupling. The highest Q-factor observed in the experiment is 2x10^7 which closely matches the upper limit of intrinsic Q-factor set by the material absorption. Given with such a high Q factor, the WGM can exert strong radiation pressure on the droplet interface, push it outward, increase the length of optical path, and produce a red-shift in WGM resonance. Our experimental results have found that the ratio of those resonance shifts and the optical power coupled into the resonator is approximately 60 fm/μW. The result closely matches to our steady-state estimation based on solving the coupled Maxwell-Navier-Stokes equation. To investigate the dynamic interplay of light and liquid, we develop a harmonic oscillator (HO) model to describe the time-domain behaviors of the coupled optofluidic system. We find a good agreement between theoretical predictions and our experimental data. The shift of WGM resonance can potentially be induced by thermal effects. To estimate the magnitude of thermal effects, we also investigate the thermally induced nonlinear behaviors of WGMs in a cylindrical fiber resonator (D~125um), where we change the mechanism of heat dissipation by changing the cladding material (e.g. air and water). For direct temperature measurements, we use a fiber optical resonator with a fiber Bragg grating (FBG) inscribed in the fiber core to observe temperature shifts induced by the high-Q WGMs. Our result shows that the temperature increase in the fiber resonator in the water is 0.13 C, whereas the fiber resonator in air shows ~4.5 C increase in temperature. Our results suggest that the relatively high thermal conductivity of water suppresses thermal nonlinearity by ~50 times, and that the red-shifts of WGMs can largely be attributed to radiation pressure effect. / Ph. D.

The nonlinear beam deflection (BD) technique is used to directly measure and time-resolve the nonlinearly-induced phase shift in a variety of materials. In this technique, a weak probe beam is spatially overlapped, while slightly displaced, with a strong excitation beam while the temporal delay is scanned. The excitation-induced index gradient, which for 3rd-order nonlinearities is proportional to the nonlinear refractive index 16 n2"> of the medium, deflects the weak probe beam. This deflection is determined using a position sensitive segmented detector after propagation to the far field. In this dissertation, we expand our previous work on BD theory to include the effects of the Gaussian spatial beam profile of the excitation, as opposed to a constant index gradient. We also explore the BD signal as we allow the spatial size and relative position of the probe with respect to the excitation beam, r, at the sample to vary to maximize the calculated signal. While the analysis requires numerical solutions, we find a simple empirical fitting function for the BD signal that allows determination of the nonlinear phase shift and thus the nonlinear refraction. We performed BD experiments at near-degenerate photon energies for various spot size ratios which resulted in very good agreement with our simulation results. In order to examine our empirical function the BD signal for various r (0.2-0.6) is measured while keeping the phase shift relatively constant. This helped us isolate the effect of spot size ratio on the BD signal. Our results showed the correct trend for the growth of BD signal as r increases, which is what is expected from our model. We also studied nondegenerate two-photon absorption (ND-2PA) in bulk silicon. We present the results of spectroscopic pump-probe measurements of ND-2PA in silicon across the indirect-gap (1.12 eV). We observed enhancement of the 2PA coefficient as the degree of degeneracy of pump and probe photon energies increased, and the dispersion compares favorably with our recently-developed semi-empirical theoretical model for the dispersion of indirect ND-2PA in silicon. Additionally, we experimentally investigated WTe2 which is a Weyl semimetal. Here, we prepared very thin flakes (10s of microns thick) of WTe2 and investigated the possibility of observing circular dichroism (CD) in pump-probe measurements, pumping at near IR and probing at mid-IR. Although we did not observe any CD, we believe this is because our pump photon energy is far from Weyl nodes and that we need to pump at mid-IR range.

<p>Detta examensarbete har utförts på Fläktwoods AB i Växjö, som producerar stora axialfläktar för olika industriapplikationer. Syftet är att öka kunskapen om fettsmorda axiella kullager genom FE analyser.</p><p>Projektet har genomförts i fem delsteg för att avgöra påverkan av en eller några få parametrar i taget. De studerade parametrarna är: elementstorlek, kontaktstyvhet, last, lagergeometri (dvs. oskulation), ickelinjär geometri och ickelinjära materialegenskaper (dvs. plasticitet).</p><p>Slutsatsen är att elementstorleken bör väljas fint nog för att ge ett jämnt resultat men grovt nog för att beräkningstiden skal vara rimlig. Kontaktstyvheten har inte stor, men tydlig, inverkan på kontakttrycket och penetrationen. Förändringar av oskulationen leder till förändringar i kontaktellipsens form medan olika laster inte påverkar formen på ellipsen, utan snarare storleken. När det handlar om plasticitet är sträckgränsen den viktigaste faktorn att beakta.</p> / <p>This thesis project was carried out at Fläktwoods AB in Växjö who produces large axial fans for different industry applications. The purpose is to increase the knowledge of grease lubricated axial ball bearings through FE analyses.</p><p>The project was executed into five sub steps to determine the influence of one or few parameters at a time. The studied parameters are: mesh density, contact stiffness, load, bearing geometry (i.e. osculation), geometrical nonlinearity and material nonlinearity (i.e. plasticity).</p><p>It is concluded that the mesh density should be selected fine enough to give a smooth result but course enough to give a reasonable calculation time. The contact stiffness has not a major, but a clear, impact on the contact pressure and penetration. Changes of the osculation lead to changes of the contact ellipse shape and applying different load level does not affect the shape of the ellipse but rather the size. When dealing with plasticity the yield strength is the most important factor to take in consideration.</p>

Detta examensarbete har utförts på Fläktwoods AB i Växjö, som producerar stora axialfläktar för olika industriapplikationer. Syftet är att öka kunskapen om fettsmorda axiella kullager genom FE analyser. Projektet har genomförts i fem delsteg för att avgöra påverkan av en eller några få parametrar i taget. De studerade parametrarna är: elementstorlek, kontaktstyvhet, last, lagergeometri (dvs. oskulation), ickelinjär geometri och ickelinjära materialegenskaper (dvs. plasticitet). Slutsatsen är att elementstorleken bör väljas fint nog för att ge ett jämnt resultat men grovt nog för att beräkningstiden skal vara rimlig. Kontaktstyvheten har inte stor, men tydlig, inverkan på kontakttrycket och penetrationen. Förändringar av oskulationen leder till förändringar i kontaktellipsens form medan olika laster inte påverkar formen på ellipsen, utan snarare storleken. När det handlar om plasticitet är sträckgränsen den viktigaste faktorn att beakta. / This thesis project was carried out at Fläktwoods AB in Växjö who produces large axial fans for different industry applications. The purpose is to increase the knowledge of grease lubricated axial ball bearings through FE analyses. The project was executed into five sub steps to determine the influence of one or few parameters at a time. The studied parameters are: mesh density, contact stiffness, load, bearing geometry (i.e. osculation), geometrical nonlinearity and material nonlinearity (i.e. plasticity). It is concluded that the mesh density should be selected fine enough to give a smooth result but course enough to give a reasonable calculation time. The contact stiffness has not a major, but a clear, impact on the contact pressure and penetration. Changes of the osculation lead to changes of the contact ellipse shape and applying different load level does not affect the shape of the ellipse but rather the size. When dealing with plasticity the yield strength is the most important factor to take in consideration.

The knowledge of geometry of molecular complexes formed via molecular association in solution through weak interactions is always important to understand the origin of stability and function of an array of molecules, supramolecular assemblies, and macromolecular networks. Simple 1:1 molecular complexes are very useful in this regard as they provide a model to understand both the nature of these interactions and their structural implications. Several weak noncovalent forces from long range (van der Waal’s, electrostatic, induction, dispersion) to short range (charge transfer) govern the geometry, that is, relative orientation of the two molecules in such a complex. On one hand, we find 1:1 electron donor acceptor (EDA) complexes such as naphthalene-tetracyanobenzene, hexamethylbenzene-chloranil etc. which stack parallel or in slipped parallel geometry in their crystals. On the other, benzene dimer has been found to stabilize in T shaped geometry in all its three physical states. In this thesis, I focus on 1:1 EDA complexes in solution. A good volume of literature is available which deals with the optical studies on the formation of such complexes. It has been suggested that the nature of the intermolecular interactions stabilizing these complexes in the gas phase or in their crystals is modified by the presence of solvent-solute interactions in solution thus bringing in difference in the solution geometry. However, the existing experimental techniques, both optical and magnetic, are unable to determine the exact geometries of 1:1 EDA complexes in solution. This opens an opportunity to probe their geometry in solution. The quadratic nonlinearity or first hyperpolarizability (β) of a molecule is a measure of the change in dipole moment (or polarization) in the second order of the applied electrical field and thus has a purely electronic origin. It is a tensorial property and can be resolved in components along the three dimensions. The number of β components and the nonlinear optical anisotropies in a typical donor-acceptor type dipolar molecule, defined as (equation) (where1, 2, 3 axes define the molecular frame, 1 being the direction along the principal axis of symmetry and pointing from the acceptor toward the donor), are determined by the symmetry /structure of the molecule. It has been shown theoretically that the 1:1 EDA complexes possess large hyperpolarizabilities. In the case of pNA dimers calculation revealed that the geometry of the dimer and its symmetry is important for obtaining the correct estimate of β from its tensorial components. Therefore, it should be possible to use the values of tensorial β components to construct the unknown geometry of such complexes. Experimentally macroscopic depolarization ratios (D and D′) in the laboratory fixed frame (XYZ, X being the direction of polarization and Z the direction of propagation of the incident light), are measured from the polarization resolved intensities of second harmonic scattering from molecules in solution using the hyper-Rayleigh scattering technique. The depolarization ratios are correlated to the anisotropy parameters, u and v through a co-ordinate transformation. In this thesis I, have first, characterized the quadratic nonlinear optical property of a variety of 1:1 electron donor acceptor complexes and used the values of u and v obtained from depolarized hyper-Rayleigh scattering to deduce their geometry in solution. Chapter 1 provides an introduction to the 1:1 electron donor acceptor complexes, their relevance to chemistry and biology. It also contains an introduction to nonlinear optical processes in molecules. The objective of the present work and scope of the investigation carried out in this thesis is presented in this chapter. Chapter 2 describes the details of the experimental polarization resolved HRS technique. The geometrical model adopted for the analysis of the HRS data has also been introduced and the method of analysis has been described in detail in this chapter. Chapter 3 presents the measurement of β values of two series of 1:1 EDA complexes of variously substituted methylbenzenes donors with tetrachloro-p-benzoquinone (CHL) and dicyanodichloro-p-benzoquinone (DDQ) acceptors at 1064 nm. In agreement with recent theoretical results we find large first hyperpolarizabilities for these complexes. The β values are greater than that of the typical push-pull molecule p-nitroaniline (pNA). We also find that in general β decreases with decrease in the donor strength. Chapter 4 presents the β values for the two series of EDA complexes of CHL and DDQ acceptors at 1907 nm. The values of β are less in magnitude at 1907 nm than that at 1064 nm which is due to the dispersion effect in β. In Chapter 5 and 6, it is described how depolarized hyper-Rayleigh scattering can be utilized to probe geometries of 1:1 complexes in solution. Chapter 5 concentrates mainly on 1:1 EDA complexes of CHL and DDQ and TCNB (tetracyanobenzene), while chapter 6 contains examples of other 1:1 molecular complexes where the noncovalent interactions are much weaker, such as in benzene-naphthalene, benzene-methoxybenzene, benzene-hexafluorobenzene and benzene-chlorobenzene pairs. We find the geometry of 1:1 EDA complexes in solution in terms of tilt angle (θ) and twist angle (ϕ) between the donor and acceptor pairs. The angle θ varies from 29°-47° for different pairs of EDA complexes, while ϕ varies within 34° and 38°. We find that the geometry of 1:1 EDA complexes in solution is different (twisted and tilted cofacial and twisted ‘V’) from those in the crystalline or gaseous states (cofacial), if known. We find that both benzene-naphthalene and benzene-chlorobenzene pairs assume twisted ‘T’ shape geometry with θ = 82° and 85°, respectively, and φ = 38°, while benzene-hexafluorobenzene assumes a twisted ‘V’ shape. A strong solvent effect is seen in the geometry of the benzene- methoxybenzene complex. The tilt angle is 55° when chloroform is used as a solvent and it is 82° without chloroform. Chapter 7 is the concluding chapter where the main work done in this thesis is summarized and future directions are presented.

This research investigates wave propagation in an elastic half-space with a quadratic nonlinearity in its stress-strain relationship. Different boundary conditions on the surface are considered that result in both one- and two-dimensional wave propagation problems. The goal of the research is to examine the generation of second-order frequency effects and static effects which may be used to determine the nonlinearity present in the material. This is accomplished by extracting the amplitudes of those effects in the frequency domain and analyzing their dependency on the third-order elastic constants (TOEC). For the one-dimensional problems, both analytical approximate solutions as well as numerical simulations are presented. For the two-dimensional problems, numerical solutions are presented whose dependency on the material's nonlinearity is compared to the one-dimensional problems. The numerical solutions are obtained by first formulating the problem as a hyperbolic system of conservation laws, which is then solved numerically using a semi-discrete central scheme. The numerical method is implemented using the package CentPack. In the one-dimensional cases, it is shown that the analytical and numerical solutions are in good agreement with each other, as well as how different boundary conditions may be used to measure the TOEC. In the two-dimensional cases, it is shown that there exist comparable dependencies of the second-order frequency effects and static effects on the TOEC. Finally, it is analytically and numerically investigated how multiple reflections in a plate can be used to simplify measurements of the material nonlinearity in an experiment.

The plane frame formed by precast elements have been increasingly used in structural systems of buildings throughout Brazil. These frames, as most of the precast structures deformable connections have a greater or lesser degree. The aim of this study is to analyze the structural behavior of portal frames composed of precast reinforced concrete taking into account the variation in the stiffness of the connections and the physical and geometric nonlinearity evaluated by approximate methods. The analysis of geometric nonlinearity was taken through the process P- and physical nonlinearity was checked by the Branson method. For the analysis of material and geometric nonlinearity simultaneously proposed a coupling of the two methods mentioned above. The proposed coupling was tested in a case study of a building with four floors, considering the loads recommended by NBR 6118 (ABNT, 2014) and evaluating the influence on the stiffness of two beam-column connection types. The results show an increase in the displaceability of the structure more deformable links analyzed five times at the most critical situation is compared with the evaluated frames with more rigid links. The structures analyzed taking into account the nonlinearities present in the internal structural strain variations of up to 162% compared with the linear analysis. Therefore it can be stated that the determination of the stiffness of the connections and the consideration of nonlinearities directly influenced the behavior of the efforts that will act in the structure, and therefore the structural analysis and design of the components of the same. / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Os pórticos planos formados por elementos pré-moldados vêm sendo cada vez mais utilizados em sistemas estruturais de edifícios em todo Brasil. Esses pórticos, como a maior parte das estruturas pré-moldadas, apresentam ligações deformáveis em maior ou menor grau. O objetivo geral deste trabalho é analisar o comportamento estrutural de pórticos planos formados por elementos pré-moldados de concreto armado levando em consideração a variação na rigidez das ligações e as não-linearidade física e geométrica avaliadas por métodos aproximados. A análise da não-linearidade geométrica foi feita através do processo P- e a não-linearidade física foi verificada através do método de Branson. Para a análise das não-linearidades física e geométrica de maneira simultânea foi proposto um acoplamento dos dois métodos citados anteriormente. O acoplamento proposto foi testado em um estudo de caso de um edifício com quatro pavimentos, considerando-se os carregamentos recomendados pelas normas brasileiras e avaliando-se a influência na rigidez de duas tipologias de ligação viga-pilar. Os resultados obtidos mostram um aumento na deslocabilidade da estrutura analisada com ligações mais deformáveis de cinco vezes na situação mais crítica se comparado com os pórticos avaliados com ligações mais rígidas . As estruturas analisadas levando em consideração as não-linearidades apresentaram variações nos esforços internos solicitantes de até 162% se comparadas com a análise linear. Portanto pode-se afirmar que a determinação da rigidez das ligações e a consideração das não-linearidades influem diretamente no comportamento dos esforços que irão atuar na estrutura e por consequência, na análise estrutural e dimensionamento dos elementos que compõem a mesma.

In the quest for ultrathin materials compatible with CMOS technology for all-optical signal processing applications in integrated photonics, graphene appears to be a promising candidate, with broadband1 optical properties and a high and broadband optical nonlinearity. However, researchers do not agree on the value of its nonlinear refractive index, and commonly used characterization methods do not provide a clear picture of the optical nonlinearity, in terms of its tensor nature or relaxation time. In the first part of this thesis, apart from the previously used Z-scan method, we have also used the ultrafast Optical Kerr Effect method coupled to Optical Heterodyne Detection (OHD-OKE) for the characterization of the third order optical nonlinearity of monolayer CVD graphene at telecom wavelengths. This method allows to separately measure the real and the imaginary part of the third-order nonlinearity, as well as their dynamics. With respect to the Z-scan method, OHD-OKE presents the major advantage of being robust against inhomogeneities of the sample. As such, we have demonstrated that graphene has a negative nonlinear refractive index, contrary to previously reported results. In addition, we have studied the real and imaginary part of graphene’s nonlinearity, when electrostatic gating is applied to change the chemical potential of graphene. Furthermore, we have proposed an enhanced version of the OHD-OKE method, together with the appropriate theoretical framework, in order to extract the tensor elements of the nonlinearity including the out-of-plane tensor elements. In particular, we have measured separately the time response of the two main tensor elements of the nonlinear susceptibility and we have experimentally verified that the out-of-plane tensor components are negligible. In the second part of this thesis, we have investigated, from an experimental point of view, the use of the nonlinear optical response of graphene for all-optical switching applications in integrated photonics. Namely, we have designed simple silicon nitride waveguide structures that constitute basic building blocks of switching devices, which were then fabricated and covered by graphene patches. Finally, we have experimentally tested the graphene-covered structures at low and high power levels and discussed the results. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

Cette thèse est consacrée à l'étude d'équations de réaction-diffusion dans un environnement périodique en temps. Ces équations modélisent l'évolution d'une tumeur cancéreuse en présence d'un traitement qui correspond à une immunothérapie dans la première partie du manuscrit, et à une chimiothérapie cytotoxique dans la suite.On considère dans un premier temps des nonlinéarités périodiques en temps pour lesquelles 0 et 1 sont des états d'équilibre linéairement stables. On étudie l'unicité, la monotonie et la stabilité de fronts pulsatoires. On exhibe également des cas d'existence et de non-existence de telles solutions. Dans la deuxième partie de la thèse, on commence par travailler sur des nonlinéarités périodiques en temps qui sont la somme d'une fonction positive traduisant la croissance de la tumeur et d'un terme de mort de cellules cancéreuses du au traitement. On s'intéresse aux états d'équilibres de telles nonlinéarités, et on va déduire de cette étude des propriétés de propagation de perturbations et l'existence de fronts pulsatoires. On raffine ensuite le modèle en considérant des nonlinéarités qui sont la somme d'une fonction asymptotiquement périodique en temps et d'un terme perturbatif. On prouve notamment que les propriétés relatives à la propagation de perturbations restent valables dans ce cadre là. Pour finir, on s'intéresse à l'influence du protocole de traitement. / This phD thesis investigates reaction-diffusion equations in a time periodic environment. These equations model the evolution of a cancerous tumor in the presence of a treatment that corresponds to an immunotherapy in the firs part of the manuscript, and to a cytotoxic chemotherapy after. We begin by considering time-periodic nonlinearities for which 0 and 1 are linearly stable equilibrium states. We study uniqueness, monotonicity and stability of pulsating fronts. We also provide some conditions for the existence and non-existence of such solutions.In the second part of the manuscript, we begin by working on time-periodic nonlinearities which are the sum of a positive function which stands for the growth of the tumor in the absence of treatment and of a death term of cancerous cells due to treatment. We are interested in equilibrium states of such nonlinearities, and we will infer from this study spreading properties and existence of pulsating fronts. We then refine the model by considering nonlinearities which are the sum of an asymptotic periodic nonlinearity and of a small perturbation. In particular we prove that the spreading properties remain valid in this case. To finish, we are interested in the influence of the protocol of the treatment.

Cette thèse apporte une nouvelle contribution à la caractérisation non destructive des matériaux en béton endommagés par la réaction alcalis-silice (RAS). À cette fin, des techniques et outils de caractérisation non linéaire ont été mises au point. Un banc de test de résonance non linéaire a été développé. Des améliorations ont été apportées au niveau de la chaine d'excitation et au niveau du traitement du signal afin d'optimiser la sensibilité du banc de test. Les essais non linéaires ont été effectues sur sept échantillons de béton endommagés par la RAS, trois échantillons de béton endommagés thermiquement, deux échantillons de béton endommagés mécaniquement et trois échantillons de béton sains. La non-linéarité comportementale des matériaux étant souvent attribuée au comportement hystérétique des micro-défauts contenus dans ces derniers, il a été montré dans un premier temps, que le béton endommagé par la RAS présente un comportement hystérétique. Cette étude a été faite à partir des essais de l'acousto-élastique. Le banc de test de résonance non linéaire a été ensuite utilisé pour la caractérisation des bétons sains et endommagés par la RAS. II a été montré que la technique non linéaire, en plus de permettre une caractérisation sans avoir l'historique de l'état du matériau, permet également de détecter de façon précoce l'endommagement du matériau réactif. L'influence de la teneur en eau sur les paramètres non linéaires a également été étudiée. II a été montré que les valeurs mesurées sur des échantillons de béton conservés dans des conditions de forte teneur en eau sont plus faibles. Dans l'objectif de trouver une particularité à l'endommagement causé par la RAS, la nature visqueuse du gel produit par la RAS a été utilisée. Une démarche, qui s'inspire des essais de fluage en statique réalisés sur des matériaux, a été utilisée pour répondre à cette question de signature de la RAS, tout en appliquant la technique de résonance non linéaire. Le modèle ressort-amortissement de Maxwell a été utilise pour l'interpretation des résultats. Ainsi, le temps de réponse au fluage a été analysé sur des échantillons endommagés par la RAS. II en ressort que le gel issu de la RAS rallonge le temps de réponse au fluage. Enfin, les limites de la technique de résonance non linéaire pour une application in situ ont été expliquées et une nouvelle technique non linéaire estimée applicable sur site a été initiée. Cette technique consiste à utiliser une source externe telle une masse pour provoquer la non-linéarité dans le matériau, pendant qu'une onde ultrasonore sonde le milieu.||The aim of this thesis is to contribute to the non-destructive characterization of concrete materials damaged by alkali-silica reaction (ASR). For this purpose, some nonlinear characterization techniques have been developed, as well as a nonlinear resonance test device. In order to optimize the sensitivity of the test device, the excitation module and signal processing have been improved. The nonlinear tests were conducted on seven samples of concrete damaged by ASR, three samples of concrete damaged by heat, three concrete samples damaged mechanically and three sound concrete samples. Since, nonlinear behaviour of the material is often attribute to its micro-defects hysteretic behaviour, it was shown at first that concrete damaged by ASR exhibits an hysteresis behaviour. To conduct this study, an acoustoelastic test was set, and then nonlinear resonance test device was used for characterizing sound concrete and concrete damaged by ASR. It was shown that the nonlinear technique can be used for characterizing the material without knowing its initial state, and also for detecting early damage in the reactive material. Studies were also carried out on the effect of moisture regarding the nonlinear parameters; they allowed understanding the low values of nonlinear parameters measured on concrete samples that were kept in high moisture conditions. In order to find a specific characteristic of damage caused by ASR, the viscosity of ASR gel was used. An approach, based on static creep analysis, performed on the material, while applying the nonlinear resonance technique. The spring-damping model of Maxwell was used for the interpretation of the results. Then, the creep time was analysed on samples damaged by ASR. It appears that the ASR gel increases the creep time. Finally, the limitations of the nonlinear resonance technique for in situ application have been explained and a new applicable nonlinear technique was initiated. This technique use an external source such as a mass for making non-linearity behaviour in the material, while an ultrasound wave is investigating the medium.