Nano scale processing of IC chips has become the prime production technique as the microelectronic industry aims towards scaling down product dimensions while increasing accuracy and performance. Accurate control of temperature and a good monitoring mechanism for thickness of the deposition layers during epitaxial growth are critical parameters influencing a good yield. The two-fold objective of this thesis is to establish the feasibility of an alternative to the current pyrometric and ellipsometric techniques to simultaneously measure temperature and thickness during wafer processing. TAP-NDE is a non-contact, non-invasive, laser-based ultrasound technique that is employed in this study to contemporarily profile the thermal and spatial characteristics of the wafer. The Gabor wavelet transform allows the wave dispersion to be unraveled and the group velocity of individual frequency components to be extracted from the experimentally acquired time waveform. The thesis illustrates the formulation of a theoretical model that is used to identify the frequencies sensitive to temperature and thickness changes. The group velocity of the corresponding frequency components is determined and their corresponding changes with respect to temperature for different thickness are analytically modeled. TAP-NDE is then used to perform an experimental analysis on Silicon wafers of different thickness to determine the maximum possible resolution of TAP-NDE towards temperature sensitivity, and to demonstrate the ability to differentiate between wafers of different deposition layer thickness at temperatures up to 600?C. Temperature resolution is demonstrated for ?10?C resolution and for ?5?C resolution; while thickness differentiation is carried out with wafers carrying 4000? and 8000? of aluminum deposition layer. The experimental group velocities of a set of selected frequency components extracted using the Gabor Wavelet time-frequency analysis as compared to their corresponding theoretical group velocities show satisfactory agreement. As a result of this work, it is seen that TAP-NDE is a suitable tool to identify and characterize thickness and temperature changes simultaneously during thermal annealing that can replace the current need for separate characterization of these two important parameters in semiconductor manufacturing.

