Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.

Isotopic data of atmospheric methane from EDML and Vostok ice cores over a full glacial cycle

The response of natural CH4 sources to climate changes will be an important factor to consider as concentrations of this potent greenhouse gas continue to increase. Polar ice cores provide the means to assess this sensitivity in the past and have shown a close connection between CH4 levels and northern hemisphere temperature variability over the last glacial cycle. However, the contribution of the various CH4 sources and sinks to these changes is still a matter of debate. Contemporaneous stable CH4 isotope records in ice cores provide additional boundary conditions for assessing changes in the CH4 sources and sinks. Here we present new ice core CH4 isotope data covering the last 160,000 years, showing a clear decoupling between CH4 loading and carbon isotopic variations over most of the record. We suggest that d13CH4 variations were not dominated by a change in the source mix but rather by climate- and CO2-related ecosystem control on the isotopic composition of the methane precursor material, especially in seasonally inundated wetlands in the tropics. In contrast, relatively stable d13CH4 intervals occurred during large CH4 loading changes concurrently with past climate changes implying that most CH4 sources (most notably tropical wetlands) responded simultaneously.

Dynamics of spatial heterogeneity in a landfill with interacting phase densities:a stochastic analysis

A landfill represents a complex and dynamically evolving structure that can be stochastically perturbed by exogenous factors. Both thermodynamic (equilibrium) and time varying (non-steady state) properties of a landfill are affected by spatially heterogenous and nonlinear subprocesses that combine with constraining initial and boundary conditions arising from the associated surroundings. While multiple approaches have been made to model landfill statistics by incorporating spatially dependent parameters on the one hand (data based approach) and continuum dynamical mass-balance equations on the other (equation based modelling), practically no attempt has been made to amalgamate these two approaches while also incorporating inherent stochastically induced fluctuations affecting the process overall. In this article, we will implement a minimalist scheme of modelling the time evolution of a realistic three dimensional landfill through a reaction-diffusion based approach, focusing on the coupled interactions of four key variables - solid mass density, hydrolysed mass density, acetogenic mass density and methanogenic mass density, that themselves are stochastically affected by fluctuations, coupled with diffusive relaxation of the individual densities, in ambient surroundings. Our results indicate that close to the linearly stable limit, the large time steady state properties, arising out of a series of complex coupled interactions between the stochastically driven variables, are scarcely affected by the biochemical growth-decay statistics. Our results clearly show that an equilibrium landfill structure is primarily determined by the solid and hydrolysed mass densities only rendering the other variables as statistically "irrelevant" in this (large time) asymptotic limit. The other major implication of incorporation of stochasticity in the landfill evolution dynamics is in the hugely reduced production times of the plants that are now approximately 20-30 years instead of the previous deterministic model predictions of 50 years and above. The predictions from this stochastic model are in conformity with available experimental observations.

Theoretical modeling and experimental validation of fluid-loaded micro rectangular plates

This paper presents a theoretical model on the vibration analysis of micro scale fluid-loaded rectangular isotropic plates, based on the Lamb's assumption of fluid-structure interaction and the Rayleigh-Ritz energy method. An analytical solution for this model is proposed, which can be applied to most cases of boundary conditions. The dynamical experimental data of a series of microfabricated silicon plates are obtained using a base-excitation dynamic testing facility. The natural frequencies and mode shapes in the experimental results are in good agreement with the theoretical simulations for the lower order modes. The presented theoretical and experimental investigations on the vibration characteristics of the micro scale plates are of particular interest in the design of microplate based biosensing devices. Copyright © 2009 by ASME.

Efficient Computation of Electromagnetic Waves in Hydrocarbon Exploration Using the Improved Numerical Mode Matching (NMM) Method

In this study, we developed and improved the numerical mode matching (NMM) method which has previously been shown to be a fast and robust semi-analytical solver to investigate the propagation of electromagnetic (EM) waves in an isotropic layered medium. The applicable models, such as cylindrical waveguide, optical fiber, and borehole with earth geological formation, are generally modeled as an axisymmetric structure which is an orthogonal-plano-cylindrically layered (OPCL) medium consisting of materials stratified planarly and layered concentrically in the orthogonal directions.

