Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solar optical performance”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Photovoltaic Cable Jackets: A Comparison of Representative Products Using Combined-Accelerated Stress Testing [Slides]

Photovoltaic (PV) cables facilitate the distribution of electricity collected from modules to the energy grid. Durable cabling enables continuous operation of PV installations, whereas cables with a lifetime less than the modules must be replaced - reducing electricity generation and adding to the operating expense. This study primarily focusses on the aging of the key cable types using the combined-accelerated stress testing (C-AST) protocol. Representative cables for utility, building, and control/auxiliary applications were examined. Cable jacket materials examined include: polyolefin, polyethylene, polyamide, poly(vinyl chloride), chlorinated polyethylene, thermoplastic elastomer, and ethylene propylene diene monomer rubber. Specimen characterizations applied include: optical microscopy, mechanical profilometry, instrumented indentation, scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDS), Fourier-transform infrared spectroscopy (FTIR), thermogravimetric analysis (TGA), and differential scanning calorimetry (DSC). A variety of performance and durability characteristics were observed, depending on the base material, polymer formulation, and jacket color. The results from C-AST are analyzed and discussed relative to a recent industry survey on electronic balance of system components in addition to a recent study where similar cables were aged using steady state ultraviolet weathering (International Electrotechnical Commission Technical Specification 62788-7-2). Recommendations are made for the screening, industry qualification, and service life prediction of PV cable jackets.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2007

GEO Population Estimates using Optical Survey Data

Optical survey data taken using the NASA Michigan Orbital Debris Survey Telescope (MODEST) gives us an opportunity to statistically sample faint object population in the Geosynchronous (GEO) and near-GEO environment. This paper will summarize the MODEST survey work that has been conducted by NASA since 2002, and will outline the techniques employed to arrive at the current population estimates in the GEO environment for dim objects difficult to detect and track using current systems in the Space Surveillance Network (SSN). Some types of orbits have a higher detection rate based on what parts of the GEO belt is being observed, a straightforward statistical technique is used to debias these observations to arrive at an estimate of the total population potentially visible to the telescope. The size and magnitude distributions of these fainter debris objects are markedly different from the catalogued population. GEO debris consists of at least two different populations, one which follows the standard breakup power law and one which has anomalously high Area-to-Mass Ratios (1 to approx. 30 square meters per kilogram; a sheet of paper = approx. 13 square meters per kilogram). The Inter-Agency Space Debris Coordination Committee (IADC) is investigating objects in GEO orbits with anomalously high Area-to-Mass Ratios (AMRs). The ESA Space Debris Telescope discovered this population and has and its properties of inclinations (0 to 30 degrees), changing eccentricities (0 and 0.6), and mean motions (approx. 1 rev), will be presented. The accepted interpretation of this orbital behavior is that solar radiation pressure drives the perturbations causing time varying inclinations and eccentricities. The orbital parameters are unstable for this population and thus difficult to predict. Their dim visual magnitudes and photometric variability make observations a challenge. The IADC has enlisted a series of observatories (participating institutions: University of Michigan/CTIO, Astronomical Institute University of Bern, Boeing LTS / AMOS, Keldysh Institute of Applied Mathematics) at different longitudes. Complete observational coverage over periods of days to months will provide a better understanding of the properties, such as solar radiation pressure effects on orbital elements, size, shape, attitude, color variations, and spectral characteristics. Results from recent observational programs will be summarized, and includes a description of the orbit elements prediction processes, a summary of the metric tracking performance, and some photometric characteristics of this class of debris.

Barker, Edwin S.↗
NASA NTRS1996

Astrophysical Adaptation of Points, the Precision Optical Interferometer in Space

POINTS (Precision Optical INTerferometer in Space) would perform microarcsecond optical astrometric measurements from space, yielding submicroarcsecond astrometric results from the mission. It comprises a pair of independent Michelson stellar interferometers and a laser metrology system that measures both the critical starlight paths and the angle between the baselines. The instrument has two baselines of 2 m, each with two subapertures of 35 cm; by articulating the angle between the baselines, it observes targets separated by 87 to 93 deg. POINTS does global astrometry, i.e., it measures widely separated targets, which yields closure calibration, numerous bright reference stars, and absolute parallax. Simplicity, stability, and the mitigation of systematic error are the central design themes. The instrument has only three moving-part mechanisms, and only one of these must move with sub-milliradian precision; the other two can tolerate a precision of several tenths of a degree. Optical surfaces preceding the beamsplitter or its fold flat are interferometrically critical; on each side of the interferometer, there are only three such. Thus, light loss and wavefront distortion are minimized. POINTS represents a minimalistic design developed ab initio for space. Since it is intended for astrometry, and therefore does not require the u-v-plane coverage of an imaging, instrument, each interferometer need have only two subapertures. The design relies on articulation of the angle between the interferometers and body pointing to select targets; the observations are restricted to the 'instrument plane.' That plane, which is fixed in the pointed instrument, is defined by the sensitive direction for the two interferometers. Thus, there is no need for siderostats and moving delay lines, which would have added many precision mechanisms with rolling and sliding parts that would be required to function throughout the mission. Further, there is no need for a third interferometer, as is required when out-of-plane observations are made. An instrument for astrometry, unlike those for imaging, can be compact and yet scientifically productive. The POINTS instrument is compact and therefore requires no deployment of precision structures, has no low-frequency (i.e., under 100 Hz) vibration modes, and is relatively easy to control thermally. Because of its small size and mass, it is easily and quickly repointed between observations. Further, because of the low mass, it can be economically launched into high Earth orbit which, in conjunction with a solar shield, yields nearly unrestricted sky coverage and a stable thermal environment.

