Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level set method”

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 955 records · Page 53

Distributed, Intelligent Edge-Sensing for a Smarter Grid

The electric grid is undergoing major transformations and developments resulting in unprecedented levels of volatility, uncertainty, and stress on grid infrastructure. Smart sensors and methods aiding in advanced visibility and situational awareness are key for tackling these issues. In this work, a decentralized architecture is proposed, where sensing, local computation and control capability are embedded in the edge devices, communicating with a set of trusted 'data mules' in a 'delay-tolerant' manner, while functioning autonomously. This system has been designed and implemented as an overall platform – called Global Asset Monitoring, Management and Analytics (GAMMA) Platform intended to provide the backbone for a global array of sensors and actuators. Further, as a building block for advanced current sensing solutions, a smart, low-cost ‘clip-on’ current sensor based on PCB-embedded Rogowski coil has been developed. The sensor hosts a novel signal conditioning stage allowing an 'auto-tuning' feature, resulting in a universal current sensor design for measuring a wide range of currents, including faults for smart grid applications. Finally, the research proposes a method to instrument and monitor key parameters for the most common electric utility asset – the pole-top distribution transformer. The work done in this research enables scalable, edge-intelligent sensing solutions for monitoring grid infrastructure, allowing utility operators to gain advanced visibility in an economical way.

Kulkarni, Shreyas Bhalchandra↗

Focus determination for the James Webb Space Telescope Science Instruments: A Survey of Methods

The James Webb Space Telescope (JWST) is a segmented deployable telescope that will require on-orbit alignment using the Near Infrared Camera as a wavefront sensor. The telescope will be aligned by adjusting seven degrees of freedom on each of 18 primary mirror segments and five degrees of freedom on the secondary mirror to optimize the performance of the telescope and camera at a wavelength of 2 microns. With the completion of these adjustments, the telescope focus is set and the optical performance of each of the other science instruments should then be optimal without making further telescope focus adjustments for each individual instrument. This alignment approach requires confocality of the instruments after integration and alignment to the composite metering structure, which will be verified during instrument level testing at Goddard Space Flight Center with a telescope optical simulator. In this paper, we present the results from a study of several analytical approaches to determine the focus for each instrument. The goal of the study is to compare the accuracies obtained for each method, and to select the most feasible for use during optical testing.

Davila, Pamela S.↗

SU-F-J-214: Dose Reduction by Spatially Optimized Image Quality Via Fluence Modulated Proton CT (FMpCT)

Purpose: Proton CT (pCT) is a promising imaging modality for reducing range uncertainty in image-guided proton therapy. Range uncertainties partially originate from X-ray CT number conversion to stopping power ratio (SPR) and are limiting the exploitation of the full potential of proton therapy. In this study we explore the concept of spatially dependent fluence modulated proton CT (FMpCT), for achieving optimal image quality in a clinical region of interest (ROI), while reducing significantly the imaging dose to the patient. Methods: The study was based on simulated ideal pCT using pencil beam (PB) scanning. A set of 250 MeV protons PBs was used to create 360 projections of a cylindrical water phantom and a head and neck cancer patient. The tomographic images were reconstructed using a filtered backprojection (FBP) as well as an iterative algorithm (ITR). Different fluence modulation levels were investigated and their impact on the image was quantified in terms of SPR accuracy as well as noise within and outside selected ROIs, as a function of imaging dose. The unmodulated image served as reference. Results: Both FBP reconstruction and ITR without total variation (TV) yielded image quality in the ROIs similar to the reference images, for modulation down to 0.1 of the full proton fluence. The average dose was reduced by 75% for the water phantom and by 40% for the patient. FMpCT does not improve the noise for ITR with TV and modulation 0.1. Conclusion: This is the first work proposing and investigating FMpCT for producing optimal image quality for treatment planning and image guidance, while simultaneously reducing imaging dose. Future work will address spatial resolution effects and the impact of FMpCT on the quality of proton treatment plans for a prototype pCT scanner capable of list mode data acquisition. Acknowledgement: DFG-MAP DFG - Munich-Centre for Advanced Photonics (MAP)

60 APPLIED LIFE SCIENCES↗

Metal Hydroxide and Oxyhydroxide Gas Database