In this report, several important improvements have been made to extend applications of this efficient solver to the anisotropic OCPL medium. The formulas for anisotropic media with three different diagonal elements in the cylindrical coordinate system are deduced to expand its application to more general materials. The perfectly matched layer (PML) is incorporated along the radial direction as an absorbing boundary condition (ABC) to make the NMM method more accurate and efficient for wave diffusion problems in unbounded media and applicable to scattering problems with lossless media. We manipulate the weak form of Maxwell's equations and impose the correct boundary conditions at the cylindrical axis to solve the singularity problem which is ignored by all previous researchers. The spectral element method (SEM) is introduced to more efficiently compute the eigenmodes of higher accuracy with less unknowns, achieving a faster mode matching procedure between different horizontal layers. We also prove the relationship of the field between opposite mode indices for different types of excitations, which can reduce the computational time by half. The formulas for computing EM fields excited by an electric or magnetic dipole located at any position with an arbitrary orientation are deduced. And the excitation are generalized to line and surface current sources which can extend the application of NMM to the simulations of controlled source electromagnetic techniques. Numerical simulations have demonstrated the efficiency and accuracy of this method.

Finally, the improved numerical mode matching (NMM) method is introduced to efficiently compute the electromagnetic response of the induction tool from orthogonal transverse hydraulic fractures in open or cased boreholes in hydrocarbon exploration. The hydraulic fracture is modeled as a slim circular disk which is symmetric with respect to the borehole axis and filled with electrically conductive or magnetic proppant. The NMM solver is first validated by comparing the normalized secondary field with experimental measurements and a commercial software. Then we analyze quantitatively the induction response sensitivity of the fracture with different parameters, such as length, conductivity and permeability of the filled proppant, to evaluate the effectiveness of the induction logging tool for fracture detection and mapping. Casings with different thicknesses, conductivities and permeabilities are modeled together with the fractures in boreholes to investigate their effects for fracture detection. It reveals that the normalized secondary field will not be weakened at low frequencies, ensuring the induction tool is still applicable for fracture detection, though the attenuation of electromagnetic field through the casing is significant. A hybrid approach combining the NMM method and BCGS-FFT solver based integral equation has been proposed to efficiently simulate the open or cased borehole with tilted fractures which is a non-axisymmetric model.

Physical Insights, Steady Aerodynamic Effects, and a Design Tool for Low-Pressure Turbine Flutter

The successful, efficient, and safe turbine design requires a thorough understanding of the underlying physical phenomena. This research investigates the physical understanding and parameters highly correlated to flutter, an aeroelastic instability prevalent among low pressure turbine (LPT) blades in both aircraft engines and power turbines. The modern way of determining whether a certain cascade of LPT blades is susceptible to flutter is through time-expensive computational fluid dynamics (CFD) codes. These codes converge to solution satisfying the Eulerian conservation equations subject to the boundary conditions of a nodal domain consisting fluid and solid wall particles. Most detailed CFD codes are accompanied by cryptic turbulence models, meticulous grid constructions, and elegant boundary condition enforcements all with one goal in mind: determine the sign (and therefore stability) of the aerodynamic damping. The main question being asked by the aeroelastician, ``is it positive or negative?'' This type of thought-process eventually gives rise to a black-box effect, leaving physical understanding behind. Therefore, the first part of this research aims to understand and reveal the physics behind LPT flutter in addition to several related topics including acoustic resonance effects. A percentage of this initial numerical investigation is completed using an influence coefficient approach to study the variation the work-per-cycle contributions of neighboring cascade blades to a reference airfoil. The second part of this research introduces new discoveries regarding the relationship between steady aerodynamic loading and negative aerodynamic damping. Using validated CFD codes as computational wind tunnels, a multitude of low-pressure turbine flutter parameters, such as reduced frequency, mode shape, and interblade phase angle, will be scrutinized across various airfoil geometries and steady operating conditions to reach new design guidelines regarding the influence of steady aerodynamic loading and LPT flutter. Many pressing topics influencing LPT flutter including shocks, their nonlinearity, and three-dimensionality are also addressed along the way. The work is concluded by introducing a useful preliminary design tool that can estimate within seconds the entire aerodynamic damping versus nodal diameter curve for a given three-dimensional cascade.

