SEARCH · Engineering Papers

Results for “multiple iterations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A deep learning-enhanced framework for multiphysics joint inversion

Joint inversion has drawn considerable attention due to the availability of multiple geophysical data sets, ever-increasing computational resources, the development of advanced algorithms, and its ability to reduce inversion uncertainty. A key issue of joint inversion is to develop effective strategies to link different geophysical data in a unified mathematical framework, in which the information obtained from different models can complement each other. We have developed a deep learning-enhanced joint inversion framework to simultaneously reconstruct different physical models by fusing different types of geophysical data. Traditionally, structure similarity constraints are pursued by joint inversion algorithms using manually crafted formulations (e.g., cross gradient). The constraint is constructed by a deep neural network (DNN) during the learning process. The framework is designed to combine the DNN and a traditional independent inversion workflow and improve the joint inversion result iteratively. The network can be easily extended to incorporate multiphysics without structural changes. Numerical experiments on the joint inversion of 2D DC resistivity data and seismic traveltime are used to validate our method. In addition, this learning-based framework demonstrates excellent generalization abilities when tested on data sets using different geologic structures. It also can handle different sensing configurations and nonconforming discretization.

Geochemistry & Geophysics↗

DOE OSTI2025

Tailoring optical response of MXene thin films

Due to their attractive optical properties, 2D MXenes have garnered interest in nanophotonic and optoelectronic applications. However, tuning their properties typically requires the iterative synthesis of MXenes with a specific set of properties, such as the absorption band position, electronic conductivity, and dielectric constant. We demonstrate how to tailor the optical properties of MXene thin films over a broad 1500-nm wavelength range by mixing different ratios of highly conductive Ti 3 C 2 T x with poorly conductive Nb 2 CT x . By changing the MXene film composition, the epsilon-near-zero (ENZ) point, where the optical properties transit from dielectric to metallic, was varied in the spectral range from 1.1 to 2.6 µm. Additionally, we observed a reduction in absorption in some compositions compared to the absorption of the pure MXene films. Compared to other methods, this approach enables simple and continual tuning of MXene optical properties without requiring multiple time-consuming synthesis steps.

2D materials↗

DOE OSTI2023

Quantifying Optical Loss of High-Voltage Degradation Modes in PV Modules Using Spectral Analysis [Slides]

The direct current bias for photovoltaic (PV) modules interconnected in series-strings may include both high voltage negative ("HV-") and positive ("HV+") polarity with respect to the electrical ground. Multiple degradation modes, resulting in quantifiable optical loss, were found to occur during HV-/HV+ sequential stress, including: corrosion of the external glass surface, encapsulant delamination (at its interfaces with the glass and the PV cell), internal haze formation (resulting from a chemical interaction between the glass and the encapsulant), corrosion and migration of the gridlines, and corrosion of the silicon nitride (SixNy) antireflective coating on the cell. The effects of these separate modes were examined using monolithic (e.g., glass or PV cell) and laminated-coupon (glass/encapsulant/glass or glass/encapsulant/cell/encapsulant/backsheet) specimens. Characterizations during and after unbiased accelerated testing at 85 degrees C/85% relative humidity included: spectrophotometry, optical microscopy, electron microscopy, and ellipsometry. For some module components (i.e., the glass and the SixNy coating), the optical performance was determined through iterative analysis of empirical measurements. Concentrating on just their spectral effect, a novel model was then developed to estimate the transfer of light to the PV cell and the return of light from the PV module with simultaneous degradation mechanisms, which was compared to a mini-module previously subjected to HV-/HV+ stress. The model suggests one-third of the current loss observed for the mini-module can be attributed to the optical degradation of the packaging materials. The dominant degradation modes include encapsulant delamination and corrosion of the SixNy coating. Recommendations are given so that the optical model may be improved relative to accelerated testing and validated relative to field aging.

14 SOLAR ENERGY↗

NASA NTRS1993

Using neural networks to assist in OPAD data analysis

