Wintner's collinear and flat solutions are nowhere dense

It is shown that in then-body problem with generalized attraction law, the sets of initial conditions which lead to Wintner's collinear, respectively flat, motion are nowhere dense relative to the set of initial conditions that define solutions inR ³.     

On the planar syzygy solutions of the 3-body problem

Relations between the rectilinear, collinear and syzygy solutions of the N-body problem are first pointed out. It is shown that, along a solution, the set of the non-collinear syzygy configuration instants is formed by isolated points. Then we restrict the study to the planar 3-body problem and prove that for Dirichlet-stable solutions, a non-syzygy solution cannot be as close as possible to a syzygy one. It is also true that, in the case of a syzygy solution, the orbit of one particle crosses the line of the other two and can not be tangent to this line in the transition point. Finally we prove that the set of initial conditions leading to non-collinear syzygy solutions is non-empty and open.     

The effect of dispersion on the stability of density-driven flows in saturated homogeneous porous media

This work is a continuation of Musuuza et al. [37] in which a stability criterion for density-driven flow in a saturated homogeneous medium was derived. The criterion predicted the stability of a system as a function of the density and viscosity contrasts, the magnitude of the flow velocity and the concentration gradients for flow aligned orthogonal to gravity. It could not accurately predict stability transition with increasing velocity, a failure we attributed to dispersion effects that were not included. Small-scale dispersion and molecular diffusion, the main stabilising mechanisms in homogeneous media can stabilise flow if the instability wavelengths are below a certain cutoff. The width of the mixing zone is also central in controlling the range of wavelengths that persist into fingers. We propose a method of quantifying the cutoff wavelength and the width of the mixing zone, which are incorporated into the earlier criterion as constituents of the dispersive part. The earlier criterion is reformulated in terms of the Rayleigh number and with the dispersive part added, we attempt to predict the number of fingers formed which is directly related to the physical stability of the system. The inclusion of the Rayleigh number and dispersion into a single stability criterion provides new insight in the way dispersion affects vertical flow systems. Stability numbers computed with the new criterion are in reasonable agreement with numerical simulations for a range of physical variables. The numerical computations are performed with the software package d³f, which uses the cell-centred finite volume and the implicit Euler methods for the spatial and temporal discretisations, respectively. The admission of the density and dispersivities as inputs into the criterion makes it usable in practical problems.     

The stability of density-driven flows in saturated heterogeneous porous media

This work concludes the investigations into the stability of haline flows in saturated porous media. In the first part [33] a stability criterion for density-driven flow in a saturated homogeneous medium was derived excluding dispersion. In the second part [34], the effects of dispersion were included. The latter criterion made reasonable predictions of the stability regimes (indicated by the number of fingers present) as a function of density and dispersivity variations. We found out that destabilising variables caused an increase in the number of fingers and vice versa. The investigation is extended here for the effects of the medium heterogeneity. The cell problem derived via homogenization theory [20] is solved and its solution used to evaluate the elements of the macrodispersion tensor as functions of time for flow aligned parallel to gravity. The longitudinal coefficient exhibits asymptotic behaviour for favourable and moderately unfavourable density contrasts while it grows indefinitely for higher density contrasts. The range of densities stabilised by medium heterogeneities can thus be estimated from the behaviour of the coefficient. The d³f software program is used for the numerical simulations. The code uses the cell-centred finite volume and the implicit Euler techniques for the spatial and temporal discretisations respectively.     

Mixed finite element discretization and Newton iteration for a reactive contaminant transport model with nonequilibrium sorption: convergence analysis and error estimates

We present a numerical scheme for reactive contaminant transport with nonequilibrium sorption in porous media. The mass conservative scheme is based on Euler implicit, mixed finite elements, and Newton method. We consider the case of a Freundlich-type sorption. In this case, the sorption isotherm is not Lipschitz but just Hölder continuous. To deal with this, we perform a regularization step. The convergence of the scheme is analyzed. An explicit order of convergence depending only on the regularization parameter, the time step, and the mesh size is derived. We give also a sufficient condition for the quadratic convergence of the Newton method. Finally, relevant numerical results are presented.     

