Search NASA⌕ Search

SEARCH · Search NASA

Results for “large deformations”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

The influence of a thoria dispersion on preferred orientation in nickel alloys.

Preferred orientation has been examined in extruded, drawn, and swaged nickel, nickel-20% chromium, and nickel-20% chromium-10% tungsten, with and without a dispersion of 2 vol % of thoria particles. Duplex 100 (in carets) + 111 (in carets) fiber textures were observed in most cases. The 100 (in carets) component was stabilized by the presence of particles and also by the application of intermediate anneals during deformation. These results and the stability of the deformed state in thoriated nickel are attributed largely to the action of the particles in generating homogeneous deformation.

Hutchinson, W. B.↗
NASA NTRS2006

The Use of RGPS Kinematic Data to Estimate Nonlinear Sea Ice Motion

In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motion In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motio

Pruis, M.↗
NASA NTRS1984

Identification and control of structures in space

The derivation of the equations of motion for the Spacecraft Control Laboratory Experiment (SCOLE) is reported and the equations of motion of a similar structure orbiting the earth are also derived. The structure is assumed to undergo large rigid-body maneuvers and small elastic deformations. A perturbation approach is proposed whereby the quantities defining the rigid-body maneuver are assumed to be relatively large, with the elastic deformations and deviations from the rigid-body maneuver being relatively small. The perturbation equations have the form of linear equations with time-dependent coefficients. An active control technique can then be formulated to permit maneuvering of the spacecraft and simultaneously suppressing the elastic vibration.

Meirovitch, L.↗
NASA NTRS1985

Maneuvering of flexible spacecraft with application to SCOLE

This paper is concerned with the derivation of the equations of motion for the Spacecraft Control Laboratory Experiment (SCOLE). For future reference, the equations of motion of a similar structure orbiting the earth are also derived. The structure is assumed to undergo large rigid-body maneuvers and small elastic deformations. A perturbation approach is presented where the quantities defining the rigid-body maneuver are assumed to be relatively large, with the elastic deformations and deviations from the rigid-body maneuver being relatively small. The perturbation equations have the form of linear, non-self-adjoint equations with time-dependent coefficients. An active control technique can then be formulated to permit maneuvering of the spacecraft and simultaneously suppressing the elastic vibration.

Meirovitch, L.↗
NASA NTRS1994

Plains tectonism on Venus: Inferences from canali longitudinal profiles

Canali-type channels on Venus show inverted profile segments; that is, large portions of the channels trend uphill. The original gradients for these channels were not horizontal, but must have progressively trended downhill. Therefore, undulation patterns imposed on the profiles have to be the result of tectonism occurring since channel formation. This implies that some of the uppermost geological units of the plains, which are coincident with canali, experienced significant postemplacement tectonic deformation. The pattern of deformation is hierarchical. Within observational limits, at least two scales of deformation are indicated. The longer scale deformation (thousands of kilometers) corresponds to large-scale basins; the shorter scale deformation (hundreds of kilometers) corresponds to ridge belts or to small-scale domes/basins. The latter features, at scales up to a few hundred kilometers, extensively deform some plain areas. Although not appearing in the topographic profiles because of resolution limitations, deformation scales of tens of kilometers or less, which mostly reflect wrinkle ridges, also overlap the longer deformation scales. The channels probably formed relatively quickly in comparison to the time scale of deformation. Canali formation is closely related to the genesis of plains, and canali profile deformation reflects tectonic processes operating at multiple scales. The lower limit rates of large-scale tectonic warping are comparable to epeirogenetic deformation rates for Earth's intraplate continental interiors.

Komatsu, Goro↗
NASA NTRS2006

Videogrammetric Model Deformation Measurement Technique for Wind Tunnel Applications

