Engineering Papers⌕ Search

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.

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

At least 91 records · Page 5

Radiation pattern synthesis of planar antennas using the iterative sampling method

A synthesis method is presented for determining an excitation of an arbitrary (but fixed) planar source configuration. The desired radiation pattern is specified over all or part of the visible region. It may have multiple and/or shaped main beams with low sidelobes. The iterative sampling method is used to find an excitation of the source which yields a radiation pattern that approximates the desired pattern to within a specified tolerance. In this paper the method is used to calculate excitations for line sources, linear arrays (equally and unequally spaced), rectangular apertures, rectangular arrays (arbitrary spacing grid), and circular apertures. Examples using these sources to form patterns with shaped main beams, multiple main beams, shaped sidelobe levels, and combinations thereof are given.

Stutzman, W. L.↗

Optical implementation of inner product neural associative memory

An optical implementation of an inner-product neural associative memory is realized with a first spatial light modulator for entering an initial two-dimensional N-tuple vector and for entering a thresholded output vector image after each iteration until convergence is reached, and a second spatial light modulator for entering M weighted vectors of inner-product scalars multiplied with each of the M stored vectors, where the inner-product scalars are produced by multiplication of the initial input vector in the first iterative cycle (and thresholded vectors in subsequent iterative cycles) with each of the M stored vectors, and the weighted vectors are produced by multiplication of the scalars with corresponding ones of the stored vectors. A Hughes liquid crystal light valve is used for the dual function of summing the weighted vectors and thresholding the sum vector. The thresholded vector is then entered through the first spatial light modulator for reiteration of the process cycle until convergence is reached.

Liu, Hua-Kuang↗

Massively parallel computation of RCS with finite elements

One of the promising combinations of finite element approaches for scattering problems uses Whitney edge elements, spherical vector wave-absorbing boundary conditions, and bi-conjugate gradient solution for the frequency-domain near field. Each of these approaches may be criticized. Low-order elements require high mesh density, but also result in fast, reliable iterative convergence. Spherical wave-absorbing boundary conditions require additional space to be meshed beyond the most minimal near-space region, but result in fully sparse, symmetric matrices which keep storage and solution times low. Iterative solution is somewhat unpredictable and unfriendly to multiple right-hand sides, yet we find it to be uniformly fast on large problems to date, given the other two approaches. Implementation of these approaches on a distributed memory, message passing machine yields huge dividends, as full scalability to the largest machines appears assured and iterative solution times are well-behaved for large problems. We present times and solutions for computed RCS for a conducting cube and composite permeability/conducting sphere on the Intel ipsc860 with up to 16 processors solving over 200,000 unknowns. We estimate problems of approximately 10 million unknowns, encompassing 1000 cubic wavelengths, may be attempted on a currently available 512 processor machine, but would be exceedingly tedious to prepare. The most severe bottlenecks are due to the slow rate of mesh generation on non-parallel machines and the large transfer time from such a machine to the parallel processor. One solution, in progress, is to create and then distribute a coarse mesh among the processors, followed by systematic refinement within each processor. Elimination of redundant node definitions at the mesh-partition surfaces, snap-to-surface post processing of the resulting mesh for good modelling of curved surfaces, and load-balancing redistribution of new elements after the refinement are auxiliary steps expected to result in a robust low i/o system for very large finite element problems.

Parker, Jay↗

On the measurement of precipitation frequencies by passive microwave radiometery

Over most of the microwave spectrum, raindrops both absorb and scatter radiation producing large changes in brightness temperatures relative to clear or cloudy conditions. Since the structure of rain varies substantially for different rain rates and climatological backgrounds, the raindrop size distribution, the rain layer thickness and the ice clouds above the rain layer are all important inputs to the model computations. The subsequent modeling involves applying the Mie theory to derive the absorption and scattering effects and the radiative transfer calculation is based upon a variational iterative approach which takes account of the multiple scattering effect of the rain layer. Results over both ocean and land backgrounds are demonstrated. It is also demonstrated that by using discrimination tests of the radiometric data, the rain/no rain decision can be made and the rainfall rate can be retrieved from a statistical inversion technique.

