In this thesis I present a new numerical scheme for 2.5-D/3-D direct current resistivity modelling in heterogeneous, anisotropic media. This method, named the ‘Gaussian quadrature grid’ (GQG) method, co-operatively combines the solution of the Variational Principle of the partial differential equation, Gaussian quadrature abscissae and local cardinal functions so that it has the main advantages of the spectral element method. The formulation shows that the GQG method is a modification of the spectral element method and does not employ the constant elements and require the mesh generator to match the earth’s surface. This makes it much easier to deal with geological models having a 2-D/3-D complex topography than using traditional numerical methods. The GQG technique can achieve a similar convergence rate to the spectral element method. It is shown that it transforms the 2.5-D/3-D resistivity modelling problem into a sparse and symmetric linear equation system, which can be solved by an iterative or matrix inversion method. Comparison with analytic solutions for homogeneous isotropic and anisotropic models shows that the error depends on the Gaussian quadrature order (abscissae number) and the sub-domain size. The higher order or smaller the subdomain size employed, the more accurate the solution. Several other synthetic examples, both homogeneous and inhomogeneous, incorporating sloping, undulating and severe topography are presented and found to yield results comparable to finite element solutions involving a dense mesh. The thesis also presents for the first time explicit expressions for the Fréchet derivatives or sensitivity functions in resistivity imaging of a heterogeneous and fully anisotropic earth. The formulation involves the Green’s functions and their gradients, and is developed both from a formal perturbation analysis and by means of a numerical (finite element) method. A critical factor in the equations is the derivative of the electrical conductivity tensor with respect to the principal conductivity values and the angles defining the axes of symmetry; these are given analytically. The Fréchet derivative expressions are given for both the 2.5-D and the 3-D problem using both constant point and constant block model parameterisations. Special cases like the isotropic earth and tilted transversely isotropic (TTI) media are shown to emerge from the general solutions. Numerical examples are presented for the various sensitivities as functions of the dip angle and strike of the plane of stratification in uniform TTI media. In addition, analytic solutions are derived for the electric potential, current density and Fréchet derivatives at any interior point within a 3-D transversely isotropic homogeneous medium having a tilted axis of symmetry. The current electrode is assumed to be on the surface of the Earth and the plane of stratification given arbitrary strike and dip. Profiles can be computed for any azimuth. The equipotentials exhibit an elliptical pattern and are not orthogonal to the current density vectors, which are strongly angle dependent. Current density reaches its maximum value in a direction parallel to the longitudinal conductivity direction. Illustrative examples of the Fréchet derivatives are given for the 2.5-D problem, in which the profile is taken perpendicular to strike. All three derivatives of the Green’s function with respect to longitudinal conductivity, transverse resistivity and dip angle of the symmetry axis (dG/dσ₁,dG/dσ₁,dG/dθ₀ ) show a strongly asymmetric pattern compared to the isotropic case. The patterns are aligned in the direction of the tilt angle. Such sensitivity patterns are useful in real time experimental design as well as in the fast inversion of resistivity data collected over an anisotropic earth. / Thesis (Ph.D.) -- University of Adelaide, School of Chemistry and Physics, 2009

Thesis (Ph. D.)--Brigham Young University. Dept. of Chemical Engineering, 2005. / Includes bibliographical references (p. 145-156).

Thesis (M.S.)--West Virginia University, 2010. / Title from document title page. Document formatted into pages; contains viii, 75 p. : ill. (some col.). Includes abstract. Includes bibliographical references (p. 72-75).

abstract: Presented is a design approach and test of a novel compact waveguide that demonstrated the outer dimensions of a rectangular waveguide through the introduction of parallel raised strips, or flanges, which run the length of the rectangular waveguide along the direction of wave propagation. A 10GHz waveguide was created with outer dimensions of a=9.0mm and b=3.6mm compared to a WR-90 rectangular waveguide with outer dimensions of a=22.86mm and b=10.16mm which the area is over 7 times the area. The first operating bandwidth for a hollow waveguide of dimensions a=9.0mm and b=3.6mm starts at 16.6GHz a 40% reduction in cutoff frequency. The prototyped and tested compact waveguide demonstrated an operating close to the predicted 2GHz with predicted vs measured injection loss generally within 0.25dB and an overall measured injection loss of approximately 4.67dB/m within the operating bandwidth. / Dissertation/Thesis / Doctoral Dissertation Electrical Engineering 2019

