Engineering topics
Gustafson, B.
Publications and source records attributed to Gustafson, B..
Dynamical evolution of dust particles in the Kuiper disk
A Kuiper belt dust disk will have a resonant structure, arising because the Plutinos are in the 3:2 mean motion resonance with Neptune. We run numerical integrations of particles originating from Plutinos to determine what percentage of particles remain in the resonance for a variety of particle and source body sizes. The dynamical evolution of the particles is followed from source to sink with Poynting-Robertson light drag, solar wind drag, radiation pressure, the Lorentz force, neutral interstellar gas drag, and the effects of planetary gravitational perturbations included. The number of particles in the 3:2 resonance increases with decreasing p for the cases where the initial source bodies are small and the percentage of particles in resonance is not significantly changed by either the addition of the Lorentz force, as long as the potential of the particles is small (U = 5 V), or the effect of neutral interstellar gas drag.
Origin and evolution of the zodiacal dust cloud
The astrophysical importance of the zodiacal cloud became more apparent. The most useful source of information on the structure of the zodiacal cloud is the Infrared Astronomical Satellite (IRAS) observations. A substantial fraction of the extensive IRAS data set was analyzed. Also, a numerical model was developed (SIMUL) that allows to calculate the distribution of night-sky brightness that would be produced by any particular distribution of dust particle orbits. This model includes the effects of orbital perturbations by both the planets and solar radiation, it reproduces the exact viewing geometry of the IRAS telescope, and allows for the eccentricity of the Earth's orbit. SIMUL now is used to model not just the solar system dust bands discovered by IRAS but the whole zodiacal cloud.
The choice of numerical boundary conditions for hyperbolic systems
Two fundamental problems for mixed initial boundary value problems with applications in fluid mechanics are discussed. First, different stability properties are discussed which are of importance for long time integrations and steady state calculations. Secondly, a new numerical technique for problems with an artificial boundary is introduced.