Videogrammetric measurement technique developments at NASA Langley were driven largely by the need to quantify model deformation at the National Transonic Facility (NTF). This paper summarizes recent wind tunnel applications and issues at the NTF and other NASA Langley facilities including the Transonic Dynamics Tunnel, 31-Inch Mach 10 Tunnel, 8-Ft high Temperature Tunnel, and the 20-Ft Vertical Spin Tunnel. In addition, several adaptations of wind tunnel techniques to non-wind tunnel applications are summarized. These applications include wing deformation measurements on vehicles in flight, determining aerodynamic loads based on optical elastic deformation measurements, measurements on ultra-lightweight and inflatable space structures, and the use of an object-to-image plane scaling technique to support NASA s Space Exploration program.

Barrows, Danny A.↗
NASA NTRS1962

The Effect on the Buckling of Perfect Cylinders of Prebuckling Deformations and Stresses Induced by Edge Support

Large deflection theory is used to compute buckling loads of simply supported initially perfect cylinders under axial compression, external hydrostatic pressure, and combinations of axial compression and internal or external pressure. Important results are obtained by taking into account prebuckling deformations and stresses induced by edge support. For example, the presence of these deformations and stresses can decrease the axial-compression buckling load of an unpressurized perfect cylinder by 50 percent or more.

CYLINDER↗
NASA NTRS1978

Internal friction Q factor measurements in lunar rocks

In order to better interpret recently reported values for the variation of seismic Q as a function of depth below the lunar surface, we have developed apparatus and made laboratory measurements of Q as a function of hydrostatic pressure, temperature and frequency. Our measurements of the Q associated with shear deformations have demonstrated that the large difference in Q between well outgassed and volatile rich rocks persists to pressures corresponding to a depth of at least 50 km. Here we report new measurements of Q as a function of temperature, on the development of techniques to measure the Q associated with extensional deformations under hydrostatic pressure, on the derivation of theoretical relations between our laboratory Q values and the attenuation coefficient of seismic waves, and on the development of a model for mechanism of adsorption.

Tittmann, B. R.↗
NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗
NASA NTRS1993

Complex structure of the Thaumasia region of Mars

The Thaumasia region was the first center of Tharsis tectonism, and it is the most complex and poorly understood. Therefore, a geologic map of the entire Thaumasia region (lat 15 deg to 50 deg S, long 50 deg to 115 deg), at 1:5,000,000 scale is being compiled. This region is mostly made up of the Thaumasia plateau, the highlands are fractured by Thaumasia, southern Claritas, Coracis, Melas, and Nectaris Fossae. Preliminary structural analysis of the most complexly faulted area in the region (the central part, at lat 30 deg to 45 deg S, long 80 deg to 100 deg) indicates that, unlike other regions of Mars, Thaumasia has undergone extensive deformation by both small- and large-scale extensional and compressional structures. These results indicate that the early (Noachian) style of tectonism commonly involved lithospheric-scale deformation, in contrast to most younger tectonism (which mainly affected the upper parts of the crust above mechanical discontinuities); this difference may be due to a weaker (and thus more readily deformable) early lithosphere in this region.

Tanaka, Kenneth L.↗
NASA NTRSDate not supplied

The Stop-and-Go Mechanism: Towards an Integrated Approach to Model Seismicity, Outgassing, Deformation, and Thermal Unrest at Active Volcanoes