Im Rahmen dieser Arbeit wurde ein besseres Verständnis der Kopplung der Troposphäre und der Stratosphäre in den mittleren und polaren Breiten der Nordhemisphäre (NH) auf Monatszeitskalen erzielt, die auf die Ausbreitung von quasi-stationären Wellen zurückzuführen ist. Der Schwerpunkt lag dabei auf den dynamisch aktiven Wintermonaten, welche die grösste Variabilität aufweisen. Die troposphärische Variabilität wird zum Grossteil durch bevorzugte Zirkulationsstrukturen, den Telekonnexionsmustern, bestimmt. Mittels einer rotierten EOF-Analyse der geopotenziellen Höhe in 500 hPa wurden die wichtigsten regionalen troposphärischen Telekonnexionsmuster der Nordhemisphäre berechnet. Diese lassen sich drei grossen geografischen Regionen zuordnen; dem nordatlantisch-europäischen Raum, Eurasien und dem pazifisch-nordamerikanischen Raum. <br><br> Da es sich um die stärksten troposphärischen Variabilitätsmuster handelt, wurden sie als grundlegende troposphärische Grössen herangezogen, um dynamische Zusammenhänge zwischen der troposphärischen und der stratosphärischen Zirkulation zu untersuchen. Dabei wurde anhand von instantanen und zeitverzögerten Korrelationsanalysen der troposphärischen Muster mit stratosphärischen Variablen erstmalig gezeigt, dass unterschiedliche regionale troposphärische Telekonnexionsmuster unterschiedliche Auswirkungen auf die stratosphärische Zirkulation haben. Es ergaben sich für die pazifisch-nordamerikanischen Muster signifikante instantane Korrelationen mit quasi-barotropen Musterstrukturen und für die nordatlantisch-europäischen Muster zonalsymmetrische Ringstrukturen ab 1978 mit signifikanten Korrelationswerten über tropischen und subtropischen Breiten und inversen Korrelationswerten über polaren Gebieten. <br><br> Bei einer Untersuchung des Einflusses der stratosphärischen Variabilität wurde gezeigt, dass sich die stärkste Kopplung von nordatlantisch-europäischen Telekonnexionsmustern mit der stratosphärischen Zirkulation bei einem in Richtung Europa verschobenen Polarwirbel ergibt, wodurch die signifikanten Korrelationen ab 1978 erklärt werden können. Eine zonal gemittelte und vor allem lokale Untersuchung der Wellenausbreitungsbedingungen während dieser stratosphärischen Situation zeigt, dass es zu schwächeren Windgeschwindigkeiten in der Stratosphäre im Bereich von Nordamerika und des westlichen Nordatlantiks kommt und sich dadurch die Wellenausbreitungsbedingungen in diesem geografischen Bereich für planetare Wellen verbessern. Durch die stärkere Wellenausbreitung kommt es zu einer stärkeren Wechselwirkung mit dem Polarjet, wobei dieser abgebremst wird. Diese Abbremsung führt zu einer Verstärkung der meridionalen Residualzirkulation. D. h., wenn es zu einer verstärkten Wellenanregung im Nordatlantik und über Europa kommt, ist die Reaktion der Residualzirkulation bei einem nach Europa verschobenem Polarwirbel besonders stark. <br><br> Die quasi-barotropen Korrelationsstrukturen, die sich bei den pazifisch-nordamerikanischen Mustern zeigen, weisen aufgrund von abnehmenden Störungsamplituden mit zunehmender Höhe, keiner Westwärtsneigung und einem negativen Brechungsindex im Pazifik auf verschwindende Wellen hin, die als Lösung der Wellengleichung bei negativem Brechungsindex auftreten. Dies wird durch den Polarjet, der im Bereich des Pazifiks stets sehr weit in Richtung Norden verlagert ist, verursacht. <br><br> Abschliessend wurde in dieser Arbeit untersucht, ob die gefundenen Zusammenhänge von nordatlantisch-europäischen Telekonnexionsmustern mit der stratosphärischen Zirkulation auch von einem Atmosphärenmodell wiedergegeben werden können. Dazu wurde ein transienter 40-Jahre-Klimalauf des ECHAM4.L39(DLR)/CHEM Modells mit möglichst realistischen Antrieben erstmalig auf die Kopplung der Troposphäre und der Stratosphäre analysiert. Dabei konnten sowohl die troposphärischen, als auch die stratosphärischen Variabilitätsmuster vom Modell simuliert werden. Allerdings zeigen sich in den stratosphärischen Mustern Phasenverschiebungen in den Wellenzahl-1-Strukturen und ihre Zeitreihen weisen keinen signifikanten Trend ab 1978 auf. Die Kopplung der nordatlantisch-europäischen Telekonnexionsmuster mit der stratosphärischen Zirkulation zeigt eine wesentlich schwächere Reaktion der meridionalen Residualzirkulation. Somit stellte sich heraus, dass insbesondere die stratosphärische Zirkulation im Modell starke Diskrepanzen zu den Beobachtungen zeigt, die wiederum Einfluss auf die Wellenausbreitungsbedingungen haben. Es wird damit deutlich, dass für eine richtige Wiedergabe der Wellenausbreitung und somit der Kopplung der Troposphäre und Stratosphäre die stratosphärische Zirkulation eine wichtige Rolle spielt. / Within the scope of this study a better understanding of the coupling of the troposphere and the stratosphere in the middle and polar latitudes (NH) on monthly timescales, caused by the propagation of quasi-stationary waves is improved. The approach was focused on the dynamical active winter months, including the largest variablity. <br><br> The tropospheric variability is strongly affected by preferred circulation patterns, the so called teleconnection patterns. The most important, regional, tropospheric teleconnection patterns in the Northern Hemisphere are determined by means of a rotated EOF-Analyses of the geopotential height at the 500 hPa level. They can be attributed to three geographical regions; North Atlantic/Europe, Eurasia and Pacific/North America. These strongest tropospheric variability patterns are taken as the basic tropospheric quantities to analyse the connections between the tropospheric and stratospheric circulation. By means of instantaneous and time-lagged correlation analyses, it has been shown for the first time that different regional, tropospheric teleconnection patterns have different effects on the stratospheric circulation. The Pacific/North American patterns reveal significant correlation values with quasi-barotropic structures and the North Atlantic/European patterns show significant correlations over tropical and subtropical latitudes and invers correlation values over the polar region. <br><br> The investigation of the stratospheric variability influence reveals that the strongest coupling of the North Atlantic/European teleconnection patterns with the stratospheric circulation appears during periods with a shift of the polar vortex towards Europe. The zonal averaged and particularly the local analyses of the wave propagation conditions show that weaker wind speed in the stratosphere over North America and the western part of the North Atlantic leads to improved wave propagation conditions in this geographical region. The stronger wave propagation produces a stronger interaction of the waves with the polar jet which results in enhanced wave breaking and an amplification of the residual circulation. In the case of a stronger wave forcing in the North Atlantic and over Europe these will be a stronger reaction of the residual circulation. The quasi-barotropic correlation structures, induced by the Pacific/North American patterns, are an indicator for evanescent waves because of the decreasing perturbations with increasing height, none westward declination and a negative refractive index in the Pacific. This is generated by the polar jet in the Pacific which is always shifted very far to the north. <br><br> Concluding, it was studied, whether Atmospheric General Circulation Models (AGCMs) can reproduce the detected connections of the North Atlantic/European teleconnection patterns with the stratospheric circulation. Therefore the transient model run of the interactively coupled chemistry-climate model ECHAM4.L39(DLR)/ CHEM is used for analysing the troposphere-stratosphere coupling, covering the period from 1960 to 1999. Both, the tropospheric and the stratospheric variability patterns have been simulated by the model. However the stratospheric patterns show a phase shift in the wave number 1 patterns and the time series of the wave number 1 structures do not offer a significant trend since 1978. The coupling of the North Atlantic/European teleconnection patterns with the stratospheric circulation shows a significantly weaker annular-like correlation structure. It turned out, that the stratospheric circulation particularly shows strong discrepancies to the observations which can influence the wave propagation conditions again. Therefore, the stratospheric circulation plays an important role for an accurate reproduction of the wave propagation and consequently for the coupling of the troposphere and the stratosphere.

