Search NASA⌕ Search

SEARCH · Search NASA

Results for “OBLATENESS”

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 145 records · Page 8

Computation of optimal low- and medium-thrust orbit transfers

This paper presents the formulation of the optimal low- and medium-thrust orbit transfer control problem, numerical methods for solution, and numerical solutions of the problem. The problem formulation is for final mass maximization and allows for second-harmonic oblateness, atmospheric drag, and 3D noncoplanar nonaligned elliptic terminal orbits. We set up examples to demonstrate the ability of two indirect methods to solve the resulting two point boundary value problems (TPBVP). The methods demonstrated are the multiple point shooting method as formulated in Oberle's (1987) subroutine BOUNDSCO, and the minimizing boundary-condition method (MBCM). We find that although both methods can converge solutions, there are tradeoffs to using either method. We present numerical solutions of planar transfers in which both the initial orbit exit and final orbit entry points have been optimized. These solutions include two- and three-burn transfers. The methods used show an ability to handle thrust down to at least T/W(o) = O(10 exp -3). They also show similar convergence abilities with or without the oblateness and drag terms. We discuss the issue of maximizing with respect to the final time and provide evidence that implies a local optimum at a maximum final time for a given number of burns.

Chuang, C.-H.↗

Tidal interaction: A possible explanation for geysers and other fluid phenomena in the Neptune-Triton system

Discovery of geyser-like plumes on the surface of Triton was a highlight of Voyager 2's passage through the Neptune planetary system. Remarkable as these observations were, they were not entirely without precedent. Considering the confirmed predictions for the 1979 Voyager Jovian passage, it was logical to consider other solar system bodies beside Io where tidal effects could be a significant factor in surface processes. It was our intuition that the Neptune-Triton gravitational bond acting at high inclination to the Neptune equator and the fact that Neptune was a fluid body was significant oblateness would produce tidal and mechanical forces that could be transformed into thermal energy vented on Triton's surface. Prior to the Voyager flyby, others have noted that capture and evolution of Triton's orbit from extreme eccentricity to near circular state today would have resulted in significant tidal heating, but these analysts disregard current day forces. Our calculations indicate that the time varying forces between Neptune-Triton fall midway between those exerted in the Earth-Moon and Jupiter-Io systems, and considering the low level of other energy inputs, this source of internal energy should not be ignored when seeking an explanation for surface activity. In each planet-satellite case, residual or steady-state eccentricity causes time-varying stresses on internal satellite strata. In the case of Jupiter the residual eccentricity is due largely to Galilean satellite interactions, particularly Io-Europa, but in the case of Neptune-Triton, it is the effect of Triton's inclined orbit about an oblate primary.

Kelly, W. D.↗

Intensity and polarization of light scattered by size distributions of randomly oriented nonspherical particles

Calculations of light scattering by small particles are important in many diverse fields of science and engineering. In many cases of practical interest, scattering particles are nonspherical and are distributed over sizes and orientations. However, accurate light scattering computations for ensembles of nonspherical particles are difficult and time-consuming, and the literature in which such calculations are reported is rather scarce. In this paper, the T-matrix approach, as extended recently to randomly oriented particles, is used to calculate rigorously light scattering by size distributions of randomly oriented axially symmetric particles. To model the variation of particle sizes in real ensembles, we use a power law distribution typical of some terrestrial aerosols. Contour plots of intensity and degree of linear polarization for polydisperse prolate and oblate spheroids of different aspect ratios and effective equivalent-sphere size parameters from 0 to 10 are calculated and compared with calculations for equivalent spheres. The angular scattering behavior of nonspherical polydispersions is found to be greatly different from that of spheres, while the scattering properties of oblate and prolate spheroids of the same aspect ratio are similar. With increasing particle size, both intensity and polarization become more shape-dependent. In general, nonspherical particles are stronger side scatterers and weaker backscatterers than equivalent spheres. With increasing aspect ratio of nonspherical particles polarization tends to be predominantly positive. Possible effects of particle nonsphericity on optical remote sensing of atmospheric aerosols are discussed.

Mishchenko, M. I.↗