In a high temperature combustion situation such as a gas turbine or rocket engine, combustion products include water vapor. This water vapor reacts with oxides in the hot-gas-path to form volatile hydroxides and oxyhydroxides. The oxides may be structural oxide components or protective oxides on alloys or non-oxide ceramics. In either case, this enhanced volatilization leads to degradation of components. Reliable thermochemical data for many metal hydroxides and oxyhydroxides has not been available. We review a database put together at NASA over many years based on experimental and ab initio computations for thermochemical data of Al, Si, Ti, Cr, Mn, Y, Zr, La, Ta, Gd, and Yb hydroxides and oxyhydroxides. Experimental methods for measuring thermodynamic quantities of metal hydroxides and oxyhydroxides are limited since the required high temperatures and high pressures preclude many conventional methods. Nonetheless, transpiration has proven a valuable tool in our laboratories for the Si-O-H, Cr-O-H, and Ti-O-H systems. Data are analyzed with second and third law methods. These data can be leveraged with quantum chemistry to provide reliable free energy functions for the third law calculation and comparison to the derived enthalpies of formation. Both wave function and DFT methods can be utilized for deriving the thermodynamic functions for the various metal hydroxides and oxyhydroxides. We have found that, in general, reasonable molecular geometries and spectroscopic data can be derived with a relatively low level of theory, i.e. DFT with the B3LYP functional. Accurate calculations of the enthalpy of formation require higher levels of theory, such as the Coupled Cluster, Singles, Doubles, and perturbative Triples CCSD(T) approach with complete basis sets (CBS). Also, corrections for hindered rotations and anharmonicity are necessary. For some symmetries, the anharmonic correction has proven to be a challenge and we have found that this is only treated properly by certain codes. The database lists the fHo(298), So(298), and Cp in the form Cp = A + BT + CT2 + D/T + E/T2, where A, B, C, D, and E are the fitting constants. Thus, it can be easily used in any of the common computational thermochemistry codes. To date, the database contains 51 gaseous species of the formula MOx(OH)y.

Metal Hydroxides↗

Spectroscopic and theoretical studies of UN and UN +

The low-energy electronic states of UN and UN+ have been examined using high-level electronic structure calculations and two-color photoionization techniques. The experimental measurements provided an accurate ionization energy for UN (IE = 50 802 +/- 5 cm -1 ). Spectra for UN+ yielded ro-vibrational constants and established that the ground state has the electronic angular momentum projection omega = 4. Ab initio calculations were carried out using the spin-orbit state interacting approach with the complete active space second-order perturbation theory method. A series of correlation consistent basis sets were used in conjunction with small-core relativistic pseudopotentials on U to extrapolate to the complete basis set limits. The results for UN correctly obtained an omega = 3.5 ground state and demonstrated a high density of configurationally related excited states with closely similar ro-vibrational constants. Similar results were obtained for UN + , with reduced complexity owing to the smaller number of outer-shell electrons. The calculated IE for UN was in excellent agreement with the measured value. Improved values for the dissociation energies of UN and UN + , as well as their heats of formation, were obtained using the Feller-Peterson-Dixon composite thermochemistry method, including corrections up through coupled cluster singles, doubles, triples and quadruples. An analysis of the ab initio results from the perspective of the ligand field theory shows that the patterns of electronic states for both UN and UN + can be understood in terms of the underlying energy level structure of the atomic metal ion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Systematic parameter inference in stochastic mesoscopic modeling

We propose a method to efficiently determine the optimal coarse-grained force field in mesoscopic stochastic simulations of Newtonian fluid and polymer melt systems modeled by dissipative particle dynamics (DPD) and energy conserving dissipative particle dynamics (eDPD). The response surfaces of various target properties (viscosity, diffusivity, pressure, etc.) with respect to model parameters are constructed based on the generalized polynomial chaos (gPC) expansion using simulation results on sampling points (e.g., individual parameter sets). To alleviate the computational cost to evaluate the target properties, we employ the compressive sensing method to compute the coefficients of the dominant gPC terms given the prior knowledge that the coefficients are “sparse”. The proposed method shows comparable accuracy with the standard probabilistic collocation method (PCM) while it imposes a much weaker restriction on the number of the simulation samples especially for systems with high dimensional parametric space. Fully access to the response surfaces within the confidence range enables us to infer the optimal force parameters given the desirable values of target properties at the macroscopic scale. Moreover, it enables us to investigate the intrinsic relationship between the model parameters, identify possible degeneracies in the parameter space, and optimize the model by eliminating model redundancies. The proposed method provides an efficient alternative approach for constructing mesoscopic models by inferring model parameters to recover target properties of the physics systems (e.g., from experimental measurements), where those force field parameters and formulation cannot be derived from the microscopic level in a straight forward way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Methods for Characterizing Fine Particulate Matter Using Satellite Remote-Sensing Data and Ground Observations: Potential Use for Environmental Public Health Surveillance