Kidder, S. Q.↗

An implicit formulation for the one-dimensional two-phase multi-interface Stefan Problem

In recent years there has been an intensified effort in the analysis of the Stefan Problem. Particular interest is directed towards determining the location of the interface between the solid and liquid phases of the substance being analyzed. A method has been developed using an implicit finite difference scheme which provides an accurate calculation of the phase change interface locations for a problem with simultaneous multiple phases. This is achieved through an iterative technique which requires the phase change energy to be correctly distributed within the nodes which are undergoing a phase transition. The new method predicts the interface locations without requiring time step size control or node size modification at the interfaces.

Williams, S. D.↗

Actuated Hybrid Mirror Telescope

The figure depicts the planned Actuated Hybrid Mirror Telescope (AHMT), which is intended to demonstrate a new approach to the design and construction of wide-aperture spaceborne telescopes for astronomy and Earth science. This technology is also appropriate for Earth-based telescopes. The new approach can be broadly summarized as using advanced lightweight mirrors that can be manufactured rapidly at relatively low cost. More specifically, it is planned to use precise replicated metallic nanolaminate mirrors to obtain the required high-quality optical finishes. Lightweight, dimensionally stable silicon carbide (SiC) structures will support the nanolaminate mirrors in the required surface figures. To enable diffraction- limited telescope performance, errors in surface figures will be corrected by use of mirror-shape-control actuators that will be energized, as needed, by a wave-front-sensing and control system. The concepts of nanolaminate materials and mirrors made from nanolaminate materials were discussed in several previous NASA Tech Briefs articles. Nanolaminates constitute a relatively new class of materials that can approach theoretical limits of stiffness and strength. Nanolaminate mirrors are synthesized by magnetron sputter deposition of metallic alloys and/or compounds on optically precise master surfaces to obtain optical-quality reflector surfaces backed by thin shell structures. As an integral part of the deposition process, a layer of gold that will constitute the reflective surface layer is deposited first, eliminating the need for a subsequent and separate reflective-coating process. The crystallographic textures of the nanolaminate will be controlled to optimize the performance of the mirror. The entire deposition process for making a nanolaminate mirror takes less than 100 hours, regardless of the mirror diameter. Each nanolaminate mirror will be bonded to its lightweight SiC supporting structure. The lightweight nanolaminate mirrors and SiC supporting structures will be fabricated from reusable master molds. The mirror-shape-control actuators will be low-power, high-capacitance lead magnesium niobate electrostrictive actuators that will be embedded in the SiC structures. The mode of operation of these actuators will be such that once power was applied, they will change in length and once power was removed, they will maintain dimensional stability to nanometer precision. This mode of operation will enable the use of low-power, minimally complex electronic control circuitry. The wave-front-sensing and control system will be designed and built according to a two-stage architecture. The first stage will be implemented by a Shack- Hartmann (SH) sensor subsystem, which will provide a large capture range. The second, higher-performance stage will be implemented by an image-based wave-front-sensing subsystem that will include a phase-retrieval camera (PRC), and will utilize phase retrieval and other techniques to measure wavefront error directly. Phase retrieval is a process in which multiple images of an unresolved object are iterated to estimate the phase of the optical system that acquired the images. The combination of SH and phase-retrieval sensors will afford the virtues of both a dynamic range of 105 and an accuracy of <10 nm.

Hickey, Gregory↗

Compact Microscope Imaging System with Intelligent Controls