Contributions to the Earth's obliquity rate, precession, and nutation

The precession and nutation of the Earth's equator arise from solar, lunar, and planetary torques on the oblate Earth. The mean lunar orbit plane is nearly coincident with the ecliptic plane. A small tilt out of the ecliptic is caused by planetary perturbations and the Earth's gravitational harmonic J(sub2). These planetary perturbations on the lunar orbit result in torques on the oblate Earth which contribute to precession, obliquity rate, and nutation while the J(sub 2) perturbations contribute to precession and nutation. Small additional contributions to the secular rates arise from tidal effects and planetary torques on the Earth's bulge. The total correction to the obliquity rate is -0.024sec/century, it is an observable motion in space (the much larger conventional obliquity rate is wholly from the motion of the ecliptic, not the equator), and it is not present in the IAU-adopted expressions for the orientation of the Earth's equator. The J(sub2) effects have generally been allowed for in past nutation theories and some procession theories. For the planetary effect, the contributions to the 18.6 yr nutation are -0.03 mas (milliarcseconds) for the in-phase Delta phi plus out-of-phase contributions of 0.14 mas in Delta phi and -0.03 mas in Delta epsilon. The latter terms demonstrate that out-of-phase contributions can arise by means other than dissipation. The sum of the contributions to the precession rate is considered and the inferred value of the moment of inertia combination (C-A)/C, which is used to scale the coefficients in the nutation series, is evaluated. Using an updated value for the precession rate, the rigid body (C-A)/C = 0.003 273 763 4 which, in combination with a satellite-derived J(sub2), gives a normalized polar moment of inertia C/MR(exp2) = 0.330 700 7. The planetary contributions to the precession and obliquity rates are not constant for long times causing accelerations in both quantities. Acceleration in precession also arises from tides and changing J(sub2). Contributions from the improved theory, masses, ecliptic motion, and measured values of the precession rate and obliquity are combined to give expressions (polynomials in time) for precession, obliquity, and Greenwich Mean Sidereal Time.

Williams, James G.↗

The intrinsic shape of NGC 3379

