Recent missions to solar system comets, such as ESA's Rosetta mission, raise interest for models and descriptions of their plasma environment. The interaction with various space phenomena such as stellar wind make the construction of an analytical description difficult. Instead, a simplified view of the comet environment is considered where the effects of magnetism and departures from radial symmetry are neglected. This is done in an effort to construct an approximation of the comet plasma behaviour later to be compared against observational accounts to find which plasma features are dependent on more complex phenomena and which plasma features arise as a result of the simpler comet view. Several attempts are made to construct an analytical description of comet plasma as based on the description within another branch of plasma physics: fusion. Previous work regarding the vacuum expansion of plasma after a stationary target is rapidly ablated via high-intensity lasers appears promising for adaptation to the comet environment. Before the comet environment can be considered the different natures of the two problems have to be considered. For example, the comet case is a stationary expansion problem as opposed to fast-ignition fusion where the expansion is treated as an initial value problem. Having accounted for the problems' inherent differences, a few methods are proposed to convert solutions of lab fusion distribution functions to the comet case. Additionally, a numerical approach to calculate the distribution function of comet electrons is presented employing ergodic invariance. Lastly, a toy-model simulation of the timescale for variations in the potential show that the error in the ergodic invariance may in practice have a faster convergent timescale dependence than theoretical bounds suggest. Optimistically, this suggest the possibility of future use in numerical attempts at modelling comet plasma.