Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coherent potential approximation”

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 37 records · Page 2

Orbit Determination Error Analysis Results for the Triana Sun-Earth L2 Libration Point Mission

Using the NASA Goddard Space Flight Center's Orbit Determination Error Analysis System (ODEAS), orbit determination error analysis results are presented for all phases of the Triana Sun-Earth L1 libration point mission and for the science data collection phase of a future Sun-Earth L2 libration point mission. The Triana spacecraft was nominally to be released by the Space Shuttle in a low Earth orbit, and this analysis focuses on that scenario. From the release orbit a transfer trajectory insertion (TTI) maneuver performed using a solid stage would increase the velocity be approximately 3.1 km/sec sending Triana on a direct trajectory to its mission orbit. The Triana mission orbit is a Sun-Earth L1 Lissajous orbit with a Sun-Earth-vehicle (SEV) angle between 4.0 and 15.0 degrees, which would be achieved after a Lissajous orbit insertion (LOI) maneuver at approximately launch plus 6 months. Because Triana was to be launched by the Space Shuttle, TTI could potentially occur over a 16 orbit range from low Earth orbit. This analysis was performed assuming TTI was performed from a low Earth orbit with an inclination of 28.5 degrees and assuming support from a combination of three Deep Space Network (DSN) stations, Goldstone, Canberra, and Madrid and four commercial Universal Space Network (USN) stations, Alaska, Hawaii, Perth, and Santiago. These ground stations would provide coherent two-way range and range rate tracking data usable for orbit determination. Larger range and range rate errors were assumed for the USN stations. Nominally, DSN support would end at TTI+144 hours assuming there were no USN problems. Post-TTI coverage for a range of TTI longitudes for a given nominal trajectory case were analyzed. The orbit determination error analysis after the first correction maneuver would be generally applicable to any libration point mission utilizing a direct trajectory.

Marr, G.↗

Stratospheric Meteorological Conditions for the 3-12 Nov. 1994 ATMOS/ATLAS-3 Measurements

During the ATLAS-3 mission (3-12 Nov 1994), United Kingdom Meteorological Office fields show that the SH vortex was still strong below about 700 K (approximately 16 hPa), with coherent vortex fragments apparent up to about 1300 K (approximately 3 hPa). The SH vortex was shifted off the pole toward 270 deg E throughout ATLAS-3, although its shape varied from day to day. SH temperatures were increasing during ATLAS-3; temperatures below 188 K were last seen 35-45 days before the mission, but were below 195 K as late as 2 Nov 1994. The NH polar vortex had developed at levels above about 550 K (approximately 35 hPa). Simulated high-resolution potential vorticity (PV) fields clearly show low latitude air being drawn up around the polar vortices in both hemispheres. These fields indicate that meteorological analyses underestimate the amount of atmospheric variability, including the strength of local PV gradients, and small-scale structure. Structure such as that in the simulated fields can result in apparent discrepancies between tracer measurements and PV fields.

Manney, G. L.↗

An optical model for composite nuclear scattering

The optical model of composite particle scattering is considered and compared to the accuracies of other models. A nonrelativistic Schroedinger equation with two-body potentials is used for the scattering of a single particle by an energy-dependent local potential. The potential for the elastic channel is composed of matrix elements of a single scattering operator taken between the ground states of the projectile and the target; the coherent amplitude is considered as dominating the scattering in the forward direction. A multiple scattering series is analytically explored and formally summed by the solution of an equivalent Schroedinger equation. Cross sections of nuclear scattering are then determined for He-4 and C-12 nuclei at 3.6 GeV/nucleus and O-16 projectiles at 2.1 GeV/nucleus, and the optical model approximations are found to be consistently lower and more accurate than approximations made by use of Glauber's theory.

Wilson, J. W.↗

Microwave bistatic reflectivity dependence on the moisture content and matric potential of bare soil