Photometric and kinematic data from the literature are combined with new dynamical models to derive the intrinsic shape of the 'standard' elliptical galaxy NGC 3379. The parameters that are best constrained are the dynamical triaxiality T (essentially the triaxiality of the total mass distribution) and the short-to long axis ratio of the light distribution c(sub L). The inferred shape is given by a Bayesian probability distribution in the (T, c(sub L) plane. Assuming a uniform prior, the most probable shape is oblate with a flattening of c(sub L) = 87. The distribution is strongly non-Gaussian, however, and the expectation values, (T) = .31 (c(sub L) = .75, imply a flatter and more triaxial figure. The 68% highest posterior density region allows more triaxial shapes as long as they are fairly round, or flatter shapes as long as they are nearly oblate. These results are essentially unchanged if the galaxy is assumed to rotate about its short axis, or if it is modeled as an S0 with a negligible-mass disk rather than as an elliptical. The suggestion of Capaccioli et al. (ApJ, 371, 535 (1991)) that NGC 3379 is a rather flat, triaxial S0 galaxy is found to be improbable at the 98% level; this conclusion is largely independent of the bulge-to-disk ratio or the relative rotation speeds of the two components.

Statler, Thomas S.↗

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗

Determination of Azimuth Angle at Burnout for Placing a Satellite Over a Selected Earth Position

Expressions are presented for relating the satellite position in the orbital plane with the projected latitude and longitude on a rotating earth surface. An expression is also presented for determining the azimuth angle at a given burnout position on the basis of a selected passage position on the earth's surface. Examples are presented of a satellite launched eastward and one launched westward, each passing over a selected position sometime after having completed three orbits. Incremental changes from the desired latitude and longitude due to the earth's oblateness are included in the iteration for obtaining the azimuth angles of the two examples. The results for both cases are then compared with those obtained from a computing program using an oblate rotating earth. Changes from the selected latitude and longitude resulting from incremental changes from the burn-out azimuth angle and latitude are also analyzed.

Skopinski, T. H.↗

Contributions to the Earth's Obliquity Rate, Precession, and Nutation

The precession and nutation of the Earth's equator arise from solar, lunar, and planetary torques on the oblate Earth. The mean lunar orbit plane is nearly coincident with the ecliptic plane. A small tilt out of the ecliptic is caused by planetary perturbations and the Earth's gravitational harmonic J(sub 2). These planetary perturbations on the lunar orbit result in torques on the oblate Earth which contribute to precession, obliquity rate, and nutation while the J(sub 2) perturbations contribute to precession and nutation. Small additional contributions to the secular rates arise from tidal effects and planetary torques on the Earth's bulge. The total correction to the obliquity rate is -0.024 sec/century, it is an observable motion in space (the much larger conventional obliquity rate is wholly from the motion of the ecliptic, not the equator), and it is not present in the IAU-adopted expressions for the orientation of the Earth's equator. The effects have generally been allowed for in past nutation theories and some precession theories. For the planetary effect, the contributions to the 18.6 yr nutation are -0.03 mas (milliarcseconds) for the in-phase Delta(psi) plus out-of-phase contributions of 0.14 mas in Delta(psi) and -0.03 mas in Delta(sub epsilon). The latter terms demonstrate that out-of-phase contributions can arise by means other than dissipation. The sum of the contributions to the precession rate is considered and the inferred value of the moment of inertia combination (C-A)/C, which is used to scale the coefficients in the nutation series, is evaluated. Using an updated value for the precession rate, the rigid body (C-A)/C =0.003 273 763 4 which, in combination with a satellite-derived J(sub 2), gives a normalized polar moment of inertia C/MR(exp 2) = 0.330 700 7. The planetary contributions to the precession and obliquity rates are not constant for long times causing accelerations in both quantities. Acceleration in precession also arises from tides and changing J(sub 2) Contributions from the improved theory, masses, ecliptic motion, and measured values of the precession rate and obliquity are combined to give expressions (polynomials in time) for precession, obliquity, and Greenwich Mean Sidereal Time.

Williams, James G.↗

Molded Magnetic Article

A molded magnetic article and fabrication method are provided. Particles of ferromagnetic material embedded in a polymer binder are molded under heat and pressure into a geometric shape. Each particle is an oblate spheroid having a radius-to-thickness aspect ratio approximately in the range of 15-30. Each oblate spheroid has flattened poles that are substantially in perpendicular alignment to a direction of the molding pressure throughout the geometric shape.

Bryant, Robert G.↗

Lunar Prospector Orbit Determination Uncertainties Using the High Resolution Lunar Gravity Models

The Lunar Prospector (LP) mission began on January 6, 1998, when the LP spacecraft was launched from Cape Canaveral, Florida. The objectives of the mission were to determine whether water ice exists at the lunar poles, generate a global compositional map of the lunar surface, detect lunar outgassing, and improve knowledge of the lunar magnetic and gravity fields. Orbit determination of LP performed at the Jet Propulsion Laboratory (JPL) is conducted as part of the principal science investigation of the lunar gravity field. This paper will describe the JPL effort in support of the LP Gravity Investigation. This support includes high precision orbit determination, gravity model validation, and data editing. A description of the mission and its trajectory will be provided first, followed by a discussion of the orbit determination estimation procedure and models. Accuracies will be examined in terms of orbit-to-orbit solution differences, as a function of oblateness model truncation, and inclination in the plane-of-sky. Long term predictions for several gravity fields will be compared to the reconstructed orbits to demonstrate the accuracy of the orbit determination and oblateness fields developed by the Principal Gravity Investigator.

Carranza, Eric↗

Coupled Model Simulation of Snowfall Events Over the Black Hills

Although many long-term simulations of snow accumulation and oblation have been made using stand-alone land surface models and surface models coupled with GCMs, less research has focused on short-term event simulations. Actually, accurate event simulations of snow-related processes are the basis for successful long-term simulation. Three advantages of event simulations of snowfall and snow melting are availability of: (1) intensive observation data from field experiments for validation; (2) more physically-realistic precipitation schemes for use in atmospheric models to simulate snowfall; and (3) a more detailed analysis of the snow melting processes. In addition to the complexities of snow related processes themselves, terrain-induced effects on snowfall/snow melting make simulations of snow events more difficult. Climatological observations indicate that terrain features such as the Black Hills of South Dakota and Wyoming can exert important effects on snow accumulation and snow oblation processes. One of the primary effects is that the orography causes forced uplift of airflow and causes atmospheric waves to form both upwind and downwind of it. Airflow often splits around the obstacle, converging on the lee side. This convergence may lead to precipitation enhancement. It also provides an elevated heat and moisture source that enhances atmospheric instability. During the period of April 5-May 5, 1999, the Upper Missouri River Basin Pilot Project (UMRBPP) made intensive observations on precipitation events occurring in the Black Hills. Two moderate snowfall events were captured during the period. The resulting high temporal and spatial resolution data provides opportunities to investigate terrain effects on snowfall amount, distribution, and melting. Successful simulation of snowfall amount, distribution, and evolution using atmospheric models is important to subsequent modeling of snow melting using snow sub-models in land surface schemes. In this paper, a coupled model system, consisting of an atmosphere model (ARPS) and a land-surface model (revised NCAR LSM), is used to simulate one of these UMRBPP cases.

Wang, Jianzhong↗

Influence of a Fluid Lunar Core on the Moons Orientation

Oblateness of and dissipation at the lunar liquid-core/solid-mantle boundary affects the precession of core and mantle. Analysis of Lunar Laser ranges gives a weak detection of oblateness and a strong determination of dissipation. Additional information is contained in the original extended abstract.

Williams, J. G.↗

Fabrication of Molded Magnetic Article

A molded magnetic article and fabrication method are provided. Particles of ferromagnetic material embedded in a polymer binder are molded under heat and pressure into a geometric shape. Each particle is an oblate spheroid having a radius-to-thickness aspect ratio approximately in the range of 15-30. Each oblate spheroid has flattened poles that are substantially in perpendicular alignment to a direction of the molding pressure throughout the geometric shape.

Robert G Bryant↗

Organization of P, S, and Fe Inclusions in a Freshwater Magnetococcus

Magnetotactic bacteria are a heterogeneous group of motile, mainly aquatic procaryotes that align and swim along geomagnetic field lines. They are of interest to astrobiologists because of the magnetite crystals found in the Mars meteorite ALH84001 which share many characteristics with the magnetite produced intracellularly by magnetotactic bacteria. These bacteria are diverse morphologically, physiologically and phylogenetically, sharing a few key characteristics: 1) the presence of intracellular membrane-bound magnetic crystals (magnetosomes), usually but not always arranged in chains; 2) motility by means of flagella; and 3) microaerophillic or anaerobic physiology. The bilophotrichous (having two flagella bundles) magnetotactic cocci (MC) are ubiquitous in aquatic habitats but have proven extremely difficult to cultivate. Because only several marine strains have been isolated and grown in axenic culture, little is known about the physiology and the biogeochemical roles of the MC. We studied the composition and distribution of intracellular structures in an uncultured MC, designated ARB-1. To do this, a combination of light microscopy, environmental scanning electron microscopy (ESEM), transmission electron microscopy (TEM) and energy dispersive X-ray spectroscopy (EDS) were used. Cells of ARB-1 were separated from sediments collected from Baldwin Lake (Los Angeles Arboretum, Arcadia, CA). They are large spherical to oblate spheroidal Gram-negative cells, ranging from 1 to 4 micrometers along the maximum dimension, which is perpendicular to the direction of swimming. Cells have two large phosphorus-containing inclusions that comprise a large percentage of the cell volume. Many smaller sulfur inclusions are located at the convex end of the cell. Most of the cellular Fe is present in the magnetosomes. These may be arranged as a clump at the concave end of the cell, near the two flagella bundles, or as chains, or as both a clump and chains. The magnetosomes were identified as magnetite (Fe3O4) by selected area electron diffraction (SAED) and high resolution TEM. We saw a trend between cell size and organization of the magnetosomes. Smaller, more spherical cells were more likely to have chains than were larger, more oblate cells. This may indicate different populations of cells, or it may be attributed to variations in cell growth cycle. The size distribution (length) of magnetosomes in chains was similar to that of magnetosomes in clusters, except that there was a larger size range for clustered magnetosomes. Magnetosomes from ARB-1 cells average 82 nm in length. If plotted on a graph of length as a function of aspect ratio, they fall within the single domain region of the plot. If compared with the size distributions of magnetite from ALH84001 and magnetosomes from the cultured magnetotactic vibrio MV-1, the magnetites produced by ARB-1 cells are, on average, larger and have a wider range of aspect ratio. ARB-1 cells have a specific organization of the P, S, and Fe inclusions. The P inclusions always occupy the majority of the cell volume and separate the S inclusions from the disorganized clumps of magnetosomes and the flagella bundles. The P inclusions may contain polyphosphate, which could play several roles in motility, adaptation to stress, growth and division, buoyancy, and energy. The S inclusions might be a way to store S, a potential energy source, when the cells move from sulfide to oxygen zones. The consistency of P, S, and Fe organization in ARB-1 cells suggests that these inclusions have some specific and interactive functions.