Reasenberg, Robert D.↗
DOE OSTI2021

Hot-carrier dynamics in InAs/AlAsSb multiple-quantum wells

A type-II InAs/AlAs 0.16 Sb 0.84 multiple-quantum well sample is investigated for the photoexcited carrier dynamics as a function of excitation photon energy and lattice temperature. Time-resolved measurements are performed using a near-infrared pump pulse, with photon energies near to and above the band gap, probed with a terahertz probe pulse. The transient terahertz absorption is characterized by a multi-rise, multi-decay function that captures long-lived decay times and a metastable state for an excess-photon energy of >100 meV. For sufficient excess-photon energy, excitation of the metastable state is followed by a transition to the long-lived states. Excitation dependence of the long-lived states map onto a nearly-direct band gap (E g ) density of states with an Urbach tail below E g . As temperature increases, the long-lived decay times increase < E g , due to the increased phonon interaction of the unintentional defect states, and by phonon stabilization of the hot carriers > E g . Additionally, Auger (and/or trap-assisted Auger) scattering above the onset of the plateau may also contribute to longer hot-carrier lifetimes. Meanwhile, the initial decay component shows strong dependence on excitation energy and temperature, reflecting the complicated initial transfer of energy between valence-band and defect states, indicating methods to further prolong hot carriers for technological applications.

14 SOLAR ENERGY↗
DOE OSTI2020

Photoinduced electron transfer from zinc meso -tetraphenylporphyrin to a one-dimensional perylenediimide aggregate: Probing anion delocalization effects

