Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process coupling”

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 1,639 records · Page 91

Parallel eigenanalysis of finite element models in a completely connected architecture

A parallel algorithm is presented for the solution of the generalized eigenproblem in linear elastic finite element analysis, (K)(phi) = (M)(phi)(omega), where (K) and (M) are of order N, and (omega) is order of q. The concurrent solution of the eigenproblem is based on the multifrontal/modified subspace method and is achieved in a completely connected parallel architecture in which each processor is allowed to communicate with all other processors. The algorithm was successfully implemented on a tightly coupled multiple-instruction multiple-data parallel processing machine, Cray X-MP. A finite element model is divided into m domains each of which is assumed to process n elements. Each domain is then assigned to a processor or to a logical processor (task) if the number of domains exceeds the number of physical processors. The macrotasking library routines are used in mapping each domain to a user task. Computational speed-up and efficiency are used to determine the effectiveness of the algorithm. The effect of the number of domains, the number of degrees-of-freedom located along the global fronts and the dimension of the subspace on the performance of the algorithm are investigated. A parallel finite element dynamic analysis program, p-feda, is documented and the performance of its subroutines in parallel environment is analyzed.

Akl, F. A.↗

Eigensolution of finite element problems in a completely connected parallel architecture

A parallel algorithm for the solution of the generalized eigenproblem in linear elastic finite element analysis, (K)(phi)=(M)(phi)(omega), where (K) and (M) are of order N, and (omega) is of order q is presented. The parallel algorithm is based on a completely connected parallel architecture in which each processor is allowed to communicate with all other processors. The algorithm has been successfully implemented on a tightly coupled multiple-instruction-multiple-data (MIMD) parallel processing computer, Cray X-MP. A finite element model is divided into m domains each of which is assumed to process n elements. Each domain is then assigned to a processor, or to a logical processor (task) if the number of domains exceeds the number of physical processors. The macro-tasking library routines are used in mapping each domain to a user task. Computational speed-up and efficiency are used to determine the effectiveness of the algorithm. The effect of the number of domains, the number of degrees-of-freedom located along the global fronts and the dimension of the subspace on the performance of the algorithm are investigated. For a 64-element rectangular plate, speed-ups of 1.86, 3.13, 3.18 and 3.61 are achieved on two, four, six and eight processors, respectively.

Akl, Fred A.↗

The abundances of elements and isotopes in the solar wind

Solar wind abundances have now been measured for eleven elements and the isotopes of the noble gases. Aside from solar wind protons and alpha particles, which have been studied extensively since the 1960's, information for heavier elements is limited. Nevertheless, two effects stand out. First is the enrichment of abundances of elements with low first ionization potential (FIP), most likely the combined result of an atom-ion separation process in the upper chromosphere, and a marginal coupling of low-charge-state heavy ions to protons and alphas during the acceleration of the solar wind. Second, there is variability in the solar wind composition over a whole range of time scales. Recent measurements carried out in the earth's magnetosheath during times that included high-speed coronal-hole-associated flows indicate a significantly lower overabundance of low FIP elements. Given the fact that the He/H ratio is remarkably constant in the coronal hole solar wind, this result suggests that both enrichment and variability are reduced in such flows.

Gloeckler, George↗

An x ray scatter approach for non-destructive chemical analysis of low atomic numbered elements

A non-destructive x-ray scatter (XRS) approach has been developed, along with a rapid atomic scatter algorithm for the detection and analysis of low atomic-numbered elements in solids, powders, and liquids. The present method of energy dispersive x-ray fluorescence spectroscopy (EDXRF) makes the analysis of light elements (i.e., less than sodium; less than 11) extremely difficult. Detection and measurement become progressively worse as atomic numbers become smaller, due to a competing process called 'Auger Emission', which reduces fluorescent intensity, coupled with the high mass absorption coefficients exhibited by low energy x-rays, the detection and determination of low atomic-numbered elements by x-ray spectrometry is limited. However, an indirect approach based on the intensity ratio of Compton and Rayleigh scattered has been used to define light element components in alloys, plastics and other materials. This XRS technique provides qualitative and quantitative information about the overall constituents of a variety of samples.

Ross, H. Richard↗