The Development of a Hybrid Optimization Algorithm for the Evaluation and Optimization of the Asynchronous Pulse Unit

The effectiveness of an optimization algorithm can be reduced to its ability to navigate an objective function’s topology. Hybrid optimization algorithms combine various optimization algorithms using a single meta-heuristic so that the hybrid algorithm is more robust, computationally efficient, and/or accurate than the individual algorithms it is made of. This thesis proposes a novel meta-heuristic that uses search vectors to select the constituent algorithm that is appropriate for a given objective function. The hybrid is shown to perform competitively against several existing hybrid and non-hybrid optimization algorithms over a set of three hundred test cases. This thesis also proposes a general framework for evaluating the effectiveness of hybrid optimization algorithms. Finally, this thesis presents an improved Method of Characteristics Code with novel boundary conditions, which better characterizes pipelines than previous codes. This code is coupled with the hybrid optimization algorithm in order to optimize the operation of real-world piston pumps.

What Does Board Capital Really Bring to the Table? Exploring the Effect of Directors’ Human and Social Capital on Effective Governance During International Expansion

What constitutes effective corporate governance? Which director characteristics render boards effective at positively influencing firm-level performance outcomes? This dissertation examines these questions by taking a multilevel, multidisciplinary approach to corporate governance. I explore the individual-, team-, and firm- level factors that enable directors to serve effectively as strategic resources during international expansion. I argue that directors’ international experience improves their ability to serve as effective strategic consultants and resource providers to firms during the complex internationalization process. However, unlike prior research, which tends to assume that directors with the potential to provide important resources uniformly do so, I acknowledge contextual factors (i.e. board cohesiveness, strategic relevance of directors’ experience) that affect their propensity to actually influence outcomes. I explore these issues in three essays: one review essay and two empirical essays. In the first empirical essay, I integrate resource dependence theory with insights from social-psychological research to explore the influence of board capital on firms’ cross-border M&A performance. Using a sample of cross-border M&As completed by S&P 500 firms from 2004-2009, I find evidence that directors’ depth of international experience is associated with superior pre-deal outcomes. This suggests that boards’ deep, market-specific knowledge is valuable during the target selection phase. I further find that directors’ breadth of international experience is associated with superior post-deal performance, suggesting that these directors’ global mindset helps firms in the post-M&A integration phase. I also find that these relationships are positively moderated by board cohesiveness, measured by boards’ internal social ties. In the second empirical essay, I explore the boundary conditions of international board capital by examining how the characteristics of firms’ internationalization strategy moderate the relationship between board capital and firm performance. Using a panel of 377 S&P 500 firms observed from 2004-2011, I find that boards’ depth of international experience and social capital are more important during early stages of internationalization, when firms tend to lack market knowledge and legitimacy in the host markets. On the other hand, I find that breadth of international experience has a stronger relationship with performance when firms’ have higher scope of internationalization, when information-processing demands are higher.

(Fig. 1D, 3E) Sea surface temperature and alkenone C37:4 for sediment core MD01-2444 during the penultimate glaciation

Past glacials can be thought of as natural experiments in which variations in boundary conditions influenced the character of climate change. However, beyond the last glacial, an integrated view of orbital- and millennial-scale changes and their relation to the record of glaciation has been lacking. Here, we present a detailed record of variations in the land-ocean system from the Portuguese margin during the penultimate glacial and place it within the framework of ice-volume changes, with particular reference to European ice-sheet dynamics. The interaction of orbital- and millennial-scale variability divides the glacial into an early part with warmer and wetter overall conditions and prominent climate oscillations, a transitional mid-part, and a late part with more subdued changes as the system entered a maximum glacial state. The most extreme event occurred in the mid-part and was associated with melting of the extensive European ice sheet and maximum discharge from the Fleuve Manche river. This led to disruption of the meridional overturning circulation, but not a major activation of the bipolar seesaw. In addition to stadial duration, magnitude of freshwater forcing, and background climate, the evidence also points to the influence of the location of freshwater discharges on the extent of interhemispheric heat transport.