Organic photovoltaics incorporating non-fullerene acceptors based on perylenediimide (PDI) now rival fullerene acceptor-based devices in performance, although the mechanisms of charge generation in PDI-based devices are not yet fully understood. Fullerene-based systems are proposed to undergo electron transfer directly from the photoexcited donor into a band of delocalized acceptor states, thus increasing charge generation efficiency. Similarly, anion delocalization has been shown to enhance the rate of electron transfer from a photoexcited donor to two electronically coupled PDI acceptors. Here we investigate how additional electron acceptors may further increase the rate of electron transfer from the donor zinc meso-tetraphenylporphyrin (ZnTPP) to an aggregate of PDI acceptors (PDI[Formula: see text]. Femtosecond transient visible and mid-infrared absorption spectroscopies show that the rate of electron transfer from 1 *ZnTPP to the PDI assembly ZnTPP 2 -PDI 3 is statistically identical to that of the previously examined ZnTPP-PDI 2 . A Marcus theory analysis indicates that the parameters governing electron transfer are nearly identical for the two molecules, suggesting that the maximum electron transfer rate enhancement has been achieved in a cofacial PDI dimer because the ZnTPP directly couples to the first two PDI acceptors whereas the coupling to the third PDI is too weak.

bio-inspired↗
NASA NTRS1994

On the decay of outbursts in dwarf novae nad X-ray novae

We perform computations using a time-dependent model for the accretion disk limit-cycle mechanism to examine the decay of the optical light following the peak of a dwarf nova outburst. We present the results of a parameter study of the physical input variables which affect the decay rate. In the model, the decay is brought about by a cooling transition front which begins at large radii in the disk and moves inward. The nature of the decay is strongly influenced by the radial dependence of the accretion disk viscosity parameter alpha. To obtain exponential decays for typical dwarf nova parameters, we require alpha proportional to r(exp epsilon(sub 0)), where epsilon(sub 0) approximately = 0.3-0.4. The exact value of epsilon(sub 0) which produces exponential decays depends on factors such as the mass of the accreting star and the inner radius of the accretion disk. Therefore, the observed ubiquity of exponential decays in two different types of systems (dwarf novae and X-ray novae) leads us to believe that alpha is an unnatural scaling for the viscosity. The physics of the cooling transition front must be self-regulating in that the timescale (-parital derivative of lnSigma(r)/partial derivative +)(exp -1) (where Sigma is the surface density) for mass extraction across the front remains constant. This may be consistent with a scaling alpha proportional to (h/r)(exp n), where h is the local disk semi-thickness and n approximately 1-2. As regards the speed of the cooling front, we find v(sub F)(r) proportional to r(exp p), where p approximately 3 at large radii, with an abrupt transition to p approximately 0 at some smaller radius. The r(exp 3) dependence is much steeper than has been found by previous workers and appears to result from the strong variation of specific heat within the cooling front when the front resides at a large radius in the disk. The outflow of disk material across the cooling front causes a significant departure of dln T(sub dff0/dln r from the standard value of -0.75 (expected from steady state accretion) within about 0.2 dex in radius of the break associated with the cooling front -- T(sub eff) aproximately 10(exp 3.9) K (r/10(exp 10 cm)) (exp -0.1). These effects should be observable with eclipse mapping. Finally, it appears that the relatively slow decay rate for the optical flux in the 1975 outburst of A0620-00 can be accounted for if the primary is a approximately 10 Solar mass black hole.

Cannizzo, John K.↗
NASA NTRS2016

3D Printing in Zero G Technology Demonstration Mission: Summary of On-Orbit Operations, Material Testing, and Future Work

Human space exploration to date has been limited to low Earth orbit and the moon. The International Space Station (ISS), an orbiting laboratory 200 miles above the earth, provides a unique and incredible opportunity for researchers to prove out the technologies that will enable humans to safely live and work in space for longer periods of time and venture farther into the solar system. The ability to manufacture parts in-space rather than launch them from earth represents a fundamental shift in the current risk and logistics paradigm for human spaceflight. In particularly, additive manufacturing (or 3D printing) techniques can potentially be deployed in the space environment to enhance crew safety (by providing an on-demand part replacement capability) and decrease launch mass by reducing the number of spare components that must be launched for missions where cargo resupply is not a near-term option. In September 2014, NASA launched the 3D Printing in Zero G technology demonstration mission to the ISS to explore the potential of additive manufacturing for in-space applications and demonstrate the capability to manufacture parts and tools on-orbit. The printer for this mission was designed and operated by the company Made In Space under a NASA SBIR (Small Business Innovation Research) phase III contract. The overarching objectives of the 3D print mission were to use ISS as a testbed to further maturation of enhancing technologies needed for long duration human exploration missions, introduce new materials and methods to fabricate structure in space, enable cost-effective manufacturing for structures and mechanisms made in low-unit production, and enable physical components to be manufactured in space on long duration missions if necessary. The 3D print unit for fused deposition modeling (FDM) of acrylonitrile butadiene styrene (ABS) was integrated into the ISS Microgravity Science Glovebox (MSG) in November 2014 and phase I printing operations took place from November through December of that year. Phase I flight operations yielded 14 unique parts (21 total specimens) that could be directly compared against ground-based prints of identical geometry manufactured using the printer prior to its launch to ISS. The 3DP unit functioned safely and produced specimens necessary to advance the understanding of the critical design and operational parameters for the FDM process as affected by the microgravity environment. From the standpoint of operations, 3DP demonstrated the ability to remove parts from the build-tray on-orbit, teleoperate the printer from the ground, perform critical maintenance functions within defined human factors limits, produce a functional tool that could be evaluated for form/fit/function, and uplink a new part file from the ground and produce it on the printer. The flight parts arrived at NASA Marshall Space Flight Center in Huntsville, Alabama in April 2015, where they underwent months of testing in the materials and processes laboratory. Ground and flight prints completed the following phases of testing: photographic/visual inspection, mass and density evaluation, structured light scanning, XRay and CT, mechanical testing, optical microscopy, scanning electron microscopy, and chemical analysis. This presentation will discuss the results of this testing as well as phase II operations for the printer, which took place in June and July of 2016. Lessons learned from the tech demo and their impacts on the design and development of the second generation 3D printer for ISS, the Additive Manufacturing Facility (AMF) by Made In Space will also be presented. In addition, progress in other elements of NASA's In Space Manufacturing (ISM) initiative such as the on-demand ISM utilization catalog, in-space Recycler ISS Technology Demonstration development, launch packaging recycling, in-space printable electronics, development of higher strength polymeric materials for 3D printing and Additive Construction by Mobile Emplacement (ACME) will also be addressed.

Prater, Tracie↗
NASA NTRS2014

Remote Sensing of Cloud Top Height from SEVIRI: Analysis of Eleven Current Retrieval Algorithms

The role of clouds remains the largest uncertainty in climate projections. They influence solar and thermal radiative transfer and the earth's water cycle. Therefore, there is an urgent need for accurate cloud observations to validate climate models and to monitor climate change. Passive satellite imagers measuring radiation at visible to thermal infrared (IR) wavelengths provide a wealth of information on cloud properties. Among others, the cloud top height (CTH) - a crucial parameter to estimate the thermal cloud radiative forcing - can be retrieved. In this paper we investigate the skill of ten current retrieval algorithms to estimate the CTH using observations from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard Meteosat Second Generation (MSG). In the first part we compare ten SEVIRI cloud top pressure (CTP) data sets with each other. The SEVIRI algorithms catch the latitudinal variation of the CTP in a similar way. The agreement is better in the extratropics than in the tropics. In the tropics multi-layer clouds and thin cirrus layers complicate the CTP retrieval, whereas a good agreement among the algorithms is found for trade wind cumulus, marine stratocumulus and the optically thick cores of the deep convective system. In the second part of the paper the SEVIRI retrievals are compared to CTH observations from the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) and Cloud Profiling Radar (CPR) instruments. It is important to note that the different measurement techniques cause differences in the retrieved CTH data. SEVIRI measures a radiatively effective CTH, while the CTH of the active instruments is derived from the return time of the emitted radar or lidar signal. Therefore, some systematic differences are expected. On average the CTHs detected by the SEVIRI algorithms are 1.0 to 2.5 kilometers lower than CALIOP observations, and the correlation coefficients between the SEVIRI and the CALIOP data sets range between 0.77 and 0.90. The average CTHs derived by the SEVIRI algorithms are closer to the CPR measurements than to CALIOP measurements. The biases between SEVIRI and CPR retrievals range from −0.8 kilometers to 0.6 kilometers. The correlation coefficients of CPR and SEVIRI observations vary between 0.82 and 0.89. To discuss the origin of the CTH deviation, we investigate three cloud categories: optically thin and thick single layer as well as multi-layer clouds. For optically thick clouds the correlation coefficients between the SEVIRI and the reference data sets are usually above 0.95. For optically thin single layer clouds the correlation coefficients are still above 0.92. For this cloud category the SEVIRI algorithms yield CTHs that are lower than CALIOP and similar to CPR observations. Most challenging are the multi-layer clouds, where the correlation coefficients are for most algorithms between 0.6 and 0.8. Finally, we evaluate the performance of the SEVIRI retrievals for boundary layer clouds. While the CTH retrieval for this cloud type is relatively accurate, there are still considerable differences between the algorithms. These are related to the uncertainties and limited vertical resolution of the assumed temperature profiles in combination with the presence of temperature inversions, which lead to ambiguities in the CTH retrieval. Alternative approaches for the CTH retrieval of low clouds are discussed.

Hamann, U.↗
NASA NTRS1999

The MECA Payload as a Dust Analysis Laboratory on the MSP 2001 Lander

In a companion abstract, the "Mars Environmental Compatibility Assessment" (MECA) payload for Mars Surveyor Program 2001 (MSP 2001) is described in terms of its capabilities for addressing exobiology on Mars. Here we describe how the same payload elements perform in terms of gathering data about surface dust on the planet. An understanding of the origin and properties of dust is important to both human exploration and planetary geology. The MECA instrument is specifically designed for soil/dust investigations: it is a multifunctional laboratory equipped to assess particulate properties with wet chemistry, camera imagery, optical microscopy (potentially with LTV fluorescence capability), atomic force microscopy (AFM; potentially with mineral-discrimination capabilities), electrometry, active & passive external materials-test panels, mineral hardness testing, and electrostatic & magnetic materials testing. Additionally, evaluation of soil chemical and physical properties as a function of depth down to about 50 cm will be facilitated by the Lander/MECA robot arm on which the camera (RAC) and electrometer are mounted. Types of data being sought for the dust include: (1) general textural and grain-size characterization of the soil as a whole --for example, is the soil essentially dust with other components or is it a clast-supported material in which dust resides only in the clast interstices, (2) size frequency distribution for dust particles in the range 0.01 to 10.00 microns, (3) particle-shape distribution of the soil components and of the fine dust fraction in particular, (4) soil fabric such as grain clustering into clods, aggregates, and cemented/indurated grain amalgamations, as well as related porosity, cohesiveness, and other mechanical soil properties, (5) cohesive relationship that dust has to certain types of rocks and minerals as a clue to which soil materials may be prime hosts for dust "piggybacking", (6) particle, aggregate, and bulk soil electrostatic properties, (7) particle hardness, (8) particle magnetic properties, (9) bulk dust geochemistry (solubility, reactivity, ionic and mineral species). All of these quantities are needed in order for the human exploration program to make assessments of hazards on Mars, and to better enable the production on earth, of soil/dust simulants that can act as realistic test materials in terms of those properties that render dust a contaminant.Such properties include the small grain size that enables penetration of space-suit joints, mechanical interfaces and bearings, seals, etc., and presents difficulty for filtration systems. Size also plays a critical role in the potential for lung disease in long-term habitats. The properties of grain shape and hardness are important parameters in determining the abrasiveness of dust as it enters mechanical systems, or bombards helmet visors and habitat windows in dust-laden winds. Adhesive electrostatic and magnetic properties of dust will be prime causes of contamination of space suits and equipment. Contamination causes mechanical malfunction, tracking of dirt into habitats, "piggybacking" of toxins on dust into habitats, changes in albedo and efficiency of solar arrays and heat exchangers, and changes in electrical conductivity of suit surfaces and other materials that may have specific safety requirements regarding electrical conductivity. Other potentially hazardous properties of dust include the possibility of high solubility of some component grains (rendering them reactive), and toxicity of some materials --grains of superoxidants and heavy metals (there is always the slim, but not inconceivable possibility of biogenic components such as spores). Because Mars has no active surface aqueous regime, volcanic emissions, meteoritic debris, weathering products, and photochemical products of Mars have nowhere to go except reside in the surface; there are few mechanical or chemical (buffering) processes to remove the accumulation of eons. From a planetology perspective, there are many enigmatic issues relating to dust and the aeolian regime in general. MECA will be able to address many questions in this area. For example, if MECA determines a particular particle size distribution (size and sorting values), it will be possible to make inferences about the origin of the dust - - is it all aeolian, or a more primitive residue of weathering, volcanic emissions, and meteoritic gardening? Trenching with the Lander/MECA robot arm will enable local stratigraphy to be determined in terms of depositional rates, amounts and cyclicity in dust storms and/or local aeolian transport. Grain shape will betray the origin of the dust fragments as being the product of recent or ancient weathering, or the comminution products of aeolian transport --the dust-silt ratio might be a measure of aeolian comminution energy. Additional information is contained in the original.

Marshall, J.↗
NASA NTRS2021

Validation Assessment for the Soil Moisture Active Passive (SMAP) Level 4 Carbon (L4_C) Data Product Version 5

The post-launch Cal/Val phase of the SMAP mission is guided by two primary objectives for each science product team: 1) to calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the SMAP Level-1 mission science requirements. Algorithm science and product maintenance activities during the SMAP extended mission phase have also involved periodic algorithm calibration and product refinements to maintain or enhance product consistency and performance as well as science utility. This report provides an assessment of the latest (Version 5) SMAP Level 4 Carbon (L4_C) product. The L4_C Version 5 (v5) global record now spans more than six years (March 2015 – present) of SMAP operations and has benefited from five major reprocessing updates to the operational product. These reprocessing events and L4_C product release updates have incorporated various algorithm refinements and calibration adjustments to account for similar refinements to the upstream GEOS land model assimilation system, SMAP brightness temperatures, and MODIS vegetation inputs used for L4_C processing. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by daily surface and root zone soil moisture information contributed from the SMAP Level 4 Soil Moisture (L4_SM) product along with optical remote sensing-based (e.g. MODIS-based) land cover and canopy fractional photosynthetic active radiation (fPAR), and other ancillary biophysical data. The carbon flux model estimates global daily net ecosystem CO2 exchange (NEE) and the component carbon fluxes, namely, vegetation gross primary production (GPP) and soil heterotrophic respiration (Rh). Other L4_C product elements include surface (~0-5 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture-related controls on GPP and ecosystem respiration (Kimball et al. 2014, Jones et al. 2017). The L4_C product addresses SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying freeze/thaw and soil moisture-related constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity. The SMAP L4_C algorithms and operational product are mature and at a CEOS Validation Stage 4 level (Jackson et al. 2012) based on extensive validation of the multi-year record against a diverse array of independent benchmarks, well characterized global performance, and systematic refinements gained from five major reprocessing events. There are no Level-1 mission science requirements for the L4_C product; however, self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications (e.g., Liu et al. 2019, Endsley et al. 2020); however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased root-mean-square error (ubRMSE, or standard deviation of the error) for NEE of 1.6 g C m-2 d-1 and 30 g C m-2 yr-1, emphasizing northern (≥45°N) boreal and arctic ecosystems; this accuracy is similar to that of tower eddy covariance measurement-based observations (Baldocchi 2008). Methods used for the latest v5 L4_C product performance and validation assessment have been established from the SMAP Cal/Val plan and previous studies (Jackson et al. 2012, Jones et al. 2017) and include: 1) consistency evaluations of the product fields against earlier product releases (version 4 or earlier); 2) comparisons of daily carbon flux estimates with independent tower eddy 4covariance measurement-based daily carbon (CO2) flux observations from core tower validation sites (CVS); and 3) consistency checks against other global carbon products, including soil carbon inventory records, global GPP records derived from tower observation upscaling methods, and satellite-based observations of canopy solar induced chlorophyll fluorescence (SIF) as a surrogate for GPP. Metrics used to evaluate relative agreement between L4_C product fields and observational benchmarks include correlation (r-value), RMSE differences, bias and model sensitivity diagnostics. Following these validation criteria, the present report provides a validation assessment of the latest L4_C product release (v5). Detailed descriptions of the L4_C algorithm and additional global product accuracy and performance results are given elsewhere (Jones et al. 2017, Endsley et al. 2020). The v5 L4_C product replaces earlier product versions and continues to show: (i) accuracy and performance levels meeting or exceeding SMAP L4_C science requirements; (ii) improvement over the previous product version (version 4); and (iii) suitability for a diversity of science applications. Example L4_C applications from the recent literature include clarifying environmental trends and controls on the northern terrestrial carbon sink (Liu et al. 2019), diagnosing drought-related impacts on ecosystem productivity (Li et al. 2020), and regional monitoring of cropland conditions for projecting annual yields (Wurster et al. 2020). 2EXPECTED L4_C ALGORITHM AND PRODUCTPERFORMANCE The L4_C algorithm performance, including variance and uncertainty estimates of model outputs, was determined during the mission pre-launch phase through spatially explicit model sensitivity studies using available model inputs similar to those currently being used for operational production and evaluating the resulting model simulations over the observed range of northern (≥45 °N) and global conditions (Kimball et al. 2012, Entekhabi et al. 2014). The L4_C algorithm options were also evaluated during the mission prelaunch phase, including deriving canopy fPAR from lower order NDVI (Normalized Difference Vegetation Index) inputs in lieu of using MODIS (MOD15) fPAR; and including an explicit model representation of boreal fire disturbance recovery impacts. These results indicated that the L4_C accuracy requirements (i.e., NEE ubRMSE ≤ 30 g C m-2 yr-1or ≤ 1.6 g C m-2 d-1) could be met from the baseline algorithms over more than 82% and 89% of global and northern vegetated land areas, respectively (Yi et al. 2013, Kimball et al. 2014). The global L4_C algorithm error budget for NEE derived during the mission prelaunch phase indicated that the estimated NEE ubRMSE uncertainty is proportional to GPP and is therefore larger in higher biomass productivity areas, including forests and croplands (Kimball et al. 2014). Likewise, NEE ubRMSE uncertainty is expected to be lower in less-productive areas, including grasslands and shrublands. Expected model NEE ubRMSE levels were also generally within targeted accuracy levels for characteristically less-productive boreal and Arctic biomes, even though relative model error as a proportion of total productivity (NEE RMSE / GPP) may be large in these areas. The estimated NEE uncertainty was lower than expected in some warmer tropical high biomass productivity areas (e.g. Amazon rainforest) because of reduced low temperature and moisture constraints to the L4_C respiration calculations so that the bulk of model uncertainty is contributed by GPP in these areas. Model NEE uncertainty in the African Congo was estimated to