The space shuttle main engine (SSME) became the subject of plume emission spectroscopy in 1986 when researchers from NASA-Marshall Space Flight Center (MSFC), Arnold Engineering Development Center (AEDC), and Rocketdyne went to the SSME test stands at the NASA-Stennis Space Center and at Rocketdyne's Santa Susan Field Laboratory to optically observe the plume. Since then, plume spectral acquisitions have recorded many nominal tests and the qualitative spectral features of the SSME plume are now well established. Significant discoveries made with both wide-band and narrow-band plume emission spectroscopy systems led MSFC to promote the Optical Plume Anomaly Detection (OPAD) program with a goal of instrumenting all SSME test stands with customized spectrometer systems. A prototype OPAD system is now installed on the SSME Technology Test Bed (TTB) at MSFC. The OPAD system instrumentation consists of a broad-band, optical multiple-channel analyzer (OMA) and a narrow-band device called a polychrometer. The OMA is a high-resolution (1.5-2.0 Angstroms) 'super-spectrometer' covering the near-ultraviolet to near-infrared waveband (2800-7400 Angstroms), providing two scans per second. The polychrometer consists of sixteen narrow-band radiometers: fourteen monitoring discrete wavelengths of health and condition monitoring elements and two dedicated to monitoring background emissions. All sixteen channels are capable of providing 500 samples per second. To date, the prototype OPAD system has been used during 43 SSME firings on the TTB, collecting well over 250 megabytes of plume spectral data. One goal of OPAD data analysis is to determine interatively with the help of a computer code, SPECTRA4, developed at AEDC. Experience has shown that iteration with SPECTRA4 is an incredibly labor intensive task and not one to be performed by band. What is really needed is the 'inverse' of SPECTRA4 but the mathematical model for this inverse mapping is tenuous at best. However, the robustness of PSECTRA4 run in the 'forward' direction means that accurate input/output mappings can be obtained. If the mappings were inverted (i.e., input becomes output and output becomes input) then an 'inverse' of SPECTRA4 would be at hand but the 'model' would be specific to the data utilized and would in no way be general. Building a generalized model based upon known input/output mappings while ignoring the details of the governing physical model is possible through the use of a neural network. The research investigation described involves the development of a neural network to provide a generalized 'inverse' of SPECTRA4. The objectives of the research were to design an appropriate neural network architecture, train the network, and then evaluate its performance.

Whitaker, Kevin W.↗

NASA NTRSDate not supplied

In-Depth Analysis of Subject Functional Performance and Subjective Feedback Data for Evaluating Argos Simulation Quality

Mechanical offloading systems such as the Active Response Gravity Offload System (ARGOS) at NASA Johnson Space Center (JSC) are used to simulate various partial gravity environments. ARGOS offloads pressurized suited subjects through a gimbaling pivot point system attached to a volumetric Portable Life Support System (PLSS). The pivot point can be configured to multiple locations with respect to the human-suit system center of gravity (CG). Previous work has tested and documented the interaction between different pivot point locations and measured CG. Observations indicated that small changes in pivot point location substantially affect the stability and functional performance (simulation quality) of a subject. This led to the development of a standardized gimbal assessment protocol to assess the functional performance of a pressurized suited subject as a function of the pivot point location, using the data obtained from a variety of planetary suit mobility tasks. As a result, an expansive repository of CG-related task metrics was generated from numerous ARGOS spacesuit test events, for different gimbal configurations and different test subjects. Based on the evaluation metrics and observed trends, several pivot point locations were iteratively identified to determine an “appropriate/optimal” configuration. This study thus describes the in-depth analysis of the CG evaluation metrics to better quantify trends and determine if specific factors are strongly associated with appropriate/optimal configurations. Quantitative and qualitative variables relating to subject task performance, subjective feedback, and anthropometry will be evaluated via probabilistic methods. The results from this study are expected to improve our understanding of optimal ARGOS gimbal settings, which will better inform the gimbal configuration identification process and improve simulation quality for extravehicular activity (EVA) training.

Joseph Yao↗

DOE OSTI2023

Quantifying optical loss of high-voltage degradation modes in photovoltaic modules using spectral analysis