Benthic foraminifera assemblages of the Gulf of Aden

Benthic and selected planktic foraminifera and stable isotope records were determined in a piston core from the Gulf of Aden, NW Arabian Sea that spans the last 530 ka. The benthic foraminifera were grouped into four principal assemblages using Q-mode Principal Component Analyses. Comparison of each of these assemblages with the fauna of the nearby regions enabled us to identify their specific environmental requirements as a function of variability in food supply and strength of the oxygen minimum zone and by that to use them as indicators of surface water productivity. The benthic foraminiferal productivity indicators coupled with the record of Globigerina bulloides, a planktic foraminifer known to be sensitive to productivity changes in the region, all indicate higher productivity during glacial intervals and productivity similar to present or even reduced during interglacial stages. This trend is opposite to the productivity pattern related to the SW summer monsoon of the Arabian Sea and indicates the role of the NE winter monsoon on the productivity of the Gulf of Aden. A period of exceptionally enhanced productivity is recognized in the Gulf of Aden region between ~60 and 13 kyr indicating the intensification of the NE winter monsoon to its maximal activity. Contemporaneous indication of increased productivity in other parts of the Arabian Sea, unexplained so far by the SW summer monsoon variability, might be related to the intensification of the NE winter monsoon. Another prominent event of high productivity, second in its extent to the last glacial productivity event is recognized between 430 and 460 kyr. These two events seem to correspond to periods of similar orbital positioning of rather low precession (and eccentricity) amplitude for a relatively long period. Glacial boundary conditions seem to control to a large extent the NE winter monsoon variability as also indicated by the dominance of the 100 ka cycle in the investigated time series. Secondary in their importance are the 23 and 41 ka cycles which seem also to contribute to the NE monsoonal variability. Following the identification of productivity events related to the NE winter monsoon in the Gulf of Aden, it is possible now to extend this observation to other parts of the Arabian Sea and consider the contribution of this monsoonal system to the productivity fluctuations preserved in the sedimentary records.

Stable carbon and oxygen isotope record of foraminifera from DSDP Hole 30-289

Sediments accumulate on the sea floor far from land with rates of a few millimetres to a few centimetres per thousand years. Sediments have been accumulating under broadly similar conditions, subject to similar controls, for the past 10 8 years and more. In principle we should be able to study the distribution of climatic variance with frequencies over the range 10**-3 to 10**-7 cycles per year with comparative ease. In fact, nearly all our data are heavily weighted towards the youngest part of the geological record. We study frequencies higher than 10**-4 cycles per year in the special case of a Pleistocene interglacial (the present one), and frequencies in the range 10**-4 to 10**-5 cycles per year in the special case of an ice-age. Although these may be of more direct interest to mankind than earlier periods, it may well be that we will understand the causes of climatic variability better if we can examine their operation over a longer time scale and under different boundary conditions. Rather than review the available data, I have collected some new data to show the feasibility of gathering a data base for examining climatic variability without this usual bias toward the recent. The most widely applicable tool for extracting climatic information from deep-sea sediments is oxygen isotope analysis of calcium carbonate microfossils. It is generally possible to select from the sediment both specimens of benthonic Foraminifera (that is, those that lived in ocean deep water at the sediment-water interface) and specimens of planktonic Foraminifera (that is, those that lived and formed their shells near the ocean surface, and fell to the sediment after death). Thus one is able to monitor conditions at the surface and at depth at simultaneous moments in the geological past. The necessity to analyse calcareous microfossils restricts investigation to calcareous sediments, but even with this restriction in sediment type there are many factors governing the rate of sediment accumulation. On a global scale, sediment accumulates so as to balance the input to the oceans from continental erosion. Even when averaged globally, long-term accumulation rates have varied by almost a factor of ten (Davies et al., 1977, doi:10.1126/science.197.4298.53). At the regional scale, surface productivity and deep-water physical and chemical conditions also affect the sediment accumulation rate. Since all these are susceptible to variation and may well vary in response to climatic change as well as other factors, it is extremely hazardous to attempt to express any climatic variable as a function of time on the basis of measurements originally made as a function of depth in sediment. Although time has been used as a basis for plotting Figs. i-8, these should be regarded as freehand sketches of climatic history rather than as time-series plots.