SMAP↗
DOE OSTI2024

Parallax-corrected VISST-derived pixel-level products from satellite GOES-16

The NASA Langley group led by William Smith produced GOES-16 satellite cloud retrievals over an approximate 10 by 10 degree region over the CACTI field campaign location. These retrievals are described here: https://www.arm.gov/capabilities/vaps/visst and are available for download here . They use algorithms historically called VISST that are now referred to as SatCORPS. More information can be found in Trepte et al. (2019), Minnis et al. (2021), and Yost et al. (2021). If using this dataset, please cite these references, the CACTI VISST dataset DOI found at the download link above, and this dataset’s DOI. The CACTI VISST pixel-level retrievals are on a 2 km spatial grid and available every 15 minutes (every 10 minutes late in the campaign), producing 21,765 files for the entire field campaign between October 2018 and April 2019. They are not corrected for parallax error, which is an offset in the actual geographical location of a cloud above the surface due to the satellite viewing the cloud partly from the side off nadir. This dataset applies a correction for parallax using the location relative to the satellite and the retrieved cloud top height above the surface, which allows the dataset to be geo-located with surface-based observations. The parallax correction for each location depends on the longitude, latitude and cloud top height above ground level (AGL) for that longitude and latitude in the original VISST files. The cloud top height AGL requires first computing the surface elevation at each VISST grid point. Data from the Advanced Spaceborne Thermal Emission and Reflection (ASTER) Global Digital Elevation Map Version 3 at 30-m resolution is projected onto the VISST grid using conservative coarsening (conserving surface elevation) in the xESMF Python package. The surface elevation is then subtracted from the VISST-retrieved cloud top height above mean sea level. These cloud top heights AGL are then combined with longitude and latitude to estimate the latitude and longitude corrections. Due to variability in cloud top height, the parallax shifts produce an irregular grid of values since higher cloud tops are shifted further than lower cloud tops. A ball tree-based neighbor search with Haversine distance is performed using the Python-based scikit-learn library to find the nearest VISST grid point to each parallax correction-shifted point. The data value of the shifted point is then assigned to that VISST grid point. In this manner, the irregular geographical shifts to correct for parallax are projected back to the rectilinear VISST grid. Because relatively higher clouds should obscure lower clouds, the variable values for the highest cloud top are preferentially chosen if two or more values are assigned to a grid point. The parallax correction should be viewed as an improved but still imperfect estimation of the cloud top locations, largely because the cloud top height is an imperfect retrieval. Please see the attached README document for further information. Users are encouraged to contact the authors with any additional questions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