The figure presents selected views of a compact microscope imaging system (CMIS) that includes a miniature video microscope, a Cartesian robot (a computer- controlled three-dimensional translation stage), and machine-vision and control subsystems. The CMIS was built from commercial off-the-shelf instrumentation, computer hardware and software, and custom machine-vision software. The machine-vision and control subsystems include adaptive neural networks that afford a measure of artificial intelligence. The CMIS can perform several automated tasks with accuracy and repeatability . tasks that, heretofore, have required the full attention of human technicians using relatively bulky conventional microscopes. In addition, the automation and control capabilities of the system inherently include a capability for remote control. Unlike human technicians, the CMIS is not at risk of becoming fatigued or distracted: theoretically, it can perform continuously at the level of the best human technicians. In its capabilities for remote control and for relieving human technicians of tedious routine tasks, the CMIS is expected to be especially useful in biomedical research, materials science, inspection of parts on industrial production lines, and space science. The CMIS can automatically focus on and scan a microscope sample, find areas of interest, record the resulting images, and analyze images from multiple samples simultaneously. Automatic focusing is an iterative process: The translation stage is used to move the microscope along its optical axis in a succession of coarse, medium, and fine steps. A fast Fourier transform (FFT) of the image is computed at each step, and the FFT is analyzed for its spatial-frequency content. The microscope position that results in the greatest dispersal of FFT content toward high spatial frequencies (indicating that the image shows the greatest amount of detail) is deemed to be the focal position.

McDowell, Mark↗

Calculation of the angular radiance distribution for a coupled atmosphere and canopy

The radiative transfer equations for a coupled atmosphere and canopy are solved numerically by an improved Gauss-Seidel iteration algorithm. The radiation field is decomposed into three components: unscattered sunlight, single scattering, and multiple scattering radiance for which the corresponding equations and boundary conditions are set up and their analytical or iterational solutions are explicitly derived. The classic Gauss-Seidel algorithm has been widely applied in atmospheric research. This is its first application for calculating the multiple scattering radiance of a coupled atmosphere and canopy. This algorithm enables us to obtain the internal radiation field as well as radiances at boundaries. Any form of bidirectional reflectance distribution function (BRDF) as a boundary condition can be easily incorporated into the iteration procedure. The hotspot effect of the canopy is accommodated by means of the modification of the extinction coefficients of upward single scattering radiation and unscattered sunlight using the formulation of Nilson and Kuusk. To reduce the computation for the case of large optical thickness, an improved iteration formula is derived to speed convergence. The upwelling radiances have been evaluated for different atmospheric conditions, leaf area index (LAI), leaf angle distribution (LAD), leaf size and so on. The formulation presented in this paper is also well suited to analyze the relative magnitude of multiple scattering radiance and single scattering radiance in both the visible and near infrared regions.

Liang, Shunlin↗

HITEMP Material and Structural Optimization Technology Transfer

The feasibility of adding viscoelasticity and the Generalized Method of Cells (GMC) for micromechanical viscoelastic behavior into the commercial HyperSizer structural analysis and optimization code was investigated. The viscoelasticity methodology was developed in four steps. First, a simplified algorithm was devised to test the iterative time stepping method for simple one-dimensional multiple ply structures. Second, GMC code was made into a callable subroutine and incorporated into the one-dimensional code to test the accuracy and usability of the code. Third, the viscoelastic time-stepping and iterative scheme was incorporated into HyperSizer for homogeneous, isotropic viscoelastic materials. Finally, the GMC was included in a version of HyperSizer. MS Windows executable files implementing each of these steps is delivered with this report, as well as source code. The findings of this research are that both viscoelasticity and GMC are feasible and valuable additions to HyperSizer and that the door is open for more advanced nonlinear capability, such as viscoplasticity.

Collier, Craig S.↗

Abstract-Reasoning Software for Coordinating Multiple Agents

A computer program for scheduling the activities of multiple agents that share limited resources has been incorporated into the Automated Scheduling and Planning Environment (ASPEN) software system, aspects of which have been reported in several previous NASA Tech Briefs articles. In the original intended application, the agents would be multiple spacecraft and/or robotic vehicles engaged in scientific exploration of distant planets. The program could also be used on Earth in such diverse settings as production lines and military maneuvers. This program includes a planning/scheduling subprogram of the iterative repair type that reasons about the activities of multiple agents at abstract levels in order to greatly improve the scheduling of their use of shared resources. The program summarizes the information about the constraints on, and resource requirements of, abstract activities on the basis of the constraints and requirements that pertain to their potential refinements (decomposition into less-abstract and ultimately to primitive activities). The advantage of reasoning about summary information is that time needed to find consistent schedules is exponentially smaller than the time that would be needed for reasoning about the same tasks at the primitive level.