Cox, Lea↗

Small-Body Extensions for the Satellite Orbit Analysis Program (SOAP)

An extension to the SOAP software allows users to work with tri-axial ellipsoid-based representations of planetary bodies, primarily for working with small, natural satellites, asteroids, and comets. SOAP is a widely used tool for the visualization and analysis of space missions. The small body extension provides the same visualization and analysis constructs for use with small bodies. These constructs allow the user to characterize satellite path and instrument cover information for small bodies in both 3D display and numerical output formats. Tri-axial ellipsoids are geometric shapes the diameters of which are different in each of three principal x, y, and z dimensions. This construct provides a better approximation than using spheres or oblate spheroids (ellipsoids comprising two common equatorial diameters as a distinct polar diameter). However, the tri-axial ellipsoid is considerably more difficult to work with from a modeling perspective. In addition, the SOAP small-body extensions allow the user to actually employ a plate model for highly irregular surfaces. Both tri-axial ellipsoids and plate models can be assigned to coordinate frames, thus allowing for the modeling of arbitrary changes to body orientation. A variety of features have been extended to support tri-axial ellipsoids, including the computation and display of the spacecraft sub-orbital point, ground trace, instrument footprints, and swathes. Displays of 3D instrument volumes can be shown interacting with the ellipsoids. Longitude/latitude grids, contour plots, and texture maps can be displayed on the ellipsoids using a variety of projections. The distance along an arbitrary line of sight can be computed between the spacecraft and the ellipsoid, and the coordinates of that intersection can be plotted as a function of time. The small-body extension supports the same visual and analytical constructs that are supported for spheres and oblate spheroids in SOAP making the implementation of the more complex algorithms largely transparent to the user.