Results are presented of an experimental program to determine the functional dependence of the microwave reflectivity of nonvegetated soil surfaces upon volumetric soil moisture and matric potential. A combination evaporation-drainage field experiment was conducted on a bare Captina slit loam with reflectivity, soil moisture content, and matric potential monitored for extended time periods. Results show that for a restricted pressure range (approximately -0.05 to -0.75 bar) there is excellent linear correlation between the log of bistatic reflectivity and both volumetric moisture content and matric potential. Layering effects due to steep moisture content (and matric potential) gradients in the profile are demonstrated to have two distinct and significant effects on the reflectivity response. At near saturation of rough surfaces a very thin dry surface layer appears to modify the effective roughness. This leads to a saturation of reflectivity at high moisture contents. As the surface proceeds to dry further, deeper layers produce coherent interference patterns in the reflectivity response, particularly at the higher frequencies.

Waite, W. P.↗

Integrated optic waveguide devices

Integrated optic waveguide circuits with a phase bias and modulator on the same chip were designed, fabricated, and tested for use in a fiber-optic rotation sensor (gyro) under development. Single mode fiber-optic pigtails were permanently coupled to the four ports of the chip. The switch format was based on coherent coupling between waveguides formed in Z-cut LiNbO3. The control of the coupling was achieved by electro-optically varying the phase propagation constants of each guide. Fiber-to-chip interfacing required the development of appropriate fixturing and manipulation techniques to achieve the close tolerance needed for high coupling efficiency between a fiber with an approximately 5 micron m core and a channel guide with a roughly 2 micron m by 5 micron m cross section. Switch and chip performance at 0.85 micron m is discussed as well as potential improvements related to insertion loss reduction, switching voltages, and suppression of Li2O out-diffusion.

Ramer, O. G.↗

Dynamical Decoupling of Crosstalk on Superconducting Qubit Devices

Current NISQ devices are prone to errors. In order to be used for practical applications or achieve fault-tolerant thresholds, strategies to suppress error rates will be needed to maximize the potential of noisy devices. Dynamical decoupling (DD) is one such strategy for suppressing — or at least alleviating — the effects of decoherence, in which sequences of pulses are applied to qubits to decouple their interaction with the environment. Through experimental runs performed on several Rigetti quantum computing units (QPUs), we first demonstrate that DD is capable of improving coherence times for isolated qubits, as well as suppressing errors caused by the ZZ coupling between pairs of qubits. Extending this framework to cycles containing2-qubit gates, we show that DD can be inserted to decouple qubits from crosstalk occurring during neighboring 2-qubit gates, and demonstrate the efficacy of this procedure on quantum approximate optimization algorithm (QAOA) circuits. We also explore the usage of tailored DD sequences for the suppression of characterized error channels. We are grateful for support from the NASA Ames Research Center and from the DARPA ONISQ program under interagency agreement IAA 8839,Annex 114. HYH is supported by the USRA Feynman QuantumAcademy funded by the NAMS R&D Student Program and a UCHellman Fellowship. JS, ZGI and ZW are supported by USRA NASAAcademic Mission Service (NNA16BD14C).

Dynamical decoupling↗

Wavespace-Based Coherent Deconvolution

Array deconvolution is commonly used in aeroacoustic analysis to remove the influence of a microphone array's point spread function from a conventional beamforming map. Unfortunately, the majority of deconvolution algorithms assume that the acoustic sources in a measurement are incoherent, which can be problematic for some aeroacoustic phenomena with coherent, spatially-distributed characteristics. While several algorithms have been proposed to handle coherent sources, some are computationally intractable for many problems while others require restrictive assumptions about the source field. Newer generalized inverse techniques hold promise, but are still under investigation for general use. An alternate coherent deconvolution method is proposed based on a wavespace transformation of the array data. Wavespace analysis offers advantages over curved-wave array processing, such as providing an explicit shift-invariance in the convolution of the array sampling function with the acoustic wave field. However, usage of the wavespace transformation assumes the acoustic wave field is accurately approximated as a superposition of plane wave fields, regardless of true wavefront curvature. The wavespace technique leverages Fourier transforms to quickly evaluate a shift-invariant convolution. The method is derived for and applied to ideal incoherent and coherent plane wave fields to demonstrate its ability to determine magnitude and relative phase of multiple coherent sources. Multi-scale processing is explored as a means of accelerating solution convergence. A case with a spherical wave front is evaluated. Finally, a trailing edge noise experiment case is considered. Results show the method successfully deconvolves incoherent, partially-coherent, and coherent plane wave fields to a degree necessary for quantitative evaluation. Curved wave front cases warrant further investigation. A potential extension to nearfield beamforming is proposed.

