Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “configurational disorder”

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.

15 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Two-dimensional O2 adsorbed on graphite

Properties of two-dimensional O2 adsorbed on graphite are calculated in the extremely low-coverage delta region and for monolayers, with use of pattern-recognition optimization and Monte Carlo techniques. Equilibrium configurations and orientations, orientational order-disorder, melting, and dissociation transitions are predicted at various O2 densities. Phase characteristics, including a plastic crystallite phase, are compared with experiment.

Etters, R. D.↗
NASA NTRS2005

Effect of electron-nuclear spin interactions for electron-spin qubits localized in InGaAs self-assembled quantum dots

The effect of electron-nuclear spin interactions on qubit operations is investigated for a qubit represented by the spin of an electron localized in an InGaAs self-assembled quantum dot. The localized electron wave function is evaluated within the atomistic tight-binding model. The electron Zeeman splitting induced by the electron-nuclear spin interaction is estimated in the presence of an inhomogeneous environment characterized by a random nuclear spin configuration, by the dot-size distribution, alloy disorder, and interface disorder. Due to these inhomogeneities, the electron Zeeman splitting varies from one qubit to another by the order of 10(-6), 10(-6), 10(-7), and 10(-9) eV, respectively. Such fluctuations cause errors in exchange operations due to the inequality of the Zeeman splitting between two qubits. However, the error can be made lower than the quantum error threshold if an exchange energy larger than 10(-4) eV is used for the operation. This result shows that the electron-nuclear spin interaction does not hinder quantum-dot based quantum computer architectures from being scalable even in the presence of inhomogeneous environments.

Whaley, K. Birgitta↗
NASA NTRS2003

Study of alloy disorder in quantum dots through multi-million atom simulations

A tight binding model which includes s, p, d, s orbitals is used to examine the electronic structures of an ensemble of dome-shaped In0.6 Ga0.4 As quantum dots. Given ensembles of identically sized quantum dots, variations in composition and configuration yield a linewidth broadening of less than 0.35 meV, much smaller than the total broadening determined from photoluminescence experiments. It is also found that the computed disorder-induced broadening is very sensitive to the applied boundary conditions, so that care must be taken to ensure proper convergence of the numerical results. Examination of local eigenenergies as functions of position shows similar convergence problems and indicates that an inaccurate resolution of the equilibrium atomic positions due to truncation of the simulation domain may be the source of the slow ground state convergence.

disorder↗
NASA NTRS1988

Radiation resistance and comparative performance of ITO/InP and n/p InP homojunction solar cells

The radiation resistance of ITO/InP cells processed by DC magnetron sputtering is compared to that of standard n/p InP and GaAs homojunction cells. After 20 MeV proton irradiations, it is found that the radiation resistance of the present ITO/InP cell is comparable to that of the n/p homojunction InP cell and that both InP cell types have radiation resistance significantly greater than GaAs. The relatively lower radiation resistance, observed at higher fluence, for the InP cell with the deepest junction depth, is attributed to losses in the cells emitter region. Diode parameters obtained from I sub sc - V sub oc plots, data from surface Raman spectroscopy, and determinations of surface conductivity types are used to investigate the configuration of the ITO/InP cells. It is concluded that thesee latter cells are n/p homojunctions, the n-region consisting of a disordered layer at the oxide semiconductor.

Weinberg, I.↗
NASA NTRS1988

Radiation resistance and comparative performance of ITO/InP and n/p InP homojunction solar cells

The radiation resistance of ITO/InP cells processed by dc magnetron sputtering is compared to that of standard n/p InP and GaAs homojunction cells. After 20 MeV proton irradiations, it is found that the radiation resistance of the present ITO/InP cell is comparable to that of the n/p homojunction InP cell and that both InP cell types have radiation resistances significantly greater than GaAs. The relatively lower radiation resistance, observed at higher fluence, for the InP cell with the deepest junction depth, is attributed to losses in the cells emitter region. Diode parameters obtained from I sub sc - V sub oc plots, data from surface Raman spectrosocpy, and determinations of surface conductivity type are used to investigate the configuration of the ITO/InP cells. It is concluded that these latter cells are n/p homojunctions, the n-region consisting of a disordered layer at the oxide semiconductor.

Weinberg, I.↗
NASA NTRS1994

Formation Energy of Disordered Alloys form the Energetics of Ordered Compounds

The cluster expansion technique is used in combination with first-principles calculations of the total energy of ordered compounds to study the energetics of the disordered phase for a number of alloy systems. The effect of short-range order, as seen in the energy differences between an alloy with the configuration corresponding to finite temperature and the perfectly random one is studied. The role of relaxation due to large size mismatch is discussed in terms of an effective cluster volume approximation. Very good agreement with measured energies of formation is obtained when short-range order and relaxation effects are taken into account.

Amador, Carlos↗
NASA NTRS1999

Centrosome and microtubule instability in aging Drosophila cells