The direct current bias for photovoltaic (PV) modules interconnected in series-strings may include both high voltage negative (“HV-”) and positive (“HV+”) polarity with respect to the electrical ground. Multiple degradation modes, resulting in quantifiable optical loss, were found to occur during HV-/HV+ sequential stress, including corrosion of the external glass surface, encapsulant delamination (at its interfaces with the glass and the PV cell), internal haze formation (resulting from a chemical interaction between the glass and the encapsulant), corrosion and migration of the gridlines, and corrosion of the silicon nitride (Si x N y ) antireflective coating on the cell. The effects of these separate modes were examined using monolithic (e.g., glass or PV cell) and laminated-coupon (glass/encapsulant/glass or glass/encapsulant/cell/encapsulant/backsheet) specimens. Characterizations during and after unbiased accelerated testing at 85°C/85% relative humidity included spectrophotometry, optical microscopy, electron microscopy, and ellipsometry. For some module components (i.e., the glass and the Si x N y coating), the optical performance was determined through iterative analysis of empirical measurements. Concentrating on just their spectral effect, a novel model was then developed to estimate the transfer of light to the PV cell and the return of light from the PV module with simultaneous degradation mechanisms, which was compared with a mini-module previously subjected to HV-/HV+ stress. Here the model suggests that one third of the current loss observed for the mini-module can be attributed to the optical degradation of the packaging materials. The dominant degradation modes include encapsulant delamination and corrosion of the Si x N y coating. Recommendations are given so that the optical model may be improved relative to accelerated testing and validated relative to field aging.

14 SOLAR ENERGY↗

DOE OSTI2024

Three-dimensional cooperative inversion of airborne magnetic and gravity gradient data using deep-learning techniques

Using multiple geophysical methods has become a prevailing approach in numerous geophysical applications to investigate subsurface structures and parameters. These multimethod-based exploration strategies have the potential to greatly diminish uncertainties and ambiguities encountered during geophysical data analysis and interpretation. One of the applications is the cooperative inversion of airborne magnetic and gravity gradient data for the interpretation of data obtained in mineral, oil and gas, and geothermal explorations. In this paper, a unified cooperative inversion framework is designed by combining the standard separate inversions with a deep neural network (DNN), which serves as the link between different types of data. A well-trained DNN takes the separately inverted susceptibility and density models as the inputs and provides improved models that will be used as the initial models of deterministic inversions. A two-round iteration strategy is adopted to guarantee the reasonability of the recovered models and overall efficiency of the inversion. In addition, this deep-learning (DL)-based framework demonstrates excellent generalization abilities when tested on models that are entirely distinct from the training data sets. The framework can easily incorporate multiphysics without necessitating any structural changes to the network. Synthetic experiments validate that our DL-based method outperforms conventional separate inversions and cross-gradient-based joint inversion in view of the accuracy of the recovered models and inversion efficiency. Successful application to field data further verifies the effectiveness of our DL-based method.

Geochemistry & Geophysics↗

NASA NTRS2017

Design and Development of the Heat Redistribution System for the Mission to Europa Spacecraft:

This paper details the conceptual design of the Heat Redistribution System (HRS) used for thermal control of the Europa Clipper spacecraft. The solar powered spacecraft will launch in the early 2020s and will study Jupiter’s Icy moon, Europa, where previous investigations have indicated the possibility of a liquid water ocean underneath its ice-encrusted surface. The spacecraft will execute multiple flybys of Europa while relying on the HRS for thermal control of the spacecraft bus. The HRS utilizes a Mechanically Pumped Fluid Loop (MPFL) to reclaim waste heat from an avionics and payload electronics compartment and redistribute it to other parts of the spacecraft, such as the propulsion module structure and the louver-radiator assembly. A “replacement heater block” provides supplemental heat when needed to maintain the spacecraft hardware above its allowable flight temperatures. Additionally, a set of thermal control valves autonomously actuate to reduce the fluid flow through the radiator during cold environment cases. The HRS was used in previous flight projects at JPL such as Mars Science Laboratory, but the uniqueness of the Europa Clipper mission has mandated changes to previous HRS designs; the most significant design drivers and design changes for this iteration of the HRS are discussed.