Clement, Bradley↗

Multiple scattering from finite inhomogeneous media

Utilizing the characteristic information concerning the apparent phase constant difference between the electric and magnetic fields propagating in an inhomogeneous medium, a theoretical analysis of the multiple scattering of electromagnetic waves in finite inhomogeneous media is presented. The solution is obtained by first approximating the coefficients of a pair of exact coupled first-order differential equations and then solving the equations by first-order iteration. The present first-order approximate solution with multiple scattering considerations is shown to be more accurate than the WKB solution. Methods to improve the accuracy of the first-order solution further are discussed. Application of the solution to slowly varying finite media with periodic properties demonstrates the validity of the solution. The same approach can be extended to frequencies in the optical region by retaining additional terms in the coefficients of the coupled differential equations.

Tang, C. C. H.↗

Performance Analysis of Distributed Object-Oriented Applications

The purpose of this research was to evaluate the efficiency of a distributed simulation architecture which creates individual modules which are made self-scheduling through the use of a message-based communication system used for requesting input data from another module which is the source of that data. To make the architecture as general as possible, the message-based communication architecture was implemented using standard remote object architectures (Common Object Request Broker Architecture (CORBA) and/or Distributed Component Object Model (DCOM)). A series of experiments were run in which different systems are distributed in a variety of ways across multiple computers and the performance evaluated. The experiments were duplicated in each case so that the overhead due to message communication and data transmission can be separated from the time required to actually perform the computational update of a module each iteration. The software used to distribute the modules across multiple computers was developed in the first year of the current grant and was modified considerably to add a message-based communication scheme supported by the DCOM distributed object architecture. The resulting performance was analyzed using a model created during the first year of this grant which predicts the overhead due to CORBA and DCOM remote procedure calls and includes the effects of data passed to and from the remote objects. A report covering the distributed simulation software and the results of the performance experiments has been submitted separately. The above report also discusses possible future work to apply the methodology to dynamically distribute the simulation modules so as to minimize overall computation time.

Schoeffler, James D.↗

Kepler Data Validation II–Transit Model Fitting and Multiple-Planet Search

This paper discusses the transit model-fitting and multiple-planet search algorithms and performance of the Kepler Science Data Processing Pipeline, developed by the Kepler Science Operations Center (SOC). Threshold crossing events (TCEs), which are transit candidate events, are generated by the Transiting Planet Search (TPS) component of the pipeline and subsequently processed in the data validation (DV) component. The transit model is used in DV to fit TCEs to characterize planetary candidates and to derive parameters that are used in various diagnostic tests to classify them. After the signature associated with the TCE is removed from the light curve of the target star, the residual light curve goes through TPS again to search for additional TCEs. The iterative process of transit model fitting and multiple-planet search continues until no TCE is generated from the residual light curve or an upper limit is reached. The transit model-fitting and multiple-planet search performance of the final release (9.3, 2016January) of the pipeline is demonstrated with the results of the processing of four years (17 quarters) of flight data from the primary Kepler Mission. The transit model-fitting results are accessible from the NASA Exoplanet Archive. The final version of the SOC codebase is available through GitHub.