Several cytoskeletal changes are associated with aging which includes alterations in muscle structure leading to muscular atrophy, and weakening of the microtubule network which affects cellular secretion and maintenance of cell shape. Weakening of the microtubule network during meiosis in aging oocytes can result in aneuploidy or trisomic zygotes with increasing maternal age. Imbalances of cytoskeletal organization can lead to disease such as Alzheimer's, muscular disorders, and cancer. Because many cytoskeletal diseases are related to age we investigated the effects of aging on microtubule organization in cell cultures of the Drosophila cell model system (Schneider S-1 and Kc23 cell lines). This cell model is increasingly being used as an alternative system to mammalian cell cultures. Drosophila cells are amenable to genetic manipulations and can be used to identify and manipulate genes which are involved in the aging processes. Immunofluorescence, scanning, and transmission electron microscopy were employed for the analysis of microtubule organizing centers (centrosomes) and microtubules at various times after subculturing cells in fresh medium. Our results reveal that centrosomes and the microtubule network becomes significantly affected in aging cells after 5 days of subculture. At 5-14 days of subculture, 1% abnormal out of 3% mitoses were noted which were clearly distinguishable from freshly subcultured control cells in which 3% of cells undergo normal mitosis with bipolar configurations. Microtubules are also affected in the midbody during cell division. The midbody in aging cells becomes up to 10 times longer when compared with midbodies in freshly subcultured cells. During interphase, microtubules are often disrupted and disorganized, which may indicate improper function related to transport of cell organelles along microtubules. These results are likely to help explain some cytoskeletal disorders and diseases related to aging.

Non-NASA Center↗
NASA NTRS1980

A model for the formation and stabilization of charged water clathrates

A model for the formation and stabilization of charged water clathrates is presented which accounts for observed anomalies in H(+)(H2O)n ion distributions. These anomalies are observed in both ion cluster and neutral expansions and are consistent with the sizes expected for clathrate ions. That the same sizes are observed in both ion cluster and neutral expansions strongly suggests that a rapid ionic process is responsible for their formation. The proposed model is based on the high mobility and bonding effects of the excess proton in water. Computer simulations suggest that excess proton movement in a water clathrate would be suitable for stabilizing the clathrate structure as well as giving it access to a large number of nearly degenerate proton configurations. The formation of clathrates in charged water clusters of proper size can be ascribed to the following: rapid excess proton movement, a strong preference of the H3O(+) for a three-coordinate bonding structure (which is compatible with hydrogen bonding), and finally, relatively slow processes leading to thermal disorder.

Holland, P. M.↗
NASA NTRSDate not supplied

Mixed-Domain Charge Transport in S-Se Alloys as a Li-S Battery Cathode Material

Lithium-sulfur batteries are emerging candidate systems for the next-generation long-range electric aircraft application owing to their potential to deliver high energy density per weight. Their cathodes are to be made primarily with sulfur, but sulfur by itself has much too low electron conductivity to serve as a useful cathode. An idea that has been suggested to overcome this bottleneck is to alloy sulfur with selenium in the hope that electron mobility, which is thought to improve charge transport at some expense of increased weight and reduced energy density. However, understanding of these alloy structures and their transport mechanism is insufficient, and electron mobility at varying degrees of selenium content is not well charted, which are critical for determining the optimum selenium content and testing whether this strategy is generally feasible. The difficulty of computationally characterizing these alloy systems is rooted in the structure. The structures of both sulfur and selenium exhibit eight-atom rings that pack together in various orientations to form their respective crystals. For one, because these crystals (and presumably their alloys as well) feature multiple polymorphs such that their structural coherence is rather unclear. Secondly, these are also relatively large systems (32 atoms per cell) with low-symmetry, which complicate computation both in terms of accuracy and efficiency. Thirdly, such semi-molecular, semi-crystalline structures lead to a combination of band-transport characteristics and hopping-transport characteristics, each of which is a domain with its own physics and set of computational, theoretical challenges. In addition, there is a general dearth of experimental transport data for sulfur-selenium alloys. In this study, we make a comprehensive attempt to tackle this problem using a wide array of first-principles methods. We generate special quasirandom structures to simulate alloy structures, and compute their first-principles Raman spectra for comparison with experimental data to ensure their structural soundness. We then proceed to use recently developed, state-of-the-art tool (AMSET) in order to efficiently compute electronic mobilities under band transport of pure sulfur, pure selenium, as well as their alloys of various compositions at a reasonable accuracy. We then sample numerous dimer configurations of nearest-neighbor eight-atom-ring-pairs and compute electronic hopping rates between them, which yields hopping mobility. A combination of these efforts lead to a general mapping of electron transport behaviors throughout the alloy range. Preliminary results show that introduction of selenium into sulfur initially damages electron mobility due to disorder but eventually improves it beyond that of pure sulfur with additional selenium content, ultimately peaking at the pure-selenium limit. Band transport dominates for holes in sulfur and electrons in selenium, but in all other cases, hopping transport is dominant. Ongoing efforts include determination of charge concentration and conductivity, as well as performing multiphysics modeling to determine the optimum selenium content for best tradeoff between conductivity and energy density for aircraft range.