Connecting the geophysical and geochemical signals recorded at and above the surface of volcanoes with source mechanisms is fundamental to understand transitions from quiescence to eruption, and to integrate the behavior of volcanoes with their regional seismotectonic and hydrological context. Some of these signals include shallow volcanic tremor, a long-lasting (minutes-to-years) ground vibration detected in volcanic areas during unrest; volcanic outgassing, which shows intriguing periodicities over multiple scales, as revealed by space-borne instruments and by ground-based, high-frequency (~1 Hz), sampling methods; inflation and deflation of volcanic edifices, commonly associated with subsurface volume changes; and large-scale (from a few to tens of km2), low-temperature, thermal anomalies, a newly-discovered signal that has been observed to emerge on volcanic flanks from months-to-years prior to gas- and magma-driven eruptions. Traditionally, these geophysical and geochemical signals have been modelled independently, which limits our interpretation of subsurface processes and thus our assessment of unrest. In this work, we show through lumped-parameter models, numerical simulations, and preliminary laboratory experiments that many of the signals recorded around volcanoes can be integrated into a common “stop-and-go” mechanism. In particular, we found that they can arise spontaneously when taking into account the feedbacks between the temporary accumulation of gas and/or magma in the crust (“stop”), and the diffusive transfer of gas and heat towards the surface (“go”). For example, we find that the “stop-and-go” mechanism can explain: (i) the emergence of monochromatic, broadband, and harmonic tremor; (ii) the periodic components commonly observed in volcanic outgassing time series; (iii) the link between large-scale thermal anomalies and shallow hydrothermal systems; (iv) the reported time lags between large-scale thermal anomalies and deformation; and (v) the breathing of potentially destructive volcanic calderas. The “stop-and-go” mechanism provides a consistent and realistic framework to link multiple geophysical and geochemical signals with the processes leading to volcanic unrest and eruption.

volcanoes↗
NASA NTRS1984

Crustal Dynamics Project - Status and plans

The progress of the NASA Crustal Dynamics Project is summarized. The Project is performed to measure regional deformation and strain accumulation due to large earthquakes in the plate boundary region of western N. America, to measure relative plate motions around the globe, to measure internal deformations of continental and oceanic lithospheric plates, to relate the rotational dynamics of the earth to earthquakes, plate motions and other geophysical phenomena and to measure regional fault motions and strain accumulation globally. VLBI and satellite laser ranging (SLR) have been employed using several base and mobile stations in N. America, Europe, S. America and the S. Pacific. Satellite data are being gathered with the LAGEOS spacecraft. Measurement campaigns in S. America and the Carribbean and Mediterranean regions will be carried out with new instrumentation in 1985.

Coates, R. J.↗
NASA NTRS1999

The Considere Condition and Rapid Stretching of Linear and Branched Polymer Melts

We analyze the onset of "necking" and subsequent filament failure during the transient uniaxial elongation of viscoelastic fluid samples in extensional rheometers. In the limit of rapid elongation (such that no molecular relaxation occurs), the external work applied is all stored elastically and the Considere criterion originally developed in solid mechanics can be used to quantitatively predict the critical Hencky strain to failure. By comparing the predictions of the Doi-Edwards model for linear homopolymer melts with those of the "Pom-Pom" model for prototypical branched melts we show that the critical strain to failure in rapid elongation of a rubbery material is intimately linked to the molecular topology of the chain, especially the degree of chain branching. The onset of necking instability is monotonically shifted to larger Hencky strains as the number of branches is increased. Numerical computations at finite Deborah numbers also show that there is an optimal range of deformation rates over which homogeneous extensions can be maintained to large strain. We also consider other rapid homogeneous stretching deformations, such as biaxial and planar stretching, and show that the degree of stabilization afforded by inclusion of material with long-chain branching is a sensitive function of the imposed mode of deformation.

McKinley, Gareth H.↗
NASA NTRS2017

Near Earth Asteroid Scout Solar Sail Thrust and Torque Model