Bahr, Christopher J.↗

Measurements of large-scale density fluctuations in the solar wind using dual-frequency phase scintillations

It is demonstrated that phase-difference scintillations measured with a coherent dual-frequency radio system such as that on Mariner 10 can be used to study the structure of density fluctuations in the solar wind covering a wider range of scale sizes than has ever been possible before. The Mariner 10 observations at solar elongations of 11.5 and 12.6 deg show that the density spectrum in the frequency range from 0.0001 to 0.5 Hz, which corresponds to the spatial wavenumber range of 2 millionths to 0.001 inverse km if the solar wind velocity is assumed to be 350 km/s, is approximately power-law and close to Kolmogorov (spectral index of 11/3). The results are consistent with direct spacecraft observations near earth and provide strong evidence that the density fluctuations are produced by turbulence. The potential and benefits of future extensive measurements are also discussed.

Woo, R.↗

Transport in Carbon Nanotubes with Defects

The linear response conductance of a defect free carbon nanotube connected to ideal contacts is approximately 4e(sup 2)/h. In this paper, we first calculate the conductance characteristics of a uniform tube in the presence of an applied bias by assuring different forms of the potential drop across the wire. It is shown that for long tubes, the conduction mechanism involve a only the four non degenerate subands when the electric field drops slowly across the the tube. We then proceed to calculate the transmission and low bias conductance of tubes with defects. Here, we show that a weak uniform disorder in the tube does not cause a dramatic change in the conductance. However the presence of certain substitutional defects, the transmission coefficient versus energy develops a gap and this can cause a dramatic change in the low bias conductance. The effect of coherence is neglected.

Anatram, M. P.↗

Chasing Shadows: Rotation of the Azimuthal Asymmetry in the TW Hya Disk

We have obtained new images of the protoplanetary disk orbiting TW Hya in visible, total intensity light with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST), using the newly commissioned BAR5 occulter. These HSTSTIS observations achieved an inner working angle of approximately 0 "2, or11.7 au, probing the system at angular radii coincident with recent images of the disk obtained by ALMA and in polarized intensity near-infrared light. By comparing our new STIS images to those taken with STIS in 2000 and with NICMOS in 1998, 2004, and 2005, we demonstrate that TW Hyas azimuthal surface brightness asymmetry moves coherently in position angle. Between 50 au and 141 au we measure a constant angular velocity in the azimuthal brightness asymmetry of 22.deg7 yr(exp -1) in a counterclockwise direction, equivalent to a period of 15.9yrassuming circular motion. Both the (short) inferred period and lack of radial dependence of the moving shadow pattern are inconsistent with Keplerian rotation at these disk radii. We hypothesize that the asymmetry arises from the fact that the disk interior to 1 au is inclined and precessing owing to a planetary companion, thus partially shadowing the outer disk. Further monitoring of this and other shadows on protoplanetary disks potentially opens anew avenue for indirectly observing the sites of planet formation.

protoplanetary disk orbiting TW Hya↗

Millimeter Wave Synthetic Aperture Imaging System with a Unique Rotary Scanning System

In recent years, millimeter wave imaging techniques, using synthetic aperture focusing and holographical approaches, have shown tremendous potential for nondestructive testing applications, involving materials and structures used in space vehicles, including the space shuttle external fuel tank spray on foam insulation and its acreage heat tiles. The ability of signals at millimeter wave frequencies (30 - 300 GHz) to easily penetrate inside of low loss dielectric materials, their relatively small wavelengths, and the possibility of detecting coherent (magnitude and phase) reflections make them suitable for high resolution synthetic aperture focused imaging the interior of such materials and structures. To accommodate imaging requirements, commonly a scanning system is employed that provides for a raster scan of the desired structure. However, most such scanners, although simple in design and construction, are inherently slow primarily due to the need to stop and start at the beginning and end of each scan line. To this end, a millimeter wave synthetic aperture focusing system including a custom-designed transceiver operating at 35 - 45 GHz (Q-band) and unique and complex rotary scanner was designed and developed. The rotary scanner is capable of scanning an area with approximately 80 cm in diameter in less than 10 minutes at step sizes of 3 mm and smaller. The transceiver is capable of producing accurate magnitude and phase of reflected signal from the structure under test. Finally, a synthetic aperture focusing algorithm was developed that translates this rotary-obtained magnitude and phase into a synthetic aperture focusing image of inspected structures. This paper presents the design of the transceiver and the rotary scanning system along with showing several images obtained with this system from various complicated structures.