Junsoo Park↗
NASA NTRS2009

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗
NASA NTRS2003

Numerical Modeling of Nanoelectronic Devices

Nanoelectronic Modeling 3-D (NEMO 3-D) is a computer program for numerical modeling of the electronic structure properties of a semiconductor device that is embodied in a crystal containing as many as 16 million atoms in an arbitrary configuration and that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantummechanical behavior of the device resolved to the atomistic level of granularity. The system of electrons in the device is represented by a sparse Hamiltonian matrix that contains hundreds of millions of terms. NEMO 3-D solves the matrix equation on a Beowulf-class cluster computer, by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. In a recent update of NEMO 3-D, a new strain treatment, parameterized for bulk material properties of GaAs and InAs, was developed for two tight-binding submodels. The utility of the NEMO 3-D was demonstrated in an atomistic analysis of the effects of disorder in alloys and, in particular, in bulk In(x)Ga(l-x)As and in In0.6Ga0.4As quantum dots.

Klimeck, Gerhard↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
NASA NTRS2019

NASA’s Galactic Cosmic Ray Simulator at Brookhaven National Laboratory: Enabling Human Exploration Missions to the Moon and Mars

With exciting new Agency plans for a sustainable return to the moon, astronauts will once again leave earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic rays (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPE’s with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever-penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include epithelial carcinogenesis and leukemia, central nervous system effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including cataracts, circulatory and heart disease, as well as, potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy ion exposures and the validation of countermeasures in a relevant space environment. NASA has developed the “GCR Simulator” at Brookhaven National Laboratory to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (~65-75%) and alpha particles (~10-20%) with heavier ions (Z≤3) contributing the remainder. The “GCR Simulator” exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the five heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires 75-90 minutes. To more closely simulate the low dose rates found in space, sequential field exposures can be divided into daily fractions over 2-4 weeks, with individual fractions as low as 0.1-0.2 mGy. In the large beam configuration (60 x 60 cm(exp 2)), 54 special housing cages can accommodate 2-3 mice each for a 70-75 min duration or ~15 individually housed rats. Emerging research results from our 2018 runs utilizing mixed heavy ion fields and protracted space exposures are forthcoming and deepen our understanding of the numerous health risks faced by our astronauts. This paper discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NASA Space Radiation Laboratory to enable future exploration missions.

Lisa C Simonsen↗
NASA NTRSDate not supplied

NASA Space Radiation Laboratory Galactic Cosmic Ray Simulator

With exciting new Agency plans for a sustainable return to the moon, astronauts will once again leave earth’s protective magnetosphere only to endure higher levels of radiation from galactic cosmic rays (GCR) and the possibility of a large solar particle event (SPE). Gateway, lunar landers, and surface habitats will be designed to protect crew against SPEs with vehicle optimization, storm shelter concepts, and/or active dosimetry; however, the ever-penetrating GCR will continue to pose the most significant health risks especially as lunar missions increase in duration and as NASA sets its aspirations on Mars. The primary risks of concern include carcinogenesis, central nervous system effects resulting in potential in-mission cognitive or behavioral impairment and/or late neurological disorders, degenerative tissue effects including cataracts, circulatory and heart disease, as well as, potential immune system decrements impacting multiple aspects of crew health. Characterization and mitigation of these risks requires a significant reduction in the large biological uncertainties of chronic (low-dose rate) heavy ion exposures and the validation of countermeasures in a relevant space environment. NASA has developed the “GCR Simulator” at Brookhaven National Laboratory to generate a spectrum of ion beams that approximates the primary and secondary GCR field experienced at human organ locations within a deep-space vehicle. The majority of the dose is delivered from protons (~65-75%) and helium ions (~10-20%) with heavier ions (Z>3) contributing the remainder. The “GCR Simulator” exposes state-of-the art cellular and animal model systems to 33 sequential beams including 4 proton energies plus degrader, 4 helium energies plus degrader, and the five heavy ions of C, O, Si, Ti, and Fe. A polyethylene degrader is used with the 100 MeV/n H and He beams to provide a nearly continuous distribution of low energy particles. A 500 mGy exposure, delivering doses from each of the 33 beams, requires 75-90 minutes. To more closely simulate the low dose rates found in space, sequential field exposures can be divided into daily fractions over 2-6 weeks, with individual beam fractions as low as 0.1-0.2 mGy. In the large beam configuration (60 x 60 cm2), 54 special housing cages can accommodate 2-3 mice each for a 70-75 min duration or ~15 individually housed rats. Emerging research results from our 2018 runs utilizing mixed heavy ion fields and protracted space exposures are forthcoming and will deepen our understanding of the numerous health risks faced by astronauts. This talk discusses NASA’s innovative technology solution for a ground-based GCR simulator at the NASA Space Radiation Laboratory to enable future exploration missions.

galactic cosmic ray simulator↗
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↗