Carnright, Robert↗

The Thermal Evolution and Internal Structure of Saturn's Midsize Icy Satellites

The Cassini-Huygens mission is returning new geophysical data for the midsize, icy satellites of Saturn (i.e., satellites with radii between 100 and 1,000 km). These data have enabled a new generation of geophysical model studies for Phoebe, Iapetus, Rhea, Mimas, Tethys, Dione, as well as Enceladus (see Spencer et al. 2009). Here we consider the new model studies that have reported significant results elucidating the evolutionary histories and internal structures of these satellites. Those results have included their age, the development of their internal structures and mineralogies, which for greatest fidelity must be done concomitantly with coupled dynamical evolutions. Surface areas, volumes, bulk densities, spin rates, orbit inclinations, eccentricities, and distance from Saturn have changed as the satellites have aged. Heat is required to power the satellites evolution, but is not overly abundant for the midsized satellites. All sources of heat must be evaluated and taken into account. This includes their intensities and when they occur and are available to facilitate evolution, both internal and dynamical. The mechanisms of heat transport must also be included. However, to model these to high fidelity the material properties of the satellite interiors must be accurately known. This is not the case. Thus, we discuss what is known about these properties and how the uncertainties affect the estimation of heat sources, transport processes, and the consequential changes in composition and evolution. Phoebe has an oblate shape that may be in equilibrium with its spin period of ~9.3 h. Its orbital properties suggest that it is not one of the regular satellites, but is a captured body. Its density is higher than that of the other satellites, consistent with formation in the solar nebula rather than from material around Saturn. Oblate shape and high density are unusual for objects in this size range, and may indicate that Phoebe was heated by Al-26 decay soon after its formation, which is consistent with some models of the origin of Kuiper belt objects. Iapetus has the shape of a hydrostatic body with a rotation period near 16 h. It subsequently despun to its current synchronous rotation state, ~79 day period. These observations are sufficient to constrain the required heating in Iapetus early history, suggesting that it formed several My after CAI condensation. Because Saturn had to be present for Iapetus to form, this date also constrains the age of Saturn and how long it took to form. Both shape and gravitational data are available for Rhea. Gravity data were obtained from the single Cassini flyby during the prime mission and within the uncertainties cannot distinguish between hydrostatic and non-hydrostatic gravitational fields. Both Dione and Tethys display evidence of smooth terrains, with Dione's appearing considerably younger. Both are conceivably linked to tidal heating in the past, but the low rock abundance within Tethys and the lack of eccentricity excitation of Tethys orbit today make explaining this satellite's geology challenging.

