Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An analytic modeling and system identification study of rotor/fuselage dynamics at hover

A combination of analytic modeling and system identification methods have been used to develop an improved dynamic model describing the response of articulated rotor helicopters to control inputs. A high-order linearized model of coupled rotor/body dynamics including flap and lag degrees of freedom and inflow dynamics with literal coefficients is compared to flight test data from single rotor helicopters in the near hover trim condition. The identification problem was formulated using the maximum likelihood function in the time domain. The dynamic model with literal coefficients was used to generate the model states, and the model was parametrized in terms of physical constants of the aircraft rather than the stability derivatives, resulting in a significant reduction in the number of quantities to be identified. The likelihood function was optimized using the genetic algorithm approach. This method proved highly effective in producing an estimated model from flight test data which included coupled fuselage/rotor dynamics. Using this approach it has been shown that blade flexibility is a significant contributing factor to the discrepancies between theory and experiment shown in previous studies. Addition of flexible modes, properly incorporating the constraint due to the lag dampers, results in excellent agreement between flight test and theory, especially in the high frequency range.

Hong, Steven W.↗
NASA NTRS2013

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗
NASA NTRS2000

Generation and Radiation of Acoustic Waves from a 2-D Shear Layer

A parallel numerical simulation of the radiation of sound from an acoustic source inside a 2-D jet is presented in this paper. This basic benchmark problem is used as a test case for scattering problems that are presently being solved by using the Impedance Mismatch Method (IMM). In this technique, a solid body in the domain is represented by setting the acoustic impedance of each medium, encountered by a wave, to a different value. This impedance discrepancy results in reflected and scattered waves with appropriate amplitudes. The great advantage of the use of this method is that no modifications to a simple Cartesian grid need to be made for complicated geometry bodies. Thus, high order finite difference schemes may be applied simply to all parts of the domain. In the IMM, the total perturbation field is split into incident and scattered fields. The incident pressure is assumed to be known and the equivalent sources for the scattered field are associated with the presence of the scattering body (through the impedance mismatch) and the propagation of the incident field through a non-uniform flow. An earlier version of the technique could only handle uniform flow in the vicinity of the source and at the outflow boundary. Scattering problems in non-uniform mean flow are of great practical importance (for example, scattering from a high lift device in a non-uniform mean flow or the effects of a fuselage boundary layer). The solution to this benchmark problem, which has an acoustic wave propagating through a non-uniform mean flow, serves as a test case for the extensions of the IMM technique.

Agarwal, Anurag↗
NASA NTRS2001

Flame Synthesis Used to Create Metal-Catalyzed Carbon Nanotubes

Metal-catalyzed carbon nanotubes are highly ordered carbon structures of nanoscale dimensions. They may be thought of as hollow cylinders whose walls are formed by single atomic layers of graphite. Such cylinders may be composed of many nested, concentric atomic layers of carbon or only a single layer, the latter forming a single-walled carbon nanotube. This article reports unique results using a flame for their synthesis. Only recently were carbon nanotubes discovered within an arc discharge and recognized as fullerene derivatives. Today metal-catalyzed carbon nanotubes are of great interest for many reasons. They can be used as supports for the metal catalysts like those found in catalytic converters. Open-ended nanotubes are highly desirable because they can be filled by other elements, metals or gases, for battery and fuel cell applications. Because of their highly crystalline structure, they are significantly stronger than the commercial carbon fibers that are currently available (10 times as strong as steel but possessing one-sixth of the weight). This property makes them highly desirable for strengthening polymer and ceramic composite materials. Current methods of synthesizing carbon nanotubes include thermal pyrolysis of organometallics, laser ablation of metal targets within hydrocarbon atmospheres at high temperatures, and arc discharges. Each of these methods is costly, and it is unclear if they can be scaled for the commercial synthesis of carbon nanotubes. In contrast, flame synthesis is an economical means of bulk synthesis of a variety of aerosol materials such as carbon black. Flame synthesis of carbon nanotubes could potentially realize an economy of scale that would enable their use in common structural materials such as car-body panels. The top figure is a transmission electron micrograph of a multiwalled carbon nanotube. The image shows a cross section of the atomic structure of the nanotube. The dark lines are individual atomic layer planes of carbon, seen here in cross section. They form a nested series of concentric cylinders, much like the growth rings on a tree. This sample was obtained by the supported catalyst method, whereby the nanoscale catalysts are dispersed on a substrate providing their support. The substrate with catalyst particles was immersed within an acetylene diffusion flame to which nitrogen had been added to eliminate soot formation. Upon removal from the flame, the nanotubes were dispersed on a holder suitable for electron microscopy. Although not seen in the figure, the tube diameter reflects that of the catalyst particle.