This study describes and demonstrates different techniques for surfacing daily environmental / hazards data of particulate matter with aerodynamic diameter less than or equal to 2.5 micrometers (PM2.5) for the purpose of integrating respiratory health and environmental data for the Centers for Disease Control and Prevention (CDC s) pilot study of Health and Environment Linked for Information Exchange (HELIX)-Atlanta. It described a methodology for estimating ground-level continuous PM2.5 concentrations using B-Spline and inverse distance weighting (IDW) surfacing techniques and leveraging National Aeronautics and Space Administration (NASA) Moderate Resolution Imaging Spectrometer (MODIS) data to complement The Environmental Protection Agency (EPA) ground observation data. The study used measurements of ambient PM2.5 from the EPA database for the year 2003 as well as PM2.5 estimates derived from NASA s satellite data. Hazard data have been processed to derive the surrogate exposure PM2.5 estimates. The paper has shown that merging MODIS remote sensing data with surface observations of PM2.5 not only provides a more complete daily representation of PM2.5 than either data set alone would allow, but it also reduces the errors in the PM2.5 estimated surfaces. The results of this paper have shown that the daily IDW PM2.5 surfaces had smaller errors, with respect to observations, than those of the B-Spline surfaces in the year studied. However the IDW mean annual composite surface had more numerical artifacts, which could be due to the interpolating nature of the IDW that assumes that the maxima and minima can occur only at the observation points. Finally, the methods discussed in this paper improve temporal and spatial resolutions and establish a foundation for environmental public health linkage and association studies for which determining the concentrations of an environmental hazard such as PM2.5 with good accuracy levels is critical.

Al-Hamdan, Mohammad Z.↗

Niobium coaxial cavities with internal quality factors exceeding 1.4 × 10 9 for circuit quantum electrodynamics

Group-V materials such as niobium and tantalum have become popular choices for extending the performance of circuit quantum electrodynamics (cQED) platforms, allowing for quantum processors and memories with reduced error rates and more modes. The complex surface chemistry of niobium, however, makes identifying the main modes of decoherence difficult at millikelvin temperatures and single-photon powers. We use niobium coaxial quarter-wave cavities to study the impact of etch chemistry, prolonged atmospheric exposure, and the significance of cavity conditions prior to and during cooldown—in particular, niobium hydride evolution—on single-photon coherence. We demonstrate cavities with quality factors Q int ≳ 1.4 × 10 9 in the single-photon regime, a 15-fold improvement over aluminum cavities of the same geometry. We rigorously quantify the sensitivity of our fabrication process to various loss mechanisms and demonstrate a two- to fourfold reduction in the two-level system loss tangent and a three- to fivefold improvement in the residual resistance over traditional buffered chemical polishing etching techniques. Finally, we demonstrate transmon integration and coherent cavity control while maintaining a cavity coherence of 11.3 ms. The accessibility of our method, which can be easily replicated in academic laboratory settings, together with the demonstration of its performance, mark an advancement in three-dimensional cQED.