Measuring Shock Stage of ltokawa Regolith Grains by Electron Back-Scattered Diffraction and Synchrotron X-Ray Diffraction

We have been analyzing Itokawa samples in order to definitively establish the degree of shock experienced by the regolith of asteroid Itokawa, and to devise a bridge between shock determinations by standard light optical petrography, crystal structures as determined by electron and X-ray diffraction techniques. We are making measurements of olivine crystal structures and using these to elucidate critical regolith impact processes. We use electron back-scattered diffraction (EBSD) and synchrotron X-ray diffraction (SXRD). We are comparing the Itokawa samples to L and LL chondrite meteorites chosen to span the shock scale experienced by Itokawa, specifically Chainpur (LL3.4, Shock Stage 1), Semarkona (LL3.00, S2), Kilabo (LL6, S3), NWA100 (L6, S4) and Chelyabinsk (LL5, S4). In SXRD we measure the line broadening of olivine reflections as a measure of shock stage. In this presentation we concentrate on the EBSD work. We employed JSC's Supra 55 variable pressure FEG-SEM and Bruker EBSD system. We are not seeking actual strain values, but rather indirect strain-related measurements such as extent of intra-grain lattice rotation, and determining whether shock state "standards" (meteorite samples of accepted shock state, and appropriate small grain size) show strain measurements that may be statistically differentiated, using a sampling of particles (number and size range) typical of asteroid regoliths. Using our system we determined that a column pressure of 9 Pa and no C-coating on the sample was optimal. We varied camera exposure time and gain to optimize mapping performance, concluding that 320x240 pattern pixilation, frame averaging of 3, 15 kV, and low extractor voltage yielded an acceptable balance of hit rate (>90%), speed (11 fps) and map quality using an exposure time of 30 ms (gain 650). We found that there was no strong effect of step size on Grain Orientation Spread (GOS) and Grain Reference Orientation Deviation angle (GROD-a) distribution; there was some effect on grain average Kernel Average Misorientation (KAM) (reduced with smaller step size for the same grain), as expected. We monitored GOS, Maximum Orientation Spread (MOS) and GROD-a differences between whole olivine grains and sub-sampled areas, and found that there were significant differences between the whole grain dataset and subsets, as well as between subsets, likely due to sampling-related "noise". Also, in general (and logically) whole grains exhibit greater degrees of cumulative lattice rotation. Sampling size affects the apparent strain character of the grain, at least as measured by GOS, MOS and GROD-a. There were differences in the distribution frequencies of GOS and MOS between shock stages, and in plots of MOS and GOS vs. grain diameter. These results are generally consistent with those reported this year. However, it is unknown whether the differences between samples of different shock states exceeds the clustering of these values to the extent that shock stage determinations can still be made with confidence. We are investigating this by examination of meteorites with higher shock stage 4 to 5. Our research will improve our understanding of how small, primitive solar system bodies formed and evolved, and improve understanding of the processes that determine the history and future of habitability of environments on other solar system bodies. The results will directly enrich the ongoing asteroid and comet exploration missions by NASA and JAXA, and broaden our understanding of the origin and evolution of small bodies in the early solar system, and elucidate the nature of asteroid and comet regolith.