Gravity and the geoid in the Nepal Himalaya

Materials within the Himalaya are rising due to convergence between India and Asia. If the rate of erosion is comparable to the rate of uplift the mean surface elevation will remain constant. Any slight imbalance in these two processes will lead to growth or attrition of the Himalaya. The process of uplift of materials within the Himalaya coupled with surface erosion is similar to the advance of a glacier into a region of melting. If the melting rate exceeds the rate of downhill motion of the glacier then the terminus of the glacier will receed up-valley despite the downhill motion of the bulk of the glacier. Thus although buried rocks, minerals and surface control points in the Himalaya are undoubtably rising, the growth or collapse of the Himalaya depends on the erosion rate which is invisible to geodetic measurements. Erosion rates are currently estimated from suspended sediment loads in rivers in the Himalaya. These typically underestimate the real erosion rate since bed-load is not measured during times of heavy flood, and it is difficult to integrate widely varying suspended load measurements over many years. An alternative way to measure erosion rate is to measure the rate of change of gravity in a region of uplift. If a control point moves vertically it should be accompanied by a reduction in gravity as the point moves away from the Earth's center of mass. There is a difference in the change of gravity between uplift with and without erosion corresponding to the difference between the free-air gradient and the gradient in the acceleration due to gravity caused by a corresponding thickness of rock. Essentially gravity should change precisely in accord with a change in elevation of the point in a free-air gradient if erosion equals uplift rate. We were funded by NASA to undertake a measurement of absolute gravity simultaneously with measurements of GPS height within the Himalaya. Since both absolute gravity and time are known in an absolute sense to 1 part in 10(exp 10) it is possible to estimate gravity with a precision of 0.1 mu gal. Known systematic errors reduce the measurement to an absolute uncertainty of 6 mu gal. The free air gradient at the point of measurement is typically about 3 mu gals/cm. At Simikot where our experiment was conducted we determined a vertical gravity gradient of 4.4 mu gals/cm.

Bilham, Roger↗

Consistent approach to describing aircraft HIRF protection

The high intensity radiated fields (HIRF) certification process as currently implemented is comprised of an inconsistent combination of factors that tend to emphasize worst case scenarios in assessing commercial airplane certification requirements. By examining these factors which include the process definition, the external HIRF environment, the aircraft coupling and corresponding internal fields, and methods of measuring equipment susceptibilities, activities leading to an approach to appraising airplane vulnerability to HIRF are proposed. This approach utilizes technically based criteria to evaluate the nature of the threat, including the probability of encountering the external HIRF environment. No single test or analytic method comprehensively addresses the full HIRF threat frequency spectrum. Additional tools such as statistical methods must be adopted to arrive at more realistic requirements to reflect commercial aircraft vulnerability to the HIRF threat. Test and analytic data are provided to support the conclusions of this report. This work was performed under NASA contract NAS1-19360, Task 52.

Rimbey, P. R.↗

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗

An Algorithm for Integrated Subsystem Embodiment and System Synthesis

Consider the statement,'A system has two coupled subsystems, one of which dominates the design process. Each subsystem consists of discrete and continuous variables, and is solved using sequential analysis and solution.' To address this type of statement in the design of complex systems, three steps are required, namely, the embodiment of the statement in terms of entities on a computer, the mathematical formulation of subsystem models, and the resulting solution and system synthesis. In complex system decomposition, the subsystems are not isolated, self-supporting entities. Information such as constraints, goals, and design variables may be shared between entities. But many times in engineering problems, full communication and cooperation does not exist, information is incomplete, or one subsystem may dominate the design. Additionally, these engineering problems give rise to mathematical models involving nonlinear functions of both discrete and continuous design variables. In this dissertation an algorithm is developed to handle these types of scenarios for the domain-independent integration of subsystem embodiment, coordination, and system synthesis using constructs from Decision-Based Design, Game Theory, and Multidisciplinary Design Optimization. Implementation of the concept in this dissertation involves testing of the hypotheses using example problems and a motivating case study involving the design of a subsonic passenger aircraft.

Lewis, Kemper↗

Performance of the Microwave Anisotropy Probe AST-201 Star Trackers