Oriani, Andrew E. [Chicago U., EFI; Chicago U.] (O↗

SU-E-T-422: Fast Analytical Beamlet Optimization for Volumetric Intensity-Modulated Arc Therapy

Purpose: To implement a fast optimization algorithm on CPU/GPU heterogeneous computing platform and to obtain an optimal fluence for a given target dose distribution from the pre-calculated beamlets in an analytical approach. Methods: The 2D target dose distribution was modeled as an n-dimensional vector and estimated by a linear combination of independent basis vectors. The basis set was composed of the pre-calculated beamlet dose distributions at every 6 degrees of gantry angle and the cost function was set as the magnitude square of the vector difference between the target and the estimated dose distribution. The optimal weighting of the basis, which corresponds to the optimal fluence, was obtained analytically by the least square method. Those basis vectors with a positive weighting were selected for entering into the next level of optimization. Totally, 7 levels of optimization were implemented in the study.Ten head-and-neck and ten prostate carcinoma cases were selected for the study and mapped to a round water phantom with a diameter of 20cm. The Matlab computation was performed in a heterogeneous programming environment with Intel i7 CPU and NVIDIA Geforce 840M GPU. Results: In all selected cases, the estimated dose distribution was in a good agreement with the given target dose distribution and their correlation coefficients were found to be in the range of 0.9992 to 0.9997. Their root-mean-square error was monotonically decreasing and converging after 7 cycles of optimization. The computation took only about 10 seconds and the optimal fluence maps at each gantry angle throughout an arc were quickly obtained. Conclusion: An analytical approach is derived for finding the optimal fluence for a given target dose distribution and a fast optimization algorithm implemented on the CPU/GPU heterogeneous computing environment greatly reduces the optimization time.

60 APPLIED LIFE SCIENCES↗

The MidAmerica Regional Microgrid Education and Training (MARMET) Consortium

The primary goal of this project was to integrate cutting-edge research and advanced instructional methods to create a flexible, evolving approach to microgrid and distributed energy resources training for all levels of student. This goal entailed the rapid development of modular material to capture and reflect the newest trends in electric power microgrid engineering and to make it available in traditional and non-traditional settings. It was our goal to develop course material that can be readily integrated into either existing courses or organized into new courses. The target audiences include technicians, contractors and consultants, inspectors, policy makers, supervisors, and engineers and undergraduate and graduate students. The main project objectives were to: 1) Perform a gap analysis to identify the research/educational/training needs and define specific implementation plans for these activities in distributed technology and Smart Grid power engineering by working with Department of Energy, the National Network Administrator (NNA) and other Distributed Technology Training Consortia (DTTC). 2) Develop learning content related to, and driven by, microgrid research, consisting of one or more multimedia-based modules. Over the life of the project MARMET provided training to 11,012 students and developed 154,432 contact hours of material (see the IREC NNA, DE-EE0006342, final project report for added details on MARMET evaluation and assessment metrics, as well as an overview of the entire GEARED network. 3) Provide a variety of delivery mechanisms and venues for face-to-face and on-line instruction in power system engineering and related smart grid operations, including university-based courses, short courses, and workshops. 4) Assemble and engage highly-qualified advisory and steering committees of experts from utilities, universities, national laboratories, and product suppliers for this DTTC, to guide and advise the activities of the DTTC, and to assist in creating student internships, coops, research projects and academic competitions for GEARED students within the DTTC and nationally. 5) Convene regular meetings of Consortia partners to disseminate research results among the participants (universities, utilities, and industry). 6) Develop recruitment mechanisms to publicize GEARED student opportunities and highlight the educational benefits of GEARED participation to assure long-term sustainability of student recruitment and retention. 7) Participate in GEARED National Network activities to identify and adopt best practices across all GEARED consortia. 8) Organize, host and manage at least one annual Student-Centered Research Conference on behalf of the GEARED Network during the period of performance; participate in each annual national Student-Centered Research Conference.

14 SOLAR ENERGY↗

A new scheme for multidimensional line transfer. II - ETLA method in one dimension with application to iron K-alpha lines

The 1D version of the radiation hydrodynamics computer code ALTAIR is discribed. It is an implementation of the equivalent two-level atom method, and also incorporates a variety of other iterative techniques to achieve a fast, accurate solution of the problem of coupled atomic kinetics and radiative transfer. All the techniques are described in sufficient detail to enable their use by others. The application of the program is illustrated by a set of calculations of the iron X-ray emission spectrum from intercloud hot gas that may exist in the broad-line region of active galactic nuclei. These calculations demonstrate that such spectra may be a rich load to mine for information about the hot gas in AGNs, and also that the concept of 'the iron K-alpha line' inadequately depicts the spectrum.

Castor, John I.↗

An Integrated Planning Representation Using Macros, Abstractions, and Cases

Planning will be an essential part of future autonomous robots and integrated intelligent systems. This paper focuses on learning problem solving knowledge in planning systems. The system is based on a common representation for macros, abstractions, and cases. Therefore, it is able to exploit both classical and case based techniques. The general operators in a successful plan derivation would be assessed for their potential usefulness, and some stored. The feasibility of this approach was studied through the implementation of a learning system for abstraction. New macros are motivated by trying to improve the operatorset. One heuristic used to improve the operator set is generating operators with more general preconditions than existing ones. This heuristic leads naturally to abstraction hierarchies. This investigation showed promising results on the towers of Hanoi problem. The paper concludes by describing methods for learning other problem solving knowledge. This knowledge can be represented by allowing operators at different levels of abstraction in a refinement.

Baltes, Jacky↗

