MPM study of water wave interaction with porous seawalls

Flood defences are becoming increasingly vital for protecting key coastal infrastructure from rising sea levels and storm surge waves. Waves attenuate rapidly as they propagate through porous media, corresponding to significant energy dissipation. The ability of porous armour layers to absorb wave energy is therefore of great interest to hydraulic engineers, as more and more natural and artificial porous structures are constructed to defend vulnerable coastlines or existing flood defences from wave attack. Design parameters for permeable flood defences include the grain size of the sediment particles that form the barrier and the width of the barrier. Previous studies have demonstrated that the runup height, which determines the risk of overtopping, is primarily dependent on the still water depth, wave amplitude and ground slope. This thesis investigates the effect of manipulating the mean grain size of a permeable barrier on the runup response to a dam-break flood wave and a solitary wave using the Material Point Method (MPM). Traditional methods for ascertaining the stability of protective barriers have used small-scale physical models, however, these are expensive and have been shown to suffer from scaling problems. Numerical methods are therefore gaining popularity in flood risk management applications. MPM is capable of handling large deformation problems within a Lagrangian framework, allowing for simple application of boundary conditions. MPM has been widely used to solve solid mechanics problems with history-dependent variables, but the application of this method to fluid mechanics has been rare. In MPM, the background mesh is only used to solve the governing equations. The material properties are stored at the material points so that issues arising from mesh distortions such as those exhibited in the classical finite element method (FEM) are easily avoided when coping with large deformations. Although it stores no permanent information, the background mesh allows for simple application of boundary conditions. Boundary conditions can be applied directly at the nodes, unlike in meshfree methods such as SPH where boundary conditions must be applied to a series of "boundary particles", which much first be identified. In single-point MPM, both the liquid and solid velocity fields are tracked by a single material point. The double-point MPM introduces two sets of material points representing the solid phase and liquid phase separately, so that it is capable of capturing the relative acceleration between the water and the soil skeleton. It is therefore capable of accurately modelling situations where the fluid moves significantly with respect to the soil, such as in wave run-up on structures. The goal of this research is to determine the effectiveness of sloped and vertical permeable barriers on preventing flooding resulting from storm surges and tsunami waves, and to ascertain the most effective permeability for reducing wave impact, to provide design guidance for coastal flood barriers. To this end, MPM is used to examine the influence of two key design parameters for porous, permeable flood defences, including the grain size of the composition particles and the width of the barrier. The term porosity is used to describe the dimensionless ratio of the volume of voids to solid material in the structure, and the term permeability is used to describe how resistant the porous structure is to flow, related to the mean grain size forming the material, in accordance with the Ergun equation used to determine the body force between the solid and liquid phases. This equation is presented in Chapter 3. A larger mean grain size results in larger voids in the material, so that there is less resistance to flow and the permeability increases. The multiphase version of the MPM package Anura3D (www.anura3d.com) has been adopted in this study. It is shown that increasing the grain size, and therefore the permeability, of the porous dam effectively reduces the overall runup height. The grain size, rather than the wall width, is shown to be the dominant parameter affecting the runup.EPSR

CFD modelling approach for dam break flow studies

Abstract. This paper presents numerical simulations of free surface flows induced by a dam break comparing the shallow water approach to fully three-dimensional simulations. The latter are based on the solution of the complete set of Reynolds-Averaged Navier-Stokes (RANS) equations coupled to the Volume of Fluid (VOF) method. The methods assessment and comparison are carried out on a dam break over a flat bed without friction, a dam break over a triangular bottom sill and a dam break flow over a 90° bend. Experimental and numerical literature data are compared to present results. The results demonstrate that the shallow water approach, even if able to sufficiently reproduce the main aspects of the fluid flows, loses some three-dimensional phenomena, due to the incorrect shallow water idealization that neglects the three-dimensional aspects related to the gravity force

GeomInt–Mechanical Integrity of Host Rocks

