Waveguide bandpass filter analysis and design using multimode parallel FDTD diakoptics

A new TLM‐based concept analogous to the Green's function approach in classical electromagnetic theory is presented. It employs the procedure known as time domain diakoptics. The response of a TLM mesh to a unit impulse excitation at selected input points is interpreted as a discrete Green's function in the time domain, and the term ‘Johns matrix’ is proposed for this characteristic response. As in classical theory, the response of the mesh to an arbitrary excitation is found by convolving the excitation with its Johns matrix. This concept extends the generalized scattering parameter concept into the time dimension, opens unprecedented possibilities for partioning time domain problems at the field level, and permits large‐scale preprocessing of substructures for computer‐aided design. It also represents and elegant way of modelling broadband absorbing boundary conditions, imperfectly conducting walls and general frequency dispersive boundaries in the time domain.

A new discrete Green's function formulation of the finite difference time domain (DGF-FDTD) method has been developed which expresses the field response as a convolution of the current sources and the impulse response (which is referred to as the discrete Green's function) of the FDTD equation system. DGF-FDTD avoids the need for computation of free space nodes and absorbing boundary conditions. The method has been demonstrated by its application to a simple Yagi-Uda antenna with an efficiency comparable to the method of moments, which most of the time domain methods lack.

This is the second paper in a set that defines the discrete Green's function (DGF). This paper focuses first on the turbulent boundary layer and presents two different methods to estimate the DGF. The long-element formulation defines the DGF with just two simple algebraic equations, but it is not quantitatively accurate for short element lengths. A short element correction is derived, but must be recalculated for each selection of flow parameters and element lengths. A similarity solution is derived that allows accurate estimates of the DGF diagonal elements for laminar boundary layers and for turbulent boundary layers discretized with short element lengths. To illustrate other methods to derive DGFs in more complex flows, a low-resolution DGF for laminar stagnation line boundary layers is determined using the skin-friction formulation combined with similarity solutions for two different thermal boundary conditions. Stagnation line flow is shown to be highly sensitive to the thermal boundary condition, and this can be analyzed effectively using the DGF.

A modal absorbing boundary condition (ABC) based on the discrete Green's function (DGF) is introduced and applied for termination of waveguides simulated by means of the finite-difference time-domain (FDTD) method. The differences between the developed approach and implementations already demonstrated in the literature are presented. By applying DGF, a consistent theoretical approach to modal ABC in the FDTD method is obtained. Developed ABC is applicable to propagating and evanescent modes yielding the reflection error below -250 dB. Because the perfect matching may be computationally expensive, results of investigations of the tradeoff between computational efficiency and accuracy of the developed method are presented.

In this article, we first introduce interpolation operator of projection type in three dimensions, from which we then derive weak estimates for tensor‐product block finite elements of degree m ≥ 1. Finally, using estimates for the discrete Green's function and the discrete derivative Green's function, we prove that both of the gradient and the function value of the finite element solution uh and the corresponding interpolant Πmu of projection type and degree m are superclose in the pointwise sense of the L∞‐norm. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009

A novel method for transient analysis of Iossy transmission lines with arbitrary nonlinear terminal networks is presented. The uniqueness of this approach is that we develop time-domain Green's functions for the multiport transmission-line systems by terminating the ports in quasi-matched loads. This ensures Green's functions of a short duration. Hence, the amount of frequency-domain data necessary to obtain time-domain Green's functions is modest. These Green's functions are then convolved with the line port voltages. With this technique one can analyze responses of multiconductor transmission lines with arbitrary nonlinear loads (even with memory) as we have at any instant of time Thevenin's equivalent of the linear portion of the system. An example is presented to illustrate the application of this technique to multiconductor nonlinearly loaded transmission lines.

The heat fluxes across the turbine tip gap are characterized by large unsteady pressure gradients and shear from the viscous effects. The classical Newton heat convection equation, based on the turbine inlet total temperature, is inadequate. Previous research from our team relied on the use of the adiabatic wall temperature. In this paper, we propose an alternative approach to predict the convective heat transfer problem across the turbine rotor tip using discrete Green's functions (DGF). The linearity of the energy equation in the solid domain with constant thermal properties can be applied with a superposition technique to measure the data extracted from flow simulations to determine the Green's function distribution. The DGF is a matrix of coefficients that relate the temperature spatial (GF) distribution with the heat flux. This methodology is first applied to a backward facing step, validated using experimental data. The final aim of this paper is to demonstrate the method in the rotor turbine tip. A turbine stage at engine-like conditions was assessed using cfd software. The heat flux pulses were applied at different locations in the rotor tip geometry, and the increment of temperature in this zone was evaluated for different clearances, with a consequent variation of the DGF coefficients. Ultimately, a detailed uncertainty analysis of the methodology was included based on the magnitude of the heat flux pulses used in the DGF coefficients calculation and the uncertainty in the experimental measurements of the wall temperature.