The Near Earth Asteroid (NEA) Scout is a solar sail mission whose objective is to scout at least one Near Earth Asteroid to help prepare for human missions to Near Earth Asteroids. NEA Scout will launch as a secondary payload on the first SLS-Orion mission. NEA Scout will perform a small trim maneuver shortly after deploy from the spent SLS upper stage using a cold gas propulsion system, but from that point on will depend entirely on the solar sail for thrust. As such, it is important to accurately characterize the thrust of the sail in order to achieve mission success. Additionally, the solar sail creates a relatively large solar disturbance torque that must be mitigated. For early mission design studies a flat plate model of the solar sail with a fixed center of pressure was adequate, but as mission concepts and the sail design matured, greater fidelity was required. Here we discuss the progress to a three-dimensional sail model that includes the effects of tension and thermal deformation that has been derived from a large structural Finite Element Model (FEM) developed by the Langley Research Center. We have found that the deformed sail membrane affects torque relatively much more than thrust; a flat plate model could potentially model thrust well enough to close mission design studies, but a three-dimensional solar sail is essential to control system design. The three-dimensional solar sail model revealed that thermal deformations of unshielded booms would create unacceptably large solar disturbance torques. The original large FEM model was used in control and mission simulations, but was resulted in simulations with prohibitive run times. This led us to adapt the Generalized Sail Model (GSM) of Rios-Reyes. A design reference sail model has been baselined for NEA Scout and has been used to design the mission and control system for the sailcraft. Additionally, since NEA Scout uses reaction wheels for attitude pointing and control, the solar torque model is essentially to successfully design the NEA Scout momentum management control system. We have also updated the estimate of diffusivity used for the aluminized sail material based on optical testing of wrinkled sail material. The model presented here represents the current state of the art of NASA's ability to model solar sail thrust and torque.

Heaton, Andy↗
NASA NTRS2019

Wrinkle Ridges on Mercury and the Moon Within and Outside of Mascons

Found on all terrestrial planets, wrinkle ridges are anticlines formed by thrust faulting and folding resulting from crustal shortening. The MErcury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER) spacecraft's orbital phase returned high resolution images and topographic data of the previously unimaged northern high latitudes of Mercury where there are large expanses of smooth plains deformed by wrinkle ridges. Concurrently, the Lunar Reconnaissance Orbiter (LRO) is obtaining high resolution images and topographic data covering lunar mare wrinkle ridges. These data allow quantitative comparison of the scale of wrinkle ridges in smooth plains volcanic units on Mercury with mare wrinkle ridges. We evaluate the topographic relief of 300 wrinkle ridges within and outside of mascon basins on the Moon and Mercury. Measured wrinkle ridges range from ~112 to 776 m in relief with a mean of ~350 m (median = ~340 m, n = 150) on Mercury and from ~47 to 678 m in relief with a mean of ~198 m (median = ~168 m, n = 150) on the Moon. Wrinkle ridges on Mercury thus are approximately twice as large in mean relief compared to their counterparts on the Moon. The larger scale of Mercury's wrinkle ridges suggests that their formation can be attributed, in part, to global contraction. As global contraction on the Moon is estimated to be an order of magnitude smaller than on Mercury, the smaller scale of lunar wrinkle ridges suggests they most likely form primarily by load induced subsidence of the mare basalt. Wrinkle ridges located in lunar mascon basins and in the Caloris mascon on Mercury are not statistically significantly different in relief than ridges in non-mascon regions, suggesting comparable levels of contractional strain. The fact that mascon basins do not host wrinkle ridges with greater structural relief relative to non-mascon units may indicate the critical role lithospheric thickness plays in controlling subsidence and contraction of thick volcanic sequences on the Moon and Mercury.

Global contraction↗
NASA NTRS2019

Wrinkle ridges on Mercury and the Moon within and outside of mascons

Found on all terrestrial planets, wrinkle ridges are anticlines formed by thrust faulting and folding resulting from crustal shortening. The MErcury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER)spacecraft's orbital phase returned high resolution images and topographic data of the previously unimaged northern high latitudes of Mercury where there are large expanses of smooth plains deformed by wrinkle ridges. Concurrently, the Lunar Reconnaissance Orbiter (LRO) is obtaining high resolution images and topographic data covering lunar mare wrinkle ridges. These data allow quantitative comparison of the scale of wrinkle ridges in smooth plains volcanic units on Mercury with mare wrinkle ridges. We evaluate the topographic relief of 300 wrinkle ridges within and outside of mascon basins on the Moon and Mercury. Measured wrinkle ridges range from ~112 to 776 min relief with a mean of ~350m (median=~340m, n=150) on Mercury and from ~47 to 678 min relief with a mean of ~198 m (median=~168m, n=150) on the Moon. Wrinkle ridges on Mercury thus are approximately twice as large in mean relief compared to their counterparts on the Moon. The larger scale of Mercury's wrinkle ridges suggests that their formation can be attributed, in part, to global contraction. As global contraction on the Moon is estimated to be an order of magnitude smaller than on Mercury, the smaller scale of lunar wrinkle ridges suggests they most likely form primarily by load induced subsidence of the mare basalt. Wrinkle ridges located in lunar mascon basins and in the Caloris mascon on Mercury are not statistically significantly different in relief than ridges in non-mascon regions, suggesting comparable levels of contractional strain. The fact that mascon basins do not host wrinkle ridges with greater structural relief relative to non-mascon units may indicate the critical role lithospheric thickness plays in controlling subsidence and contraction of thick volcanic sequences on the Moon and Mercury.

