Simulating wave propagation in microstrip lines with Gyrotropic magnetic substrate is considered in this thesis. Since the static internal field distribution has an important effect on the device behavior, accurate determination of the internal fields are considered as well. To avoid the losses at microwave frequencies it is assumed that the magnetic substrate is saturated in the direction of local internal field. An iterative method to obtain the magnetization distribution has been developed. It is applied to a variety of nonlinear nonuniform magnetic material configurations that one may encounter in the design stage, subject to a nonuniform applied field. One of the main characteristics of the proposed iterative method to obtain the static internal field is that the results are supported by a uniqueness theorem in magnetostatics. The series of solutions Mn,Hn, where n is the iteration number, minimize the free Gibbs energy G(M) in sequence. They also satisfy the constitutive equation M = χH at the end of each iteration better than the previous one. Therefore based on the given uniqueness theorem, the unique stable equilibrium state M is determined. To simulate wave propagation in the Gyrotropic magnetic media a new FDTD formulation is proposed. The proposed formulation considers the static part of the electromagnetic field, obtained by using the iterative approach, as parameters and updates the dynamic parts in time. It solves the Landau-Lifshitz-Gilbert equation in consistency with Maxwell’s equations in time domain. The stability of the initial static field distribution ensures that the superposition of the time varying parts due to the propagating wave will not destabilize the code. Resonances in a cavity filled with YIG are obtained. Wave propagation through a microstrip line with YIG substrate is simulated. Magnetization oscillations around local internal field are visualized. It is proved that the excitation of magnetization precession which is accompanied by the excitation of magnetostatic waves is responsible for the gap in the scattering parameter S12. Key characteristics of the wide microstrip lines are verified in a full wave FDTD simulation. These characteristics are utilized in a variety of nonreciprocal devices like edgemode isolators and phase shifters.

In the literature, many methods are proposed to value American options. However, due to computational difficulty, there are only approximate solution or numerical method to evaluate American options. It is not easy for general investors either to understand nor to apply. In this thesis, we build up an option pricing and virtual asset model system, which provides a friendly environment for general public to calculate early exercise boundary of an American option. This system modularize the well-handled pricing models to provide the investors an easy way to value American options without learning difficult financial theories. The system consists two parts: the first one is an option pricing system, the other one is an asset model simulation system. The option pricing system provides various option pricing methods to the users; the virtual asset model system generates virtual asset prices for different underlying models.

With the trends of fast edge rates, high clock frequencies, and low voltage levels for the high-speed digital computer systems, the ground bounce noise (GBN) or simultaneously switching noise (SSN) on the power/ground planes is becoming one of the major challenges for designing the high-speed circuits. In order to analyze the impact of the GBN on signal integrity (SI) and electromagnetic interference (EMI), an accurate and efficient modeling approach that considers the active devices and passive interconnects is required. This thesis focuses on two points. One is developing modeling approaches for analyzing the GBN effects, and the other is proposing solutions to reduce it. First, based on the FDTD algorithm several efficient modeling approaches including equivalent current-source method (ECSM), Kirchoff surface integral representation (KSIR), and slot-corrected 2D-FDTD are developed. After that, a power/ground-planes design for efficiently eliminating the GBN in high-speed digital circuits is proposed by using low-period coplanar electromagnetic bandgap (LPC-EBG) structure. Its extinctive behaviors of low radiation and broadband suppression of the GBN is demonstrated numerically and experimentally. Good agreements are seen.

The dissertation focuses on the analysis and design for the new fiber-optic passive components based on the photonic-crystal-fiber (PCF). The vector boundary element method (VBEM) and the finite-difference time-domain (FDTD) method are employed to the propagation characteristics of PCF components. A novel octagonal microstructured fiber (OMF) with eight air-holes in the first ring has been proposed. The OMF has significantly wider wavelength range for single-mode operation, more circular-like field distribution, and less confinement loss. In addition, a novel compact polarization beam splitter (PBS) based on the twin-elliptical-core PCF (TEC-PCF) has also been proposed. It behaves with high extinction ration and broad bandwidth with significantly short splitter length. The design concept and the coupling mechanism are presented in this dissertation based on the normal-mode coupling theory and VBEM.

First, one of the purposes of this thesis is to estimate the complete crosstalk effects including the package and the pads on the surface acoustic wave (SAW) substrate. A new approach based on finite-difference time-domain (FDTD) with equivalent current source method is applied. Two kinds of patterns of one-port SAW resonators with the same package structure and inter-digital transducer (IDT) design are studied. Verification with the measurement results shows that our method is able to obtain good agreement and be used to observe the influence from the SAW pattern. Second, the equivalent current source method is extended to model the excitation of human-body¡¦s electrostatic discharge (ESD) situations. The efficiencies of sacrificial electrodes are also discussed. Finally, a novel sacrificial electrode with fractal to protect SAW devices from ESD break is proposed. Comparing with traditional electrode, the simulation results show that fractal can improve the protective efficiency greatly. Finally, a novel analysis model that can be used to analyze and optimize the impedance of an RFID transponder integrated circuit (IC) which uses backscatter encoding based on simultaneously maintaining the BER of the reader and maximizing the received power of the transponder IC is proposed. The analysis method utilizes mapping from signal constellation of the backscattered signal to the Smith chart to relate the two parameters. Given the system specification and characteristics of the reader and transponder antennas, the optimum impedances of transponder IC for binary communication system can be easily designed by using this model.

