Gravitational Collapse and Singularity Removal in Rastall Theory

In the present work, we study spherically symmetric gravitational collapse of a homogeneous fluid in the framework of Rastall gravity. Considering a nonlinear equation of state (EoS) for the fluid profiles, we search for a class of nonsingular collapse solutions and the possibility of singularity removal. We find that depending on the model parameters, the collapse scenario halts at a minimum value of the scale factor at which a bounce occurs. The collapse process then enters an expanding phase in the postbounce regime, and consequently the formation of a spacetime singularity is prevented. We also find that, in comparison to the singular case where the apparent horizon forms to cover the singularity, the formation of apparent horizon can be delayed allowing thus the bounce to be causally connected to the external universe. The nonsingular solutions we obtain satisfy the weak energy condition (WEC) which is crucial for physical validity of the model.

Comparison of occurrence-bias-adjusting methods for hydrological impact modelling

Over the past decade, various methods for bias adjustment of precipitation occurrence or intensity have been proposed. However, the performance of combined methods has not yet been thoroughly evaluated, especially in a hydrological and climate change context. In this study, four occurrence-bias-adjusting methods are combined with one univariate and one multivariate intensity-bias-adjusting method. The occurrence-bias-adjusting methods include thresholding, Stochastic Singularity Removal, Triangular Distribution Adjustment, and are compared with the intensity-bias-adjusting methods without specific adjustment as a baseline. These combined methods are compared with respect to precipitation amount, precipitation occurrence and discharge. This comparison, summarized in terms of the residual bias relative to both the observations and the model bias,shows significant differences in performance. Occurrence-bias-adjusting methods that add stochasticity perform worse, an effect that is reinforced by multivariate intensity-bias-adjusting methods. The use of simpler methods is thus advised until the uncertainty caused by combining methods is better understood.

A singularity removal method for coupled 1D–3D flow models

In reservoir simulations, the radius of a well is inevitably going to be small compared to the horizontal length scale of the reservoir. For this reason, wells are typically modelled as lower-dimensional sources. In this work, we consider a coupled 1D–3D flow model, in which the well is modelled as a line source in the reservoir domain and endowed with its own 1D flow equation. The flow between well and reservoir can then be modelled in a fully coupled manner by applying a linear filtration law. The line source induces a logarithmic-type singularity in the reservoir pressure that is difficult to resolve numerically. We present here a singularity removal method for the model equations, resulting in a reformulated coupled 1D–3D flow model in which all variables are smooth. The singularity removal is based on a solution splitting of the reservoir pressure, where it is decomposed into two terms: an explicitly given, lower-regularity term capturing the solution singularity and some smooth background pressure. The singularities can then be removed from the system by subtracting them from the governing equations. Finally, the coupled 1D–3D flow equations can be reformulated so they are given in terms of the well pressure and the background reservoir pressure. As these variables are both smooth (i.e. non-singular), the reformulated model has the advantage that it can be approximated using any standard numerical method. The reformulation itself resembles a Peaceman well correction performed at the continuous level.

The vortex-entrainment sheet in an inviscid fluid: theory and separation at a sharp edge

In this paper a model for viscous boundary and shear layers in three dimensions is introduced and termed a vortex-entrainment sheet. The vorticity in the layer is accounted for by a conventional vortex sheet. The mass and momentum in the layer are represented by a two-dimensional surface having its own internal tangential flow. Namely, the sheet has a mass density per-unit-area making it dynamically distinct from the surrounding outer fluid and allowing the sheet to support a pressure jump. The mechanism of entrainment is represented by a discontinuity in the normal component of the velocity across the sheet. The velocity field induced by the vortex-entrainment sheet is given by a generalized Birkhoff–Rott equation with a complex sheet strength. The model was applied to the case of separation at a sharp edge. No supplementary Kutta condition in the form of a singularity removal is required as the flow remains bounded through an appropriate balance of normal momentum with the pressure jump across the sheet. A pressure jump at the edge results in the generation of new vorticity. The shedding angle is dictated by the normal impulse of the intrinsic flow inside the bound sheets as they merge to form the free sheet. When there is zero entrainment everywhere the model reduces to the conventional vortex sheet with no mass. Consequently, the pressure jump must be zero and the shedding angle must be tangential so that the sheet simply convects off the wedge face. Lastly, the vortex-entrainment sheet model is demonstrated on several example problems.