Modelling and simulation of atmosphere-breathing electric propulsion intakes via direct simulation Monte Carlo

The Air-Breathing Ion Engine (ABIE) is an electric propulsion system capable of compensating for drag at low altitudes by ingesting the surrounding atmospheric particles to be utilized as propellant. The current state of the art of the ABIE performance is evaluated via Direct Simulation Monte Carlo (DSMC), due to the rarefied nature of the atmosphere in Very-Low Earth Orbit (VLEO). Nevertheless, the scarce availability of relevant simulation methodologies in the literature limits the repeatability of such numerical studies. Therefore, this paper proposes an independent methodology applicable to the modelling and simulation of Atmosphere-Breathing Electric Propulsion (ABEP) intakes that aims to validate the ABIE DSMC results retrieved from the literature. This is achieved by investigating the ABIE intake collection efficiency and compression ratio through the open-source solver dsmcFoam+ and by assessing the results against the available RARAC-3D DSMC data. First, the variation of grid transparency is discussed and compared between both solvers, yielding a mean percentage error of $$2.97\%$$ 2.97 % for the compression ratio and $$2.06\%$$ 2.06 % for the collection efficiency. Second, the absence of intermolecular collisions is verified for which errors of $$1.61\%$$ 1.61 % for collection efficiency and $$3.49\%$$ 3.49 % for compression ratio are observed. Then, the variation of flow incidence angle is simulated between $$0^{\circ }$$ 0 ∘ and $$15^{\circ }$$ 15 ∘ , yielding differences lower than $$1.80\%$$ 1.80 % . Consecutively, the intake aspect ratio is varied between 10 and 40, for which a maximum discrepancy of $$1.83\%$$ 1.83 % is measured and, finally, the drag coefficient of the intake is obtained in dsmcFoam+ to define the power density requirements.

Numerical simulation of pulsed planar evaporation into background gas based on direct Monte Carlo simulation and solution of the BGK model kinetic equation

A numerical study of the planar gas expansion under pulsed evaporation into the background gas is carried out. The chosen conditions are typical for nanosecond laser deposition of thin films and nanostructure synthesis, with the saturated gas pressure at the surface of 5.4 MPa and the background pressure of 50 and 500 Pa. The problem is solved based on the direct simulation Monte Carlo method and direct numerical solution of the BGK model kinetic equation. A generally good agreement was obtained for all computed macroscopic quantities, with the exception of the higher density peak in the compressed layer and a wider shock front in the background gas for the BGK model.

Improved theory for shock waves using the OBurnett equations

The main goal of the present study is to thoroughly test the recently derived OBurnett equations for the normal shock wave flow problem for a wide range of Mach number ( $3 \leq Ma \leq 9$ ). A dilute gas system composed of hard-sphere molecules is considered and the numerical results of the OBurnett equations are validated against in-house results from the direct simulation Monte Carlo method. The primary focus is to study the orbital structures in the phase space (velocity–temperature plane) and the variation of hydrodynamic fields across the shock. From the orbital structures, we observe that the heteroclinic trajectory exists for the OBurnett equations for all the Mach numbers considered, unlike the conventional Burnett equations. The thermodynamic consistency of the equations is also established by showing positive entropy generation across the shock. Further, the equations give smooth shock structures at all Mach numbers and significantly improve upon the results of the Navier–Stokes equations. With no tweaking of the equations in any way, the present work makes two important contributions by putting forward an improved theory of shock waves and establishing the validity of the OBurnett equations for solving complex flow problems.