Results of Closed-Flow Column Experiments (zip-archive 2.3 MB)

The identification of transport parameters by inverse modeling often suffers from equifinality or parameter correlation when models are fitted to observations of the solute breakthrough in column outflow experiments. This parameters uncertainty can be approached by the application of multiple experimental designs such as column experiments in open-flow mode and the recently proposed closed-flow mode. Latter are characterized by the recirculation of the column effluent into the solution supply vessel that feeds the inflow. Depending on the experimental conditions, the solute concentration in the solution supply vessel and the effluent follows a damped sinusoidal oscillation. As a result, the closed-flow experiment provides additional observables in the breakthrough curve. The evaluation of these emergent features allows intrinsic control over boundary conditions and impacts the uncertainty of parameters in inverse modeling. We present a comprehensive sensitivity analysis to illustrate the potential application of closed-flow experiments. We show that the sensitivity with respect to the apparent dispersion can be controlled by the experimenter leading to a decrease in parameter uncertainty as compared to classical experiments by an order of magnitude for optimal settings. With these finding we are also able to reduce the equifinality found for situations, where rate-limited interactions impede a proper determination of the apparent dispersion and rate coefficients. Furthermore, we show the expected breakthrough curve for equilibrium and kinetic sorption, the latter showing strong similarities to the behavior found for completely mixed batch reactor experiments. This renders the closed-flow mode a useful complementary approach to classical column experiments.

Tight-Binding Approximations in 1D and 2D Coupled-Cavity Photonic Crystal Structures

Light confinement and controlling an optical field has numerous applications in the field of telecommunications for optical signals processing. When the wavelength of the electromagnetic field is on the order of the period of a photonic microstructure, the field undergoes reflection, refraction, and coherent scattering. This produces photonic bandgaps, forbidden frequency regions or spectral stop bands where light cannot exist. Dielectric perturbations that break the perfect periodicity of these structures produce what is analogous to an impurity state in the bandgap of a semiconductor. The defect modes that exist at discrete frequencies within the photonic bandgap are spatially localized about the cavity-defects in the photonic crystal. In this thesis the properties of two tight-binding approximations (TBAs) are investigated in one-dimensional and two-dimensional coupled-cavity photonic crystal structures We require an efficient and simple approach that ensures the continuity of the electromagnetic field across dielectric interfaces in complex structures. In this thesis we develop \textrm{E} -- and \textrm{D} --TBAs to calculate the modes in finite 1D and 2D two-defect coupled-cavity photonic crystal structures. In the \textrm{E} -- and \textrm{D} --TBAs we expand the coupled-cavity \overrightarrow{E} --modes in terms of the individual \overrightarrow{E} -- and \overrightarrow{D} --modes, respectively. We investigate the dependence of the defect modes, their frequencies and quality factors on the relative placement of the defects in the photonic crystal structures. We then elucidate the differences between the two TBA formulations, and describe the conditions under which these formulations may be more robust when encountering a dielectric perturbation. Our 1D analysis showed that the 1D modes were sensitive to the structure geometry. The antisymmetric \textrm{D} mode amplitudes show that the \textrm{D} --TBA did not capture the correct (tangential \overrightarrow{E} --field) boundary conditions. However, the \textrm{D} --TBA did not yield significantly poorer results compared to the \textrm{E} --TBA. Our 2D analysis reveals that the \textrm{E} -- and \textrm{D} --TBAs produced nearly identical mode profiles for every structure. Plots of the relative difference between the \textrm{E} and \textrm{D} mode amplitudes show that the \textrm{D} --TBA did capture the correct (normal \overrightarrow{E} --field) boundary conditions. We found that the 2D TBA CC mode calculations were 125-150 times faster than an FDTD calculation for the same two-defect PCS. Notwithstanding this efficiency, the appropriateness of either TBA was found to depend on the geometry of the structure and the mode(s), i.e. whether or not the mode has a large normal or tangential component.