Global Pattern of Potential Evaporation Calculated from the Penman-Monteith Equation Using Satellite and Assimilated Data

Potential evaporation (E(0)) has been found to be useful in many practical applications and in research for setting a reference level for actual evaporation. All previous estimates of regional or global E(0) are based upon empirical formulae using climatologic meteorologic measurements at isolated stations (i.e., point data). However, the Penman-Monteith equation provides a physically based approach for computing E(0), and by comparing 20 different methods of estimating E(0), Jensen et al. (1990) showed that the Penman-Monteith equation provides the most accurate estimate of monthly E(0) from well-watered grass or alfalfa. In the present study, monthly total E(0) for 24 months (January 1987 to December 1988) was calculated from the Penman-Monteith equation, with prescribed albedo of 0.23 and surface resistance of 70 s/m, which are considered to be representative of actively growing well-watered grass covering the ground. These calculations have been done using spatially representative data derived from satellite observations and data assimilation results. Satellite observations were used to obtain solar radiation, fractional cloud cover, air temperature, and vapor pressure, while four-dimensional data assimilation results were used to calculate the aerodynamic resistance. Meteorologic data derived from satellite observations were compared with the surface measurements to provide a measure of accuracy. The accuracy of the calculated E(0) values was assessed by comparing with lysimeter observations for evaporation from well-watered grass at 35 widely distributed locations, while recognizing that the period of present calculations was not concurrent with the lysimeter measurements and the spatial scales of these measurements and calculations are vastly different. These comparisons suggest that the error in the calculated E(0) values may not be exceeded, on average, 20% for any month or location, but are more likely to be about 15%. These uncertainties are difficult to quantify for mountainous areas or locations close to extensive water bodies. The difference between the calculated and observed E(0) is about 5% when all month and locations were considered. Errors are expected to be less than 15% for averages of E(0) over large areas or several months. Further comparisons with lysimeter observations could provide a better appraisal of the calculated values. Global pattern of E(0) was presented, together with zonal average values.

Choudhury, Bhaskar J.↗

Reactions of SiCl2 and SiHCl with H and Cl Atoms

Calculations have been carried out for the reaction of SiCl2 and SiHCl with H and Cl atoms. In each case, the stationary point geometries and harmonic frequencies were characterized using CASSCF/derivative methods and the cc-pVDZ basis set. Accurate energetics were obtained by combining the CCSD(T) results using the a-cc-pVTZ basis set with an extrapolation to the basis set limit using the a-cc-pVDZ, a-cc-pVTZ, and a-cc-pVQZ basis sets at the MP2 level. The geometries, energetics, and harmonic frequencies were used to obtain rate constants using conventional transition state theory or a Gorin-like model. In each case we find direct abstraction pathways compete with an addition elimination pathway. In the case of SiClH + H the two direct pathways are H abstraction which is barrierless and Cl abstraction with a barrier of 13.5 kcal/mol, while the addition elimination process has a barrier of 26.9 kcal/mol. In the case of SiCl2 + H the direct pathway is Cl abstraction with a barrier of 16.4 kcal/mol, while the addition elimination pathway has a barrier of 29.6 kcal/mol. In the case of SiClH + Cl the direct pathway is H abstraction which is barrierless and the addition elimination pathway has a barrier of 2.0 kcal/mol.

Walch, Stephen P.↗

Microbiological Contamination of Spacecraft

The International Space Station (ISS) Phase1 Program resulted in seven US astronauts residing aboard the Russian Space Station Mir between March 1995 and May 1998. Collaboration between U.S. and Russian scientists consisted of collection and analyses of samples from the crewmembers and the Mir and Shuttle environments before, during, and after missions that lasted from 75 to 209 days in duration. The effects of long-duration space flight on the microbial characteristics of closed life support systems and the interactions of microbes with the spacecraft environment and crewmembers were investigated. Air samples were collected using a Russian or U.S.-supplied sampler (SAS, RCS, or Burkard,) while surface samples were collected using contact slides (Hycon) or swabs. Mir recycled condensate and stored potable water sources were analyzed using the U.S.-supplied Water Experiment Kit. In-flight analysis consisted of enumeration of levels of bacteria and fungi. Amounts of microorganisms seen in the air and on surfaces were mostly within acceptability lin1its; observed temporal fluctuations in levels of microbes probably reflect changes in environmental conditions (e.g., humidity). All Mir galley hot water samples were within the standards set for Mir and the ISS. Microbial isolates were returned to Earth for identification of bacterial and fungal isolates. Crew samples (nose, throat, skin, urine, and feces) were analyzed using methods approved for the medical evaluations of Shuttle flight crews. No significant changes in crew microbiota were found during space flight or upon return relative to preflight results. Dissemination of microbes between the crew and environment was demonstrated by D A fingerprinting. Some biodegradation of spacecraft materials was observed. Accumulation of condensate allowed for the recovery of a wide range of bacteria and fungi as well as some protozoa and dust mites.

