Adaptive Time Stepping for Transient Network Flow Simulation in Rocket Propulsion Systems

Fluid and thermal transients found in rocket propulsion systems such as propellant feedline system is a complex process involving fast phases followed by slow phases. Therefore their time accurate computation requires use of short time step initially followed by the use of much larger time step. Yet there are instances that involve fast-slow-fast phases. In this paper, we present a feedback control based adaptive time stepping algorithm, and discuss its use in network flow simulation of fluid and thermal transients. The time step is automatically controlled during the simulation by monitoring changes in certain key variables and by feedback. In order to demonstrate the viability of time adaptivity for engineering problems, we applied it to simulate water hammer and cryogenic chill down in pipelines. Our comparison and validation demonstrate the accuracy and efficiency of this adaptive strategy

A novel finite element technique for the solution of engineering flow problems

A new technique known as the bubble function method is developed for the modelling of fluid flow problems. The main motivation for this work has been the desire to resolve difficulties that traditional methods show in dealing with multi-scale behaviour in flow regimes. All of the traditionally used methods require excessive mesh refinement in the simulations of systems that combine different scale of behaviour in one domain. The present bubble function method avoids such crude remedies and instead of using an elegant mathematical technique for the conjunctive approximation of fine and coarse scale phenomena. Using numerical experiments it is shown that the implementation of the bubble function method generates accurate and stable solutions for a wide range of problems. This range includes convection and reaction dominated transport phenomena, and various types of porous flow systems, it can also be extended to transient flow simulations. To demonstrate the applicability of the present technique it has been used to solve a realistic problem, namely solute dispersion in an estuary. The results of this simulation show good agreement with field survey data

Robust stabilised finite element solvers for generalised Newtonian fluid flows

Various materials and solid-fluid mixtures of engineering and biomedical interest can be modelled as generalised Newtonian fluids, as their apparent viscosity depends locally on the flow field. Despite the particular features of such models, it is common practice to combine them with numerical techniques originally conceived for Newtonian fluids, which can bring several issues such as spurious pressure boundary layers, unsuitable natural boundary conditions and coupling terms spoiling the efficiency of nonlinear solvers and preconditioners. In this work, we present a finite element framework dealing with such issues while maintaining low computational cost and simple implementation. The building blocks of our algorithm are (i) an equal-order stabilisation method preserving consistency even for lowest-order discretisations, (ii) robust extrapolation of velocities in the time-dependent case to decouple the rheological law from the overall system, (iii) adaptive time step selection and (iv) a fast physics-based preconditioned Krylov subspace solver, to tackle the relevant range of discretisation parameters including highly varying viscosity. Selected numerical experiments are provided demonstrating the potential of our approach in terms of robustness, accuracy and efficiency for problems of practical interest

Python framework for HP adaptive discontinuous Galerkin methods for two phase flow in porous media

In this paper we present a framework for solving two-phase flow problems in porous media. The discretization is based on a Discontinuous Galerkin method and includes local grid adaptivity and local choice of polynomial degree. The method is implemented using the new Python frontend Dune-FemPy to the open source framework Dune. The code used for the simulations is made available as Jupyter notebook and can be used through a Docker container. We present a number of time stepping approaches ranging from a classical IMPES method to a fully coupled implicit scheme. The implementation of the discretization is very flexible allowing to test different formulations of the two-phase flow model and adaptation strategies

Adaptive time-stepping for incompressible flow part I: scalar advection-diffusion

Even the simplest advection-diffusion problems can exhibit multiple time scales. This means that robust variable step time integrators are a prerequisite if such problems are to be efficiently solved computationally. The performance of the second order trapezoid rule using an explicit Adams–Bashforth method for error control is assessed in this work. This combination is particularly well suited to long time integration of advection-dominated problems. Herein it is shown that a stabilized implementation of the trapezoid rule leads to a very effective integrator in other situations: specifically diffusion problems with rough initial data; and general advection-diffusion problems with different physical time scales governing the system evolution

Computational fluid dynamics of coupled free/porous regimes: a specialised case of pleated cartridge filter

The multidisciplinary project AEROFIL has been defined and coordinated with the idea of developing novel filter designs to be employed in aeronautic hydraulic systems. The cartridge filters would be constructed using eco-friendly filtration media supported by unconventional disposable or reusable solid components. My main contribution to this project is the development of a robust and cost-effective design and analysis tool for simulating the hydrodynamics in these pleated cartridge filters. The coupled free and porous flow regimes are generally observed in filtration processes. These processes have been the subject of intense investigation for researchers over the decades who are striving hard to resolve some of the critical issues related to the free/porous interfacial constraints and their mathematical representations concerning its industrial applications. [Continues.