Nitrous oxide is often used in the space industry, as an oxidiser or monopropellant, mostly in self-pressurised configurations. It has potential for growth in use due to the recent rising interest in green propellants. At the same time, modelling the behaviour of a self-pressurising nitrous oxide tank is a challenging task, and few accurate numerical models are currently available. Two-phase flow, heat transfer and rapid changes of mass and temperature in the investigated system all increase the difficulty of accurately predicting this process. To get a get better understanding of the emptying of a self-pressurised nitrous oxide tank, two models were developed: a phase equilibrium model (single node equilibrium), treating the control volume as a single node in equilibrium state, and a phase interface model, featuring a moving interface between parts of the investigated medium. The single node equilibrium model is a variation of equilibrium model previously described in the literature, while the phase interface model involves a novel approach. The results show that the models are able to capture general trends in the main parameters, such as pressure or temperature. The phase interface model predicts nitrous oxide as a liquid, a two-phase mixture, and vapour in the lower part of the tank, which is reflected in the dynamics of changes in pressure and mass flow rate. The models developed for self-pressurisation, while created for predicting nitrous oxide behaviour, could be adapted for other media in conditions near vapour– liquid equilibrium by adding appropriate state equations.

Shell and tube heat exchangers are commonly used in a wide range of practical engineering. The key issue in such a system is the heat exchange between the hot and cold working media. An increased cost of production of these devices has forced all manufacturing companies to reduce the total amount of used materials by better optimizing their construction. Numerous studies on the heat exchanger design codes have been carried out, basically focusing on the use of fully time-dependent partial differential equations for mass, momentum, and energy balance. They are very complex and time-consuming, especially when the designers want to have full information in a full 3D system. The paper presents the 1D mathematical model for analysis of the thermal performance of the counter-current heat exchanger comprised of mixed time-dependent and time-independent equations, solved by the upwind numerical solution method, which allows for a reduction in the CPU time for obtaining the proper solution. The comparison of numerical results obtained from an in-house program called Upwind Heat Exchanger Solver written in a Fortran code, with those derived using commercial software package ASPEN, and those obtained experimentally, shows very good agreement in terms of the temperature and pressure distribution predictions. The proposed method for fast designing calculations appears beneficial for other tube shapes and types of heat exchangers.