Ghasr, M. T.↗

Diurnal Cycle of Convection in the East Pacific ITCZ during EPIC-2001

During the last three weeks of September 2001, the EPIC-2001 intensive field campaign focused on studies of deep convection in the ITCZ over the Mexican warm pool region (10N, 95W) of the East Pacific. This study focuses on the pronounced observed diurnal cycle of environmental and convective parameters within the experiment domain. Data from three primary sources are examined: the R/V Ronald H. Brown C-band weather radar, 4-hourly soundings from the Brown and the Global Atmospherics, Inc. National Lightning Detection Network (long range product). Satellite data from TRMM, GOES and OV-1 are also used. The domain boundary layer shows a robust daily evolution of moist enthalpy (as reflect by equivalent potential temperature, theta-e, or wet bulb potential temperature, theta-w), with contributions from changes in both dry and moist entropy. Peak theta-w is found after local nightfall; the average diurnal range of theta-w is approximately 1 deg C. A composite diurnal cycle of convective properties was derived from the C-band volume scans, sampled continuously through the experiment at 10 minute updates. Products derived from the volumetric data include a surface PPI, 15 and 30 dBZ echo top height, vertically integrated liquid, and 6 km (mixed phase region) reflectivity CAPPIs. For almost all products, the parameter means showed virtually no diurnal cycle. However, for the upper-level products, the parameter spectra showed a clear peak in the occurrence of deep/vigorous convection (the "tail end of the distribution") between 7-9 UTC (1-3 AM local), while overall frequency of occurrence peaked later, from 12-15 UTC (6-9 AM local). This represents a daily "outbreak" of isolated deep cells a couple of hours after sunset and subsequent growth, organization and decay through the nighttime hours. The coherence of the diurnal cycle of the convective spectrum is impressive given the wide variety of convective organization observed during the experiment, and given the modulation by passage of 3-5 day easterly waves. While earlier satellite OLR composites suggested an offshore coastal migration of storms into the domain at night, examination of the 150 km and 300 km range radar products showed little evidence of such organization; almost all convection developed "in-place" within the analysis domain. Consistent with the diurnal thermodynamic and microphysical evolution, a clear cycle in cloud-to-ground (CG) lightning occurrence was observed. The local CG diurnal cycle is significantly stronger than the satellite-derived tropical ocean diurnal cycle of total (IC+CG) lightning. Flash rates of 3-4 fl/min were often visually observed after nightfall; these are fairly 'healthy' flash rates for tropical ocean storms, and the domain was qualitatively noted to be unusually lightning-productive by the R/V Brown crew (also consistent with satellite-based climatologies).

Boccippio, Dennis J.↗

Long-Term Space Astrophysics Program

This is the final report for our Long-Term Space Astrophysics Program (NRA 94-OSS-12) grant NAG 5-3225. The proposal is entitled 'Spectral and Temporal Properties of Black Hole Candidates', and began funding in May 1995, and ran through 31 Aug 2000. The project summary from the original proposal was as follows: 'We will study the spectral and temporal properties of black hole candidates (BHC) by using data from archival sources (e.g., EXOSAT, Ginga, ROSAT) and proposed follow-up observations with modern instruments (e.g., ASCA, XTE). Our spectral studies will focus on identifying the basic characteristics and luminosities of the emission components in the various 'states' of BHC. We hope to understand and quantify the global energetics of these states. Our temporal studies will focus on expanding and classifying our knowledge of BHC variability properties in each state. We will explore the nature of quasi-periodic oscillations in BHC. We will combine our spectral and temporal studies by analyzing time lags and variability coherence between energy channels. In addition, we will investigate ways of correlating observed variability behavior with specific emission components.' We have accomplished many of these goals laid out within the original proposal. As originally proposed, we have utilized both archival and proprietary satellite data. In terms of archival data, we have utilized data from the Advanced Satellite for Cosmology and Astrophysics (ASCA), ROSAT, and the Rossi X-ray Timing Explorer (RXTE). We also obtained proprietary data from ASCA, RXTE, and the Extreme Ultraviolet Explorer (EUVE). In terms of sources, we have examined a wide variety of both galactic black hole candidates and extra-galactic black holes. For the galactic black holes we have observed and analyzed both the low/hard state and the high/soft state. We have performed both spectral and timing analyses on all of these objects. In addition, we have also examined a number of neutron stars or potential neutron stars. All of our research on the above mentioned objects has resulted in one or more publications in peer-reviewed journals. Attached is a list of refereed publications of research results which have been funded by this grant over approximately the past five and a half years. In addition, we have included a list of conference proceedings and other similar reports that have been associated with this grant.

