Higher Order Bases in a 2D Hybrid BEM/FEM Formulation

The advantages of using higher order, interpolatory basis functions are examined in the analysis of transverse electric (TE) plane wave scattering by homogeneous, dielectric cylinders. A boundary-element/finite-element (BEM/FEM) hybrid formulation is employed in which the interior dielectric region is modeled with the vector Helmholtz equation, and a radiation boundary condition is supplied by an Electric Field Integral Equation (EFIE). An efficient method of handling the singular self-term arising in the EFIE is presented. The iterative solution of the partially dense system of equations is obtained using the Quasi-Minimal Residual (QMR) algorithm with an Incomplete LU Threshold (ILUT) preconditioner. Numerical results are shown for the case of an incident wave impinging upon a square dielectric cylinder. The convergence of the solution is shown versus the number of unknowns as a function of the completeness order of the basis functions

Evaluating the Gradient of the Thin Wire Kernel

Recently, a formulation for evaluating the thin wire kernel was developed that employed a change of variable to smooth the kernel integrand, canceling the singularity in the integrand. Hence, the typical expansion of the wire kernel in a series for use in the potential integrals is avoided. The new expression for the kernel is exact and may be used directly to determine the gradient of the wire kernel, which consists of components that are parallel and radial to the wire axis

Integrating the Gradient of the Thin Wire Kernel

A formulation for integrating the gradient of the thin wire kernel is presented. This approach employs a new expression for the gradient of the thin wire kernel derived from a recent technique for numerically evaluating the exact thin wire kernel. This approach should provide essentially arbitrary accuracy and may be used with higher-order elements and basis functions using the procedure described in [4].When the source and observation points are close, the potential integrals over wire segments involving the wire kernel are split into parts to handle the singular behavior of the integrand [1]. The singularity characteristics of the gradient of the wire kernel are different than those of the wire kernel, and the axial and radial components have different singularities. The characteristics of the gradient of the wire kernel are discussed in [2]. To evaluate the near electric and magnetic fields of a wire, the integration of the gradient of the wire kernel needs to be calculated over the source wire. Since the vector bases for current have constant direction on linear wire segments, these integrals reduce to integrals of the for

Refinement of Methods for Evaluation of Near-Hypersingular Integrals in BEM Formulations

In this paper, we present advances in singularity cancellation techniques applied to integrals in BEM formulations that are nearly hypersingular. Significant advances have been made recently in singularity cancellation techniques applied to 1 R type kernels [M. Khayat, D. Wilton, IEEE Trans. Antennas and Prop., 53, pp. 3180-3190, 2005], as well as to the gradients of these kernels [P. Fink, D. Wilton, and M. Khayat, Proc. ICEAA, pp. 861-864, Torino, Italy, 2005] on curved subdomains. In these approaches, the source triangle is divided into three tangent subtriangles with a common vertex at the normal projection of the observation point onto the source element or the extended surface containing it. The geometry of a typical tangent subtriangle and its local rectangular coordinate system with origin at the projected observation point is shown in Fig. 1. Whereas singularity cancellation techniques for 1 R type kernels are now nearing maturity, the efficient handling of near-hypersingular kernels still needs attention. For example, in the gradient reference above, techniques are presented for computing the normal component of the gradient relative to the plane containing the tangent subtriangle. These techniques, summarized in the transformations in Table 1, are applied at the sub-triangle level and correspond particularly to the case in which the normal projection of the observation point lies within the boundary of the source element. They are found to be highly efficient as z approaches zero. Here, we extend the approach to cover two instances not previously addressed. First, we consider the case in which the normal projection of the observation point lies external to the source element. For such cases, we find that simple modifications to the transformations of Table 1 permit significant savings in computational cost. Second, we present techniques that permit accurate computation of the tangential components of the gradient; i.e., tangent to the plane containing the source element

Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals

Recently, significant progress has been made in the handling of singular and nearly-singular potential integrals that commonly arise in the Boundary Element Method (BEM). To facilitate object-oriented programming and handling of higher order basis functions, cancellation techniques are favored over techniques involving singularity subtraction. However, gradients of the Newton-type potentials, which produce hypersingular kernels, are also frequently required in BEM formulations. As is the case with the potentials, treatment of the near-hypersingular integrals has proven more challenging than treating the limiting case in which the observation point approaches the surface. Historically, numerical evaluation of these near-hypersingularities has often involved a two-step procedure: a singularity subtraction to reduce the order of the singularity, followed by a boundary contour integral evaluation of the extracted part. Since this evaluation necessarily links basis function, Green s function, and the integration domain (element shape), the approach ill fits object-oriented programming concepts. Thus, there is a need for cancellation-type techniques for efficient numerical evaluation of the gradient of the potential. Progress in the development of efficient cancellation-type procedures for the gradient potentials was recently presented. To the extent possible, a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. However, since the gradient kernel involves singularities of different orders, we also require that the transformation leaves remaining terms that are analytic. The terms "normal" and "tangential" are used herein with reference to the source element. Also, since computational formulations often involve the numerical evaluation of both potentials and their gradients, it is highly desirable that a single integration procedure efficiently handles both