This open access book summarizes the results of the collaborative project “GeomInt: Geomechanical integrity of host and barrier rocks - experiment, modeling and analysis of discontinuities” within the Program: Geo Research for Sustainability (GEO: N) of the Federal Ministry of Education and Research (BMBF). The use of geosystems as a source of resources, a storage space, for installing underground municipal or traffic infrastructure has become much more intensive and diverse in recent years. Increasing utilization of the geological environment requires careful analyses of the rock–fluid systems as well as assessments of the feasibility, efficiency and environmental impacts of the technologies under consideration. The establishment of safe, economic and ecological operation of underground geosystems requires a comprehensive understanding of the physical, (geo)chemical and microbiological processes on all relevant time and length scales. This understanding can only be deepened on the basis of intensive laboratory and in-situ experiments in conjunction with reliable studies on the modeling and simulation (numerical experiments) of the corresponding multi-physical/chemical processes. The present work provides a unique handbook for experimentalists, modelers, analysts and even decision makers concerning the characterization of various types of host rocks (salt, clay, crystalline formations) for various geotechnical applications

A particle finite element method for fluid-related problems in civil engineering

The work presented in this Thesis is a set of developments focused on the Particle Finite Element Method (PFEM) and its applicability in several fields in Civil Engineering. The PFEM had already been proven to be a powerful tool for the free surface flows with large deformation and domain separation, but the application to actual engineering problems requires many more advances. The interaction between the fluid and many solids contacting with each other, the erosion of soils and the transport of small particles are some of these advances, which are main topics addressed in this document. Apart from them, other developments related with the fluid solution are included, which are intended to get deeper than ever before in the practical use of PFEM.El treball que es presenta en aquesta Tesi és un conjunt de desenvolupaments centrats en el Particle Finite Element Method (PFEM) i en la seva aplicació a diversos camps de l'enginyeria civil. El PFEM ja havia demostrat ser una eina potent pels fluxes amb superfície lliure amb grans deformacions i separació de dominis, però l'aplicació a problemes d'enginyeria reals requereix molts més avenços. La interacció entre el fluid i molts sòlids que contacten els uns amb els altres, l'erosió de sòls i el transport de partícules petites són alguns d'aquests avenços, que són els principals temes tractats en aquesta Tesi. Apart d'aquests, s'inclouen altres desenvolupaments relacionats amb la solució del fluid, que miren d'arribar més profunditat que mai abans en l'ús pràctic del PFEM i la seva implementació. Primer es presenta el PFEM, es descriuen els desenvolupaments de l'autor millorant la solució de la dinàmica de fluids i altres capacitats simples que s'hi han afegit. Després, tres capítols principals es centren en a) l'algoritme d'interacció fluid-estructura amb contacte b) l'erosió de sòls c) el transport de partícules. A continuació altres aplicacions del mètode s'expliquen, així com una llista dels projectes d'investigació amb els quals aquesta tesi ha tingut vincle.Postprint (published version

Multi-scale modeling of inertial flows through propped fractures

Non-Darcy flows are expected to be ubiquitous in near wellbore regions, completions, and in hydraulic fractures of high productivity gas wells. Further, the prevailing dynamic effective stress in the near wellbore region is expected to be an influencing factor for the completion conductivity and non-Darcy flow behavior in it. In other words, the properties (fracture permeability and β-factor) can vary with the time and location in the reservoir (especially in regions close to the wellbore). Using constant values based on empirical correlations for reservoirs/completions properties can lead to erroneous cumulative productivity predictions. With the recent advances in the imaging technology, it is now possible to reconstruct pore geometries of the proppant packs under different stress conditions. With further advances in powerful computing platforms, it is possible to handle large amount of computations such as Lattice Boltzmann (LB) simulations faster and more efficiently. Calculated properties of the proppant pack at different confining stresses show reasonable agreement with the reported values for both permeability and β-factor. These predicted stress-dependent permeability and β-factors corresponding to the effective stress fields around the hydraulic fractured completions is included in a 2D gas reservoir simulator to calculate the productivity index. In image-based flow simulations, spatial resolution of the digital images used for modeling is critical not only because it dictates the scale of features that can be resolved, but also because for most techniques there is at least some relationship between voxel size in the image data and numerical resolution applied to the computational simulations. In this work we investigate this relationship using a computer-generated consolidated porous medium, which was digitized at voxel resolutions in the range 2-10 microns. These images are then used to compute permeability and tortuosity using lattice Boltzmann (LB) and compared against finite elements methods (FEM)simulation results. Results show how changes in computed permeability are affected by image resolution (which dictates how well the pore geometry is approximated) versus grid or mesh resolution (which changes numerical accuracy). For LB, the image and grid resolution are usually taken to be the same; we show at least one case where effects of grid and image resolution appear to counteract one another, giving the mistaken appearance of resolution-independent results. For FEM, meshing can provide certain attributes (such as better conformance to surfaces), but it also adds an extra step for error or approximation to be introduced in the workflow

Numerical Modeling of River Diversions in the Lower Mississippi River

The presence of man-made levees along the Lower Mississippi River (MR) has significantly reduced the River sediment input to the wetlands and much of the River\u27s sediment is now lost to the Gulf of Mexico. The sediment load in the River has also been decreased by dams and river revetments along the Upper MR. Freshwater and sediment diversions are possible options to help combat land loss. Numerical modeling of hydrodynamics and sediment transport of the MR is a useful tool to evaluate restoration projects and to improve our understanding of the resulting River response. The emphasis of this study is on the fate of sand in the river and the distributaries. A 3-D unsteady flow mobile-bed model (ECOMSED; HydroQual 2002) of the Lower MR reach between Belle Chasse (RM 76) and downstream of Main Pass (RM 3) was calibrated using field sediment data from 2008 – 2010 (Nittrouer et al. 2008; Allison, 2010). The model was used to simulate River currents, diversion sand capture efficiency, erosional and depositional patterns with and without diversions over a short period of time (weeks). The introduction of new diversions at different locations, e.g., Myrtle Grove (RM 59) and Belair (RM 65), with different geometries and with different outflows was studied. A 1-D unsteady flow mobile-bed model (CHARIMA; Holly et al. 1990) was used to model the same Lower MR reach. This model was used for longer term simulations (months). The simulated diversions varied from 28 m3/s (1, 000 cfs) to 5, 700 m3/s (200, 000 cfs) for river flows up to 35, 000 m3/s (1.2x106 cfs). The model showed that the smaller diversions had little impact on the downstream sand transport. However, the larger diversions had the following effects: 1) reduction in the slope of the hydraulic grade line downstream of the diversion; 2) reduction in the available energy for transport of sand along distributary channels; 3) reduced sand transport capacity in the main channel downstream of the diversion; 4) increased shoaling downstream of the diversion; and 5) a tendency for erosion and possible head-cutting upstream of the diversion

Proposal for Numerical Benchmarking of Fluid-Structure Interaction in Cerebral Aneurysms

Computational fluid dynamics is intensively used to deepen the understanding of aneurysm growth and rupture in the attempt to support physicians during therapy planning. Numerous studies have assumed fully-rigid vessel walls in their simulations, whose sole hemodynamics may fail to provide a satisfactory criterion for rupture risk assessment. Moreover, direct in-vivo observations of intracranial aneurysm pulsation have been recently reported, encouraging the development of fluid-structure interaction for their modelling and for new assessments. In this work, we describe a new fluid-structure interaction benchmark setting for the careful evaluation of different aneurysm shapes. The studied configurations consist of three real aneurysm domes positioned on a toroidal channel. All geometric features, meshing characteristics, flow quantities, comparisons with a rigid-wall model and corresponding plots are provided. Reported results emphasize the alteration of flow patterns and hemodynamic descriptors when moving from the rigid-wall model to the complete fluid-structure interaction framework, thereby underlining the importance of the coupling between hemodynamics and the surrounding vessel tissue.Comment: 23 pages, 14 figure