The sound generated by high speed trains can be exacerbated by the presence of trackside structures. Tunnels are the principal structures that have a strong influence on the noise produced by trains. A train entering a tunnel causes air to flow in and out of the tunnel portal, forming a monopole source of low frequency sound ["infrasound"] whose wavelength is large compared to the tunnel diameter. For the compact case, when the tunnel diameter is small, incompressible flow theory can be used to compute the Green's function that determines the monopole sound. However, when the infrasound is "shielded" from the far field by a large "flange" at the tunnel portal, the problem of calculating the sound produced in the far field is more complex. In this case, the monopole contribution can be calculated in a first approximation in terms of a modified Compact Green's function, whose properties are determined by the value at the center of a. disk (modelling the flange) of a diffracted potential produced by a thin circular disk. In this thesis this potential is calculated numerically. The scattering of sound by a thin circular disk is investigated using the Finite Difference Method applied to the three dimensional Helmholtz equation subject to appropriate boundary conditions on the disk. The solution is also used to examine the unsteady force acting on the disk.

The inverse source problem for the Helmholtz equation is studied. An unknown source is to be identified from the knowledge of its radiated wave. The focus is placed on the effect that multi-frequency data has on establishing uniqueness. In particular, we prove that data obtained from finitely many frequencies is not sufficient. On the other hand, if the frequency varies within an open interval of the positive real line, then the source is determined uniquely. An algorithm is based on an incomplete Fourier transform of the measured data and we establish an error estimate under certain regularity assumptions on the source function. We conclude that multi-frequency data not only leads to uniqueness for the inverse source problem, but in fact it contributes with a stability result for the reconstruction of an unknown source.

This thesis presents the Differential Quadrature Method (DQM) for solving Helmholtz, modified Helmholtz and Helmholtz eigenvalue-eigenvector equations. The equations are discretized by using Polynomial-based and Fourier-based differential quadrature technique wich use basically polynomial interpolation for the solution of differential equation.

We study two finite element methods for solving time-harmonic electromagnetic and acoustic problems: the discontinuous Petrov-Galerkin (DPG) method and the hybrid discontinuous Galerkin (HDG) method. The DPG method for the Helmholtz equation is studied using a test space normed by a modified graph norm. The modification scales one of the terms in the graph norm by an arbitrary positive scaling parameter. We find that, as the parameter approaches zero, better results are obtained, under some circumstances. A dispersion analysis on the multiple interacting stencils that form the DPG method shows that the discrete wavenumbers of the method are complex, explaining the numerically observed artificial dissipation in the computed wave approximations. Since the DPG method is a nonstandard least-squares Galerkin method, its performance is compared with a standard least-squares method having a similar stencil. We study the HDG method for complex wavenumber cases and show how the HDG stabilization parameter must be chosen in relation to the wavenumber. We show that the commonly chosen HDG stabilization parameter values can give rise to singular systems for some complex wavenumbers. However, this failure is remedied if the real part of the stabilization parameter has the opposite sign of the imaginary part of the wavenumber. For real wavenumbers, results from a dispersion analysis for the Helmholtz case are presented. An asymptotic expansion of the dispersion relation, as the number of mesh elements per wave increase, reveal values of the stabilization parameter that asymptotically minimize the HDG wavenumber errors. Finally, a dispersion analysis of the mixed hybrid Raviart-Thomas method shows that its wavenumber errors are an order smaller than those of the HDG method. We conclude by presenting some contributions to the development of software tools for using the DPG method and their application to a terahertz photonic structure. We attempt to simulate field enhancements recently observed in a novel arrangement of annular nanogaps.