Bertagne, Christopher↗

DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗

NASA NTRS2018

Formulation and Implementation of Inflow/Outflow Boundary Conditions to Simulate Propulsive Effects

Boundary conditions appropriate for simulating flow entering or exiting the computational domain to mimic propulsion effects have been implemented in an adaptive Cartesian simulation package. A robust iterative algorithm to control mass flow rate through an outflow boundary surface is presented, along with a formulation to explicitly specify mass flow rate through an inflow boundary surface. The boundary conditions have been applied within a mesh adaptation framework based on the method of adjoint-weighted residuals. This allows for proper adaptive mesh refinement when modeling propulsion systems. The new boundary conditions are demonstrated on several notional propulsion systems operating in flow regimes ranging from low subsonic to hypersonic. The examples show that the prescribed boundary state is more properly imposed as the mesh is refined. The mass-flowrate steering algorithm is shown to be an efficient approach in each example. To demonstrate the boundary conditions on a realistic complex aircraft geometry, two of the new boundary conditions are also applied to a modern low-boom supersonic demonstrator design with multiple flow inlets and outlets.

Boundary↗

NASA NTRSDate not supplied

Integrated In-Situ Resource Utilization Modeling of a Lunar Water Processing System

A key element of achieving a sustained surface presence, such as defined in NASA’s Artemis plan, is In-Situ Resource Utilization (ISRU). ISRU is the practice of using local resources to provide mission consumables that reduce system launch mass requirements, and regenerate resources (chiefly, water and oxygen) for propulsion and life support supporting both Lunar and Martian missions. ISRU systems require multiple complex processes, such as excavation, chemical reactors, and electrolysis subsystems that must operate in harmony to optimize the overall system process from beginning to end. The Mission Analysis and Integration Tool (MAIT) connects individual subsystem models into a customized, flexible framework for the purpose of technology downselect, optimization, and end-to-end process planning. MATLAB was favorable to use as the main software integration tool due to its ability to communicate with a vast number of other programming languages and makes up the backbone of data flow between inputs and outputs to the subsystem models. MAIT initially evaluated a suite of technologies, including but not limited to, the water processing Lunar Auger Dryer for ISRU (LADI) system and an oxygen extraction carbothermal reduction process. With individual models consolidated, the MAIT tool generated over 60,000 cases during its parametric sweeps; these system iterations produced valuable insight into the optimal LADI geometry for minimizing heater energy demands, estimating carbothermal reactor and radiator mass relationships, and calculated the power dynamics of the electrolysis unit. Efforts to update the MAIT tool are ongoing to accommodate and scale ISRU technologies supporting the Space Technology Mission Directorate’s (STMD) commercialization strategy, and increase the MAIT software capability to handle a wide array of ISRU system models beyond the Lunar environment, e.g. production of propellant for a Martian lander.

Avery Carlson↗

NASA NTRS2016

SLS Flight Software Testing: Using a Modified Agile Software Testing Approach

NASA's Space Launch System (SLS) is an advanced launch vehicle for a new era of exploration beyond earth's orbit (BEO). The world's most powerful rocket, SLS, will launch crews of up to four astronauts in the agency's Orion spacecraft on missions to explore multiple deep-space destinations. Boeing is developing the SLS core stage, including the avionics that will control vehicle during flight. The core stage will be built at NASA's Michoud Assembly Facility (MAF) in New Orleans, LA using state-of-the-art manufacturing equipment. At the same time, the rocket's avionics computer software is being developed here at Marshall Space Flight Center in Huntsville, AL. At Marshall, the Flight and Ground Software division provides comprehensive engineering expertise for development of flight and ground software. Within that division, the Software Systems Engineering Branch's test and verification (T&V) team uses an agile test approach in testing and verification of software. The agile software test method opens the door for regular short sprint release cycles. The idea or basic premise behind the concept of agile software development and testing is that it is iterative and developed incrementally. Agile testing has an iterative development methodology where requirements and solutions evolve through collaboration between cross-functional teams. With testing and development done incrementally, this allows for increased features and enhanced value for releases. This value can be seen throughout the T&V team processes that are documented in various work instructions within the branch. The T&V team produces procedural test results at a higher rate, resolves issues found in software with designers at an earlier stage versus at a later release, and team members gain increased knowledge of the system architecture by interfacing with designers. SLS Flight Software teams want to continue uncovering better ways of developing software in an efficient and project beneficial manner. Through agile testing, there has been increased value through individuals and interactions over processes and tools, improved customer collaboration, and improved responsiveness to changes through controlled planning. The presentation will describe agile testing methodology as taken with the SLS FSW Test and Verification team at Marshall Space Flight Center.