Zolensky, Michael↗
DOE OSTI2025

Diagnosis of PV Cell Antireflective Coating Degradation Resulting From Hot-Humid High-Voltage Potential Aging

Corrosion of the antireflective coating on a photovoltaic cell ("ARc corrosion") has previously been observed in studies using hot-humid test conditions with external high-voltage (HV) bias. This study primarily focuses on known vulnerable legacy aluminum back surface field cells in mini-modules (MiMos) put through comparative stepped stress tests. Each cell type had MiMos at +1500 V, -1500 V, or unbiased ("Voc") potential, which were sequentially subjected to test conditions of 60 degrees C/60% relative humidity (RH) for 96 h, as in International Electrotechnical Commission Technical Specification 62804-1; 70 degrees C/70% RH for 200 h; and 85 degrees C/85% RH for 200 h. Characterizations at each step included visual camera and electroluminescence (EL) imaging, colorimetry, and current-voltage curve tracing. Final characterizations included: Suns-Voc, spatial mapping of external quantum efficiency, high-resolution photoluminescence, EL, and dark lock-in thermography imaging. Forensics were performed on extracted cores, including scanning electron microscopy (SEM) with energy-dispersive X-ray spectroscopy (EDS), X-ray photoelectron spectroscopy, and scanning Auger microscopy (SAM). Forensics were also conducted on MiMos from previous studies that underwent stepped HV aging and separate outdoor aged full-sized modules. ARc corrosion was specifically seen for the glass/encapsulant/cell side of the +1500 V (HV+) stressed MiMos and modules. Appearance, color, and reflectance were the most distinguishing characteristics relative to glass corrosion, gridline corrosion and delamination, and other concurrent degradation modes. SEM/EDS and SAM identified the conversion of silicon nitride to hydrated silica, hydrous silica, or hydrated amorphous silica, which preferentially occurred at the edges and tips of the pyramidal textured cell surface.