The increased knowledge about the complexity of the physiological processes increases the demand on the analytical techniques employed to explore them. A comprehensive analysis of the entire sample content is today the most common approach to investigate the molecular interplay behind a physiological deviation. For this purpose a method that offers a number of important properties, such as speed and simplicity, high resolution and sensitivity, minimal sample volume requirements, cost efficiency and robustness, possibility of automation, high-throughput and wide application range of analysis is requested. Capillary electrophoresis (CE) coupled to mass spectrometry (MS) has a great potential and fulfils many of these criteria. However, further developments and improvements of these techniques and their combination are required to meet the challenges of complex biological samples. Protein analysis using CE is a challenging task due to protein adsorption to the negatively charged fused-silica capillary wall. This is especially emphasised with increased basicity and size of proteins and peptides. In this thesis, the adsorption problem was addressed by using an in-house developed physically adsorbed polyamine coating, named PolyE-323. The coating procedure is fast and simple that generates a coating stable over a wide pH range, 2-11. By coupling PolyE-323 modified capillaries to MS, either using electrospray ionisation (ESI) or matrix-assisted laser desorption/ionisation (MALDI), successful analysis of peptides, proteins and complex samples, such as protein digests and crude human body fluids were obtained. The possibilities of using CE-MALDI-MS/MS as a proteomic tool, combined with a proper sample preparation, are further demonstrated by applying high-abundant protein depletion in combination with a peptide derivatisation step or isoelectric focusing (IEF). These approaches were applied in profiling of the proteomes of human cerebrospinal fluid (CSF) and human follicular fluid (hFF), respectively. Finally, a multiplexed quantitative proteomic analysis was performed on a set of ventricular cerebrospinal fluid (vCSF) samples from a patient with traumatic brain injury (TBI) to follow relative changes in protein patterns during the recovery process. The results presented in this thesis confirm the potential of CE, in combination with MS, as a valuable choice in the analysis of complex biological samples and clinical applications.

This study investigates wave propagation in an elastic half-space with quadratic nonlinearity due to a line load on the surface. The consideration of this problem is one of the well known Lamb problems. Even since Lamb's original solution, numerous investigators have obtained solutions to many different variants of the Lamb problem. However, most of the solutions existing in the current literature are limited to wave propagation in a linear elastic half-space. In this work, the Lamb problem in an elastic half-space with quadratic nonlinearity is considered. For this, the problem is first formulated 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. The accuracy of the numerical method is first studied by comparing the numerical solution with the analytical solution for a half-space with linear response (the original Lamb's problem). The numerical results for the half-space with quadratic nonlinearity are than studied using signal-processing tools such as the fast Fourier transform (FFT) in order to analyze and interpret any nonlinear effects. This in particular gives the possibility to evaluate the excitation of higher order harmonics whose amplitude is used to infer material properties. To quantify and compare the nonlinearity of different materials, two parameters are introduced; these parameters are similar to the acoustical nonlinearity parameter for plane waves.