Nowark, Michael A.↗

Remote Sensing of Aerosol Backscatter and Earth Surface Targets By Use of An Airborne Focused Continuous Wave CO2 Doppler Lidar Over Western North America

Airborne lidar systems are used to determine wind velocity and to measure aerosol or cloud backscatter variability. Atmospheric aerosols, being affected by local and regional sources, show tremendous variability. Continuous wave (cw) lidar can obtain detailed aerosol loading with unprecedented high resolution (3 sec) and sensitivity (1 mg/cubic meter) as was done during the 1995 NASA Multi-center Airborne Coherent Atmospheric Wind Sensor (MACAWS) mission over western North America and the Pacific Ocean. Backscatter variability was measured at a 9.1 micron wavelength cw focused CO2 Doppler lidar for approximately 52 flight hours, covering an equivalent horizontal distance of approximately 30,000 km in the troposphere. Some quasi-vertical backscatter profiles were also obtained during various ascents and descents at altitudes that ranged from approximately 0.1 to 12 km. Similarities and differences for aerosol loading over land and ocean were observed. Mid-tropospheric aerosol backscatter background mode was approximately 6 x 10(exp -11)/ms/r, consistent with previous lidar datasets. While these atmospheric measurements were made, the lidar also retrieved a distinct backscatter signal from the Earth's surface from the unfocused part of the focused cw lidar beam during aircraft rolls. Atmospheric backscatter can be highly variable both spatially and temporally, whereas, Earth-surface backscatter is relatively much less variant and can be quite predictable. Therefore, routine atmospheric backscatter measurements by an airborne lidar also give Earth surface backscatter which can allow for investigating the Earth terrain. In the case where the Earth's surface backscatter is coming from a well-known and fairly uniform region, then it can potentially offer lidar calibration opportunities during flight. These Earth surface measurements over varying Californian terrain during the mission were compared with laboratory backscatter measurements using the same lidar of various Earth surfaces giving good agreement, suggesting that the lidar efficiency, and thus a lidar calibration factor for detection, can be estimated fairly well using Earth's surface signal.

Jarzembski, Maurice A.↗

Investigation of Unsteady Pressure-Sensitive Paint (uPSP) and a Dynamic Loads Balance to Predict Launch Vehicle Buffet Environments

This NESC assessment examined the accuracy of estimating buffet loads on in-line launch vehicles without booster attachments using sparse unsteady pressure measurements. The buffet loads computed using sparse sensor data were compared with estimates derived using measurements with much higher spatial resolution. The current method for estimating launch vehicle buffet loads is through wind tunnel testing of models with approximately 400 unsteady pressure transducers. Even with this relatively large number of sensors, the coverage can be insufficient to provide reliable integrated unsteady loads on vehicles. In general, sparse sensor spacing requires the use of coherence-length-based corrections in the azimuthal and axial directions to integrate the unsteady pressures and obtain reasonable estimates of the buffet loads. Coherence corrections have been used to estimate buffet loads for a variety of launch vehicles with the assumption methodology results in reasonably conservative loads. For the Space Launch System (SLS), the first estimates of buffet loads exceeded the limits of the vehicle structure, so additional tests with higher sensor density were conducted to better define the buffet loads and possibly avoid expensive modifications to the vehicle design. Without the additional tests and improvements to the coherence-length analysis methods, there would have been significant impacts to the vehicle weight, cost, and schedule. If the load estimates turn out to be too low, there is significant risk of structural failure of the vehicle. This assessment used a combination of unsteady pressure-sensitive paint (uPSP), unsteady pressure transducers, and a dynamic force and moment balance to investigate the integration schemes used with limited unsteady pressure data by comparing them with direct integration of extremely dense fluctuating pressure measurements. An outfall of the assessment was to evaluate the potential of using the emerging uPSP technique in a production test environment for future launch vehicles. The results show that modifications to the current technique can improve the accuracy of buffet estimates. More importantly, the uPSP worked remarkably well and, with improvements to the frequency response, sensitivity, and productivity, will provide an enhanced method for measuring wind tunnel buffet forcing functions (BFFs).