Experimental and RANS-CFD Study of Nacelle Drag

The drag on a nacelle model was investigated experimentally and computationally to provide guidance and insight into the capabilities of RANS-based CFD. The research goal was to determine whether industry constrained CFD could participate in the aerodynamic design of nacelle bodies. Grid refinement level, turbulence model and near wall treatment settings, to predict drag to the highest accuracy, were key deliverables. Cold flow low-speed wind tunnel experiments were conducted at a Reynolds number of 6∙〖10〗^5, 293 K and a Mach number of 0.1. Total drag force was measured by a six-component force balance. Detailed wake analysis, using a seven-hole pressure probe traverse, allowed for drag decomposition via the far-field method. Drag decomposition was performed through a range of angles of attack between 0o and 45o. Both methods agreed on total drag within their respective uncertainties. Reversed flow at the measurement plane and saturation of the load cell caused discrepancies at high angles of attack. A parallel CFD study was conducted using commercial software, ICEM 15.0 and FLUENT 15.0. Simulating a similar nacelle geometry operating under inlet boundary conditions obtained through wind tunnel characterization allowed for direct comparisons with experiment. It was determined that the Realizable k-ϵ was best suited for drag prediction of this geometry. This model predicted the axial momentum loss and secondary flow in the wake, as well as the integrated surface forces, within experimental error up to 20o angle of attack. SST k-ω required additional surface grid resolution on the nacelle suction side, resulting in 15% more elements, due to separation point prediction sensitivity. It was further recommended to apply enhanced wall treatment to more accurately capture the viscous drag and separated flow structures. Overall, total drag was predicted within 5% at 0o angle of attack and 10% at 20o, each within experimental uncertainty. What is more, the form and induced drag predicted by CFD and measured by the wake traverse shared good agreement. Which indicated CFD captured the key flow features accurately despite simplification of the nacelle interior geometry.

Graphene-based waveguide resonators for submillimeter-wave applications

Utilization of graphene covered waveguide inserts to form tunable waveguide resonators is theoretically explained and rigorously investigated by means of full-wave numerical electromagnetic simulations. Instead of using graphene-based switching elements, the concept we propose incorporates graphene sheets as parts of a resonator. Electrostatic tuning of the graphene surface conductivity leads to changes in the electromagnetic field boundary conditions at the resonator edges and surfaces, thus producing an effect similar to varying the electrical length of a resonator. The presented outline of the theoretical background serves to give phenomenological insight into the resonator behavior, but it can also be used to develop customized software tools for design and optimization of graphene-based resonators and filters. Due to the linear dependence of the imaginary part of the graphene surface impedance on frequency, the proposed concept was expected to become effective for frequencies above 100 GHz, which is confirmed by the numerical simulations. A frequency range from 100 GHz up to 1100 GHz, where the rectangular waveguides are used, is considered. Simple, all-graphene-based resonators are analyzed first, to assess the achievable tunability and to check the performance throughout the considered frequency range. Graphene–metal combined waveguide resonators are proposed in order to preserve the excellent quality factors typical for the type of waveguide discontinuities used. Dependence of resonator properties on key design parameters is studied in detail. Dependence of resonator properties throughout the frequency range of interest is studied using eight different waveguide sections appropriate for different frequency intervals. Proposed resonators are aimed at applications in the submillimeter-wave spectral region, serving as the compact tunable components for the design of bandpass filters and other devices.