VanderWal, Randy L.↗
NASA NTRS2019

Effect of Spatial Filtering in Implicit Large-Eddy Simulations of Separated Flows

The relatively high Reynolds number of turbulent flows encountered in various applications puts these problems well beyond the reach of direct numerical simulation (DNS) at present. Meanwhile, lower-fidelity Reynolds-averaged Navier-Stokes (RANS) calculations are known to be not accurate enough in complex problems, such as smooth-body flow separation and other flows involving highly-unsteady phenomena. Hence, given the current infeasibility of DNS and the unsatisfactory performance of RANS, intermediate techniques such as large-eddy simulation (LES) and hybrid RANS-LES, whose fidelity lie between RANS and DNS, have received much attention for application to various problems of practical importance. Modeling of the effect of missing scales on resolved scales, also known as subgrid-scale (SGS) modeling, is an important subject for LES. SGS models can be broadly categorized as explicit or implicit approaches. The explicit approach is based on an SGS model that explicitly appears in the governing equations expressed in the form of so-called “filtered Navier-Stokes equations”, which describe the evolution of the turbulence scales resolved by the LES grid. The effect of the scales unresolved by the grid is represented by the SGS model. The implicit modeling approach, on the other hand, does not employ an explicit model but instead treats the intrinsic dissipation of the numerical discretization scheme as an implicit SGS model. An LES without an explicit SGS model is commonly termed as an implicit LES (ILES). The relative merits of one SGS modeling approach over another is a subject of ongoing debate. We have opted to employ an ILES methodology, based on high-order compact finite-difference and spatial filtering schemes, in our recent investigations of separated flow problems [1, 2]. Further discussion of our preference of ILES over explicit LES is provided in Uzun and Malik [2]. The spatial filtering operation, described in the next section, is treated as an implicit SGS model for the ILES. Some observations made during the course of our recent investigations, which pointed out to excessive numerical dissipation in certain parts of the flowfield, prompted us to take a closer look at the potential effect of the spatial filter on ILES predictions. This technical note is therefore devoted to spatial filter effects in the context of a high Reynolds number, transonic shock-induced separated flow.

Ali Uzun↗
NASA NTRS2016

The Space-Time Conservative Schemes for Large-Scale, Time-Accurate Flow Simulations with Tetrahedral Meshes

Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.

Venkatachari, Balaji Shankar↗
NASA NTRS1993

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗
NASA NTRS2010

Trajectory Optimization: OTIS 4

The latest release of the Optimal Trajectories by Implicit Simulation (OTIS4) allows users to simulate and optimize aerospace vehicle trajectories. With OTIS4, one can seamlessly generate optimal trajectories and parametric vehicle designs simultaneously. New features also allow OTIS4 to solve non-aerospace continuous time optimal control problems. The inputs and outputs of OTIS4 have been updated extensively from previous versions. Inputs now make use of objectoriented constructs, including one called a metastring. Metastrings use a greatly improved calculator and common nomenclature to reduce the user s workload. They allow for more flexibility in specifying vehicle physical models, boundary conditions, and path constraints. The OTIS4 calculator supports common mathematical functions, Boolean operations, and conditional statements. This allows users to define their own variables for use as outputs, constraints, or objective functions. The user-defined outputs can directly interface with other programs, such as spreadsheets, plotting packages, and visualization programs. Internally, OTIS4 has more explicit and implicit integration procedures, including high-order collocation methods, the pseudo-spectral method, and several variations of multiple shooting. Users may switch easily between the various methods. Several unique numerical techniques such as automated variable scaling and implicit integration grid refinement, support the integration methods. OTIS4 is also significantly more user friendly than previous versions. The installation process is nearly identical on various platforms, including Microsoft Windows, Apple OS X, and Linux operating systems. Cross-platform scripts also help make the execution of OTIS and post-processing of data easier. OTIS4 is supplied free by NASA and is subject to ITAR (International Traffic in Arms Regulations) restrictions. Users must have a Fortran compiler, and a Python interpreter is highly recommended.