The Microwave Anisotropy Probe (MAP) was launched to create a full-sky map of the cosmic microwave background. MAP incorporates two modified Lockheed Martin AST-201 (Autonomous Star Tracker) star trackers. The AST-201 employs an eight element radiation hardened lens assembly which is used to focus an image on a charge coupled device (CCD). The CCD image is then processed by a star identification algorithm which outputs a three-axis attitude. A CCD-shift algorithm called Time Delayed Integration (TDI) was also included in each star tracker. In order to provide some radiation effect filtering during MAP's three to five phasing loop passes through the Van Allen radiation belts, a simple pixel filtering scheme was implemented, rather than using a more complex, but more robust windowing algorithm. The trackers also include a fiber optic data interface. This paper details the ground testing that was accomplished on the MAP trackers.

Ward, David K.↗

Particle Generation And Evolution In Silane (SiH4)/Acetylene (C2H2) Flames In Microgravity

The objective of this experimental program is to advance the understanding of the coupling of particle formation with gas phase combustion processes. The work utilizes the unique SiH4/C2H2 combustion system which generates particulate products ranging from high purity, white SiC to carbonaceous soot depending on equivalence ratio (Ref. 1). A goal of this work is to identify gas phase or particle formation processes that provide the enthalpy release needed to drive the combustion wave, and to locate the steps of the particle formation process that determine SiC stoichiometry and crystallinity. In a real sense, these SiH4/C2H2 flames act like highly sooty hydrocarbon flames, but with simpler chemistry. This simplification is expected to allow them to be used as surrogates to advance understanding of soot formation in such rich hydrocarbon flames. It is also expected that this improved understanding of SiC particle generation and evolution in these self-sustaining flames will advance the commercial potential of the flame process for the generation of high purity SiC powders.

Keil, D. G.↗

Aqueous Alteration of Carbonaceous Chondrites: New Insights from Comparative Studies of Two Unbrecciated CM2 Chondrites, Y 791198 and ALH 81002

Carbonaceous chondrites are an important resource for understanding the physical and chemical conditions in the early solar system. In particular, a long-standing question concerns the role of water in the cosmochemical evolution of carbonaceous chondrites. It is well established that extensive hydration of primary nebular phases occurred in the CM and CI chondrites, but the location where this alteration occurred remains controversial. In the CM2 chondrites, hydration formed secondary phases such as serpentine, tochilinite, pentlandite, carbonate and PCP. There are several textural observations which suggest that alteration occurred before the accretion of the final CM parent asteroid, i.e. preaccretionary alteration. Conversely, there is a significant body of evidence that supports parent-body alteration. In order to test these two competing hypotheses further, we studied two CM chondrites, Y-791198 and ALH81002, two meteorites that exhibit widely differing degrees of aqueous alteration. In addition, both meteorites have primary accretionary textures, i.e. experienced minimal asteroidal brecciation. Brecciation significantly complicates the task of unraveling alteration histories, mixing components that have been altered to different degrees from different locations on the same asteroidal parent body. Alteration in Y-791198 is mostly confined to chondrule mesostases, FeNi metal and fine-grained matrix and rims. In comparison, the primary chondrule silicates in ALH81002 have undergone extensive replacement by secondary hydrous phases. This study focuses on compositional and textural relationships between chondrule mesostasis and the associated rim materials. Our hypothesis is: both these components are highly susceptible to aqueous alteration and should be sensitive recorders of the alteration process. For parent body alteration, we expect systematic coupled mineralogical and compositional changes in rims and altered mesostasis, as elemental exchange between these components occurs. Conversely, for preaccretionary alteration, there should be no clear relationships between the rims and mesostases.

Chizmadia, L. J.↗

An Arabidopsis Ran-binding protein, AtRanBP1c, is a co-activator of Ran GTPase-activating protein and requires the C-terminus for its cytoplasmic localization