Controls on the transport of oceanic heat to Kangerdlugssuaq Glacier, East Greenland

Greenland's marine-terminating glaciers may be sensitive to oceanic heat, but the fjord processes controlling delivery of this heat to glacier termini remain poorly constrained. Here we use a three-dimensional numerical model of Kangerdlugssuaq Fjord, East Greenland, to examine controls on fjord/shelf exchange. We find that shelf-forced intermediary circulation can replace up to ~25% of the fjord volume with shelf waters within 10 d, while buoyancy-driven circulation (forced by subglacial runoff from marine-terminating glaciers) exchanges ~10% of the fjord volume over a 10 d period under typical summer conditions. However, while the intermediary circulation generates higher exchange rates between the fjord and shelf, the buoyancy-driven circulation is consistent over time hence more efficient at transporting water along the full length of the fjord. We thus find that buoyancy-driven circulation is the primary conveyor of oceanic heat to glaciers during the melt season. Intermediary circulation will however dominate during winter unless there is sufficient input of fresh water from subglacial melting. Our findings suggest that increasing shelf water temperatures and stronger buoyancy-driven circulation caused the heat available for melting at Kangerdlugssuaq Glacier to increase by ~50% between 1993–2001 and 2002–11, broadly coincident with the onset of rapid retreat at this glacier

Computation of layered mixed potentials for the accurate and efficient analysis of periodic printed structures

International audienceOriginal acceleration procedures are proposed for the efficient calculation of the vertical components of the dyadic and scalar mixed-potential layered-media periodic Green's functions for various types of periodic structures. The extraction of suitable asymptotic terms, i.e., quasi-static images, is performed in order to speed up the convergence of the relevant spectral series. The extracted terms can be expressed as potentials for array of half-plane and half-line sources, depending on the type of the considered periodic Green's function. The relevant numerical results show the remarkable improvements in the efficiency of the approach

An investigation of conformable antennas for the astronaut backpack communication system

During periods of extravehicular activity it is obviously important that communication and telemetry systems continue to function independently of the astronaut. A system of antennas must therefore be designed that will provide the necessary isotropic coverage using circular polarization over both the transmit and receive frequency bands. To avoid the inherent physical limitations to motion that would be incurred with any sort of protruding antenna, it is necessary that the radiator be essentially flush-mounted or conformable to the structure on which it is attached. Several individual antenna elements are needed for the desired coverage. Both the particular elements chosen and their location determine the ultimate radiation pattern of the overall system. For these reasons a two-fold research plan was undertaken. First, individual elements were investigated and designed. Then various mounting locations were considered and the radiation patterns were predicted taking into account the effects of the astronaut's backpack

Evaluation of Static Potential Integrals on Triangular Domains

Static potential integrals for constant and linear sources on triangles are derived in a straightforward way. The new representations, as presented, are robust with respect to machine evaluation in important limiting cases. The potential integrals comprise up to six functions, each dependent on the relative position and orientation (with respect to an observation point) of a vertex and edge, respectively, of the source triangle. Gradients of the potentials are derived by differentiation, thus preserving relations between the representations. Each such vertex function reveals any anomalous functional behavior near its associated vertex or edge, which is useful information for devising test integral schemes. Potential plots in the source plane of an equilateral triangle illustrate such behavior, as do similar plots for each vertex function and gradient components near their associated edge and vertex

Pennsylvania Folklife Vol. 34, No. 4

• Coverlets • Sign Painting • Reverse Painting on Glass • Kites • Snake Lore • Horncraft • Weathervanes and Country Signs • Festival Focus • Sheep Shearing & Natural Knits • Bread Baking Among the Pennsylvania Dutch • The Craft of Rushing • Toy Soldier Casting • Pennsylvania Dutch Humor • Fireside Brooms and Whirligigs • Springerlehttps://digitalcommons.ursinus.edu/pafolklifemag/1108/thumbnail.jp