Bolton, Albanie T.↗

NASA NTRS2014

Development of a Route Crossing Tool for Shared Airspace Environments

In current-day Terminal Radar Approach Control (TRACON) operations, departure and arrival controllers maintain separate and dedicated airspace for their respective traffic flows. Although this practice has obvious safety features, it also leads to inefficiencies; for example, departure aircraft may be routinely capped beneath arrival airspace. With the right decision-support and coordination tools, departures could continue to climb through arrival airspace when sufficient gaps exist. Previous studies of ‘shared airspace’ have examined pre-arranged coordination procedures, as well as tools that gave feedback to the controllers on where gaps between arrivals were located and whether the departure aircraft could be scheduled to fly through those gaps [1, 2, 3, 4]. Since then, the Route Crossing Tool (RCT) has been developed to allow controllers to assess multiple pre-defined route options at points where the arrivals and departures cross, thereby increasing the possibility of climbing a departure through an arrival gap.The RCT aids in ensuring lateral separation between departure and arrival aircraft that pass through the same altitude. Since the RCT can be applied tactically, it can enable aircraft to fly through arrival flows even if these aircraft depart outside scheduled times. The RCT makes use of a set of predefined parallel departure routes crossing the arrival flow at equidistant intersecting points on the arrival route. The RCT uses the Estimated Time of Arrival (ETA) of the departure aircraft at each intersecting point to calculate the lateral separation with the neighboring arrivals when it crosses that point; this information is graphically displayed to the controller. Additionally, the RCT incorporates forecast winds in its ETA predictions.Multiple prototypes of the RCT have been iteratively developed with feedback from Subject Matter Experts (SMEs). This paper presents the final design, the design process, and lessons learned. Initial results from a simulation suggest that the tool was successful in helping controllers to safely climb more aircraft. Controller feedback on the tool was also positive.

shared airspace↗

DOE OSTI2021

Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

We present a topology-based method for mesh-partitioning in three-dimensional discrete fracture network (DFN) simulations that takes advantage of the intrinsic multi-level nature of a DFN. DFN models are used to simulate flow and transport through low-permeability fractured media in the subsurface by explicitly representing fractures as discrete entities. The governing equations for flow and transport are numerically integrated on computational meshes generated on the interconnected fracture networks. Modern high-fidelity DFN simulations require high-performance computing on multiple processors where performance and scalability depends partially on obtaining a high-quality partition of the mesh to balance work-loads and minimize communication across all processors. The discrete structure of a DFN naturally lends itself to various graph representations, which can be thought of as coarse-scale representations of the computational mesh. Using this concept, we develop two applications of the multilevel graph partitioning algorithm to partition the mesh of a DFN. In the first, we project a partition of the graph based on the DFN topology onto the mesh of the DFN and in the second, this DFN-based projection is used as the initial condition for further partitioning refinement of the mesh. We compare the performance of these methods with standard multi-level graph partitioning using graph-based metrics (cut, imbalance, partitioning time), computational-based metrics (FLOPS, iterations, solver time), and total run time. The DFN-based and the mesh-based partitioning methods are comparable in terms of the graph-based metrics, but the time required to obtain the partition is several orders of magnitude faster using the DFN-based partitions. The computation-based metrics show comparable performance between both methods so, in combination, the DFN-based partitions are several orders of magnitude faster than the mesh-based partition. Furthermore, the method which uses the DFN-partition solution as the initial condition of the mesh partition provided cut and imbalance values that were close to the mesh-based partition but in a fraction of the time. In turn, this hybrid method outperformed both of the other methods in terms of the total run time.