14 SOLAR ENERGY↗
DOE OSTI2025

Diagnosis of PV Cell Antireflective Coating Degradation Resulting from Hot-Humid, High Voltage Potential Aging

Corrosion of the antireflective coating on a photovoltaic cell ("ARc corrosion") has previously been observed in studies using hot-humid test conditions with external high-voltage (HV) bias. This study primarily focuses on known vulnerable, legacy aluminum back surface field (Al-BSF) cells in mini-modules (MiMos) put through comparative stepped stress tests. Each cell type had MiMos at +1500 V, -1500 V, or unbiased ("Voc") potential, which were sequentially subjected to test conditions of 60deg C/60% relative humidity (RH) for 96 h, as in International Electrotechnical Commission (IEC) Technical Specification 62804-1; 70 deg C/70% RH for 200 h; and 85 deg C/85% RH for 200 h. Characterizations at each step included visual camera and electroluminescence (EL) imaging, colorimetry, and current-voltage (I-V) curve tracing. Final characterizations included: Suns-Voc; spatial mapping of external quantum efficiency (EQE); high-resolution photoluminescence (PL), EL, and dark lock-in thermography (DLIT) imaging. Forensics were performed on extracted cores, including scanning electron microscopy (SEM) with energy-dispersive X-ray spectroscopy (EDS), X-ray photoelectron spectroscopy (XPS), and scanning Auger microscopy (SAM). Forensics were also conducted on MiMos from previous studies that underwent stepped HV aging and separate outdoor aged full-sized modules. ARc corrosion was specifically seen for the glass/encapsulant/cell side of the +1500V (HV+) stressed MiMos and modules. Appearance, color, and reflectance were the most distinguishing characteristics relative to glass corrosion, gridline-corrosion and -delamination as well as other concurrent degradation modes. SEM/EDS and SAM identified conversion of silicon nitride to hydrated silica, hydrous silica, or hydrated amorphous silica, which preferentially occurred at the edges and tips of the pyramidal textured cell surface.

14 SOLAR ENERGY↗
NASA NTRS2007

Broadband Achromatic Telecentric Lens

A new type of lens design features broadband achromatic performance as well as telecentricity, using a minimum number of spherical elements. With appropriate modifications, the lens design form can be tailored to cover the range of response of the focal-plane array, from Si (400-1,000 nm) to InGaAs (400-1,700 or 2,100 nm) or InSb/HgCdTe reaching to 2,500 nm. For reference, lenses typically are achromatized over the visible wavelength range of 480-650 nm. In remote sensing applications, there is a need for broadband achromatic telescopes, normally satisfied with mirror-based systems. However, mirror systems are not always feasible due to size or geometry restrictions. They also require expensive aspheric surfaces. Non-obscured mirror systems can be difficult to align and have a limited (essentially one-dimensional) field of view. Centrally obscured types have a two-dimensional but very limited field in addition to the obscuration. Telecentricity is a highly desirable property for matching typical spectrometer types, as well as for reducing the variation of the angle of incidence and cross-talk on the detector for simple camera types. This rotationally symmetric telescope with no obscuration and using spherical surfaces and selected glass types fills a need in the range of short focal lengths. It can be used as a compact front unit for a matched spectrometer, as an ultra-broadband camera objective lens, or as the optics of an integrated camera/spectrometer in which the wavelength information is obtained by the use of strip or linear variable filters on the focal plane array. This kind of camera and spectrometer system can find applications in remote sensing, as well as in-situ applications for geological mapping and characterization of minerals, ecological studies, and target detection and identification through spectral signatures. Commercially, the lens can be used in quality-control applications via spectral analysis. The lens design is based on the rear landscape lens with the aperture stop in front of all elements. This allows sufficient room for telecentricity in addition to making the stop easily accessible. The crucial design features are the use of a doublet with an ultra-low dispersion glass (fluorite or S-FPL53), and the use of a strong negative element, which enables flat field and telecentricity in conjunction with the last (field lens) element. The field lens also can be designed to be in contact with the array, a feature that is desirable in some applications. The lens has a 20deg field of view, for a 50-mm focal length, and is corrected over the range of wavelengths of 450-2,300 nm. Transverse color, which is the most pernicious aberration for spectroscopic work, is controlled at the level of 1 m or below at 0.7 m field and 5 m at full field. The maximum chief ray angle is less than 1.7 , providing good telecentricity. An additional feature of this lens is that it is made exclusively with glasses that provide good transmission up to 2,300 nm and even some transmission to 2,500 nm; thus, the lens can be used in applications that cover the entire solar-reflected spectrum. Alternative realizations are possible that provide enhanced resolution and even less transverse color over a narrower wavelength range.