Riehl, John P.↗
NASA NTRS2006

Fibonacci Grids

Recent years have seen a resurgence of interest in a variety of non-standard computational grids for global numerical prediction. The motivation has been to reduce problems associated with the converging meridians and the polar singularities of conventional regular latitude-longitude grids. A further impetus has come from the adoption of massively parallel computers, for which it is necessary to distribute work equitably across the processors; this is more practicable for some non-standard grids. Desirable attributes of a grid for high-order spatial finite differencing are: (i) geometrical regularity; (ii) a homogeneous and approximately isotropic spatial resolution; (iii) a low proportion of the grid points where the numerical procedures require special customization (such as near coordinate singularities or grid edges). One family of grid arrangements which, to our knowledge, has never before been applied to numerical weather prediction, but which appears to offer several technical advantages, are what we shall refer to as "Fibonacci grids". They can be thought of as mathematically ideal generalizations of the patterns occurring naturally in the spiral arrangements of seeds and fruit found in sunflower heads and pineapples (to give two of the many botanical examples). These grids possess virtually uniform and highly isotropic resolution, with an equal area for each grid point. There are only two compact singular regions on a sphere that require customized numerics. We demonstrate the practicality of these grids in shallow water simulations, and discuss the prospects for efficiently using these frameworks in three-dimensional semi-implicit and semi-Lagrangian weather prediction or climate models.

Swinbank, Richard↗
NASA NTRS2018

Optical Relay for Future NASA Geosynchronous Orbiting Satellite for High Data Rate Links to NASA User Missions

NASA is exploring options for its Next Generation Relay (NGR) architecture while the current Tracking Data Relay Satellite System (TDRSS) completes its mission. The plan is to start implementation of the NGR beginning around 2025. The new system of proposed relay satellites will greatly increase the data rates between low Earth orbiting (LEO) satellite missions and the NASA TDRSS relay satellites. This increase in data rates will allow an unprecedented increase in data throughput from the LEO satellite missions back to the principal investigators (PI). This can be accomplished at Ka-band frequencies with high order modulation or at optical frequencies using Differential Phase Shift Keying (DPSK). The first satellite in the next set of relay satellites will have to be backward compatible with current technology to support ongoing and planned missions. The new set of satellites will be launched over a 10-year period with design lifetimes of at least 15 years. To meet these requirements, we analyzed various architectures and designed both the communication payloads on the relay satellite and candidate payloads on the user spacecraft by utilizing optical heads already designed. From this analysis, a demonstration optical satellite named “the Next Generation Optical Relay Pathfinder” with Ka-band capabilities was proposed to be built and launched with the purpose of evaluating an integrated high-speed optical and Ka-band communication system. Given a cost limit for the demonstration satellite, various satellite configurations were developed by varying the number of optical communication payloads. The communication payload on the relay satellite consisted of three major sub-systems: 1) Optical communication payload, 2) Ka-band communication payload, 3) Digital processing and routing of signals. The size, mass (weight), and power (SWaP) of the communication payload and other sub-systems of the satellite were obtained. The NASA Glenn Research Center COMPASS team designed the Pathfinder satellite and performed a cost analysis for its build and launch. In this paper, we first describe the needs, drivers, and the associated challenges for the Next Generation Optical Relay Pathfinder to be capable of connecting multiple LEO and GEO satellites at high data rates. Second, we detail the concept of operations (ConOps) and the system architecture, including the satellite configurations considered, their attributes and limitations, and the size of the satellite needed for each configuration. Third, we provide a summary of the Next Generation Optical Relay Pathfinder satellite design trades and its key elements. Finally, we present the path needed for implementation and operations.

Warner, Joseph H. D.↗
NASA NTRS2006

Imaging Spectroscopy of Sunspots using IBIS