Threshold crossing events (TCEs↗

Computation of optimal output-feedback compensators for linear time-invariant systems

The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.

Platzman, L. K.↗

A Multilevel Approach For SolvingLarge-Scale QUBO Problems With Noisy Hybrid Quantum Approximate Optimization

Quantum approximate optimization is one ofthe promising candidates for useful quantum computation,particularly in the context of finding approximate solutionsto Quadratic Unconstrained Binary Optimization (QUBO)problems. However, the existing quantum processing units(QPUs) are of relatively small size, and canonical mappingsof QUBO via the Ising model require one qubit per vari-able, rendering direct large-scale optimization infeasible.In classical optimization, a general strategy for addressingmany large-scale problems is via multilevel/multigrid meth-ods, where the large target problem is iteratively coarsenedand the global solution is constructed from multiple small-scale optimization runs. In this work, we experimentallytest how existing QPUs perform when used as a sub-solverwithin such a multilevel strategy. To this aim, we com-bine and extend (via additional classical processing steps)the recently proposed Noise-Directed Adaptive Remapping(NDAR) and Quantum Relax&Round (QRR) algorithms.We first demonstrate the effectiveness of our heuristicextensions on Rigetti’s superconducting transmon deviceAnkaa-2. We find approximate solutions to10instances offully connected82-qubit Sherrington-Kirkpatrick graphswith random integer-valued coefficients obtaining normal-ized approximation ratios (ARs) in the range∼0.98−1.0,and the same class with real-valued coefficients (ARs∼0.94−1.0). Then, we implement the extended NDAR andQRR algorithms as subsolvers in the multilevel algorithmfor6large-scale graphs with at most∼27,000variables.In practice, the QPU (with classical post-processing steps)is used to find approximate solutions to dozens of at most82-qubit problems, which are iteratively used to constructthe global solution. We observe that quantum optimizationresults are competitive in terms of the quality of solutionswhen compared to classical heuristics used as subsolverswithin the multilevel approach.Reproducibility: source code and data are available at[TBA upon acceptance]

quantum computing↗

Compensator improvement for multivariable control systems

A theory and the associated numerical technique are developed for an iterative design improvement of the compensation for linear, time-invariant control systems with multiple inputs and multiple outputs. A strict constraint algorithm is used in obtaining a solution of the specified constraints of the control design. The result of the research effort is the multiple input, multiple output Compensator Improvement Program (CIP). The objective of the Compensator Improvement Program is to modify in an iterative manner the free parameters of the dynamic compensation matrix so that the system satisfies frequency domain specifications. In this exposition, the underlying principles of the multivariable CIP algorithm are presented and the practical utility of the program is illustrated with space vehicle related examples.

Mitchell, J. R.↗

Research in computer science

Synopses are given for NASA supported work in computer science at the University of Virginia. Some areas of research include: error seeding as a testing method; knowledge representation for engineering design; analysis of faults in a multi-version software experiment; implementation of a parallel programming environment; two computer graphics systems for visualization of pressure distribution and convective density particles; task decomposition for multiple robot arms; vectorized incomplete conjugate gradient; and iterative methods for solving linear equations on the Flex/32.

Ortega, J. M.↗

Accurate orbit determination strategies for the tracking and data relay satellites

The National Aeronautics and Space Administration (NASA) has developed the Tracking and Data Relay Satellite (TDRS) System (TDRSS) for tracking and communications support of low Earth-orbiting satellites. TDRSS has the operational capability of providing 85% coverage for TDRSS-user spacecraft. TDRSS currently consists of five geosynchronous spacecraft and the White Sands Complex (WSC) at White Sands, New Mexico. The Bilateration Ranging Transponder System (BRTS) provides range and Doppler measurements for each TDRS. The ground-based BRTS transponders are tracked as if they were TDRSS-user spacecraft. Since the positions of the BRTS transponders are known, their radiometric tracking measurements can be used to provide a well-determined ephemeris for the TDRS spacecraft. For high-accuracy orbit determination of a TDRSS user, such as the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft, high-accuracy TDRS orbits are required. This paper reports on successive refinements in improved techniques and procedures leading to more accurate TDRS orbit determination strategies using the Goddard Trajectory Determination System (GTDS). These strategies range from the standard operational solution using only the BRTS tracking measurements to a sophisticated iterative process involving several successive simultaneous solutions for multiple TDRSs and a TDRSS-user spacecraft. Results are presented for GTDS-generated TDRS ephemerides produced in simultaneous solutions with the TOPEX/Poseidon spacecraft. Strategies with different user spacecraft, as well as schemes for recovering accurate TDRS orbits following a TDRS maneuver, are also presented. In addition, a comprehensive assessment and evaluation of alternative strategies for TDRS orbit determination, excluding BRTS tracking measurements, are presented.

Oza, D. H.↗