Ran-binding proteins (RanBPs) are a group of proteins that bind to Ran (Ras-related nuclear small GTP-binding protein), and thus either control the GTP/GDP-bound states of Ran or help couple the Ran GTPase cycle to a cellular process. AtRanBP1c is a Ran-binding protein from Arabidopsis thaliana (L.) Heynh. that was recently shown to be critically involved in the regulation of auxin-induced mitotic progression [S.-H. Kim et al. (2001) Plant Cell 13:2619-2630]. Here we report that AtRanBP1c inhibits the EDTA-induced release of GTP from Ran and serves as a co-activator of Ran-GTPase-activating protein (RanGAP) in vitro. Transient expression of AtRanBP1c fused to a beta-glucuronidase (GUS) reporter reveals that the protein localizes primarily to the cytosol. Neither the N- nor C-terminus of AtRanBP1c, which flank the Ran-binding domain (RanBD), is necessary for the binding of PsRan1-GTP to the protein, but both are needed for the cytosolic localization of GUS-fused AtRanBP1c. These findings, together with a previous report that AtRanBP1c is critically involved in root growth and development, imply that the promotion of GTP hydrolysis by the Ran/RanGAP/AtRanBP1c complex in the cytoplasm, and the resulting concentration gradient of Ran-GDP to Ran-GTP across the nuclear membrane could be important in the regulation of auxin-induced mitotic progression in root tips of A. thaliana.

NASA Discipline Plant Biology↗

Detection, Identification, Location, and Remote Sensing using SAW RFID Sensor Tags

In this presentation, we will consider the problem of simultaneous detection, identification, location estimation, and remote sensing for multiple objects. In particular, we will describe the design and testing of a wireless system capable of simultaneously detecting the presence of multiple objects, identifying each object, and acquiring both a low-resolution estimate of location and a high-resolution estimate of temperature for each object based on wireless interrogation of passive surface acoustic wave (SAW) radiofrequency identification (RFID) sensor tags affixed to each object. The system is being studied for application on the lunar surface as well as for terrestrial remote sensing applications such as pre-launch monitoring and testing of spacecraft on the launch pad and monitoring of test facilities. The system utilizes a digitally beam-formed planar receiving antenna array to extend range and provide direction-of-arrival information coupled with an approximate maximum-likelihood signal processing algorithm to provide near-optimal estimation of both range and temperature. The system is capable of forming a large number of beams within the field of view and resolving the information from several tags within each beam. The combination of both spatial and waveform discrimination provides the capability to track and monitor telemetry from a large number of objects appearing simultaneously within the field of view of the receiving array. In the presentation, we will summarize the system design and illustrate several aspects of the operational characteristics and signal structure. We will examine the theoretical performance characteristics of the system and compare the theoretical results with results obtained from experiments in both controlled laboratory environments and in the field.

Barton, Richard J.↗

The Evolution and Development of the Lunar Regolith and Implications for Lunar Surface Operations and Construction

The lunar regolith consists of about 90% submillimeter particles traditionally termed lunar soil. The remainder consists of larger particles ranging up to boulder size rocks. At the lower size end, soil particles in the 10s of nanometer sizes are present in all soil samples. Lunar regolith overlies bedrock which consists of either lava flows in mare regions or impact-produced megaregolith in highland regions. Lunar regolith has been produced over billions of years by a combination of breaking and communition of bedrock by meteorite bombardment coupled with a variety of complex space weathering processes including solar wind implantation, solar flare and cosmic ray bombardment with attendant radiation damage, melting, vaporization, and vapor condensation driven by impact, and gardening and turnover of the resultant soil. Lunar regolith is poorly sorted compared to most terrestrial soils, and has interesting engineering properties including strong grain adhesion, over-compacted soil density, an abundance of agglutinates with sharp corners, and a variety of properties related to soil maturity. The NASA program has supported a variety of engineering test research projects, the production of bricks by solar or microwave sintering, the production of concrete, the in situ sintering and glazing of regolith by microwave, and the extraction of useful resources such as oxygen, hydrogen, iron, aluminum, silicon and other products. Future requirements for a lunar surface base or outpost will include construction of protective berms, construction of paved roadways, construction of shelters, movement and emplacement of regolith for radiation shielding and thermal control, and extraction of useful products. One early need is for light weight but powerful digging, trenching, and regolith-moving equipment.

McKay, David↗

Detection, Identification, Location, and Remote Sensing Using SAW RFID Sensor Tags