In this thesis, we consider four different stochastic partial differential equations. Firstly, we study stochastic Helmholtz equation driven by an additive white noise, in a bounded convex domain with smooth boundary of Rd (d = 2, 3). And then with the help of the perfectly matched layers technique, we also consider the stochastic scattering problem of Helmholtz type. The second part of this thesis is to investigate the time harmonic case for stochastic Maxwell's equations driven by an color noise in a simple medium, and then we expand the results to the stochastic Maxwell's equations in case of dispersive media in Rd (d = 2, 3). Thirdly, we study stochastic parabolic partial differential equation driven by space-time color noise, where the domain O is a bounded domain in R2 with boundary ∂O of class C2+alpha for 0 < alpha < 1/2. In the last part, we discuss the stochastic wave equation (SWE) driven by nonlinear noise in 1D case, where the noise 626x6t W(x, t) is the space-time mixed second-order derivative of the Brownian sheet. / Many physical and engineering phenomena are modeled by partial differential equations which often contain some levels of uncertainty. The advantage of modeling using so-called stochastic partial differential equations (SPDEs) is that SPDEs are able to more fully capture interesting phenomena; it also means that the corresponding numerical analysis of the model will require new tools to model the systems, produce the solutions, and analyze the information stored within the solutions. / One of the goals of this thesis is to derive error estimates for numerical solutions of the above four kinds SPDEs. The difficulty in the error analysis in finite element methods and general numerical approximations for a SPDE is the lack of regularity of its solution. To overcome such a difficulty, we follow the approach of [4] by first discretizing the noise and then applying standard finite element methods and discontinuous Galerkin methods to the stochastic Helmholtz equation and Maxwell equations with discretized noise; standard finite element method to the stochastic parabolic equation with discretized color noise; Galerkin method to the stochastic wave equation with discretized white noise, and we obtain error estimates are comparable to the error estimates of finite difference schemes. / We shall focus on some SPDEs where randomness only affects the right-hand sides of the equations. To solve the four types of SPDEs using, for example, the Monte Carlo method, one needs many solvers for the deterministic problem with multiple right-hand sides. We present several efficient deterministic solvers such as flexible CG method and block flexible GMRES method, which are absolutely essential in computing statistical quantities. / Zhang, Kai. / Adviser: Zou Jun. / Source: Dissertation Abstracts International, Volume: 70-06, Section: B, page: 3552. / Thesis (Ph.D.)--Chinese University of Hong Kong, 2008. / Includes bibliographical references (leaves 144-155). / Electronic reproduction. Hong Kong : Chinese University of Hong Kong, [2012] System requirements: Adobe Acrobat Reader. Available via World Wide Web. / Electronic reproduction. [Ann Arbor, MI] : ProQuest Information and Learning, [200-] System requirements: Adobe Acrobat Reader. Available via World Wide Web. / Abstracts in English and Chinese. / School code: 1307.

The main focus of this thesis is on fast numerical methods, where adaptivity is an important mechanism to lowering the methods' complexity. The application of the methods are in the areas of wireless communication, antenna design, radar signature computation, noise prediction, medical ultrasonography, crystal growth, flame propagation, wave propagation, seismology, geometrical optics and image processing.   We first consider high frequency wave propagation problems with a variable speed function in one dimension, modeled by the Helmholtz equation. One significant difficulty of standard numerical methods for such problems is that the wave length is very short compared to the computational domain and many discretization points are needed to resolve the solution. The computational cost, thus grows algebraically with the frequency w. For scattering problems with impenetrable scatterer in homogeneous media, new methods have recently been derived with a provably lower cost in terms of w. In this thesis, we suggest and analyze a fast numerical method for the one dimensional Helmholtz equation with variable speed function (variable media) that is based on wave-splitting. The Helmholtz equation is split into two one-way wave equations which are then solved iteratively for a given tolerance. We show rigorously that the algorithm is convergent, and that the computational cost depends only weakly on the frequency for fixed accuracy.  We next consider interface tracking problems where the interface moves by a velocity field that does not depend on the interface itself. We derive fast adaptive  numerical methods for such problems. Adaptivity makes methods robust in the sense that they can handle a large class of problems, including problems with expanding interface and problems where the interface has corners. They are based on a multiresolution representation of the interface, i.e. the interface is represented hierarchically by wavelet vectors corresponding to increasingly detailed meshes. The complexity of standard numerical methods for interface tracking, where the interface is described by marker points, is O(N/dt), where N is the number of marker points on the interface and dt is the time step. The methods that we develop in this thesis have O(dt^(-1)log N) computational cost for the same order of accuracy in dt. In the adaptive version, the cost is O(tol^(-1/p)log N), where tol is some given tolerance and p is the order of the numerical method for ordinary differential equations that is used for time advection of the interface.   Finally, we consider time-dependent Hamilton-Jacobi equations with convex Hamiltonians. We suggest a numerical method that is computationally efficient and accurate. It is based on a reformulation of the equation as a front tracking problem, which is solved with the fast interface tracking methods together with a post-processing step.  The complexity of standard numerical methods for such problems is O(dt^(-(d+1))) in d dimensions, where dt is the time step. The complexity of our method is reduced to O(dt^(-d)|log dt|) or even to O(dt^(-d)). / <p>QC 20121116</p>