Experimentation and simulation of meteoroid ablation
Abstract/Contents
- Abstract
- Numerous satellite failures over the years have shown the various dangers posed to man made objects by the space environment. Meteoroids historically represent one of the foremost external phenomena to cause catastrophic loss of spacecrafts. With velocities between 11~km/s and 72~km/s, meteoroids with masses above 1~mg are easily capable of penetrating the hulls of satellites. Such an encounter is likely to cause irreparable mechanical damage and possibly electrical damage to the spacecraft and cause end of mission. Although these type of events happen infrequently, they are a major concern that needs to be taken into account when deploying payload in Lower Earth Orbit (LEO). In addition to presenting engineering risks to our space assets, meteoroids are the primary source of extra-terrestrial material to Earth. Meteoroid ablation deposits between 5 and 300 tons of material in our atmosphere every day. This daily input is responsible for multiple phenomena in the Mesospheric and Lower Thermospheric region (MLT) such as the creation of sporadic E layers, the creation of Noctilucent clouds, and the formation of stratospheric aerosol. A fundamental parameter to understand how meteoroids are responsible for the phenomena described above is their mass. Unfortunately, the only way we can observe meteoroids is via the plasma and light formation --- called a meteor --- they create as they descend through the atmosphere. Over the years, the scientific community has primarily employed two instruments to understand and estimate the Earth-bound meteoroid mass flux from the observation of meteors, namely radars and (CCD) cameras. Mass estimation from optical devices is based on the idea of energy conservation: the energy lost during the entry process of a meteoroid is converted to light. This relation is mediated via a parameter called the luminous efficiency, $\tau$. Radar mass (RM) relies instead on a more refined theoretical basis which, similarly to optical masses (OM), depends on an equivalent empirical value: the ionization coefficient $\beta$. The existing formulations for $\tau$ and $\beta$ have come from many different sources: theoretical models, ground experiments, \textit{in-situ} experiments, and ground observations. The results of these studies have produced a multitude of empirical values, each in partial disagreement with the others. The discrepancy between different modeling parameters then adversely affects our ability to correctly estimate a meteoroid's mass. This work aims to bridge these discrepancies in two major way: first, by cross-estimating the value of $\tau$ and $\beta$ from concurrent observation via both radar and optical device. Second, we provide a numerical framework to simulate the entire entry process and mass loss for millimeter-sized meteoroids. For the analysis on the cross calibration of $\tau$ and $\beta$ we used the Poker Flat Inchoerent Scatter Radar --- located in Fairbanks, Alaska --- in combination with an Andor camera. We found 28 common events over a period of 10 hours of observation over two nights. Our analysis shows that out of the existing $\beta-\tau$ values, the combination consisting of $\beta$ by \cite{Bronshten1983} with $\tau$ by \cite{Campbell2012} provides the more consistent results based off an error ratio metric. We conclude the study by providing an updated value of $\tau$ based off of our gathered data and the newest experimental results for $\beta$ conducted at the University of Colorado's dust accelerator facility \citep{thomas2016, deluca2017}. This new $\beta-\tau$ provides an error ratio value approximately 10\% lower than the previously computed best $\beta-\tau$ pair. For our numerical results we implemented a Direct Simulation Monte Carlo flow solver loosely coupled with an ablation routine. This coupled system allow us to provide the first full 2D simulation of the atmospheric entry of a sub-centimeter meteoroid into Earth's atmosphere. Our simulations show three important results: first, the existing zero order models make for a reasonable approximation of the mass loss for sub-millimeter objects, albeit the mass deposition profile is different between the two approximations. The small scales of the meteoroid considered, coupled with its high conductivity makes for a quasi-uniform meteoric temperature. Second, the neutral density surrounding the meteoroid significantly increases from it's background value. At the stagnation point it increases between two to three order of magnitude, depending on the initial velocity of the meteoroid, its ablation rate and the original neutral density. This density increase is due to a combination of the ablated mass and the compression of the impinging fluid. Lastly, existing zero order models under-estimate the total heat flux received by the meteoroid. Heat shielding for meteoroids in free molecular flow is negligible even when taking into account increases in neutral densities caused by the evaporated mass. Our simulations show meteoroids receive up to 95\% of the available heat flux from the flow. Shielding effects become significant as meteoroid size increases ($> $1~cm), which is connected with lower altitudes and a different flow regime. In addition to analyzing the meteoroid mass loss process, we investigated the possibility of using meteoroids to probe our atmosphere. Specifically we tried and succeeded in computing the neutral density from meteoroid mass loss measurements. Comparing our results with the popular MSIS-90 model, we obtain averaged values 20-40\% lower than those computed by the semi-empirical global model. The discrepancies found between our methodology and MSIS-90 are consistent with the differences obtained between in-situ measurements from previous experiments and the semi-empirical global model. In addition, we are capable of observing short term changes in density of a magnitude consistent with experimental results. This result suggests that the presented methodology is able to estimate the thermospheric neutral density. Nonetheless, the high error bar associated with our computation, as well as the diversity of the meteoroid flux, presents some limitation to the application of the proposed technique.
Description
| Type of resource | text |
|---|---|
| Form | electronic resource; remote; computer; online resource |
| Extent | 1 online resource. |
| Place | California |
| Place | [Stanford, California] |
| Publisher | [Stanford University] |
| Copyright date | 2018; ©2018 |
| Publication date | 2018; 2018 |
| Issuance | monographic |
| Language | English |
Creators/Contributors
| Author | Limonta, Lorenzo |
|---|---|
| Degree supervisor | Close, Sigrid |
| Thesis advisor | Close, Sigrid |
| Thesis advisor | Alonso, Juan |
| Thesis advisor | Senesky, Debbie |
| Degree committee member | Alonso, Juan |
| Degree committee member | Senesky, Debbie |
| Associated with | Stanford University, Department of Aeronautics and Astronautics. |
Subjects
| Genre | Theses |
|---|---|
| Genre | Text |
Bibliographic information
| Statement of responsibility | Lorenzo Limonta. |
|---|---|
| Note | Submitted to the Department of Aeronautics and Astronautics. |
| Thesis | Thesis Ph.D. Stanford University 2019. |
| Location | electronic resource |
Access conditions
- Copyright
- © 2018 by Lorenzo Limonta
- License
- This work is licensed under a Creative Commons Attribution Non Commercial 3.0 Unported license (CC BY-NC).
Also listed in
Loading usage metrics...