Lisa S. Schleicher↗
NASA NTRSDate not supplied

A Solar Sail Shape Modeling Approach for Attitude Control Design and Analysis

Solar sails operating in the space environment experience deformations in sail shape that result in relatively large disturbance torques which dictate the required performance of the spacecraft attitude control and momentum management systems. These deformations are driven by thermal loads on the booms (due to uneven solar heating), manufacturing and assembly tolerances, and variations in membrane tension. The Solar Cruiser spacecraft utilizes a four-quadrant sail design with four 30-meter length booms and four triangular sail membranes, creating a square sail structure of >1600 m2. Medium-fidelity mesh models were developed based on a characteristic deformed shape. A series of parametric studies were conducted using this shape paradigm to determine worst-case deformed sail shapes which produce bounding disturbance torques. A large database of shapes was produced, and the forces and moments induced by each individual shape were calculated using a Rios-Reyes reduced order generalized sail model. Two were selected as reference worst-case shapes for the Solar Cruiser mission: one which produced the highest pitch/yaw root-sum-squared (RSS) torque, and one which produced the highest roll torque. The results showed that the worst-case shapes at high solar incidence angles induce significantly higher (2-10x) disturbance torques than an ideal, flat-plate sail. Even with considerable safety margins, assuming an ideal sail is unlikely to sufficiently bound the disturbances, which is critical when designing the attitude control system and sizing actuators. Accurate sail shape modeling methodologies should therefore be employed on future solar sail missions.

Solar Sail↗
NASA NTRSDate not supplied

A Solar Sail Shape Modeling Approach for Attitude Control Design and Analysis

Solar sails operating in the space environment experience deformations in sail shape that result in relatively large disturbance torques which dictate the required performance of the spacecraft attitude control and momentum management systems. These deformations are driven by thermal loads on the booms (due to uneven solar heating), manufacturing and assembly tolerances, and variations in membrane tension. The Solar Cruiser spacecraft utilizes a four-quadrant sail design with four 30-meter length booms and four triangular sail membranes, creating a square sail structure of >1600 m2. Medium-fidelity mesh models were developed based on a characteristic deformed shape. A series of parametric studies were conducted using this shape paradigm to determine worst-case deformed sail shapes which produce bounding disturbance torques. A large database of shapes was produced, and the forces and moments induced by each individual shape were calculated using a Rios-Reyes reduced order generalized sail model [1]. Two were selected as reference worst-case shapes for the Solar Cruiser mission: one which produced the highest pitch/yaw root-sum-squared (RSS) torque, and one which produced the highest roll torque. The results showed that the worst-case shapes at high solar incidence angles induce significantly higher (2-10x) disturbance torques than an ideal, flat-plate sail. Even with considerable safety margins, assuming an ideal sail is unlikely to sufficiently bound the disturbances, which is critical when designing the attitude control system and sizing actuators. Accurate sail shape modeling methodologies should therefore be employed on future solar sail missions. References [1] L. Rios-Reyes and D. J. Scheeres, “Generalized Model for Solar Sails,” Journal of Spacecraft and Rockets, Vol. 42, Jan. 2005, pp. 182–185, 10.2514/1.9054.

Solar Sail↗