This thesis focuses on the modeling and solutions of the simultaneous switching noise (SSN) problems in the power delivery networks (PDN) of high-speed digital circuit packages. An efficient numerical approach based on two-dimension (2D) finite-difference time-domain (FDTD) method combined with the lumped circuit model of the interconnection is proposed to model the PDN of a package and PCB. Based on this approach, the mechanism of noise coupling between package and PCB can be analyzed. In addition, a novel photonic crystal power layer (PCPL) design for the PDN of the package or PCB is proposed to suppress the SSN. The periodic High-Dk material is embedded into the substrate layer between the power and ground planes. Both modeling and measurement demonstrate the PCPL can form a wide stopband well with excellent suppression of the SSN propagation in the substrate and the corresponding electromagnetic interference (EMI).

Contactless electrical impedance imaging technique is based upon the measurement of secondary electromagnetic fields caused by induced currents inside the body. In this study, a circular single-coil is used as a transmitter and a receiver. The purpose of this study is twofold: (1) to solve the induced current density distribution inside the realistic head model resulting from a sinusoidal excitation, (2) to calculate the impedance change of the same coil from the induced current distribution inside the head model. The Finite Difference Method is used to solve the induced current density in the head. The realistic head model is formed by seven tissues with a 1 mm resolution. The electrical properties of the model are assigned as a function of frequency. The quasi-stationary assumptions, especially for head tissues, are explored. It is shown that, numerical solution of only the scalar potential is sufficient to obtain the induced current density in the head below 10 MHz operating frequency. This simplification not only reduce the excessive size of the solution domain, but also reduces the number of unknowns by a factor of 4. For higher frequencies (depending on the application) induction and propagation effects become important. Additionally it is observed that dynamic monitoring of hemorrhage at any frequency seems feasible. It is concluded that the methodology provides useful information about the electrical properties of the human head via contactless measurements and has a potent as a new imaging modality for different clinical applications.

In this thesis, the Generalized Finite Difference (GFD) method is applied to analyze the elastoplastic deformation behavior of a long functionally graded (FGM) tube subjected to internal pressure. First, the method is explained in detail by considering the elastic response of a rotating FGM tube. Then, the pressurized tube problem is treated. A long FGM tube with fixed ends (axially constrained ends) is taken into consideration. The two cases in which the modulus of elasticity only and both the modulus of elasticity and the yield limit are graded properties are analyzed. The plastic model here is based on incremental theory of plasticity, Tresca&#039 / s yield criterion and its associated flow rule. The numerical results are compared to those of analytical ones. Furthermore, the elastic response of an FGM tube with free ends is studied considering graded modulus of elasticity and Poisson&#039 / s ratio. The results of these computations are compared to those of Shooting solutions. In the light of analyses and comparisons stated above, the applicability of the GFD method to the solution of similar problems is discussed. It is observed that, in purely elastic deformations the accuracy of the method is sufficient. However, in case of elastic-plastic deformations, the discrepancies between numerical and analytical results may increase in determining plastic displacements. It is also noteworthy that the predictions for tubes with two graded properties, i. e. the modulus of elasticity and the yield limit, turn out to be better than those with one graded property in this regard.

Accuracy of the finite difference sensitivity calculations are improved by calculating the optimum finite difference interval sizes. In an aerodynamic inverse design algorithm, a compressor cascade geometry is perturbed by shape functions and finite differences sensitivity derivatives of the flow variables are calculated with respect to the base geometry flow variables. Sensitivity derivatives are used in an optimization code and a new airfoil is designed verifying given design characteristics. Accurate sensitivities are needed for optimization process. In order to find the optimum finite difference interval size, a method is investigated. Convergence error estimation techniques in iterative solutions and second derivative estimations are investigated to facilitate this method. For validation of the method, analytical sensitivity calculations of Euler equations are used and several applications are performed. Efficiency of the finite difference sensitivity calculations is improved by parallel computing. Finite difference sensitivity calculations are independent tasks in an inverse aerodynamic design algorithm and can be computed separately. Sensitivity calculations are performed on parallel processors and computing time is decreased.

Slope stability is an important aspect of geotechnical engineering. Input parameters for the analysis are the governing factors and they must be determined accurately and precisely. Field investigations, laboratory testing and back analyses are vital instruments for the input parameters. This study presents the results of slope stability analysis for the soil slopes at Elbistan-&Ccedil / &ouml / llolar lignite mine. After executing the drilling programme, samples taken from the drilling work, delivered to soil mechanics laboratory for testing. The basic input parameters, namely cohesion and friction angle determined at soil mechanics laboratory were compared to the parameters obtained from back analysis of a large scale slope failure. Input parameters for the analysis are determined by this way. After determining the input parameters, slope stability analyses were carried out both for the permanent and temporary slopes in AfSin-Elbistan lignite basin, &Ccedil / &ouml / llolar sector. The effect of ground water on the stability of slopes was investigated in detail and maximum safe slope angles were determined for different water levels. For limit equilibrium analysis, Rocscience SLIDE software, for finite difference analysis in 3-D, Itasca FLAC3D was used. In the limit equilibrium analyses both circular and composite failures were considered. Shear strength reduction method is used for the finite difference method. The results between limit equilibrium and 3-D finite difference methods were compared. When the failure surfaces obtained from the finite difference analyses were imposed to limit equilibrium analysis, computations are resulted in lower factor of safety values for limit equilibrium analysis.