We use the Interferometric BIdimensional Spectrometer (IBIS) of the INAF/Arcetri Astrophysical Observatory and installed at the National Solar Observatory (NSO) Dunn Solar Telescope, to understand the structure of sunspots. These high resolution observations were acquired on 2004 July 30-31, of active region NOAA 10654, using the high order NSO adaptive optics system. We map the spatio-temporal variation of the penumbral Doppler signatures in three spectral lines, FeI 6301.5\AA, FeII 7224.4\AA, and CaII 8542.6\AA, from the photosphere to the chromosphere. From a 70-minute temporal average of individual 32-second cadence Doppler observations we find that the averaged velocities decrease with height, about 3.5 times larger in the deeper photosphere (FeII 7224.4\AA; height-of-formation $\approx$\50 km) than in the upper photosphere FeI 6301.5\AA; height-of-formation $\approx$\350 km), There is a remarkable coherence of Doppler signals over the height difference of 300 km. From a highspeed animation of the Doppler sequence we find evidence for what appears to be ejection of high speed gas concentrations from edges of penumbral filaments into the surrounding granular photosphere. The Evershed flow persists a few arcseconds beyond the traditionally demarcated penumbra-granulation boundary. We present these and other results and discuss the implications of these measurements for sunspot models

Balasubramaniam, K. S.↗
NASA NTRS1992

Practical aspects of spatially high accurate methods

The computational qualities of high order spatially accurate methods for the finite volume solution of the Euler equations are presented. Two dimensional essentially non-oscillatory (ENO), k-exact, and 'dimension by dimension' ENO reconstruction operators are discussed and compared in terms of reconstruction and solution accuracy, computational cost and oscillatory behavior in supersonic flows with shocks. Inherent steady state convergence difficulties are demonstrated for adaptive stencil algorithms. An exact solution to the heat equation is used to determine reconstruction error, and the computational intensity is reflected in operation counts. Standard MUSCL differencing is included for comparison. Numerical experiments presented include the Ringleb flow for numerical accuracy and a shock reflection problem. A vortex-shock interaction demonstrates the ability of the ENO scheme to excel in simulating unsteady high-frequency flow physics.

Godfrey, Andrew G.↗
NASA NTRS1994

Theoretical study of the incompressible Navier-Stokes equations by the least-squares method

Usually the theoretical analysis of the Navier-Stokes equations is conducted via the Galerkin method which leads to difficult saddle-point problems. This paper demonstrates that the least-squares method is a useful alternative tool for the theoretical study of partial differential equations since it leads to minimization problems which can often be treated by an elementary technique. The principal part of the Navier-Stokes equations in the first-order velocity-pressure-vorticity formulation consists of two div-curl systems, so the three-dimensional div-curl system is thoroughly studied at first. By introducing a dummy variable and by using the least-squares method, this paper shows that the div-curl system is properly determined and elliptic, and has a unique solution. The same technique then is employed to prove that the Stokes equations are properly determined and elliptic, and that four boundary conditions on a fixed boundary are required for three-dimensional problems. This paper also shows that under four combinations of non-standard boundary conditions the solution of the Stokes equations is unique. This paper emphasizes the application of the least-squares method and the div-curl method to derive a high-order version of differential equations and additional boundary conditions. In this paper, an elementary method (integration by parts) is used to prove Friedrichs' inequalities related to the div and curl operators which play an essential role in the analysis.

Jiang, Bo-Nan↗
NASA NTRS1982

Observation of additional low-degree 5-min modes of solar oscillation

High-order solar oscillations with degrees l=3, 4, and 5 could be detected. The observations were made by measuring the difference between the shifts in the Fe 5,124 spectrum line from light integrated from a central circular portion of the solar disk and from an annular portion exterior to it. The frequencies of the octupole modes agree well with the values obtained from whole-disk measurements at the South Pole. A least-squares fit of the observed frequencies to values interpolated between and extrapolated from the predictions of a sequence of solar models with different chemical compositions selects two models. One, a helium-rich solution, agrees with that of similar analyses of whole-disk data. The extrapolated solution has a relatively deep convection zone, and is thus consistent with analyses of 5-min oscillations of high degree.

Scherrer, P. H.↗
NASA NTRS1996

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗
NASA NTRS1993

An analysis of nearfield normal mode amplitude anomalies of the Landers earthquake