This thesis analyses the governing equation of elastic wave propagation by the finite difference method , and considered absorbing boundary condition and the material damping to simulate behavior of wave propagation. Otherwise, we combined with the mechanics of the soil pore water pressure raised by shear stress effected repeatedly and the soil property is changed by water pressure effected to simulate physical phenomenon in half-space, and probe into the soil liquefaction process during different force types. Using the developed numerical wave propagation model probe into reducing vibration by dug trench and filler trench, and analyzed data by 1/3 octave band method. This thesis discuss with reducing vibration effect by different trench disposed¡Bdifferent filler material property, complex filler, and extending the force source pile length.

Nano scale processing of IC chips has become the prime production technique as the microelectronic industry aims towards scaling down product dimensions while increasing accuracy and performance. Accurate control of temperature and a good monitoring mechanism for thickness of the deposition layers during epitaxial growth are critical parameters influencing a good yield. The two-fold objective of this thesis is to establish the feasibility of an alternative to the current pyrometric and ellipsometric techniques to simultaneously measure temperature and thickness during wafer processing. TAP-NDE is a non-contact, non-invasive, laser-based ultrasound technique that is employed in this study to contemporarily profile the thermal and spatial characteristics of the wafer. The Gabor wavelet transform allows the wave dispersion to be unraveled and the group velocity of individual frequency components to be extracted from the experimentally acquired time waveform. The thesis illustrates the formulation of a theoretical model that is used to identify the frequencies sensitive to temperature and thickness changes. The group velocity of the corresponding frequency components is determined and their corresponding changes with respect to temperature for different thickness are analytically modeled. TAP-NDE is then used to perform an experimental analysis on Silicon wafers of different thickness to determine the maximum possible resolution of TAP-NDE towards temperature sensitivity, and to demonstrate the ability to differentiate between wafers of different deposition layer thickness at temperatures up to 600?C. Temperature resolution is demonstrated for ?10?C resolution and for ?5?C resolution; while thickness differentiation is carried out with wafers carrying 4000? and 8000? of aluminum deposition layer. The experimental group velocities of a set of selected frequency components extracted using the Gabor Wavelet time-frequency analysis as compared to their corresponding theoretical group velocities show satisfactory agreement. As a result of this work, it is seen that TAP-NDE is a suitable tool to identify and characterize thickness and temperature changes simultaneously during thermal annealing that can replace the current need for separate characterization of these two important parameters in semiconductor manufacturing.

Presented in this thesis are the following theoretical investigations carried out on the non-invasive microwave hyperthermia of malignant tumours in the human body: Fundamental concepts of electromagnetic wave propagation through a biomass and its interaction with it, are discussed. Various types of applicators used for producing hyperthermia in a biomass, are also discussed. Propagation of a uniform plane electromagnetic wave through a human body is investigated for the general case of oblique incidence. Various models used for the human body have been discussed and the planar multilayer model has been chosen for this study. Reflection and transmission coefficients for both the parallel and perpendicular linear polarisations of the wave, have been determined. For normal incidence, power transfer ratio at the muscle has been defined and calculated at 433, 915 and 2450 MHz (ISM frequencies). Efects of skin thickness and also of fat thickness, on the power transfer ratio at muscle, have been studied. Effects of the thickness and dielectric constant of a bolus, and also of the dielectric constant of an initial layer, on the power transfer ratio, have been studied and their optimum values obtained at the ISM frequencies. For microwave hyperthermia, 915 MHz is recommended as the frequency of operation. Steady-state solution of the bioheat transfer equation has been obtained, assuming the biomass to be a semi-infinite homogeneous medium. Effects of various physical parameters on the temperature profile in the biomass, have been studied. Also studied is the effect of the surface temperature on the magnitude, location and the width of the temperature peak attained in the biomass. A method to determine the microwave power and the surface temperature required to produce a prescribed temperature profile in the biomass, has been developed. The transient-state solution of the bioheat transfer equation has been obtained to study the building up of the temperature profile. Procedures for the design of an open-ended rectangular metal waveguide applicator and for estimating the total microwave power requirement to produce hyperthermia in the human body, have been developed. Performance of the applicators employing linear as well as planar arrays of open-ended rectangular metal waveguide antennas, has also been studied. In order to reduce the overall physical size of the applicators, filling up of the feed waveguide with a high dielectric constant but low loss material is suggested. A simple method of obtaining the elements of the array by partitioning a large aperture by using metal walls has been adopted. Calculation of the total microwave power required by various applicators for producing hyperthermia at various depths in a biomas, have been made and a comparison of the performance of various applicators, has been presented.