58 GEOSCIENCES↗

NASA NTRS2024

Comparison of Multisine Peak Factor Minimization Algorithms for Aircraft System Identification

Two phase-optimized multisine peak factor minimization algorithms are presented and evaluated. The first algorithm minimizes peak factor by iteratively clipping the peaks of generated multisine signals. The second algorithm optimizes peak factor indirectly through minimization of an approximation of the infinity norm of the multisine. Algorithm performance was evaluated as a function of different signal properties, including the number of harmonics, harmonic spacing, and number of snow harmonics (extra harmonics included for further reduction of the peak factor). The two algorithms are compared against results obtained by minimizing peak factor directly using a simplex algorithm, which has been a common approach when designing phase-optimized multisines for system identification flight tests. Sample results show that the clipping and infinity norm algorithms produced multisine signals with comparable peak factors that were lower than that of the simplex algorithm. However, the clipping algorithm runs an order of magnitude faster than the other two algorithms, which also makes it practical to repeat the algorithm multiple times to achieve even lower peak factors.

system identification↗

NASA NTRSDate not supplied

Comparison of Multisine Peak Factor Minimization Algorithms for Aircraft System Identification(Presentation)

flight test↗

NASA NTRS2014

Explaining Earths Energy Budget: CERES-Based NASA Resources for K-12 Education and Public Outreach

Among atmospheric scientists, the importance of the Earth radiation budget concept is well understood. Papers have addressed the topic for over 100 years, and the large Clouds and the Earth's Radiant Energy System (CERES) science team (among others), with its multiple on-orbit instruments, is working hard to quantify the details of its various parts. In education, Earth's energy budget is a concept that generally appears in middle school and Earth science curricula, but its treatment in textbooks leaves much to be desired. Students and the public hold many misconceptions, and very few people have an appreciation for the importance of this energy balance to the conditions on Earth. More importantly, few have a correct mental model that allows them to make predictions and understand the effect of changes such as increasing greenhouse gas concentrations. As an outreach element of the core CERES team at NASA Langley, a multi-disciplinary group of scientists, educators, graphic artists, writers, and web developers has been developing and refining graphics and resources to explain the Earth's Energy budget over the last few decades. Resources have developed through an iterative process involving ongoing use in front of a variety of audiences, including students and teachers from 3rd to 12th grade as well as public audiences.

Chambers, L. H.↗

NASA NTRS1995

Development of Fast Algorithms Using Recursion, Nesting and Iterations for Computational Electromagnetics

In the first phase of our work, we have concentrated on laying the foundation to develop fast algorithms, including the use of recursive structure like the recursive aggregate interaction matrix algorithm (RAIMA), the nested equivalence principle algorithm (NEPAL), the ray-propagation fast multipole algorithm (RPFMA), and the multi-level fast multipole algorithm (MLFMA). We have also investigated the use of curvilinear patches to build a basic method of moments code where these acceleration techniques can be used later. In the second phase, which is mainly reported on here, we have concentrated on implementing three-dimensional NEPAL on a massively parallel machine, the Connection Machine CM-5, and have been able to obtain some 3D scattering results. In order to understand the parallelization of codes on the Connection Machine, we have also studied the parallelization of 3D finite-difference time-domain (FDTD) code with PML material absorbing boundary condition (ABC). We found that simple algorithms like the FDTD with material ABC can be parallelized very well allowing us to solve within a minute a problem of over a million nodes. In addition, we have studied the use of the fast multipole method and the ray-propagation fast multipole algorithm to expedite matrix-vector multiplication in a conjugate-gradient solution to integral equations of scattering. We find that these methods are faster than LU decomposition for one incident angle, but are slower than LU decomposition when many incident angles are needed as in the monostatic RCS calculations.

Chew, W. C.↗