Pierson, D. L.↗

Space Mission Concept Development Using Concept Maturity Levels

Over the past five years, pre-project formulation experts at the Jet Propulsion Laboratory (JPL) has developed and implemented a method for measuring and communicating the maturity of space mission concepts. Mission concept development teams use this method, and associated tools, prior to concepts entering their Formulation Phases (Phase A/B). The organizing structure is Concept Maturity Level (CML), which is a classification system for characterizing the various levels of a concept's maturity. The key strength of CMLs is the ability to evolve mission concepts guided by an incremental set of assessment needs. The CML definitions have been expanded into a matrix form to identify the breadth and depth of analysis needed for a concept to reach a specific level of maturity. This matrix enables improved assessment and communication by addressing the fundamental dimensions (e.g., science objectives, mission design, technical risk, project organization, cost, export compliance, etc.) associated with mission concept evolution. JPL's collaborative engineering, dedicated concept development, and proposal teams all use these and other CML-appropriate design tools to advance their mission concept designs. This paper focuses on mission concept's early Pre-Phase A represented by CMLs 1- 4. The scope was limited due to the fact that CMLs 5 and 6 are already well defined based on the requirements documented in specific Announcement of Opportunities (AO) and Concept Study Report (CSR) guidelines, respectively, for competitive missions; and by NASA's Procedural Requirements NPR 7120.5E document for Projects in their Formulation Phase.

Concept Maturity Levels (CMLs)↗

A new database of building-space-specific internal loads and load schedules for performance based code compliance modeling of commercial buildings

Building-level loads and load profiles prescribed by current modeling rules save modelers time and avoid gaming during whole building performance modeling. However, recent studies show that they sometimes insufficiently capture the entire building performance due to the varied loads and load profiles for different space types. As a solution to this issue, this paper develops a database of building-space-specific loads and load profiles used in code compliance modeling. The existing sets of loads and load profiles are reviewed and the challenges behind using them for specific research topics are discussed. Then, the proposed method to develop the building-space-specific loads and load profiles is introduced. After that, the database for these building-space-specific loads and load profiles is presented. In addition, one case is studied to demonstrate the applications of these loads and load profiles. In this case study, three methods are used to develop building energy models: space-specific (using knowledge of the distribution and location of space types and applying the space-specific data in the developed database), building-level (assuming a lack of knowledge of the space types and using the building-level data in the developed database), and calculated-ratio (assuming knowledge of the distribution of space types but not their locations and calculating weighted average values based on the space-specific data in the developed database). Finally, the energy results simulated by using these three methods are compared, which show building-level methods can produce energy results up to 20% different than the space-specific methods. Finally, this paper discusses the application scope and maintenance of this new database.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Application of p-Multigrid to Discontinuous Galerkin Formulations of the Poisson Equation

We investigate p-multigrid as a solution method for several different discontinuous Galerkin (DG) formulations of the Poisson equation. Different combinations of relaxation schemes and basis sets have been combined with the DG formulations to find the best performing combination. The damping factors of the schemes have been determined using Fourier analysis for both one and two-dimensional problems. One important finding is that when using DG formulations, the standard approach of forming the coarse p matrices separately for each level of multigrid is often unstable. To ensure stability the coarse p matrices must be constructed from the fine grid matrices using algebraic multigrid techniques. Of the relaxation schemes, we find that the combination of Jacobi relaxation with the spectral element basis is fairly effective. The results using this combination are p sensitive in both one and two dimensions, but reasonable convergence rates can still be achieved for moderate values of p and isotropic meshes. A competitive alternative is a block Gauss-Seidel relaxation. This actually out performs a more expensive line relaxation when the mesh is isotropic. When the mesh becomes highly anisotropic, the implicit line method and the Gauss-Seidel implicit line method are the only effective schemes. Adding the Gauss-Seidel terms to the implicit line method gives a significant improvement over the line relaxation method.

Helenbrook, B. T.↗