The 1992 Landers earthquake (M(sub w) = 7.3) occurred in the middle of the TERRAscope network. Long-period Rayleigh waves recorded at the TERRAscope stations (delta less than or = 3 deg) after traveling around the Earth show large amplitude anomalies, one order of magnitude larger than spherical Earth predictions up to a period of about 600 s. The ground motions over the epicentral region at and after the arrival of R4-5 are in phase at all stations. These observations are inconsistent with the nearly vertical strike slip mechanism of the Landers earthquake. Synthetic seismograms for a rotating, elliptic, and laterally heterogeneous Earth model calculated by the variational method agree well with the observed waveforms. Calculations for various 3D Earth models demonstrate that the amplitudes are very sensitive to the large scale aspherical structure in the crust and the mantle. The anomalies for modes shorter than 300 s period can be explained by lateral heterogeneity shallower than the upper mantle. Rotation of the Earth and lower mantle heterogeneity are required to explain mode amplitudes at longer periods. Current whole mantle seismic tomographic models can fully explain the observed amplitudes longer than 300 s. To assess the effect of the high order lateral heterogeneity in the mantle more precise estimate of the crustal correction is required.

Watada, Shingo↗
NASA NTRS2018

High-Fidelity Computational Aerodynamics of Multi-Rotor Unmanned Aerial Vehicles

High-fidelity Computational Fluid Dynamics (CFD) simulations have been carried out for several multi-rotor Unmanned Aerial Vehicles (UAVs). Three vehicles have been studied: the classic quadcopter DJI Phantom 3, an unconventional quadcopter specialized for forward flight, the SUI Endurance, and an innovative concept for Urban Air Mobility (UAM), the Elytron 4S UAV. The three-dimensional unsteady Navier-Stokes equations are solved on overset grids using high-order accurate schemes, dual-time stepping, and a hybrid turbulence model. The DJI Phantom 3 is simulated with different rotors and with both a simplified airframe and the real airframe including landing gear and a camera. The effects of weather are studied for the DJI Phantom 3 quadcopter in hover. The SUI En- durance original design is compared in forward flight to a new configuration conceived by the authors, the hybrid configuration, which gives a large improvement in forward thrust. The Elytron 4S UAV is simulated in helicopter mode and in airplane mode. Understanding the complex flows in multi-rotor vehicles will help design quieter, safer, and more efficient future drones and UAM vehicles.

High-Fidelity↗
NASA NTRS2013

Planar Superconducting Millimeter-Wave/Terahertz Channelizing Filter

This innovation is a compact, superconducting, channelizing bandpass filter on a single-crystal (0.45 m thick) silicon substrate, which operates from 300 to 600 GHz. This device consists of four channels with center frequencies of 310, 380, 460, and 550 GHz, with approximately 50-GHz bandwidth per channel. The filter concept is inspired by the mammalian cochlea, which is a channelizing filter that covers three decades of bandwidth and 3,000 channels in a very small physical space. By using a simplified physical cochlear model, and its electrical analog of a channelizing filter covering multiple octaves bandwidth, a large number of output channels with high inter-channel isolation and high-order upper stopband response can be designed. A channelizing filter is a critical component used in spectrometer instruments that measure the intensity of light at various frequencies. This embodiment was designed for MicroSpec in order to increase the resolution of the instrument (with four channels, the resolution will be increased by a factor of four). MicroSpec is a revolutionary wafer-scale spectrometer that is intended for the SPICA (Space Infrared Telescope for Cosmology and Astrophysics) Mission. In addition to being a vital component of MicroSpec, the channelizing filter itself is a low-resolution spectrometer when integrated with only an antenna at its input, and a detector at each channel s output. During the design process for this filter, the available characteristic impedances, possible lumped element ranges, and fabrication tolerances were identified for design on a very thin silicon substrate. Iterations between full-wave and lumped-element circuit simulations were performed. Each channel s circuit was designed based on the availability of characteristic impedances and lumped element ranges. This design was based on a tabular type bandpass filter with no spurious harmonic response. Extensive electromagnetic modeling for each channel was performed. Four channels, with 50-GHz bandwidth, were designed, each using multiple transmission line media such as microstrip, coplanar waveguide, and quasi-lumped components on 0.45- m thick silicon. In the design process, modeling issues had to be overcome. Due to the extremely high frequencies, very thin Si substrate, and the superconducting metal layers, most commercially available software fails in various ways. These issues were mitigated by using alternative software that was capable of handling them at the expense of greater simulation time. The design of on-chip components for the filter characterization, such as a broadband antenna, Wilkinson power dividers, attenuators, detectors, and transitions has been completed.

Ehsan, Negar↗