Schuster, David M.↗

The Detection of Change in the Arctic Using Satellite and Buoy Data

The decade of the 1990s is the warmest decade of the last century while the year 1998 is the warmest year ever observed by modern techniques with 9 out of 12 months of the year being the warmest month. Since the Arctic is expected to provide early signals of a possible warming scenario, detailed examination of changes in the Arctic environment is important. In this study, we examined available satellite ice cover and surface temperature data, wind and pressure data, and ocean hydrographic data to gain insights into the warming phenomenon. The areas of open water in both western and eastern regions of the Arctic were found to follow a cyclical pattern with approximately decadal period but with a lag of about three years between the two regions. The pattern was interrupted by unusually large anomalies in open water area in the western region in 1993 and 1998 and in the eastern region in 1995. The big 1998 open water anomaly occurred at the same time when a large surface temperature anomaly was also occurring in the area and adjacent regions. The infrared temperature data show for the first time the complete spatial scope of the warming anomalies and it is apparent that despite the magnitude of the 1998 anomaly, it is basically confined to North America and the Western Arctic. The large increases in open water areas in the Western Sector form 1996 to 1998 were observed to be coherent with changing wind directions which was predominantly cyclonic in 1996 and anti-cyclonic in 1997 and 1998. Detailed hydrography measurements up to 500 m depth over the same general area in April 1996 and April 1997 also indicate significant freshening and warming in the upper part of the mixed layer suggesting increases in ice melt. Continuous ocean temperature and salinity data from ocean buoys confirm this result and show significant seasonal changes from 1996 to 1998, at depths of 8 m, 45 m, and 75 m. Long data records of temperature and hydrography were also examined and the potential impact of a warming, freshening, and the presence of abnormally large open areas on the state of the Arctic climate system are discussed.

Comiso, Josefino C.↗

Hypervelocity impact survivability experiments for carbonaceous impactors, part 2

Hypervelocity impact experiments were performed to further test the survivability of carbonaceous impactors and to determine potential products that may have been synthesized during impact. Diamonds were launched by the Ames two-stage light gas gun into Al plate at velocities of 2.75 and 3.1 km sec(exp -1). FESEM imagery confirms that diamond fragments survived in both experiments. Earlier experiments found that diamonds were destroyed on impact above 4.3 km sec(exp -1). Thus, the upper stability limit for diamond on impact into Al, as determined from our experimental conditions, is between 3.1 and 4.3 km sec(exp -1). Particles of the carbonaceous chondrite Nogoya were also launched into Al at a velocity of 6.2 km sec (exp -1). Laser desorption (L (exp 2) MS) analyses of the impactor residues indicate that the lowest and highest mass polycyclic aromatic hydrocarbons (PAH's) were largely destroyed on impact; those of intermediate mass (202-220 amu) remained at the same level or increased in abundance. In addition, alkyl-substituted homologs of the most abundant pre-impacted PAH's were synthesized during impact. These results suggest that an unknown fraction of some organic compounds can survive low to moderate impact velocities and that synthesized products can be expected to form up to velocities of, at least, 6.5 km sec(exp -1). We also present examples of craters formed by a unique microparticle accelerator that could launch micron-sized particles of almost any coherent material at velocities up to approximately 15 km sec(exp -1). Many of the experiments have a direct bearing on the interpretation of LDEF craters.

Bunch, T. E.↗

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