This thesis endeavors towards developing various concepts employed in analysis and design of acoustic filters for varied applications ranging from combination mufflers for automobiles to complex networks of gas carrying ducts to multiply connected complex automotive silencing devices to the noise control coatings for underwater applications. A two-dimensional wave modeling approach has been proposed to evaluate sound attenuation characteristics of dissipative mufflers of finite length with/without extended inlet and outlet tubes including very large mufflers. The correctness of the method has been validated through comparison with experimental results from literature. Two other frequently used approximate schemes have been discussed briefly with reference to the available literature. These three approaches have then been weighed against each other to show the effectiveness and limitations of each one. A thorough comparison study has been performed to investigate each one’s extent of applicability. A parametric study with different parameters suggests some useful design guidelines that can be put to use while designing such mufflers. Benefits and drawbacks of reactive and dissipative mufflers have been discussed with an intention of striking a compromise between them to achieve a better transmission quality over a broad frequency range. This has been accomplished by combining these two types of mufflers/filters explicitly. These combination mufflers are analyzed using a transfer matrix based approach by extending the aforesaid concept of two-dimensional wave modeling for finite dissipative ducts. The present approach has been used to analyze axi-symmetric circular lined plenum chambers also. The effectiveness of the bulk reaction assumption to model absorptive lining is illustrated. A parametric study has been carried out to investigate the effects of different thicknesses and placements of the absorptive lining. The contributions of reflective and absorptive portion of the combination mufflerto overall attenuation performance have been investigated from the designer’s point of view A generalized algorithm has been developed for studying the plane sound wave propa- gation in a system of interconnected rigid-walled acoustic filter elements. Interconnection between various elements is represented by a connectivity matrix. Equations of volume velocity continuity and pressure equilibrium at the interconnections are generated using this connectivity matrix and are solved using the Gauss-Jordan elimination scheme to get the overall transfer matrix of the system. The algorithm used for generalized labeling of the network and computation of Transmission Loss has also been discussed. The algorithm has been applied to investigate a multiply connected automobile mufflers as a network of acoustic elements which guides the way to a specialized application discussed next. Results for some configurations have been compared with those from the FEM analysis and experiments. A parametric study with respect to some geometric variables is carried out. The acoustical similarity between apparently different networks is discussed. The approach is flexible to incorporate any other acoustic elements, provided the acoustic variables at the junctions of the element can be related by a transfer matrix a priori. Commercial automotive mufflers are often too complex to be broken into a cascade of one dimensional elements with predetermined transfer matrices. The one dimensional (1-D) scheme presented here is based on an algorithm that uses user friendly visual volume elements to generate the system equations which are then solved using a Gauss-Jordan elimination scheme to derive the overall transfer matrix of the muffler. This work attempts and succeeds to a great extent in exploiting the speed of the one dimensional analysis with the flexibility, generality and user friendliness of three dimensional analysis using geometric modeling. A code based on the developed algorithm has been employed to demonstrate the generality of the proposed method in analyzing commercial muffers by considering three very diverse classes of mufflers with different kinds of combinations of reactive, perforated and absorptive elements. Though the examples presented in the thesis are not very complex for they are meant to be just representative cases of certain classes of mufflers, yet the algorithm can handle a large domain of commercial mufflers of high degree of complexity. Results from the present algorithm have been validated through comparisons with both the analytical and the more general, three-dimensional FEM based results. The forte of the proposed method is its power to construct the system matrix consistent with the boundary conditions from the geometrical model to evaluate the four pole parameters of the entire muffer and thence its transmission loss,etc. Thus, the algorithm can be used in conjunction with the transfer matrix based muffler programs to analyze the entire exhaust system of an automobile. A different kind of acoustic filter than the above mentioned cases is then taken up for investigation. These refer to the specialized underwater acoustic filters laid as linings on submerged bodies. These kind of underwater noise control linings have three different types of objectives, namely, Echo Reduction, Transmission Reduction (TL maximization) and a combination thereof. These coatings have been shown to be behaving very differently with different shape, size and number of air channels present in the layer. In this regard, a finite element model based methodology has been followed. An hybrid type finite element based on the Pian and Tong formulation has been modified and used so as to make the computational efforts less demanding as compared to the original one. The developed finite element has been shown to be immune to the difficulties that arise due to the near incompressible characteristics of the viscoelastic materials used and the high distortion of the elements of the FE mesh. The adequacy of this formulation has been shown by comparing its results with the analytical, FE based, and experimental results. Then, this methodology has been used to analyze and generate design curves to control various geometrical parameters for proper designing of these linings. Different unit cell representations for different types of distributions of air cavities on the linings have been discussed. Four different types of layers have been introduced and analyzed to address different objectives mentioned above. They have been termed as the Anechoic layer, Insulation layer and Combination Layer of coupled and decoupled type in this thesis. The first two layers have been designed to achieve very dissimilar characteristics and the next two layers have been designed to balance their disparities. A thorough parametric study has been carried out on the geometrical parameters of all the layers to come up with the design guidelines. For anechoic and the insulation layers, different distributions have been analyzed with different unit cell geometries and their usability in specific situations has been outlined. Effect of static pressure has also been studied by using an approximate finite element method. This method can be used to simulate deep-sea testing environment.