The Electromagnetic Systems Branch (EV4) of the Avionic Systems Division at NASA Johnson Space Center in Houston, TX is studying the utility of surface acoustic wave (SAW) radiofrequency identification (RFID) tags for multiple wireless applications including detection, identification, tracking, and remote sensing of objects on the lunar surface, monitoring of environmental test facilities, structural shape and health monitoring, and nondestructive test and evaluation of assets. For all of these applications, it is anticipated that the system utilized to interrogate the SAW RFID tags may need to operate at fairly long range and in the presence of considerable multipath and multiple-access interference. Towards that end, EV4 is developing a prototype SAW RFID wireless interrogation system for use in such environments called the Passive Adaptive RFID Sensor Equipment (PARSED) system. The system utilizes a digitally beam-formed planar receiving antenna array to extend range and provide direction-of-arrival information coupled with an approximate maximum-likelihood signal processing algorithm to provide near-optimal estimation of both range and temperature. The system is capable of forming a large number of beams within the field of view and resolving the information from several tags within each beam. The combination of both spatial and waveform discrimination provides the capability to track and monitor telemetry from a large number of objects appearing simultaneously within the field of view of the receiving array. In this paper, we will consider the application of the PARSEQ system to the problem of simultaneous detection, identification, localization, and temperature estimation for multiple objects. We will summarize the overall design of the PARSEQ system and present a detailed description of the design and performance of the signal detection and estimation algorithms incorporated in the system. The system is currently configured only to measure temperature (jointly with range and tag ID), but future versions will be revised to measure parameters other than temperature as SAW tags capable of interfacing with external sensors become available. It is anticipated that the estimation of arbitrary parameters measured using SAW-based sensors will be based on techniques very similar to the joint range and temperature estimation techniques described in this paper.

Barton, Richard J.↗

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗

Reducing Centroid Error Through Model-Based Noise Reduction

A method of processing the digitized output of a charge-coupled device (CCD) image detector has been devised to enable reduction of the error in computed centroid of the image of a point source of light. The method involves model-based estimation of, and correction for, the contributions of bias and noise to the image data. The method could be used to advantage in any of a variety of applications in which there are requirements for measuring precise locations of, and/or precisely aiming optical instruments toward, point light sources. In the present method, prior to normal operations of the CCD, one measures the point-spread function (PSF) of the telescope or other optical system used to project images on the CCD. The PSF is used to construct a database of spot models representing the nominal CCD pixel outputs for a point light source projected onto the CCD at various positions incremented by small fractions of a pixel.

Lee, Shinhak↗

NASA Tech Briefs, March 2005

Topics covered include: Scheme for Entering Binary Data Into a Quantum Computer; Encryption for Remote Control via Internet or Intranet; Coupled Receiver/Decoders for Low-Rate Turbo Codes; Processing GPS Occultation Data To Characterize Atmosphere; Displacing Unpredictable Nulls in Antenna Radiation Patterns; Integrated Pointing and Signal Detector for Optical Receiver; Adaptive Thresholding and Parameter Estimation for PPM; Data-Driven Software Framework for Web-Based ISS Telescience; Software for Secondary-School Learning About Robotics; Fuzzy Logic Engine; Telephone-Directory Program; Simulating a Direction-Finder Search for an ELT; Formulating Precursors for Coating Metals and Ceramics; Making Macroscopic Assemblies of Aligned Carbon Nanotubes; Ball Bearings Equipped for In Situ Lubrication on Demand; Synthetic Bursae for Robots; Robot Forearm and Dexterous Hand; Making a Metal-Lined Composite-Overwrapped Pressure Vessel; Ex Vivo Growth of Bioengineered Ligaments and Other Tissues; Stroboscopic Goggles for Reduction of Motion Sickness; Articulating Support for Horizontal Resistive Exercise; Modified Penning-Malmberg Trap for Storing Antiprotons; Tumbleweed Rovers; Two-Photon Fluorescence Microscope for Microgravity Research; Biased Randomized Algorithm for Fast Model-Based Diagnosis; Fast Algorithms for Model-Based Diagnosis; Simulations of Evaporating Multicomponent Fuel Drops; Formation Flying of Tethered and Nontethered Spacecraft; and Two Methods for Efficient Solution of the Hitting- Set Problem.

Source record↗