The possibility of a subdominant component of dark matter dissipating energy could lead to dramatic new phenomenology such as the formation of a dark disk. One rigorous way to assess this possibility and settle the debate on its feasibility is to include the dissipative dark component in a numerical hydrodynamical simulation. A necessary input to such a simulation is a prescription including energy dissipation rates of different processes and rates of processes that change the number densities of dark ions and atoms. In this article, we study the simplest dissipative dark sector which consists of a dark electron and proton, both charged under a dark gauged U(1). We present approximate analytic formulas for energy loss rates due to Compton scattering, bremsstrahlung, recombination, collisional ionization and collisional excitation as well as the rates of number density change. We also include the heating rate due to photoionization. The work serves as the first step to realize a numerical simulation including a dissipative dark sector, which hopefully can shed more light on the formation and properties of a dark disk originating from dark matter self-interactions.

The composition of Jupiter and the primordial distribution of the heavy elements are determined by its formation history. As a result, in order to constrain the primordial internal structure of Jupiter the growth of the core and the deposition and settling of accreted planetesimals must be followed in detail. In this paper we determine the distribution of the heavy elements in proto-Jupiter and determine the mass and composition of the core. We find that while the outer envelope of proto-Jupiter is typically convective and has an homogeneous composition, the innermost regions have compositional gradients. In addition, the existence of heavy elements in the envelope leads to much higher internal temperatures (several times 10$^4$ K) than in the case of a hydrogen-helium envelope. The derived core mass depends on the actual definition of the core: if the core is defined as the region in which the heavy-element mass fraction is above some limit (say 0.5), then it can be much more massive ($\sim$ 15 M$_{\oplus}$) and more extended (10\% of the planet's radius) than in the case where the core is just the region with 100\% heavy elements. In the former case Jupiter's core also consists of hydrogen and helium. Our results should be taken into account when constructing internal structure models of Jupiter and when interpreting the upcoming data from the {\it Juno} (NASA) mission.

A 3-D numerical model of comet nuclei is presented. An implicit numerical scheme was developed for the thermal evolution of a spherical nucleus composed of a mixture of ice and dust. The model was tested against analytical solutions, simplified numerical solutions, and 1-D thermal evolution codes. The 3-D code was applied to comet 67P/Churyumov-Gerasimenko; surface temperature maps and the internal thermal structure was obtained as function of depth, longitude and hour angle. The effect of the spin axis tilt on the surface temperature distribution was studied in detail. It was found that for small tilt angles, relatively low temperatures may prevail on near-pole areas, despite lateral heat conduction. A high-resolution run for a comet model of 67P/Churyumov-Gerasimenko with low tilt angle, allowing for crystallization of amorphous ice, showed that the amorphous/crystalline ice boundary varies significantly with depth as a function of cometary latitude.