Assessment of seismic risk scenarios for Bucharest,Romania

In this paper, seismic risk scenarios for Bucharest, the capital city of Romania, are proposed and assessed. Bucharest has one of the highest seismic risk levels in Europe, and this is due to a combination of relatively high seismic hazard and a building stock built mainly before the devastating Vrancea 1977 earthquake. In this study, the seismic risk of Bucharest is assessed using the most recent information regarding the characteristics of the residential building stock. The ground motion amplitudes are evaluated starting from random fields obtained by coupling a ground motion model derived for the Vrancea intermediate-depth seismic source with a spatial correlation model. The seismic risk evaluation method applied in this study is based on the well-known macroseismic method. For several structural typologies, the vulnerability parameters are evaluated based on a damage survey performed on 18,000 buildings in Bucharest after the March 1977 earthquake. Subsequently, the risk metrics are compared with those from other studies in the literature that apply a different risk assessment methodology in order to gain a better view of the uncertainties associated with a seismic risk study at city level. Finally, the impact of several Vrancea intermediate-depth earthquake scenarios is evaluated and the results show that the earthquake which has the closest epicenter to Bucharest appears to be the most damaging.     

Investigation of the variability of strong ground motions from Vrancea earthquakes

A systematic investigation of the applicability of several ground motion prediction models for Vrancea intermediate-depth seismic source is conducted in this research. Two ground motion prediction models recommended by previous evaluations (Vacareanu et al. in Bull Earthq Eng 11(6):1867–1884, 2013a; Pavel et al. in Earthq Struct 6(1):1–18, 2014), as well as two new state-of-the-art ground motion prediction equations (Vacareanu et al. in J Earthq Eng, 2013b; Earthq Struct 6(2):141–161, 2014) are tested using an increased strong ground motion database consisting of 150 recordings from Vrancea subcrustal earthquakes. The evaluation is performed by using several goodness-of-fit parameters from the literature. Moreover, the applicability of the single-station sigma method is also investigated by using the same strong ground motion database recorded in 30 seismic stations from southern and eastern Romania. The influence of the soil conditions on the numerical results obtained in this study is investigated and discussed using the results provided by the analysis of variance method. The impact of the single-station standard deviation on the levels of seismic hazard is also assessed in this study, and the results show, in the analyzed cases, significant reductions of the hazard levels.     

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.     

The formation and closure of the Big Breach of Sacalin spit associated with extreme shoreline retreat and shoreface erosion

We document a case of exceptionally large natural breaching of a sandy spit (Sacalin spit, Danube Delta) using multiannual to seasonal surveys of topography and bathymetry on successive cross-barrier and shoreface profiles, LiDAR data, satellite imagery, and wind and wave data. The large breach, which quickly reached 3.4 km in May 2014, is attributed to morphological preconditioning of the narrow (50–150 m) barrier, which was susceptible to breaching even during moderate storm conditions. The event switched the barrier's decadal evolution from low cross-shore transport to high cross-shore transport over the barrier, which is an order of magnitude larger than during the non-breach period. Upper shoreface erosion, as indicated by the extensive erosion down to −4 m, indicates that this zone is a significant source for the rebuilding of the barrier. Barrier recovery and widening trigger a negative feedback which limits the back-barrier sediment transfer. As a result, back-barrier deposition decreases whilst the barrier aggradation through overwash becomes more frequent. The Big Breach (TBB) closed naturally in three years. The very high deposition rate of sediment in the breach is a testimony of the high sediment volumes supplied by the longshore transport and the high sediment released through shoreface retreat, and resulted in widening the barrier to a maximum of 1 km. Since the newly-formed barrier shoreline retreated 500 m, this reveals that barrier breaching is an important mechanism which significantly accelerates the landward migration of the barrier system and is a proof of the highly non-linear morphodynamics involved in the barrier island translation. © 2019 John Wiley & Sons, Ltd.