Periodic arrays of radiating or scattering elements play an important role in phased array antennas, frequency selective surfaces and related applications. To gain an understanding of the sparsely explored time domain (TD) behavior of such structures, we have initiated a systematic investigation of relevant canonical TD dipole-excited Green's functions (GF), which so far include those for infinite and truncated line periodic arrays, parameterized in terms of TD-Floquet waves (FW) and truncation-induced TD-FW-modulated tip diffractions. Such waves on semi-infinite and finite square arrays of dipoles have been investigated in the frequency domain, and shown to be useful in practical array applications. This article extends our TD studies to an infinite periodic sequentially pulsed planar array. Like the predecessor GFs, this canonical prototype is simple enough to admit closed form exact solutions, whose interpretation is discussed phenomenologically aided by asymptotic parameterization in terms of instantaneous frequencies. Preliminary numerical results demonstrate the efficiency of the TD-FW algorithms.

In this work, an FFT-accelerated MOT solver is used to analyze broadband scattering phenomena in a half-space environment. The multilevel space-time FFT scheme accelerates both the matrix-fill and marching times by exploiting the translational invariance of the underlying integral kernel through the use of uniform spatial meshes. All time domain Green's functions are computed by inverse Fourier transforming their frequency domain counterparts, which, in turn, are computed by adaptive evaluation of Sommerfeld integrals. In what follows, the proposed method is applied to the analysis of transient scattering from PEC bodies residing above a dielectric half-space. The method is flexible enough to be used in the analysis of layered media with more complicated Green's functions and can be extended to nonuniformly meshed scatterers using the time domain adaptive integral method.

In this note a representation of the discrete Green's function of a compact discretization of a two point boundary value problem of order n 2 is given which among other things allows a direct proff of the convergence (and divergence) properties.

Time domain responses of porous media have been studied by some authors, but generally the possible descriptions have been given in the frequency domain. The aim of this paper, limited to materials with rigid skeleton considered as equivalent fluids, is to compare the descriptions by Johnson-Allard ($JA$%) as well as by Pride-Lafarge ($PL$) with i) some analytical, approximate formulas, based upon asymptotic high frequency expansion ; ii) the exact formula by Zwikker and Kosten for the case of cylindrical pores. The paper starts with a short summary of the statement of the different general full frequency models ($JA$ and $PL).$ The Green function in the time domain is shown to exhibit interesting properties of materials. In particular the maximum response depends on one dimensionless parameter only, which is denoted $\xi $ and is the ratio of the travelled distance to the product of the \textquotedblleft frozen\textquotedblright\ sound speed and a characteristic viscous relaxation time. The distance $\xi $ is related to a time domain Stokes number. The numerical computation of the Green function is done by FFT, with some precautions, because of the importance of the higher frequencies on the response shape. The $PL$ description is shown to be the best full frequency general model, but some discrepancies with the exact model appear at short times or short distances. When the distance $\xi $ increases from zero, the asymptotic expansion shows that the maximum of the Green function decreases first as $1/\xi ^{2}$, then exponentially.

In the use of the time-domain integral equation (TDIE) method for the analysis of layered media, it is important to have the time-domain layered medium Green's function computed for many source-to-field distances /spl rho/ and time instants t. In this paper, a numerical method is used that computes the mixed potential Green's functions G/sub v/(/spl rho/,t) and G/sub A/(/spl rho/,t) for a multilayered medium for many /spl rho/'s and t's simultaneously. The method is applicable to multilayered media and for lossless or lossy dispersive media. Salient features of the method are: 1) the use of complex /spl omega/ so that the surface wave poles are lifted off the real k/sub /spl rho// axis such that pole extractions are not required; 2) the use of half-space extraction so that the integrand for the Sommerfeld integral decays exponentially along the k/sub /spl rho// axis to obtain fast convergence of the integral; and 3) the use of the fast Hankel transform so that the Green's function is calculated for many values of /spl rho/ simultaneously. For a four-layer medium, we illustrate the numerical results by a three-dimensional plot of /spl rho/G/sub v/(/spl rho/,t) versus /spl rho/ and t and demonstrate the space-time evolution of these Green's functions. For a maximum frequency range of 8 GHz, the method requires only a few CPU minutes to compute a table of 100 (points in /spl rho/) /spl times/ 168 (points in t) uniformly spaced values of G/sub v/(/spl rho/,t) on an 867-MHz Pentium PC.

Multiscale modeling of failure in metal matrix composites

A speedy computation of the time-domain Green's function for microstrip structures is performed. The spatial Green's function in the frequency domain is written as the sum of real- and complex-image terms, which are respectively converted into the time domain analytically and numerically via fast Fourier transforms. The coefficients of the complex-image terms are both space and frequency independent and are computed only once, which results in a speedy algorithm.

Summary form only given. Classes of asymptotic methods for the direct computation of discrete Green's functions are considered.