Mouroulis, Pantazis↗
DOE OSTI2025

Selection of high-redshift Lyman-Break Galaxies from broadband and wide photometric surveys

In this paper, we investigate the possibility of selecting high-redshift Lyman-Break Galaxies (LBG) using current and future broadband wide photometric surveys, such as the Ultraviolet Near Infrared Optical Northern Survey (UNIONS) or the Vera C. Rubin Legacy Survey of Space and Time (LSST), using a Random Forest algorithm. This work is conducted in the context of future large-scale structure spectroscopic surveys like DESI-II, the next phase of the Dark Energy Spectroscopic Instrument (DESI), which will start around 2029.We use deep imaging data from the Hyper Suprime Camera (HSC) and the Canada-France-Hawaii Telescope Large Area U-band Deep Survey (CLAUDS) on the COSMOS and XMM-LSS fields. To predict the selection performance of LBGs with image quality similar to UNIONS, we degrade the u,g,r,i and z bands to UNIONS depth.The Random Forest algorithm is trained with the u,g,r,i and z bands to classify LBGs in the 2.5 < z < 3.5 range.We find that fixing a target density budget of 1,100 deg$^{-2}$, the Random Forest approach gives a density of z > 2 targets of 873 deg$^{-2}$, and a density of 493 deg$^{-2}$ of confirmed LBGs after spectroscopic confirmation with DESI. This UNIONS-like selection was tested in a dedicated spectroscopic observation campaign of 1,000 targets with DESI on the COSMOS field, providing a safe spectroscopic sample with a mean redshift of 3. This sample is used to derive forecasts for DESI-II, assuming a sky coverage of 5,000 deg$^{2}$. We predict uncertainties on Alcock-Paczynski parameters α$_{⊥}$ and α$_{∥}$ to be 0.7% and 1% for 2.6 < z < 3.2, resulting in a potential 2% measurement of the dark energy fraction at high redshift. Additionally, we estimate the uncertainty in local non-Gaussianity and predict σ$_{fNL}$ ≈ 7, which would be comparable to the current best precision achieved by Planck. The latter forecast suggests that achieving the precision required to place stringent constraints on inflationary models (σ$_{fNL}$ ≈ 1) using spectroscopic galaxy surveys necessitates the development of a next-generation (Stage V) spectroscopic survey.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

The 2.5-12 micrometers spectrum of comet Halley from the IKS-VEGA experiment

The infrared instrument IKS flown on board the VEGA space probes was designed for the detection of emission bands of parent molecules, and for a measurement of the size and temperature of the thermal emitting nuclear region. The instrument had three channels with cooled detectors: an "imaging channel" designed to modulate the signal of the nucleus and two spectroscopic channels operating at 2.5-5 and 6-12 micrometers, respectively, equipped with circular variable filters of resolving power approximately 50. This paper presents and discusses the results from the spectral channels. On VEGA 1, usable spectra were obtained at distances D from the comet nucleus ranging from 250,000 to 40,000 km corresponding to fields of view 4000 and 700 km in diameter, respectively. The important internal background signal caused by the instrument itself, which could not be cooled, had to be eliminated. Since no sky chopping was performed, we obtain difference spectra between the current spectrum and a reference spectrum with little or no cometary signal taken at the beginning of the observing sequence (D approximately 200,000 km). Final discrimination between cometary signal and instrumental background is achieved using their different time evolution, since the instrumental background is proportional to the slow temperature drift of the instrument, and the cometary signal due to parent molecules or dust grains is expected to vary in first order as D-1. The 2.5-5 micrometers IKS spectra definitely show strong narrow signals at 2.7 and 4.25 micrometers, attributed to the nu 3 vibrational bands of H2O and CO2, respectively, and a broader signal in the region 3.2-3.5 micrometers, which may be attributed to CH-bearing molecules. All these signals present the expected D-1 intensity variation. Weaker emission features at 3.6 and 4.7 micrometers could correspond to the nu 1 and nu 5 bands of H2CO and the (1 - 0) band of CO, respectively. Molecular production rates are derived from the observed emissions, assuming that they are due to resonance fluorescence excited by the Sun's infrared radiation. For the strong bands of H2O and CO2, the rovibrational lines are optically thick, and radiative transfer is taken into account. We derive production rates, at the moment of the VEGA 1 flyby, of approximately 10(30) sec-1 for H2O, approximately 2.7 x 10(28) sec-1 for CO2, approximately 5 x 10(28) sec-1 for CO, and 4 x 10(28) sec-1 for H2CO, if attributions to CO and H2CO are correct. The production rate of carbon atoms in CH-bearing molecules is approximately 9 x 10(29) sec-1 assuming fluorescence of molecules in the gas phase, but could be much less if the 3.2-3.5 micrometers emission is attributed to C-H stretch in polycyclic aromatic hydrocarbons or small organic grains. In addition, marginal features are present at 4.85 and 4.45 micrometers, tentatively attributed to OCS and molecules with the CN group, respectively. Broad absorption at 2.8-3.0 micrometers, as well as a narrow emission at 3.15 micrometers, which follow well the D-1 intensity variation, might be due to water ice. Emission at 2.8 micrometers is also possibly present, and might be due to OH created in vibrationally excited states after water photodissociation. The 6-12 micrometers spectrum does not show any molecular emission, nor emission in the 7.5-micrometers region. The spectrum is dominated by silicate emission showing a double structure with maxima at 9.0 and 11.2 micrometers, which suggests the presence of olivine.

NASA Discipline Number 52-10↗
NASA NTRS2013

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