Wave propagation is a common phenomenon in aircraft structures resulting from high velocity transient loadings like bird hit, gust etc. Apart from understanding the behavior of structures under such loading, wave propagation analysis is also important to gain knowledge about their high frequency characteristics, which have several applications. The applications include structural health monitoring using diagnostic waves and control of wave transmission for reduction of noise and vibration. Transient loadings with high frequency content are associated with wave propagation. As a result, the higher modes of the structure participate in the response. Finite element (FE) modeling for such problem requires very fine mesh to capture these higher modes. This leads to large system size and hence large computational cost. Wave propagation problems are usually solved in frequency domain using fast Fourier transform (FFT) and spectral finite element method is one such technique which follows FE procedure in the transformed frequency domain. In this thesis, a novel wavelet based spectral finite element (WSFE) is developed for wave propagation analysis in finite dimension structures. In WSFE for 1-D waveguides, the partial differential wave equations are reduced to a set of ODEs using orthogonal compactly supported Daubechies scaling functions for temporal approximation. The localized nature of the Daubechies basis functions allows finite domain analysis and imposition of the boundary conditions. The reduced ODEs are usually solved exactly, the solution of which gives the dynamic shape functions. The interpolating functions used here are exact solution of the governing differential equation and hence, the exact elemental dynamic stiffness matrix is derived. Thus, In the absence of any discontinuities, one element is sufficient to model 1-D waveguide of any length. This elemental stiffness matrix can be assembled to obtain the global matrix as in FE and after solution, the time domain responses are obtained using the inverse wavelet transform. The developed technique circumvents several serious limitations of the conventional FFT based Spectral Finite Element (FSFE). In FSFE, the wave equations are reduced to ODEs using FFT for time approximation. The remaining part of the formulation is quite similar to that of WSFE. The required assumption of periodicity in FSFE, however, does not allow modeling of finite length structures. It results in “wrap around” problem, which distorts the response simulated using FSFE and a semi-infinite (“throw-off”) element is required for imparting artificial damping. This artificial damping occurs as the “throw off” element allows leakage of energy. In some cases, a very high damping can also be considered instead of “throw off” element to remove wrap around effects. In either cases, the damping introduced is much larger than any inherent damping that may be present in the structure. It should also be mentioned that even in presence of the artificial damping, a larger time window is required for removing the distortions completely. The developed WSFE method is completely free from such problems and can efficiently handle undamped finite length structures irrespective of the time window considered. Apart from this, FSFE allows imposition of only zero initial condition and in contrary any initial conditions can be used in WSFE. Though FSFE has problem in modeling finite length undamped structures for time domain analysis, it is well suited for performing frequency domain study of wave characteristics, namely, the determination of spectrum and dispersion relations. WSFE is also capable of extracting these frequency dependent wave properties, however only up to a certain fraction of the Nyquist frequency. This constraint results from the loss in frequency resolution due to the increase in time resolution in wavelet analysis, where the basis functions are bounded both in time and frequency. A price has to be paid in frequency domain in order to obtain a bound in the time domain. The consequence of this analysis is to impose a constraint on the time sampling rate for the simulation with WSFE, to avoid spurious dispersion. WSFE for 2-D waveguides are formulated using Daubechies scaling functions for both temporal and spatial approximations. The initial and boundary conditions, however, are imposed using two different methods, which are wavelet extrapolation technique and periodic extension or restraint matrix respectively. The 2-D WSFE is bounded in both the spatial directions unlike 2-D FSFE, which is essentially unbounded in one spatial direction. Apart from this, 2-D WSFE is also free from “wrap around” problem similar to 1-D WSFE due to the localized nature of the basis functions used for temporal approximation. In this thesis, WSFE is developed for isotropic 1-D and 2-D waveguides for time and frequency domain analysis. These include elementary rod, Euler-Bernoulli and Timoshenko beams in 1-D modeling, and plates and axisymmetric cylinders in 2-D modeling. The wave propagation responses simulated using WSFE for these waveguides are validated using FE results. The advantages of the proposed technique over the corresponding FSFE method are also highlighted all through the numerical examples. Next part of the thesis involves the extension of the developed WSFE technique for modeling composite and nano-composite structures to study their wave propagation behavior. Due to their anisotropic nature, analysis of composite structures, particularly high frequency transient analysis is much more complicated compared to the corresponding metallic structures. This is due to the presence of stiffness coupling in these structures. Superior mechanical properties of composites, however, are making them integral parts of an aircraft and thus they often experience such short duration, high velocity impact Loadings. Very few literatures report the response of composite structures subjected to such high frequency excitations. Here, WSFE is formulated for a higher order composite beam with axial, flexural, shear and contractional degrees of freedom. WSFE is also formulated for composite plates using classical laminated plate theory with axial and flexural degrees of freedom. Simulations performed using these WSFE models are used to study the higher order and elastic coupling effects on the wave propagation responses. Carbon nanotubes (CNTs) and their composites are attracting a great deal of experimental and theoretical research world-wide. The recent trend in the literature shows a great interest in the dynamic and wave characteristics of CNTs and nano-composites because of their several applications. In most of these applications, CNTs are used in the embedded form as it does not requires precise alignment of the nano-tubes. In addition, the extraordinary mechanical properties of CNTs are being exploited to achieve high strength nano-composite. Apart from the experimental studies and atomistic simulation to study the mechanical properties of CNTs and nano-composites, continuum modeling is also receiving much attention, mainly due to its computational viability. In this thesis, a 1-D WSFE is formulated for multi-wall carbon nanotube (MWNT) embedded composite modeled as beam using higher order layer-wise theory. This theory allows to model partial interfacial shear stress transfer, which normally occurs due to improper dispersion of CNTs in nano-composites. The effects of different matrix materials and fraction of shear stress transfer on the wave characteristics are studied. The responses obtained using other beam theories are also compared. The beam modeling does not allow capturing the radial motions of the CNT, which are important for several applications. These can be effectively captured by modeling the CNT using a 2-D axisymmetric model. Hence, a 2-D WSFE model is constructed to capture the high frequency characteristics of single-walled carbon nanotubes (SWNTs). The response of SWNT simulated using the developed model is validated with experimental and atomistic simulation results reported in the literature. The comparison are done for dispersion relation and also radial breathing mode frequencies. The effects of geometrical parameters, namely the radius and the wall thickness of the SWNT on the higher radial, longitudinal and coupled radial-longitudinal vibrational modes are analyzed. These behaviors are studied in both time and frequency domains. Such time domain analyses of finite length SWNT are not possible with the Fourier transform based techniques reported in literature, although, such analyses are important particularly for sensor applications of SWNT. Spectral finite element method is very much suited for solution of inverse problems like force reconstruction from the measured wave response. This is because the technique is based on the concept of transfer function between the displacements (output) and applied forces (input). In the present work, WSFE is implemented for identification of impact force from the wave propagation responses simulated with FE and used as surrogate experimental results. The results show that WSFE can accurately reconstruct the impulse load applied to 1-D waveguides which include rod, Euler-Bernoulli beam and connected 2-D frame, even with highly truncated response. This is unlike FSFE, where the accuracy of the identified force depends largely on the time window of the measured responses. The detection of damage from the wave propagation analysis is another class of inverse problems considered in this thesis and is of utmost importance in the area of aircraft structural health monitoring. Here, the detection scheme is based on arrival time of the waves reflected from the damage. A novel detection technique based on wavelet filtering is proposed here and it is shown to work efficiently even in the presence of noise in the measured wave responses. Detection of damage requires an efficient damage model to simulate the mode of structural failure. In this regard, two spectrally formulated wavelet elements are proposed, one to model isotropic beam with through-width notch and the second to model composite beam with embedded de-lamination. In the first case, the response of the damaged beam is considered as the perturbation of the undamaged response and the linear perturbation analysis leads to a completely new set of dynamic stiffness matrix. In the second case, the delamination is modeled by subdividing the de-laminated region into separate waveguides and full damage model is established by imposing the kinematics. These models help to simulate wave propagation in such damaged beams to study the effect of damage on the wave response. Noise and vibration are often transmitted from the source to the other parts of the structure in the form of wave propagation. Thus, control of such wave transmission is essential for reduction of noise and vibration, which are the main cause of discomfort and in many cases cause failure of structure. Here, techniques for both passive and active controls of wave are proposed. For active control, a closed loop system is modeled using WSFE with magnetostrictive actuator for control of axial and flexural wave propagations in connected isotropic 1-D waveguides. The feedback is negative velocity and/or acceleration measured at different sensor points. A very new application of CNT reinforced composite for passive control of vibration and wave response is explored in this thesis. For this, a novel concept of nano-composite inserts is proposed. This insert can be made from CNTs dispersed in polymer. The high stiffness of the inserts helps to regulate the power flow in the form of wave propagation from the point of application of the loads to other parts of the structures. The length of the insert, volume fraction of CNTs and position are changed to achieve the required reduction in wave amplitudes. The entire thesis is split up into eight chapters. Chapter 1 presents a brief introduction, the motivation and objective of the thesis. Chapters 2 and 3 give a detail account of wavelet spectral finite element formulation for 1-D and 2-D isotropic waveguides, while Chapter 4 gives the same for composite waveguides. Chapter 5 brings out essential wave characteristics in carbon nanotubes and nano-composite structures, while Chapters 6 and 7 exclusively deal with application of WSFE to some real world problems. The thesis ends with summary and directions of future research. In summary, the thesis has brought out several new aspects of wave propagation in isotropic, composite and nano-composite structures. In addition to establishing wavelet spectral finite element as a useful tool for wave propagation analysis, several new techniques are presented, several new algorithm are proposed and several new concepts are explored.