Matson, Dennis L.↗

Polarization of Directly Imaged Young Giant Planets as a Probe of Mass, Rotation, and Clouds

Young, hot gas giant planets at large separations from their primaries have been directly imaged around several nearby stars. More such planets will likely be detected by ongoing and new imaging surveys with instruments such as the Gemini Planet Imager (GPI). Efforts continue to model the spectra of these planets in order to constrain their masses, effective temperatures, composition, and cloud structure. One potential tool for analyzing these objects, which has received relatively less attention, is polarization. Linear polarization of gas giant exoplanets can arise from the combined influences of light scattering by atmospheric dust and a rotationally distorted shape. The oblateness of gas giant planet increases of course with rotation rate and for fixed rotation also rises with decreasing gravity. Thus young, lower mass gas giant planets with youthful inflated radii could easily have oblateness greater than that of Saturn s 10%. We find that polarizations of over 1% may easily be produced in the near-infrared in such cases. This magnitude of polarization may be measurable by GPI and other instruments. Thus if detected, polarization of a young Jupiter places constraints on the combination of its gravity, rotation rate, and degree of cloudiness. We will present results of our multiple scattering analysis coupled with a self-consistent dusty atmospheric models to demonstrate the range of polarizations that might be expected from resolved exoplanets and the range of parameter space that such observations may inform.

Marley, Mark Scott↗

Retrieval of Microphysical Characteristics of Particles in Atmospheres of Distant Comets from Ground-Based Polarimetry

We summarize unique aperture data on the degree of linear polarization observed for distant comets C/2010 S1, C/2010 R1, C/2011 KP36, C/2012 J1, C/2013 V4, and C/2014 A4 with heliocentric distances exceeding 3 AU. Observations have been carried out at the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences (Nizhnij Arkhyz, Russia) during the period from 2011 to 2016. The measured negative polarization proves to be significantly larger in absolute value than what is typically observed for comets close to the Sun. We compare the new observational data with the results of numerical modeling performed with the T-matrix and superposition T-matrix methods. In our computer simulations, we assume the cometary coma to be an optically thin cloud containing particles in the form of spheroids, fractal aggregates composed of spherical monomers, and mixtures of spheroids and aggregate particles. We obtain a good semi-quantitative agreement between all polarimetric data for the observed distant comets and the results of numerical modeling for the following models of the cometary dust: (i) a mixture of submicrometer water-ice oblate spheroids with aggregates composed of submicrometer silicate monomers; and (ii) a mixture of submicrometer water-ice oblate spheroids and aggregates consisting of both silicate and organic monomers. The microphysical parameters of these models are presented and discussed.

Numerical modeling↗