We are concerned with the high-fidelity subsurface imaging of the soil, which commonly arises in geotechnical site characterization and geophysical explorations. Specifically, we attempt to image the spatial distribution of the Lame parameters in semi-infinite, three-dimensional, arbitrarily heterogeneous formations, using surficial measurements of the soil's response to probing elastic waves. We use the complete waveforms of the medium's response to drive the inverse problem. Specifically, we use a partial-differential-equation (PDE)-constrained optimization approach, directly in the time-domain, to minimize the misfit between the observed response of the medium at select measurement locations, and a computed response corresponding to a trial distribution of the Lame parameters. We discuss strategies that lend algorithmic robustness to the proposed inversion schemes. To limit the computational domain to the size of interest, we employ perfectly-matched-layers (PMLs). The PML is a buffer zone that surrounds the domain of interest, and enforces the decay of outgoing waves. In order to resolve the forward problem, we present a hybrid finite element approach, where a displacement-stress formulation for the PML is coupled to a standard displacement-only formulation for the interior domain, thus leading to a computationally cost-efficient scheme. We discuss several time-integration schemes, including an explicit Runge-Kutta scheme, which is well-suited for large-scale problems on parallel computers. We report numerical results demonstrating stability and efficacy of the forward wave solver, and also provide examples attesting to the successful reconstruction of the two Lame parameters for both smooth and sharp profiles, using synthetic records. We also report the details of two field experiments, whose records we subsequently used to drive the developed inversion algorithms in order to characterize the sites where the field experiments took place. We contrast the full-waveform-based inverted site profile against a profile obtained using the Spectral-Analysis-of-Surface-Waves (SASW) method, in an attempt to compare our methodology against a widely used concurrent inversion approach. We also compare the inverted profiles, at select locations, with the results of independently performed, invasive, Cone Penetrometer Tests (CPTs). Overall, whether exercised by synthetic or by physical data, the full-waveform inversion method we discuss herein appears quite promising for the robust subsurface imaging of near-surface deposits in support of geotechnical site characterization investigations.