Documents published in 2020 : Scipedia

Experimental study of transition and skimming flows on stepped spillways in rcc dams: qualitative analysis and pressure measurements
J. Dolz, M. Juny

Journal of Hydraulic Research (2005). Vol. 43 (5), pp. 540-548

Abstract
This paper is a presentation of an experimental facility used for analyzing the flow over the stepped spillway of an RCC dam. The technique of digitizing video sequences is applied in order to obtain a qualitative analysis of the hydraulic behavior of these chutes. The onset of skimming flow is studied using two different criteria: a visual criterion and another that makes use of the pressure pattern on the horizontal faces of the steps. The evolution along the length of the spillway of pressures measured on the centers of symmetry of the horizontal faces of the steps is also analyzed to obtain further conclusions regarding hydraulic behavior during a skimming flow regime.
Abstract
This paper is a presentation of an experimental facility used for analyzing the flow over the stepped spillway of an RCC dam. The technique of digitizing video sequences is [...]
- 18
- read
Experimental study of the dividing flow in steep street crossings
L. Nanía, M. Gómez, J. Dolz

Journal of Hydraulic Research (2004). Vol. 42 (4), pp. 406-412

Abstract
A first approach to the experimental study of the dividing flow in steep street crossings is presented. A scale model of a street crossing has been built and a series of tests has been conducted for studying the flow pattern in the crossing when the flow in the streets is supercritical. It has been found that the inflow power ratio at the entrance of the crossing can be used as a non-dimensional parameter for predicting the flow distribution as well as the flow pattern. On the basis of these relationships a one-dimensional formulation can be developed in order to predict the dividing flow in crossings.
Abstract
A first approach to the experimental study of the dividing flow in steep street crossings is presented. A scale model of a street crossing has been built and a series of tests [...]
- 3
- read
Incidencia de los embalses en el régimen térmico del río. El caso del tramo final del río Ebro
R. Val, D. Niñerola, J. Armengol, J. Dolz

Limnética (2003). Vol. 23 (1), pp. 85-92

Abstract
El agua del río mantiene un intercambio de calor con la atmósfera y con el fondo; este proceso puede lograr su equilibrio a lo largo de un tramo de río con características geológicas y climáticas similares. Su comportamiento puede ser modificado por efectos artificiales; tales como el uso del agua del río para refrigeración de centrales termoeléctricas o por la existencia de embalses. En el caso de regiones con clima mediterráneo, donde se encuentran los embalses aquí estudiados, los efectos de las alteraciones en el régimen térmico producido por éstos, provocan una tendencia a la constancia térmica estacional, un aumento de la uniformidad diaria, una elevación de las temperaturas invernales y la disminución de las temperaturas estivales. La temperatura del agua es un factor importante en la ecología del río, ya que condiciona en gran medida la vida de las comunidades fluviales. Aguas abajo de las termoeléctricas o presas, la temperatura del agua evoluciona de tal forma que tiende a lograr el equilibrio con el medio ambiente.
Abstract
El agua del río mantiene un intercambio de calor con la atmósfera y con el fondo; este proceso puede lograr su equilibrio a lo largo de un tramo de río [...]
- 3
- read
Hydraulic behaviour of the diversion conduit gates at Caruachi
M. Juny, D. Ninyerola, J. Pomares, F. Abadía, E. Egusquiza, J. Dolz, F. Quintero

International journal on hydropower and dams (2002). Vol. 9 (6), pp. 37-41

Abstract
Discharge and aeration capacity of the diversion conduits, the pressures in the ducts and vibration and downpull forces on the gates were analyzed. Model studies were carried out to study the hydraulic behavior of the gates used for closing operation. The model was built to a scale of 1:20 and was based on Froude's similarity criterion. The tests consisted of picking up the vibratory vertical response to impulsive stimulations, using accelerometers. Results indicated that on the current model, no interference was caused by the uncertainties related to the modelled areas up- and downstream
Abstract
Discharge and aeration capacity of the diversion conduits, the pressures in the ducts and vibration and downpull forces on the gates were analyzed. Model studies were carried [...]
- 7
- read
Presa de Caruachi. Estudio del comportmiento hidráulico de las compuertas de los conductos de desvío
M. Sánchez-Juny, D. Niñerola, J. Pomares, E. Egusquiza, J. Dolz, F. Quintero, F. Abadía

Revista de Obras Públicas (2002). N. 3423, pp. 12-26

Abstract
Se presenta el estudio desarrollado en el Laboratorio de Modelos Reducidos del Dpto. de Ing. Hidráulica, Marítima y Ambiental de la UPC, para el análisis del comportamiento hidráulico de las compuertas a utilizar en la operación de cierre de la presa de Caruachi en el río Caroní (Venezuela). Se han analizado la capacidad de desagüe, la aireación de los conductos de desvío, las presiones en los mismos durante la operación, así como las vibraciones y fuerza vertical descendente (down-pull) sobre las compuertas. The present article describes the study carried out at the Reduced Model Laboratory (Hydraulic, Maritime and Environmental Engineering Department, UPC) into the hydraulic behaviour of the gates employed in the closing operation of the Caruachi dam on the River Caroní (Venezuela). An analysis is made of the discharge and aeration capacity of the diversión conduits, the operation pressures in the same and the vibration and down-pull on the gates.
Abstract
Se presenta el estudio desarrollado en el Laboratorio de Modelos Reducidos del Dpto. de Ing. Hidráulica, Marítima y Ambiental de la UPC, para el análisis [...]
- 2
- read
El episodio pluviométrico del 10 de junio en Cataluña: un primer estudio hidrometeorológico
D. Sempere-Torres, C. Corral, M. Berenguer, R. Sánchez-Diezma, J. Dolz

Ingeniería del agua (2000). Vol. 7 (2), pp. 105-115

Abstract
En este artículo se estudia el episodio pluviométrico registrado el día 10 de junio de 2000 en Cataluña. En particular se analiza la lluvia caída en la cuenca de 97 km2 correspondiente a la riera de Magarola en su cruce con la autovía N-II en Esparreguera (Barcelona). Para ello se utiliza información de campo obtenida a partir de radar meteorológico, pluviómetros y pluviógrafos. Sobre la mayor parte de la cuenca se recogieron más de 150 mm de lluvia durante el episodio. También prácticamente sobre toda la cuenca cayeron más de 100 mm en dos horas. El cálculo del período de retorno asociado a la lluvia de 24 horas evidencia la naturaleza excepcional del fenómeno. No obstante la información histórica disponible es insuficiente y solo permite realizar estimaciones con una gran incertidumbre.
Abstract
En este artículo se estudia el episodio pluviométrico registrado el día 10 de junio de 2000 en Cataluña. En particular se analiza la lluvia caída [...]
- 8
- read
Las aguas superficiales y la marisma del Parque Nacional de Doñana
B. Bayán, J. Dolz

Revista de obras públicas (1995). Vol. 142 (3340), pp. 17-29

Abstract
La Marisma del Parque Nacional de Doñana ocupa aproximadamente 27.000 Ha. y constituye uno de los tres grandes ecosistemas existentes en el parque. En el presente siglo se llevaron a cabo actuaciones que redujeron notablemente su extensión y alteraron sustancialmente sus características hidrológicas. A lo largo de las últimas décadas se han tomado decisiones y realizado obras tendentes a la recuperación de la funcionalidad hidrológica original de la marisma y de este modo poder asegurar su conservación. En el presente artículo se analiza el origen, evolución histórica y estado actual de la hidrología de la marisma. También se describen las obras realizadas para su regeneración hídrica y finalmente se hacen unas consideraciones respecto la necesidad de mejorar el conocimiento cuantitativo del funcionamiento hidráulico e hidrológico de este espacio de características que, sin exageración alguna, pueden calificarse de únicas. The marsh area of Doñana, one of the three large eco-systems of the Park, covers some 65.000 acres. It was considerably reduced in size during this century and its hydrological capacity gravely affected. Work has been carried in recent years to restore the original hydrological function of the marshland and to ensure its survival. The origin, the history and the present state of the marshland are examined in this paper, and an account is given of the work of regeneration. Attention is drawn to the need for more information about the quantitative hydraulic and hydrological features of this unique area.
Abstract
La Marisma del Parque Nacional de Doñana ocupa aproximadamente 27.000 Ha. y constituye uno de los tres grandes ecosistemas existentes en el parque. En el presente siglo [...]
- 3
- read
Consecuencias hidrologicas de los incendios forestales
J. Dolz, D. Sempere-Torres

Ingeniería del agua (1994). Vol. 1 (4), pp. 33-48

Abstract
La alteración del terreno y la destrucción de la vegetación debida a los incendios provoca cambios importantes en las características hidrológicas de las cuencas afectadas. El presente estudio aborda una cuantificación de los efectos inducidos por el fuego en los caudales punta y en la modificación de la respuesta hidrológica de una pequeña cuenca mediterránea francesa gracias a la disponibilidad de una serie de 23 años de datos antes del incendio y de 3 años después. La relevancia de las consecuencias puestas de manifiesto resulta especialmente crítica en la región mediterránea, donde una climatología propicia convierte a la secuencia fuego-inundaciones-erosión en el principal agente del proceso de desertificación.
Abstract
La alteración del terreno y la destrucción de la vegetación debida a los incendios provoca cambios importantes en las características hidrológicas [...]
- 12
- read
Problematica del drenaje de aguas pluviales en zonas urbanas y del estudio hidraulico de las redes de colectores
J. Dolz, M. Gómez

Ingeniería del agua (1994). Vol. 1 (1), pp. 55-66

Abstract
Se analiza la problemática del drenaje de aguas pluviales en zonas urbanas, en particular en áreas de rápido y reciente desarrollo urbano como es el caso del litoral mediterráneo español. Se estudia la repercusión que tiene sobre el drenaje un proceso urbanizador no respetuoso con la hidrología de las cuencas naturales preexistentes. Asimismo se analiza de una forma conceptual la problemática que presenta la modelación numérica de los diferentes procesos involucrados en el drenaje urbano, especialmente el comportamiento hidráulico de las redes de colectores.
Abstract
Se analiza la problemática del drenaje de aguas pluviales en zonas urbanas, en particular en áreas de rápido y reciente desarrollo urbano como es el caso [...]
- 25
- read
Erosion fluvial y proteccion en pilas de puente: modelo reducido del puente nuevo de mora del Ebro
J. Dolz, J. Martí-Vide

Revista de Obras Públicas (1994). Vol. 141 (3328), pp. 31-40

Abstract
Este trabajo se propone llamar la atención sobre la acción del agua en las pilas de los puentes, origen del fenómeno de la socavación local, responsable de un gran número de fallos de puentes. En primer lugar, se presenta un breve resumen del conocimiento actual sobre la erosión local y las medidas de protección con mantos de escollera. El núcleo del artículo es el caso del puente nuevo de Mora sobre el río Ebro, estudiado en modelo reducido debido a la complejidad geométrica de sus pilas y estribos, con el objetivo de dimensionar unas protecciones de escollera. Se trata de un modelo a escala 65, localmente erosionable, ensayado con aguas claras y con semejanza de Froude. Finalmente se comparan los resultados del modelo con las predicciones sobre erosión local y tamaño de escollera por algunas de las fórmulas más conocidas, quedando de manifiesto el carácter aproximado de estas fórmulas.
Abstract
Este trabajo se propone llamar la atención sobre la acción del agua en las pilas de los puentes, origen del fenómeno de la socavación local, responsable [...]
- 12
- read
Kinematics of the moving hydraulic jump
J. Martín-Vide, J. Dolz, J. Estal

Journal of Hydraulic Research (1993). Vol. 31 (2), pp. 171-186

Abstract
The moving hydraulic jump, considered as an example of rapidly varied unsteady flow, is theoretically revised. Experimental work is conducted to study the kinematic characteristics of the jump propagation. The position, velocity and conjugate depths of the jump during its propagation either upstream or downstream are investigated and some semi-empirical equations are derived. From this work, some differences are observed between the moving hydraulic jump propagating upstream and downstream. A phenomenon called extinction of the moving hydraulic jump is reported in the latter case. This phenomenon is explained on the grounds of previous theoretical analysis.
Abstract
The moving hydraulic jump, considered as an example of rapidly varied unsteady flow, is theoretically revised. Experimental work is conducted to study the kinematic characteristics [...]
- 8
- read
R-adaptivity in limit analysis
J. Muñoz, J. Hambleton, S. Sloan

(2017). "Advances in direct methods for materials and structures". Berlín: Springer, p. 73-84.

Abstract
Direct methods aim to find the maximum load factor that a domain made of a plastic material can sustain before undergoing full collapse. Its analytical solution may be posed as a constrained maximisation problem, which is computationally solved by resorting to appropriate discretisation of the relevant fields such as the stress or velocity fields. The actual discrete solution is though strongly dependent on such discretisation, which is defined by a set of nodes, elements, and the type of interpolation. We here resort to an adaptive strategy that aims to perturb the positions of the nodes in order to improve the solution of the discrete maximisation problem. When the positions of the nodes are taken into account, the optimisation problem becomes highly non-linear. We approximate this problem as two staggered linear problems, one written in terms of the stress variable (lower bound problem) or velocity variables (upper bound problem), and another with respect to the nodal positions. In this manner, we show that for some simple problems, the computed load factor may be further improved while keeping a constant number of elements.
Abstract
Direct methods aim to find the maximum load factor that a domain made of a plastic material can sustain before undergoing full collapse. Its analytical solution [...]
- 2
- read
Computation of bounds for anchor problems in limit analysis and decomposition techniques
J. Muñoz, R. Syednima, A. Lyamin, A. Huerta

(2013). Direct methods for limit states in structures and materials". Springer, p. 79-100.

Abstract
Numerical techniques for the computation of strict bounds in limit analyses have been developed for more than thirty years. The efficiency of these techniques have been substantially improved in the last ten years, and have been successfully applied to academic problems, foundations and excavations. We here extend the theoretical background to problems with anchors, interface conditions, and joints. Those extensions are relevant for the analysis of retaining and anchored walls, which we study in this work. The analysis of three-dimensional domains remains as yet very scarce. From the computational standpoint, the memory requirements and CPU time are exceedingly prohibitive when mesh adaptivity is employed. For this reason, we also present here the application of decomposition techniques to the optimisation problem of limit analysis. We discuss the performance of different methodologies adopted in the literature for general optimisation problems, such as primal and dual decomposition, and suggest some strategies that are suitable for the parallelisation of large three-dimensional problems. The propo sed decomposition techniques are tested against representative problems.
Abstract
Numerical techniques for the computation of strict bounds in limit analyses have been developed for more than thirty years. The efficiency of these techniques have been [...]
- 1
- read
Stress dependent morphogenesis continuum mechanics and system of trusses
J. Muñoz, V. Conte, M. Miodownik

(2010). "Continuum mechanics". Nova Publishers, p. 223-243.

Abstract
The mechanical analysis of soft tissues in biomechanics has undergone an increasing progress during the last decade. Part of this success is due to the development and application of some techniques of continuum mechanics, in particular, the decomposition of the deformation gradient, and the introduction of mass, density or volume changes in the reference configuration. Resorting to the common terminology employed in the literature, the changes in biomechanical processes may be classified as growth (change of mass), remodelling (change of density or other material properties such as fibre orientation) or morphogenesis (change of shape). Although the use of those concepts in bone and cardiovascular analysis is well extended, their use in morphogenesis during embryo development has been far less studied. The reasons of this fact may be found in the large shape changes encountered during this process, or the complexity of the material changes involved. In this chapter we develop a general framework for the modelling of morphogenesis by introducing a growth process in the structural elements of the cell, which in turn depends on the stress state of the tissue. Some experimental observations suggest this feedback mechanism during embryo development, and only very recently this behaviour has started to be simulated. We here derive the necessary equilibrium equations of a stress controlled growth mechanism in the context of continuum mechanics. In these derivations we assume a free energy source which is responsible of the active forces during the elongation process, and a passive hyperelastic response of the material. In addition, we write the necessary conditions that the active elongation law must satisfy in order to be thermodynamically consistent. We particularise these equations and conditions for the relevant elements of the cytoskeleton, namely, microfilaments and microtubules. We apply themodel to simulate the shape changes observed during embryomorphogenesis in truss element. As a salient result, themodel reveals that by imposing boundary stress conditions, unbounded elongation would be obtained. Therefore, either prescribed displacements or cross-links between fibres are necessary to reach a homeostatic state.
Abstract
The mechanical analysis of soft tissues in biomechanics has undergone an increasing progress during the last decade. Part of this success is due to the development and application [...]
- 0
- read
Planar differential growth rates initiate precise fold positions in complex epithelia
M. Tozlouglu, M. Duda, N. Kirkland, R. Barrientos, J. Burden, J. Muñoz, Y. Mao

Developmental cell (2019). Vol. 51 (3), pp. 299-312

Abstract
Tissue folding is a fundamental process that shapes epithelia into complex 3D organs. The initial positioning of folds is the foundation for the emergence of correct tissue morphology. Mechanisms forming individual folds have been studied, but the precise positioning of folds in complex, multi-folded epithelia is less well-understood. We present a computational model of morphogenesis, encompassing local differential growth and tissue mechanics, to investigate tissue fold positioning. We use the Drosophila wing disc as our model system and show that there is spatial-temporal heterogeneity in its planar growth rates. This differential growth, especially at the early stages of development, is the main driver for fold positioning. Increased apical layer stiffness and confinement by the basement membrane drive fold formation but influence positioning to a lesser degree. The model successfully predicts the in vivo morphology of overgrowth clones and wingless mutants via perturbations solely on planar differential growth in silico.
Abstract
Tissue folding is a fundamental process that shapes epithelia into complex 3D organs. The initial positioning of folds is the foundation for the emergence of correct tissue [...]
- 24
- read
Computational modeling of epithelial wound healing: Short and long term chemo-mechanical mechanisms
L. Roldán, J. Muñoz, P. Sáez

Comput. Methods Appl. Mech. Engrg., (2019). Vol. 350, pp. 28-56

Abstract
During the lifetime of all living multicellular organisms, wounds in their tissues are frequently observed. The capability of closing those gaps is fundamental for a healthy development. If done deficiently, many diseases may occur from simple inflammation to tumor formation. The wound healing process in epithelial tissue occurs in three different stages. The first one is the assembly of a supra-cellular actomyosin cable and its migration towards the wound edge, triggered by biochemical processes in which calcium plays a distinctive role. How this process is orchestrated following damage remains unclear. Later, after its positioning, the cable contracts driving the tissue towards the gap and reducing the wound area. Finally, cell migration towards the interior of the wound ends up sealing the tissue. In this work, we make use of a mechanical continuum model for the first two stages in order to developed and 2D finite element simulations within a monolithically fully implicit implementation. The model for the actomyosin cable formation involves the coupling of transient calcium ions transport, with actin fibers and myosin motors recruitment and non-linear mechanical response of the tissue. The contraction stage, the active deformation of the previously formed actomyosin cable is taken into account. The relative motion of the myosin motors over the actin filaments is modeled so there exists an active tissue contraction in the direction of those fibers. Upon implementation, the model is capable of performing a wide range of biophysical situations reported experimentally, as we demonstrate in our numerical results. We have been able to rationalize through computational mechanics the firing of calcium in the wound right after damage infliction as well as the consequent formation of actin ring, reproducing nicely what has been reported in biological literature. Thereafter, the numerical model of acto-myosin contraction, fully integrated with the non-linear mechanics of the problem, correlates with the mechanics of wound closure at the actin-ring contraction stage. More importantly, the approach is the first of its kind in the modeling of epithelial and embryonic cell layers, where a wide number of complex mechanics has been integrated and solved though computational methods in engineering. We believe that the simulations will help to unravel new insights in open questions of developmental biology.
Abstract
During the lifetime of all living multicellular organisms, wounds in their tissues are frequently observed. The capability of closing those gaps is fundamental [...]
- 13
- read
Stress relaxation in epithelial monolayers is controlled by the actomyosin cortex
N. Khalilgharibi, J. Fouchard, N. Asadipour, R. Barrientos, M. Duda, A. Bonfanti, A. Yonis, A. Harris, P. Mosaffa, Y. Fujita, A. Kabla, Y. Mao, B. Baum, J. Muñoz, M. Miodownik, G. Charras

Nature Physics (2019). Vol. 15, pp. 839-847 (Preprint)

Abstract
Epithelial monolayers are one-cell-thick tissue sheets that line most of the body surfaces, separating internal and external environments. As part of their function, they must withstand extrinsic mechanical stresses applied at high strain rates. However, little is known about how monolayers respond to mechanical deformations. Here, by subjecting suspended epithelial monolayers to stretch, we find that they dissipate stresses on a minute timescale and that relaxation can be described by a power law with an exponential cut-off at timescales larger than about 10 s. This process involves an increase in monolayer length, pointing to active remodelling of cellular biopolymers at the molecular scale during relaxation. Strikingly, monolayers consisting of tens of thousands of cells relax stress with similar dynamics to single rounded cells, and both respond similarly to perturbations of the actomyosin cytoskeleton. By contrast, cell–cell junctional complexes and intermediate filaments do not relax tissue stress, but form stable connections between cells, allowing monolayers to behave rheologically as single cells. Taken together, our data show that actomyosin dynamics governs the rheological properties of epithelial monolayers, dissipating applied stresses and enabling changes in monolayer length.
Abstract
Epithelial monolayers are one-cell-thick tissue sheets that line most of the body surfaces, separating internal and external environments. As part of their function, they [...]
- 19
- read
Mechanical oscillations in biological tissues as a result of delayed rest-length changes
J. Muñoz, M. Dingle, M. Wenzel

Physical review E: statistical, nonlinear and soft matter physics (2018). Vol. 98 (5), pp. 1-15 (Preprint)

Abstract
Oscillatory behavior in tissue biology is ubiquitous and may be observed in the form of either chemical or mechanical signals. We present here and solve a mechanical cell model that exhibits oscillations emerging from delayed viscoelastic rheological laws. These include a time delay between the mechanical response and the rest-length changes, which evolve proportionally to the delayed cell deformation and use a remodeling rate parameter. We show that different regimes (no oscillatory response, sustained oscillations, and unstable oscillations) are obtained for different values of the delay or the remodeling rate. The results are analytically demonstrated in a one-dimensional problem with one and two cells that are represented by simple line elements. Oscillations of the cell deformations are obtained whenever different delays coexist, or when the delay is size-dependent. We also extend our results to a multicellular two-dimensional vertex model that includes the same rheological law, and which inherits the presence of critical values of the delay or remodeling rate. We numerically show that indeed the size-dependent rest-length changes induce oscillations in the cell shapes and areas.
Abstract
Oscillatory behavior in tissue biology is ubiquitous and may be observed in the form of either chemical or mechanical signals. We present here and solve a mechanical cell [...]
- 5
- read
Computation of forces from deformed visco-elastic biological tissues
J. Muñoz, D. Amat, V. Conte

Inverse Problems (2018). Vol. 34 (4), pp. 1-27 (Preprint)

Abstract
We present a least-squares based inverse analysis of visco-elastic biological tissues. The proposed method computes the set of contractile forces (dipoles) at the cell boundaries that induce the observed and quantified deformations. We show that the computation of these forces requires the regularisation of the problem functional for some load configurations that we study here. The functional measures the error of the dynamic problem being discretised in time with a second-order implicit time-stepping and in space with standard finite elements. We analyse the uniqueness of the inverse problem and estimate the regularisation parameter by means of an L-curved criterion. We apply the methodology to a simple toy problem and to an in vivo set of morphogenetic deformations of the Drosophila embryo.
Abstract
We present a least-squares based inverse analysis of visco-elastic biological tissues. The proposed method computes the set of contractile forces (dipoles) at the cell boundaries [...]
- 3
- read
Hybrid cell-centred/vertex model for multicellular systems with equilibrium-preserving remodelling
P. Mosaffa, A. Rodríguez-Ferran, J. Muñoz

Int. J. Numer. Meth. Bio. Engng. (2018). Vol. 34 (3), pp. 1-33 (Preprint)

Abstract
We present a hybrid cell‐centred/vertex model for mechanically simulating planar cellular monolayers undergoing cell reorganisation. Cell centres are represented by a triangular nodal network, while the cell boundaries are formed by an associated vertex network. The two networks are coupled through a kinematic constraint which we allow to relax progressively. Special attention is paid to the change of cell‐cell connectivity due to cell reorganisation or remodelling events. We handle these situations by using a variable resting length and applying an Equilibrium‐Preserving Mapping on the new connectivity, which computes a new set of resting lengths that preserve nodal and vertex equilibrium. We illustrate the properties of the model by simulating monolayers subjected to imposed extension and during a wound healing process. The evolution of forces and the Equilibrium‐Preserving Mapping are analysed during the remodelling events. As a by‐product, the proposed technique enables to recover fully vertex or fully cell‐centred models in a seamless manner by modifying a numerical parameter of the model.
Abstract
We present a hybrid cell‐centred/vertex model for mechanically simulating planar cellular monolayers undergoing cell reorganisation. Cell centres are represented by a triangular [...]
- 4
- read
Preface
J. Font-Llagunes, J. Muñoz, J. Ambrósio

Multibody system dynamics (2017). Vol. 39 (1), pp. 1-2
- 4
- read
Porous-based rheological model for tissue fluidisation
N. Asadipour, X. Trepat, J. Muñoz

Journal of the Mechanics and Physics of Solids (2016). Vol. 96, pp. 535-549 (Preprint)

Abstract
It has been experimentally observed that cells exhibit a fluidisation process when subjected to a transient stretch, with an eventual recovery of the mechanical properties upon removal of the applied deformation. This fluidisation process is characterised by a decrease of the storage modulus and an increase of the phase angle. We propose a rheological model which is able to reproduce this combined mechanical response. The model is described in the context of continua and adapted to a cell-centred particle system that simulates cell–cell interactions. Mechanical equilibrium is coupled with two evolution laws: (i) one for the reference configuration, and (ii) another for the porosity or polymer density. The first law depends on the actual strain of the tissue, while the second assumes different remodelling rates during porosity increase and decrease. The theory is implemented on a particle based model and tested on a stretching experiment. The numerical results agree with the experimental measurements for different stretching magnitudes
Abstract
It has been experimentally observed that cells exhibit a fluidisation process when subjected to a transient stretch, with an eventual recovery of the mechanical [...]
- 4
- read
Collective cell durotaxis emerges from long-range intercellular force transmission
R. Sunyer, V. Conte, J. Escribano, A. Elosegui-Artola, A. Labernadie, L. Valon, D. Navajas, J. García-Aznar, J. Muñoz, P. Roca-Cusachs, X. Trepat

Science (2016). Vol. 353 (6304), pp. 1157-1161 (Preprint)

Abstract
The ability of cells to follow gradients of extracellular matrix stiffness—durotaxis—has been implicated in development, fibrosis, and cancer. Here, we found multicellular clusters that exhibited durotaxis even if isolated constituent cells did not. This emergent mode of directed collective cell migration applied to a variety of epithelial cell types, required the action of myosin motors, and originated from supracellular transmission of contractile physical forces. To explain the observed phenomenology, we developed a generalized clutch model in which local stick-slip dynamics of cell-matrix adhesions was integrated to the tissue level through cell-cell junctions. Collective durotaxis is far more efficient than single-cell durotaxis; it thus emerges as a robust mechanism to direct cell migration during development, wound healing, and collective cancer cell invasion.
Abstract
The ability of cells to follow gradients of extracellular matrix stiffness—durotaxis—has been implicated in development, fibrosis, and cancer. Here, we found multicellular [...]
- 4
- read
Non-regularised inverse finite element analysis for 3D traction force microscopy
J. Muñoz

Int. J. Numer. Anal. Model. (2016). Vol. 13 (5), pp. 763-781 (Preprint)

Abstract
The tractions that cells exert on a gel substrate from the observed displacements is an increasingly attractive and valuable information in biomedical experiments. The computation of these tractions requires in general the solution of an inverse problem. Here, we resort to the discretisation with finite elements of the associated direct variational formulation, and solve the inverse analysis using a least square approach. This strategy requires the minimisation of an error functional, which is usually regularised in order to obtain a stable system of equations with a unique solution. In this paper we show that for many common three-dimensional geometries, meshes and loading conditions, this regularisation is unnecessary. In these cases, the computational cost of the inverse problem becomes equivalent to a direct finite element problem. For the non-regularised functional, we deduce the necessary and sufficient conditions that the dimensions of the interpolated displacement and traction fields must preserve in order to exactly satisfy or yield a unique solution of the discrete equilibrium equations. We apply the theoretical results to some illustrative examples and to real experimental data. Due to the relevance of the results for biologists and modellers, the article concludes with some practical rules that the finite element discretisation must satisfy.
Abstract
The tractions that cells exert on a gel substrate from the observed displacements is an increasingly attractive and valuable information in biomedical experiments. The computation [...]
- 4
- read
Control of cell-cell forces and collective cell dynamics by the intercellular adhesome
E. Bazellières, V. Conte, A. Elosegui-Artola, X. Serra-Picamal, M. Bintanel-Morcillo, P. Roca-Cusachs, J. Muñoz, M. Sales-Pardo, R. Guimerà, X. Trepat

Nature Cell Biology (2015). Vol. 17 (4), pp. 409-420

Abstract
Dynamics of epithelial tissues determine key processes in development, tissue healing and cancer invasion. These processes are critically influenced by cell–cell adhesion forces. However, the identity of the proteins that resist and transmit forces at cell–cell junctions remains unclear, and how these proteins control tissue dynamics is largely unknown. Here we provide a systematic study of the interplay between cell–cell adhesion proteins, intercellular forces and epithelial tissue dynamics. We show that collective cellular responses to selective perturbations of the intercellular adhesome conform to three mechanical phenotypes. These phenotypes are controlled by different molecular modules and characterized by distinct relationships between cellular kinematics and intercellular forces. We show that these forces and their rates can be predicted by the concentrations of cadherins and catenins. Unexpectedly, we identified different mechanical roles for P-cadherin and E-cadherin; whereas P-cadherin predicts levels of intercellular force, E-cadherin predicts the rate at which intercellular force builds up.
Abstract
Dynamics of epithelial tissues determine key processes in development, tissue healing and cancer invasion. These processes are critically influenced by cell–cell [...]
- 6
- read
Mapping forces and kinematics during collective cell migration
X. Serra-Picamal, V. Conte, R. Sunyer, J. Muñoz, X. Trepat

Methods in Cell Biology (2015). Vol. 125, pp. 309-330

Abstract
Fundamental biological processes including morphogenesis and tissue repair require cells to migrate collectively. In these processes, epithelial or endothelial cells move in a cooperative manner coupled by intercellular junctions. Ultimately, the movement of these multicellular systems occurs through the generation of cellular forces, exerted either on the substrate via focal adhesions (cell–substrate forces) or on neighboring cells through cell–cell junctions (cell–cell forces). Quantitative measurements of multicellular forces and kinematics with cellular or subcellular resolution have become possible only in recent years. In this chapter, we describe some of these techniques, which include particle image velocimetry to map cell velocities, traction force microscopy to map forces exerted by cells on the substrate, and monolayer stress microscopy to map forces within and between cells. We also describe experimental protocols to perform these measurements. The combination of these techniques with high-resolution imaging tools and molecular perturbations will lead to a better understanding of the mechanisms underlying collective cell migration in health and disease.
Abstract
Fundamental biological processes including morphogenesis and tissue repair require cells to migrate collectively. In these processes, epithelial or endothelial [...]
- 11
- read
AAR-based decomposition method for lower-bound limit analysis
J. Muñoz, R. Syednima

Proceedings of the ICE-Engineering and Computational Mechanics (2015). Vol. 168 (4), pp. 169-177 (Preprint)

Abstract
Despite the recent progress in optimisation techniques, finite-element stability analysis of realistic three-dimensional problems is still hampered by the size of the resulting optimisation problem. Current solvers may take a prohibitive computational time, if they give a solution at all. The possible remedies to this are the design of adaptive de-remeshing techniques, decomposition of the system of equations or of the optimisation problem. This paper concentrates on the last approach, and presents an algorithm especially suited for limit analysis. Optimisation problems in limit analysis are in general convex but non-linear. This fact renders the design of decomposition techniques specially challenging. The efficiency of general approaches such as Benders or Dantzig–Wolfe is not always satisfactory, and strongly depends on the structure of the optimisation problem. This work presents a new method that is based on rewriting the feasibility region of the global optimisation problem as the intersection of two subsets. By resorting to the averaged alternating reflections (AAR) method in order to find the distance between the sets, the optimisation problem is successfully solved in a decomposed manner. Some representative examples illustrate the application of the method and its efficiency with respect to other well-known decomposition algorithms
Abstract
Despite the recent progress in optimisation techniques, finite-element stability analysis of realistic three-dimensional problems is still hampered by the size of the resulting [...]
- 12
- read
Cell-centred model for the simulation of curved cellular monolayers
P. Mosaffa, N. Asadipour, D. Millán, A. Rodríguez-Ferran, J. Muñoz

Computational Particle Mechanics (2015). Vol. 2 (4), pp. 359-370 (Preprint)

Abstract
This paper presents a cell-centred model for the simulation of planar and curved multicellular soft tissues. We propose a computational model that includes stress relaxation due to cell reorganisation (intercellular connectivity changes) and cytoskeleton remodelling (intracellular changes). Cells are represented by their cell centres, and their mechanical interaction is modelled through active non-linear elastic laws with a dynamically changing resting length. Special attention is paid to the handling of connectivity changes between cells, and the relaxation that the tissues exhibit under these topological changes. Cell–cell connectivity is computed by resorting to a Delaunay triangulation, which is combined with a mapping technique in order to obtain triangulations on curved manifolds. Our numerical results show that even a linear elastic cell–cell interaction model may induce a global non-linear response due to the reorganisation of the cell connectivity. This plastic-like behaviour is combined with a non-linear rheological law where the resting length depends on the elastic strain, mimicking the global visco-elastic response of tissues. The model is applied to simulate the elongation of planar and curved monolayers.
Abstract
This paper presents a cell-centred model for the simulation of planar and curved multicellular soft tissues. We propose a computational model that includes stress relaxation [...]
- 5
- read
Forces driving epithelial wound healing
A. Brugués, E. Anon, V. Conte, J. Veldhuis, M. Gupta, J. Colombelli, J. Muñoz, G. Brodland, B. Ladoux, X. Trepat

Natura Physics (2014). Vol. 10 (9), pp. 683-690 (Preprint)

Abstract
A fundamental feature of multicellular organisms is their ability to self-repair wounds through the movement of epithelial cells into the damaged area. This collective cellular movement is commonly attributed to a combination of cell crawling and ‘purse-string’ contraction of a supracellular actomyosin ring. Here we show by direct experimental measurement that these two mechanisms are insufficient to explain force patterns observed during wound closure. At early stages of the process, leading actin protrusions generate traction forces that point away from the wound, showing that wound closure is initially driven by cell crawling. At later stages, we observed unanticipated patterns of traction forces pointing towards the wound. Such patterns have strong force components that are both radial and tangential to the wound. We show that these force components arise from tensions transmitted by a heterogeneous actomyosin ring to the underlying substrate through focal adhesions. The structural and mechanical organization reported here provides cells with a mechanism to close the wound by cooperatively compressing the underlying substrate.
Abstract
A fundamental feature of multicellular organisms is their ability to self-repair wounds through the movement of epithelial cells into the damaged area. This collective cellular [...]
- 13
- read
A truss element for modelling reversible softening in living tissues
J. Muñoz, V. Conte, N. Asadipour, M. Miodownik

Mechanics Research Communications (2013). Vol. 49, pp. 44-49 (Preprint)

Abstract
We resort to non-linear viscoelasticity to develop a truss element able to model reversible softening in lung epithelial tissues undergoing transient stretch. Such a Maxwell truss element is built by resorting to a three-noded element whose mid-node is kinematically constrained to remain on the line connecting the end-nodes. The whole mechanical system undergoes an additive decomposition of the strains along the truss direction where the total contribution of the mid-node is accounted for by using a null-space projection and static condensation techniques. Assembling of such line-elements in 3D networks allows us to model extended regions of living tissues as well as their anisotropies.
Abstract
We resort to non-linear viscoelasticity to develop a truss element able to model reversible softening in lung epithelial tissues undergoing transient stretch. Such a Maxwell [...]
- 2
- read
Modelling of mixed damage on fibre reinforced composite laminates subjected to low velocity impact
J. Curiel, S. Phaneendra, J. Muñoz

International journal of damage mechanics (2013). Vol. 22 (3), pp. 356-374 (Preprint)

Abstract
This article presents a numerical technique for the computation of damage in fibre-reinforced laminated composites which is embedded into an explicit finite element method. The technique is composed of two main components. First, a novel characterisation of the directions in which the various modes of damage propagate, which produces an effective damage localisation. This is conducted taking into consideration the strain-rate dependence. Second, a new mapping between the strain and stress spaces for the computation of the damage surfaces whereby time-stepping convergence is enhanced. Additionally, new damage initiation criteria in terms of strain damage surfaces are presented. Details of the in-house code developed are presented as well as the programming features. The capabilities of the technique are shown by means of tests on single fibre-reinforced element and low velocity impact on the laminate. It is shown that delamination is located in the expected regions by gradual progression of internal damage variables.
Abstract
This article presents a numerical technique for the computation of damage in fibre-reinforced laminated composites which is embedded into an explicit finite element method. [...]
- 11
- read
Physiology-based model of cell viscoelasticity
J. Muñoz, S. Albo

Physical review E: statistical, nonlinear and soft matter physics (2013). Vol. 88 (1), pp. 1-8 (Preprint)

Abstract
The measured viscoelastic properties of biological tissues is the result of the passive and active response of the cells. We propose an evolution law of the remodeling process in the cytoskeleton which is able to mimic the viscous properties of biological cellular tissues. Our model is based on dynamical changes of the resting length. We show that under the small strain regime, the linear rheology models are recovered, with the relaxation time being replaced by the cell resistance to remodel. We implement the one-dimensional model into network systems of two and three dimensions, and show that the same conclusions may be drawn for those systems.
Abstract
The measured viscoelastic properties of biological tissues is the result of the passive and active response of the cells. We propose an evolution law of the remodeling process [...]
- 8
- read
A biomechanical analysis of ventral furrow formation in the Drosophila melanogaster embryo
V. Conte, F. Ulrich, B. Baum, J. Muñoz, J. Veldhuis, W. Brodland, M. Miodownik

PloS One (2012). Vol. 7 (4), pp. 1-17

Abstract
The article provides a biomechanical analysis of ventral furrow formation in the Drosophila melanogaster embryo. Ventral furrow formation is the first large-scale morphogenetic movement in the fly embryo. It involves deformation of a uniform cellular monolayer formed following cellularisation, and has therefore long been used as a simple system in which to explore the role of mechanics in force generation. Here we use a quantitative framework to carry out a systematic perturbation analysis to determine the role of each of the active forces observed. The analysis confirms that ventral furrow invagination arises from a combination of apical constriction and apical–basal shortening forces in the mesoderm, together with a combination of ectodermal forces. We show that the mesodermal forces are crucial for invagination: the loss of apical constriction leads to a loss of the furrow, while the mesodermal radial shortening forces are the primary cause of the internalisation of the future mesoderm as the furrow rises. Ectodermal forces play a minor but significant role in furrow formation: without ectodermal forces the furrow is slower to form, does not close properly and has an aberrant morphology. Nevertheless, despite changes in the active mesodermal and ectodermal forces lead to changes in the timing and extent of furrow, invagination is eventually achieved in most cases, implying that the system is robust to perturbation and therefore over-determined.
Abstract
The article provides a biomechanical analysis of ventral furrow formation in the Drosophila melanogaster embryo. Ventral furrow formation is the first large-scale [...]
- 12
- read
Diffusion-reaction model for Drosophila embryo development
R. Allena, J. Muñoz, D. Aubry

Comput. Methods Appl. Mech. Engrg., (2011). Vol. 16 (3), pp. 235-248

Abstract
During the early stages of gastrulation in Drosophila embryo, the epithelial cells composing the single tissue layer of the egg undergo large strains and displacements. These movements have been usually modelled by decomposing the total deformation gradient in an (imposed or strain/stress dependent) active part and a passive response. Although the influence of the chemical and genetic activity in the mechanical response of the cell has been experimentally observed, the effects of the mechanical deformation on the latter have been far less studied, and much less modelled. Here, we propose a model that couples morphogen transport and the cell mechanics during embryogenesis. A diffusion-reaction equation is introduced as an additional mechanical regulator of morphogenesis. Consequently, the active deformations are not directly imposed in the analytical formulation, but they rather depend on the morphogen concentration, which is introduced as a new variable. In this study, we show that strain patterns similar to those observed during biological experiments can be reproduced by properly combining the two phenomena. In addition, we use a novel technique to parameterise the embryo geometry by solving two Laplace problems with specific boundary conditions. We apply the method to two morphogenetic movements: ventral furrow invagination and germ band extension. The matching between our results and the observed experimental deformations confirms that diffusion-reaction of morphogens can actually be controlling large morphogenetic movements.
Abstract
During the early stages of gastrulation in Drosophila embryo, the epithelial cells composing the single tissue layer of the egg undergo large strains and displacements. [...]
- 6
- read
Stress-dependent morphogenesis: continuum mechanics and truss systems
J. Muñoz, V. Conte, M. Miodownik

Biomechanics and Modeling in Mechanobiology (2010). Vol. 9 (4), pp. 451-467

Abstract
A set of equilibrium equations is derived for the stress-controlled shape change of cells due to the remodelling and growth of their internal architecture. The approach involves the decomposition of the deformation gradient into an active and a passive component; the former is allowed to include a growth process, while the latter is assumed to be hyperelastic and mass-preserving. The two components are coupled with a control function that provides the required feedback mechanism. The balance equations for general continua are derived and, using a variational approach, we deduce the equilibrium equations and study the effects of the control function on these equations. The results are applied to a truss system whose function is to simulate the cytoskeletal network constituted by myosin microfilaments and microtubules, which are found experimentally to control shape change in cells. Special attention is paid to the conditions that a thermodynamically consistent formulation should satisfy. The model is used to simulate the multicellular shape changes observed during ventral furrow invagination of the Drosophila melanogaster embryo. The results confirm that ventral furrow invagination can be achieved through stress control alone, without the need for other regulatory or signalling mechanisms. The model also reveals that the yolk plays a distinct role in the process, which is different to its role during invagination with externally imposed strains. In stress control, the incompressibility constraint of the yolk leads, via feedback, to the generation of a pressure in the ventral zone of the epithelium that eventually eases its rise and internalisation.
Abstract
A set of equilibrium equations is derived for the stress-controlled shape change of cells due to the remodelling and growth of their internal architecture. The approach involves [...]
- 11
- read
− 2gy) (1 + y )dx La soluc´íon es, y = tag(α)x −g 2u (1 + tag (α))x Esta es la solución que parte del origen de coordenadas con una velocidad , u , yángulo α y pasa por un punto,B(x
. A, A. Lechuga

Antonio Lechuga's personal collection (2020). 2
- 2
- read
Robust mechanisms of ventral furrow invagination require the combination of cellular shape changes
V. Conte, J. Muñoz, B. Baum, M. Miodownik

Physical Biology (2009). Vol. 6 (1), pp. 16010-16021

Abstract
Ventral furrow formation in Drosophila is the first large-scale morphogenetic movement during the life of the embryo, and is driven by co-ordinated changes in the shape of individual epithelial cells within the cellular blastoderm. Although many of the genes involved have been identified, the details of the mechanical processes that convert local changes in gene expression into whole-scale changes in embryonic form remain to be fully understood. Biologists have identified two main cell deformation modes responsible for ventral furrow invagination: constriction of the apical ends of the cells (apical wedging) and deformation along their apical–basal axes (radial lengthening/shortening). In this work, we used a computer 2D finite element model of ventral furrow formation to investigate the ability of different combinations of three plausible elementary active cell shape changes to bring about epithelial invagination: ectodermal apical–basal shortening, mesodermal apical–basal lengthening/shortening and mesodermal apical constriction. We undertook a systems analysis of the biomechanical system, which revealed many different combinations of active forces (invagination mechanisms) were able to generate a ventral furrow. Two important general features were revealed. First that combinations of shape changes are the most robust to environmental and mutational perturbation, in particular those combining ectodermal pushing and mesodermal wedging. Second, that ectodermal pushing plays a big part in all of the robust mechanisms (mesodermal forces alone do not close the furrow), and this provides evidence that it may be an important element in the mechanics of invagination in Drosophila.
Abstract
Ventral furrow formation in Drosophila is the first large-scale morphogenetic movement during the life of the embryo, and is driven by co-ordinated changes in the [...]
- 19
- read
A 3D finite element model of ventral furrow invagination in the Drosophila melanogaster embryo
V. Conte, J. Muñoz, M. Miodownik

Journal of the Mechanical Behavior of Biomedical Materials (2008). Vol. 1 (2), pp. 188-198

Abstract
The paper describes a mechanical model of epithelial tissue development in Drosophila embryos to investigate a buckling phenomenon called invagination. The finite element method is used to model this ventral furrow formation in 3D by decomposing the total deformation into two parts: an imposed active deformation, and an elastic passive deformation superimposed onto the latter. The model imposes as boundary conditions (i) a constant yolk volume and (ii) a sliding contact condition of the cells against the vitelline membrane, which is interpolated as a B-Spline surface. The active deformation simulates the effects of apical constriction and apico-basal elongation of cells. This set of local cellular mechanisms leads to global shape changes of the embryo which are associated with known gene expressions. Using the model we have tested different plausible hypotheses postulated to account for the mechanical behaviour of epithelial tissues. In particular, we conclude that only certain combinations of local cell shape change can successfully reproduce the invagination process. We have quantitatively compared the model with a 2D model and shown that it exhibits a more robust invagination phenomenon. The 3D model has also revealed that invagination causes a yolk flow from the central region to the anterior and posterior ends of the embryo, causing an accordion-like global compression and expansion wave to move through the embryo. Such a phenomenon cannot be described by 2D models.
Abstract
The paper describes a mechanical model of epithelial tissue development in Drosophila embryos to investigate a buckling phenomenon called invagination. The finite [...]
- 18
- read
On the modelling of incompressibility in linear and non-linear elasticity with the master-slave approach
J. Muñoz

International Journal for Numerical Methods in Engineering (2008). Vol. 74 (2), pp. 269-293

Abstract
The master–slave approach is adapted to model the kinematic constraints encountered in incompressibility. The method presented here allows us to obtain discrete displacement and pressure fields for arbitrary finite element formulations that have discontinuous pressure interpolations. The resulting displacements satisfy exactly the incompressibility constraints in a weak sense, and are obtained by solving a system of equations with the minimum (independent) degrees of freedom. In linear analysis, the method reproduces the well‐known stability results for inf–sup compliant elements, and permits to compute the pressure modes (physical or spurious) when they exist. By rewriting the equilibrium equations of a hyperelastic material, the method is extended to non‐linear elasticity, while retaining the exact fulfilment of the incompressibility constraints in a weak sense. Problems with analytical solution in two and three dimensions are tested and compared with other solution methods.
Abstract
The master–slave approach is adapted to model the kinematic constraints encountered in incompressibility. The method presented here allows us to obtain discrete displacement [...]
- 0
- read
Modelling unilateral frictionless contact using the null-space method and cubic B-Spline interpolation
J. Muñoz

Comput. Methods Appl. Mech. Engrg., (2008). Vol. 197 (9-12), pp. 979-993

Abstract
The analysis of unilateral sliding contact in elasticity is equivalent to a minimisation problem subjected to a set of inequality constraints. However, the presence of boundary discontinuities, such as those stemming from the spatial discretisation, appears as a major problem to determine the set of active constraints. This work introduces a smoothing technique of the master surface resorting to cubic B-Spline interpolation, which is C1 continuous in contact situations between elastic and rigid bodies, and G1 continuous in elastic–elastic contact problems. The technique is applied in conjunction with the null-space method, where the solution is searched in an unconstrained manifold. The resulting formulation eases the contact transition along the master surface, and recovers the quadratic convergence of the iterative Newton–Raphson process. The robustness of the method is demonstrated using 2D and 3D examples.
Abstract
The analysis of unilateral sliding contact in elasticity is equivalent to a minimisation problem subjected to a set of inequality constraints. However, the presence of boundary [...]
- 11
- read
A deformation gradient decomposition method for the analysis of the mechanics of morphogenesis
J. Muñoz, K. Barrett, M. Miodownik

Journal of Biomechanics (2007). Vol. 40 (6), pp. 1372-1380

Abstract
A new finite element model is proposed for the analysis of the mechanical aspects of morphogenesis and tested on the biologically well studied gastrulation phenomenon, in particular ventral furrow invagination of the Drosophila melanogaster embryo. A set of mechanisms are introduced in the numerical model, which lead to the observed deformed shapes. We split the total deformation into two parts: an imposed active deformation, and an elastic deformation superimposed onto the latter. The active deformation simulates the effects of apical constriction and apico-basal elongation. These mechanisms are associated with known gene expressions and so in this way we attempt to bridge the well explored signalling pathways, and their associated phenotypes in a mechanical model. While the former have been studied in depth, much less can be said about the forces they produce and the mechanisms involved. From the numerical results, we are able to test different plausible mechanical hypotheses that generate the necessary folding observed in the invagination process. In particular, we conclude that only certain ratios between both modes (apical constriction and apico-basal elongation) can successfully reproduce the invagination process. The model also supports the idea that this invagination requires the contribution of several mechanisms, and that their redundancy provides the necessary robustness.
Abstract
A new finite element model is proposed for the analysis of the mechanical aspects of morphogenesis and tested on the biologically well studied gastrulation phenomenon, in [...]
- 0
- read
On the numerical simulation of fatigue driven delamination with interface elements
J. Muñoz, U. Galvanetto, P. Robinson

International Journal of Fatigue (2006). Vol. 28 (10), pp. 1136-1146

Abstract
This paper investigates robustness and efficiency of a method for the numerical simulation of fatigue driven delamination growth in laminated composites. The method is based on the fatigue degradation of a bilinear interface element and the paper studies the dependence of the predicted behaviour on the interface element size, Δl, and the number of cycles per increment, ΔN. It is concluded that the value of ΔN must be small enough to ensure that at least two increments are required to advance the crack by Δl. It is also shown that Δl must be sufficiently small so that at least two interface elements lie in the cohesive zone ahead of the crack tip. The paper also proposes extensions of the model to include a non-zero minimum cyclic load value and an improved algorithm for mixed-mode fatigue driven delamination growth.
Abstract
This paper investigates robustness and efficiency of a method for the numerical simulation of fatigue driven delamination growth in laminated composites. The method [...]
- 15
- read
Sliding joints in 3D beams: conserving algorithms using the master-slave approach
J. Muñoz, G. Jelenic

Moltibody System Dynamics (2006). Vol. 16 (3), pp. 137-261

Abstract
This paper proposes two time-integration algorithms for motion of geometrically exact 3D beams under sliding contact conditions. The algorithms are derived using the so-called master–slave approach, in which constraint equations and the related time-integration of a system of differential and algebraic equations are eliminated by design. Specifically, we study conservation of energy and momenta when the sliding conditions on beams are imposed and discuss their algorithmic viability. Situations where the contact jumps to adjacent finite elements are analysed in detail and the results are tested on two representative numerical examples. It is concluded that an algorithmic preservation of kinematic constraint conditions is of utmost importance.
Abstract
This paper proposes two time-integration algorithms for motion of geometrically exact 3D beams under sliding contact conditions. The algorithms are derived using the so-called [...]
- 5
- read
Sliding contact conditions using the master–slave approach with application on geometrically non-linear beams
J. Muñoz, G. Jelenic

International Journal of Solid and Structures (2004). Vol. 41 (24-25), pp. 6963-6992

Abstract
Frictionless sliding conditions between two bodies are usually defined using either the method of Lagrangian multipliers or by prescribing an artificial (penalty) stiffness which resists the penetration at the contact point. Both of these methods impose the condition that the contact force should be normal to the contact surface, with the Lagrangian multiplier or the penalty parameter serving as a measure of this force. In this work, an alternative approach is undertaken: the frictionless sliding condition is defined through a relationship between nodal parameters of the virtual displacements of a discretised principle of virtual work. This method, which does not involve additional force parameters or degrees of freedom, is known as the master–slave or the minimum-set method and is particularly convenient for displacement-based finite-element implementation. The method is analysed in detail in context of bilateral sliding constraints characteristic of prismatic and cylindrical joints in flexible beam assemblies undergoing large overall motion. Two numerical examples are presented and assessed against the results in the literature.
Abstract
Frictionless sliding conditions between two bodies are usually defined using either the method of Lagrangian multipliers or by prescribing an artificial (penalty) stiffness [...]
- 8
- read
Master-slave approach for the modelling of joints with dependent degrees of freedom in flexible mechanisms
J. Muñoz, G. Jelenic, M. Crisfield

Commun. Numer. Meth. Engng (2003). Vol. 19 (9), pp. 689-702

Abstract
The analysis of multibody systems requires an exact description of the kinematics of the joints involved. In the present work the master–slave approach is employed and endowed with the possibility of including several more complex types of joints. We present the formulation for joints where some relation between the different released degrees of freedom exists such as the screw joint, the rack‐and‐pinion joint or the cam joint. These joints are implemented in conjunction with geometrically exact beams and an energy‐momentum conserving time‐stepping algorithm.
Abstract
The analysis of multibody systems requires an exact description of the kinematics of the joints involved. In the present work the master–slave approach is employed and [...]
- 10
- read
An effective aggregate signature scheme suitable for VANETs and its implementation using OPNET IT Guru
D. Wang

Rev. int. métodos numér. cálc. diseño ing. (2020). Vol. 36, (3), 39

Abstract
Vehicular Ad hoc NETwork (VANET) is a multi-hop and high-speed mobile wireless communication network. In order to realize the security authentication of the information transmission between vehicle nodes in vehicular ad hoc networks, a certificateless aggregate signature scheme is designed. The proposed scheme uses certificateless cryptography, which eliminates the complex maintenance cost of certificate and solves the problem of key escrow. Communicating through pseudonyms and nodes around the roadside units generated, the conditional privacy protection is achieved for vehicle users. Then, the efficiency of the scheme is analyzed, and the relationship between the traffic density in VANETs environment and the time delay of message verification is simulated. The results show that the scheme satisfies the message authentication, anonymity, unforgeability and traceability, as well as the higher communication efficiency and the shorter delay of message verification, which is more suitable for dynamic vehicular ad hoc network environment.
Abstract
Vehicular Ad hoc NETwork (VANET) is a multi-hop and high-speed mobile wireless communication network. In order to realize the security authentication of the information transmission [...]
- 288
- read
Thermic Model to Predict Biogas Production in Unheated Fixed-Dome Digesters Buried in the Ground
G. Terradas, C. Pham, J. Triolo, J. Martí-Herrero, S. Sommer

Environmental science and technology (2014). Vol. 48 (6), pp. 3253-3262

Abstract
In many developing countries simple biogas digesters are used to produce energy for domestic purposes from anaerobic digestion of animal manure. We developed a simple, one-dimensional (1-D), thermal model with easily-available input data for unheated, unstirred, uninsulated, fixed-dome digesters buried in the soil to study heat transfer between biogas digester and its surroundings. The predicted temperatures in the dome, biogas and slurry inside the digester and the resulting biogas production are presented and validated. The model was well able to estimate digester temperature (linear slope nearly 1, R2 0.96). The results show that a suitable soil temperature profile renders the use of a full thermal model superfluous for buried, unheated, unstirred and uninsulated digesters. Model validation for methane production gave root-mean-square (RMSE) and relative-root-mean-square errors (rRMSEP(%)) of 54.4 L CH4 digester-1 day-1 and 35.4%, respectively. The validation result was considerably improved if only using winter data (RMSE = 26.1 L CH4 digester-1 day-1; rRMSEP(%) = 17.7%). The model performed satisfactorily in light of the uncertainties attached to it. Since unheated digesters suffer critically low methane production during the winter, the model could be particularly useful for assessing methane production and for improving the ability of unheated digesters to provide sufficient energy during cold periods.
Abstract
In many developing countries simple biogas digesters are used to produce energy for domestic purposes from anaerobic digestion of animal manure. We developed a simple, one-dimensional [...]
- 13
- read
Low cost tubular digesters as appropriate technology for widespread application: Results and lessons learned from Bolivia
J. Martí-Herrero, M. Chipana, C. Cuevas, G. Paco, V. Serrano, B. Zymla, K. Heising, J. Sologuren, A. Gamarra

Renewable Energy (2014). Vol. 71, pp. 156-165

Abstract
This paper presents the results and lessons learned from four and a half years of implementing low cost tubular digesters in Bolivia. The selection of this technology is justified in comparison with other popular technologies such as fixed dome or floating drum digesters. The highlighted weakness of the tubular model (its short life expectancy), is transformed into a strength, making the low cost tubular digester an appropriate technology for widespread application. The experiences in Bolivia show that the success of biogas programs depend more on socio-economic factors than on the validated technology selected, suggesting that local circumstances are a critical, and often underestimated, factor to be taken into consideration in the praxis. Finally, some testimonies of the use of biol (bio-slurry or effluent) are reported, identifying the high potential of this anaerobic digestion product that provides a food sovereignty approach, reduced expansion of the agricultural frontier, increased agricultural productivity and hence family income, that other household energizing systems do not have. A brief report of lessons learned is also included
Abstract
This paper presents the results and lessons learned from four and a half years of implementing low cost tubular digesters in Bolivia. The selection of this technology is justified [...]
- 2
- read
Improvement through low cost biofilm carrier in anaerobic tubular digestion in cold climate regions
J. Martí-Herrero, R. Alvarez, R. Rojas, L. Aliaga, R. Cespedes, J. Carbonell

Bioresource Technology (2014). Vol. 167, pp. 87-93

Abstract
The aim of this research is to evaluate the increase of biogas production with low cost tubular digesters in cold climates using PET rings inside the reactor. Two similar digesters have been operated and monitored in cold weather conditions and have been fed with cow manure. Digester 1 was filled with PET – rings as a biofilm carrier, Digester 2 was kept as a reference. Through the PET – rings the functional surface could be increased by a factor 4.2. The results show that 44% more biogas per KgSV has been produced with the biofilm carrier in use (0.33 m3/kgSV) (reference digester −0.23 m3/kgSV), at an organic load rate of 0.26 kgSV/m3/d. The thermal performance shows that with an adaptation of the low cost tubular digester the slurry temperature can be raised up to 16.6 °C (average) by surrounding temperature of 6.1 °C (average) without using any active heating system.
Abstract
The aim of this research is to evaluate the increase of biogas production with low cost tubular digesters in cold climates using PET rings inside the reactor. Two similar [...]
- 7
- read
Cow, sheep and llama manure at psychrophilic anaerobic co-digestion with low cost tubular digesters in cold climate and high altitude
J. Martí-Herrero, R. Alvarez, R. Cespedes, M. Rojas, V. Conde, L. Aliaga, M. Balboa, S. Danov

Bioresource Technology (2015). Vol. 181, pp. 238-246

Abstract
The aim of this research is to evaluate the co-digestion of cow and llama manure combined with sheep manure, in psychrophilic conditions and real field low cost tubular digesters adapted to cold climate. Four digesters were monitored in cold climate conditions; one fed with cow manure, a second one with llama manure, the third one with co-digestion of cow–sheep manure and the fourth one was fed with llama–sheep manure. The slurry had a mean temperature of 16.6 °C, the organic load rate was 0.44 kgvs m−3 d−1 and the hydraulic retention time was 80 days. After one hundred days biogas production was stable, as was the methane content and the pH of the effluent. The co-digestion of cow–sheep manure results in a biogas production increase of 100% compared to the mono-digestion of cow manure, while co-digestion of llama–sheep manure results in a decrease of 50% in biogas production with respect to mono-digestion of llama manure
Abstract
The aim of this research is to evaluate the co-digestion of cow and llama manure combined with sheep manure, in psychrophilic conditions and real field low cost tubular digesters [...]
- 0
- read
The influence of users' behavior on biogas production from low cost tubular digesters: A technical and socio-cultural field analysis
J. Martí-Herrero, M. Ceron, R. García, L. Pracejus, R. Alvarez, X. Cipriano

Energy for Sustainable Development (2015). Vol. 27, pp. 73-83

Abstract
The aim of this paper is to understand the influence of the user behavior on tubular digesters performance, through a technical and a social approach in the Bolivian context. Fifteen domestic digesters were evaluated, from which 6 were installed in the Altiplano and other 6 in the Andean Valleys. Data about slurry temperature, feedstock and biogas quality were collected from these 12 digesters, while daily biogas production and feeding pattern were also monitored from further three digesters in the valleys. Because of changes in user behavior along the monitoring period and particular characteristics of the digesters monitored, 5 complete patterns of biogas production and digester management were established. Furthermore, the results of a socio-cultural study with Andean families about the perception of poverty, their needs and the role played by digesters in their expectations in improving life quality, are correlated to the obtained technical data. The technical evaluation shows how the digester management seems to have a seasonal performance throughout the year according to the agricultural calendar. This means that families are more interested in using bioslurry in crops and agricultural improvements than in the use of biogas. The Bolivian government subsidy on liquefied gas seems to be one of the key issues to understand these results. Finally, data also reveals how the thermal behavior of tubular digesters adapted to cold climate that use a passive solar design, is similar to the thermal behavior of valley digesters, and therefore intends to add the same passive heating techniques for warm and tropical climates, to increase slurry temperatures and achieve a higher biogas production.
Abstract
The aim of this paper is to understand the influence of the user behavior on tubular digesters performance, through a technical and a social approach in the Bolivian context. [...]
- 27
- read
How to report biogas production when monitoring small-scale digesters in field
J. Martí-Herrero, T. Flores, R. Alvarez, D. Pérez

Biomass and Bioenergy (2016). Vol. 84, pp. 31-36

Abstract
The aim of this research was to evaluate the error originated when biogas production from field monitoring digesters, influenced by the diurnal temperature cycle, was normalized to standard conditions for pressure and temperature (273.15 K and 100 kPa) from local conditions. The biogas production data is often reported without indicating if done under local conditions, whether these conditions have been standardized and, if they have actually been standardized, the standard temperature and pressure is not indicated. In this research ambient and biogas temperature, as well as biogas production were monitored with a 30 min frequency during three consecutive days, in three different tubular digesters. Normalization was realized using the high frequency data collected as reference values, and also using daily biogas production with mean daily biogas, ambient and nearby meteorological station temperatures. The outcome of this research shows that normalization of biogas production can be obtained using daily biogas production and the daily mean ambient temperature with an overestimation by no more than 1.5%, in comparison to the normalization achieved by using high frequency data from biogas temperature and production. Using mean daily ambient temperature or mean daily biogas temperature results in the same overestimation, while using mean daily ambient temperature from a nearby airport weather station pushes the overestimation up to 2.7%. So, if ambient temperature and altitude is identified, biogas production reported in local conditions can be normalized.
Abstract
The aim of this research was to evaluate the error originated when biogas production from field monitoring digesters, influenced by the diurnal temperature [...]
- 9
- read
Household anaerobic digesters for biogas production in Latin America: A review
M. Garfí, J. Martí-Herrero, A. Garwood, I. Ferrer

Renewable and Sustainable Energy Reviews (2016). Vol. 60, pp. 599-614

Abstract
This review aims to provide an overview of household biogas digester implementation in rural areas of Latin America. It considers the history of household digesters in Latin America, including technical, environmental, social and economic aspects. Several successful experiences have been promoted during the last decade, including the creation of the Network for Biodigesters in Latin America and the Caribbean (RedBioLAC) that provides a forum to coordinate implementation and research programmes throughout the continent. Although the potential of this technology is well demonstrated, some barriers are identified, such as the need for technical improvements, lack of social acceptance and high investment costs. Thus, further efforts should be undertaken to overcome these barriers and improve the technical performance, social acceptance, economic benefits and environmental impact in order to enhance its wide-spread dissemination in energy poor communities.
Abstract
This review aims to provide an overview of household biogas digester implementation in rural areas of Latin America. It considers the history of household digesters in Latin [...]
- 41
- read
Analysis of energy, CO2 emissions and economy of the technological migration for clean cooking in Ecuador
J. Martínez, J. Martí-Herrero, S. Villacís, D. Vaca, A. Riofrio

Energy Policy (2017). Vol. 107, pp. 182-187

Abstract
The objective of this study is to analyze the CO2 emissions and economic impacts of the implementation of the National Efficient Cooking Program (NECP) in Ecuador, which aims to migrate the population from Liquefied Petroleum Gas (LPG)-based stoves to electric induction stoves. This program is rooted in the current effort to change Ecuador's energy balance, with hydroelectric power expected to generate 83.61% of national electricity by 2022, ending the need for subsidized LPG. For this analysis, the 2014 baseline situation has been compared with two future scenarios for 2022: a business-as-usual scenario and an NECP-success scenario. This study demonstrates the viability of migration from imported fossil fuels to locally-produced renewable energy as the basis for an efficient cooking facility. The new policies scenario would save US$ 1.162 billion in annual government expenditure on cooking subsidies, and reducing CO2 emissions associated to energy for cooking in 1.8 tCO2/y.
Abstract
The objective of this study is to analyze the CO2 emissions and economic impacts of the implementation of the National Efficient Cooking Program (NECP) in Ecuador, which [...]
- 15
- read
Monitoreo a escala de real de un digestor anaeróbico de bajo costo instalado en una planta de sacrificio animal
L. Castro, H. Escalante, M. Carballo, J. Martí-Herrero

Red de Biodigestores para Latina América y el Caribe (2017). Vol. 1, pp. 37-44

Abstract
Slaughterhouses generate an average 6.5 kg of manure/animal. Day. The animals has a weight between 350 and 400 kg. Manure causes economic, environmental and legal problems for the slaughter companies. The aim of this work, was to stabilize the 5,5 % of cattle manure of FOGASA Slaughterhouse by means of low cost digester. The results show organic matter reduce of 60% volatile solid and an average biogas production of 0,9 Nm/d. During anaerobic digestion process efficiencies were reach in terms of biogas production rate, specific biogas production and reduction of the pathogen content of 0,21 m3biogás/m3digestor d y 0,20 m3biogás/KgSV y 90% respectively. The digestate contained nitrogen, phosphorous and potassium, which is an important feature of fertilizers
Abstract
Slaughterhouses generate an average 6.5 kg of manure/animal. Day. The animals has a weight between 350 and 400 kg. Manure causes economic, environmental and legal problems [...]
- 6
- read
Evaluation of the low technology tubular digesters in the production of biogas from slaughterhouse wastewater treatment: Low cost tubular digesters
J. Martí-Herrero, R. Alvarez, T. Flores

Journal of Cleaner Production (2018). Vol. 199, pp. 633-642

Abstract
The aim of this research is to identify the operational conditions and performance of low cost tubular digesters producing biogas from the slaughterhouse wastewater treatment. Three full scale sequenced tubular digesters (of more than 8 m3 each), operating in a slaughterhouse in Cochabamba (Bolivia), have been monitored during one year. The Organic Load Rate (OLR) varied from 0.04 to 1.13 kgSV/m3d (0.04–1.07 kgCOD/m3d), and the Hydraulic Retention Time (HRT) from 3.2 to 87.4 days. The passive solar heating design considered drove the slurry temperature to 27 °C, 8 °C more than mean ambient temperature. The first digester has been filled with soda bottle plastic rings in order to increase the internal effective surface and evaluate, obtaining that an increment in 67% of the internal effective surface, more than doubles the biogas production and enables higher OLRs. A Biogas Production Rate (BPR) peak of 0.2 m3/m3d was observed for OLR = 0.37 kgSV/m3d, HRT = 9.7d and Specific Biogas Production (SBP) = 0.55 m3/kgSV. COD removal achieves values above 70% from HRT>19 d. A comparison with other anaerobic digestion technologies demonstrates that low cost tubular digesters with biofilm carriers, are competitive for OLR<0.5 kg COD/m3d of biogas production, and for 0.25 kgCOD/m3d, when the focus is on COD removal.
Abstract
The aim of this research is to identify the operational conditions and performance of low cost tubular digesters producing biogas from the slaughterhouse wastewater [...]
- 7
- read
Biogas from a full scale digester operated in psychrophilic conditions and fed only with fruit and vegetable waste
J. Martí-Herrero, G. Soria-Castellón, A. Díaz-de-Basurto, R. Alvarez, D. Chemisana

Renewable Energy (2019). Vol. 133, pp. 676-684

Abstract
The aim of this work was to evaluate the feasibility of treating fruit and vegetable waste from a municipal market in a full-scale anaerobic digester with minimum implementation and operational costs, i.e. without pre-treatment, clean water consumption, active heating or mixing. For this purpose, a 13.9 m3 digester that forces the submersion of solids, gains heat through solar radiation and recirculates effluent was monitored during one year of real operational and weather conditions in Bolivia. The digester was initially loaded with cow rumen (inoculum source), and after eight weeks with only fruit and vegetable waste. The digester operated at an organic loading rate (OLR) of 1 kgVS/(m3d). The normalized specific biogas production averaged 0.26 m3/kgVS and increased to 0.34 m3/kgVS when OLR was halved to 0.57 kgVS/(m3d), while the passive solar design allowed an average digester temperature of 21 °C (4 °C above the average ambient temperature).
Abstract
The aim of this work was to evaluate the feasibility of treating fruit and vegetable waste from a municipal market in a full-scale anaerobic digester with minimum implementation [...]
- 2
- read
Energy use of Jatropha oil extraction wastes: Pellets from biochar and Jatropha shell blends
V. Ramírez, J. Martí-Herrero, M. Romero, D. Rivadeneira

Journal of Cleaner Production (2019). Vol. 215, pp. 1095-1102

Abstract
The aim of this paper is to evaluate the use of biochar (SCB) produced through the pyrolysis of Jatropha seed cake in combination with Jatropha shell waste (JSh) to form pellet combustion fuel. Mixtures using different proportions of SCB and JSh (0, 25, 50, 75 and 100% each), were prepared, with each proportion tested using two particle sizes and two different quantities of water added. Combustion characteristics and mechanical durability were analyzed to identify the best blend, for which a thermogravimetric analysis was made. Pellets composed of 50% JSh and 50% SCB, with 25% of additional water and 4 mm particle size, presented the greatest mechanical durability (96.83%), and higher heating value (22.14 MJ kg−1). Combustion of pellets made with these characteristics in an industrial burner showed average temperatures ranging from 300 to 350 °C, with a maximum temperature of 460 °C. The production cost of the selected pellet is estimated at 0.11 USD kg−1, making the energy cost of the pellet [0.005 USD MJ−1] approximately 28% of the energy cost of LPG, demonstrating the commercial viability of the selected pellet as feasible.
Abstract
The aim of this paper is to evaluate the use of biochar (SCB) produced through the pyrolysis of Jatropha seed cake in combination with Jatropha shell waste [...]
- 7
- read
Future oil extraction in Ecuador using a Hubbert approach
V. Sebastián, J. Fontalvo, J. Martí-Herrero, P. Ramírez, I. Capellán-Pérez

Energy (2019). Vol. 182, pp. 520-534

Abstract
Hubbert based models to project future oil extraction in Ecuador were developed. Two values of ultimate recoverable resources (URR) (7860–10,700 million barrels (MMbbl)) are applied to 16 models, considering symmetric and asymmetric Hubbert models and one and two cycles under top-down and bottom-up approaches. Models are discussed based on the best fit to historical data, and year and value of maximum extraction. The peak oil extraction obtained ranges between 196 and 215 MMbbl and would be reached in the years 2014–2025. An analysis of the implications of extraction models in a Business as Usual and Alternative oil demand scenarios up to 2035 was performed. Ecuador could become a net oil importer between 2024 and 2035, depending on the model and demand scenario. Economic oil trade balance could be seriously affected, decreasing from a current positive value of around 2 billion USD to incur deficits of 0.6–16.7 billion USD in 2035. Current and future oil dependence for Ecuador would increase vulnerability and compromise the country in terms of energy security and trade balance. It is critical for Ecuador to consider more ambitious policies focused on energy efficiency, renewables and diversification of the productive structure over the next few years.
Abstract
Hubbert based models to project future oil extraction in Ecuador were developed. Two values of ultimate recoverable resources (URR) (7860–10,700 million barrels (MMbbl)) [...]
- 17
- read
Management and design of biogas digesters: A non-calibrated heat transfer model
S. Vilms, J. Martí-Herrero, A. Singh, S. Sommer, S. Hafner

Bioresource Technology (2020). Vol. 296, p. 122264

Abstract
A thermal balance modeling framework is developed, based on heat transfer-resistance networks. The heat transfer model accounts for effects of digester- design, location and operation, including effects of solar irradiance, external heating and ambient climate. We demonstrate extendibility of the framework by using the model in dynamic simulations of substrate temperature for digesters comprising two very different designs. Digester designs modeled include fixed-dome, buried, uninsulated and unheated household digesters in Hanoi, Vietnam, and an industrial-scale anaerobic digester located at a wastewater treatment plant in Esbjerg, Denmark. The modeled temperature profiles were evaluated against measured substrate temperatures over long periods, from 7 months and up. For the two Hanoi digesters, root-mean-square-error were 1.43 °C and 0.92 °C, with Nash-Sutcliffe model efficiency coefficients (NS-C) of 0.87 and 0.93 respectively. For the industrial digester in Esbjerg root-mean-square-error was 0.48 °C with an NS-C of 0.94. The model was not calibrated prior simulation, suggesting good predictive performance.
Abstract
A thermal balance modeling framework is developed, based on heat transfer-resistance networks. The heat transfer model accounts for effects of digester- design, location and [...]
- 113
- read
Performance Analysis of a Small-Scale Biogas-Based Trigeneration Plant: An Absorption Refrigeration System Integrated to an Externally Fired Microturbine
J. Villarroel-Schneider, A. Malmquist, A. Araoz, J. Martí-Herrero, A. Martin

Energies (2019). Vol. 12 (20), pp. 3830

Abstract
Trigeneration or combined cooling, heat and power (CCHP) systems fueled by raw biogas can be an interesting alternative for supplying electricity and thermal services in remote rural areas where biogas can be produced without requiring sophisticated equipment. In this sense, this study considers a performance analysis of a novel small-scale CCHP system where a biogas-fired, 5 kWel externally fired microturbine (EFMT), an absorption refrigeration system (ARS) and heat exchangers are integrated for supplying electricity, refrigeration and hot water demanded by Bolivian small dairy farms. The CCHP solution presents two cases, current and nominal states, in which experimental and design data of the EFMT performance were considered, respectively. The primary energy/exergy rate was used as a performance indicator. The proposed cases show better energy performances than those of reference fossil fuel-based energy solutions (where energy services are produced separately) allowing savings in primary energy utilization of up to 31%. Furthermore, improvements in electric efficiency of the EFMT and coefficient of performance (COP) of the ARS, identified as key variables of the system, allow primary energy savings of up to 37%. However, to achieve these values in real conditions, more research and development of the technologies involved is required, especially for the EFMT
Abstract
Trigeneration or combined cooling, heat and power (CCHP) systems fueled by raw biogas can be an interesting alternative for supplying electricity and thermal services in remote [...]
- 13
- read
Biogas based polygeneration plant options utilizing dairy farms waste: A Bolivian case
J. Villarroel-Schneider, B. Mainali, J. Martí-Herrero, A. Malmquist, A. Martin, L. Alejo

Sustainable Energy Technologies and Assessments (2020). Vol. 37, 100571

Abstract
This study presents a comparative techno-economic feasibility analysis for two polygeneration plant solutions, applied to low-income dairy farms in Bolivia. The first option considers an internally fired microturbine (IFMT) and, the second, an internal combustion engine (ICE). They are integrated with an absorption refrigeration system and a fertilizer dryer. Biogas, produced with farms waste, fuels these power generators. The levelized costs of biogas for cooking, electricity, cooling and fertilizers were determined. The cost of biogas, for both options, was found to be 0.020 USD/kWh, which is lower than the subsidized price of LPG. The most competitive cost of electricity was determined for the ICE plant option; it was found to be 0.082 USD/kWh and is lower than the subsidized cost of fossil fuel-based electricity. The cost of cooling was found to be around 0.082 USD/kWh, which is slightly higher than the cost of cooling supplied by using grid electricity. In a realistic scenario, the shorter payback period was found to be 4.4 years for the ICE plant option. From this, the ICE-based plant was found as the most feasible option. Additionally, if no subsidies are applied to the fossil fuel-based services, the proposed polygeneration systems are a highly competitive alternative.
Abstract
This study presents a comparative techno-economic feasibility analysis for two polygeneration plant solutions, applied to low-income dairy farms in Bolivia. The first option [...]
- 11
- read
Techno-Economic Study of a Biogas-Based Polygeneration Plant for Small Dairy Farms in Central Bolivia
J. Martí-Herrero, A. Malmquist, A. Martin, J. Villarroel-Schneider, B. Mainali, L. Alejo

(2019). Renewable Energy and Sustainable Buildings. Innovative Renewable Energy. Springer, Cham, pp. 675-687

Abstract
This paper presents a techno-economic feasibility study of a polygeneration plant proposed for an association of dairy farmers in Bolivia. The systems have been designed in an integrated approach and are based on the resource recovery utilizing available cow dung. The production of biogas is proposed using “low cost” tubular digesters. The biogas is used for the production of electricity and heat in a 40 kWel internal combustion engine. Heat from the exhaust gas is recovered to drive an absorption cooling system for milk refrigeration and for a bio-slurry drying system. The final services are biogas, electricity, cooling and fertilizers. The techno-economic analysis focuses on determining the levelized cost of the services and compare them with the subsidized and non-subsidized prices of other competitive services in the market. Sensitivity analyses for the services prices are carried out to see the influence of the feedstock handling cost and the subsidies on investment capital. Once the selling prices of the services are defined, the payback period of investment capital is determined. The results show that the biogas and electricity costs were found lower than the subsidized prices of similar services in the market while the cost of cooling was found slightly higher. From the sensitivity analysis we have that the feedstock handling cost can be increased from 10 to 18 USD/ton while maintaining the biogas price below the subsidized LPG price, and applying subsidies to the investment capital allows reducing the rest of the services costs. The payback period is around 3.5 years when a subsidy of 15% is applied in the investment capital and the services are sold at the subsidized prices of the conventional services. Our study has shown that under the conditions of the Bolivian market it is feasible the installation of a polygeneration plant. On the other hand, a market without or with reduced subsidies for the conventional services makes the proposed system more feasible and competitive.
Abstract
This paper presents a techno-economic feasibility study of a polygeneration plant proposed for an association of dairy farmers in Bolivia. The systems have been designed in [...]
- 17
- read
Hosting solutions and other services for congresses
Abstract
Hosting solutions and other services for congresses
- 193
- read
Combinación de métodos lagrangianos de elementos finitos y partículas (PFEM) con métodos eulerianos para el análisis del comportamiento de presas de escollera durante un sobrevertido
A. LARESE, R. Rossi

(2009). "Congreso de Métodos Numéricos en Ingeniería 2009". Barcelona: Sociedad Española de Métodos Numéricos en Ingeniería (SEMNI), p. 466. ISBN9788496736665

Abstract
El objetivo de este trabajo es desarrollar una herramienta de aplicabilidad general que ayude en el estudio de la estabilidad y seguridad de una presa de escollera en una avenida extraordinaria que origina un sobrevertido. se combinan métodos eulerianos para el cálculo de flujo de agua fuera y dentro de la presa con métodos lagrangianos de elementos finitos y partículas para analizar el comportamiento tenso-deformacional de la presa y su posible rotura y eventual desmoronamiento bajo la acción del agua. Finalmente se presenta un método para la proyección de datos entre mallas no coincidentes que permita el acomplamiento fluido-estructura.
Abstract
El objetivo de este trabajo es desarrollar una herramienta de aplicabilidad general que ayude en el estudio de la estabilidad y seguridad de una presa de escollera en una [...]
- 5
- read
Un método de cálculo para el estudio de la seguridad de presas de escollera durante un sobrevertido
A. LARESE, R. Rossi

(2008). "VIII Jornadas Españolas de Presas", p. 1-10. ISBN978-84-380-0406-7

Abstract
El objetivo de este trabajo es desarrollar y validar un nuevo método de cálculo de aplicabilidad general para el estudio de la estabilidad y seguridad de una presa de escollera durante un sobrevertido. Esta nueva herramienta permite analizar con detalle el estado tenso-deformacional de la presa y su posible rotura y eventual desmoronamiento bajo la acción del flujo de agua y combina métodos de elementos finitos avanzados con métodos de partículas. Esta técnica, denominada PFEM (particle finite element method), permite tener en cuenta de las complejas no-linealidades geométricas y mecánicas en problemas de interacción fluido-estructura, incluyendo fenómenos de superficie libre y de erosión. En el trabajo se describen las características esenciales del PFEM y se detalla la aplicación de este método para modelar el sobrevertido en presas de escollera y su eventual desplome. Este tabajo numérico se combina con una parte experimental que se desarrolla en el CEDEX y en la ETS de Ingenieros de Caminos de Madrid en el marco del proyecto XPRES.
Abstract
El objetivo de este trabajo es desarrollar y validar un nuevo método de cálculo de aplicabilidad general para el estudio de la estabilidad y seguridad de una [...]
- 1
- read
Análisis numérico del proceso de filtración en presas de escollera
C. Tapia, E. Oñate, A. LARESE

(2009). Research Report, IT-591

Abstract
Las presas de materiales sueltos funcionan como presas de gravedad, el peso propio de los materiales que la componen aportan la resistencia necesaria para soportar el empuje horizontal y la flexión que producen las aguas en el trasdós. El inconveniente principal de éstas presas está en el sobrevertido. Si se produce una avenida y el agua vierte sobre la presa, ésta se comporta de manera inestable y puede colapsarse y desmoronarse. Se sabe que este fenómeno del sobrevertido es una de las principales razones por las que colapsan este tipo de presas, pero se trata de un fenómeno poco conocido y poco desarrollado. Por ello existe un interés cada vez mayor en el estudio del comportamiento de las presas de materiales sueltos cuando se lleva a cabo el llenado, y el sobrevertido de agua sobre el talud aguas abajo. Este estudio está limitado por la ausencia de un método de cálculo, por lo que es esencial llegar a constituir un método numérico que modelice el comportamiento físico de estos procesos en dichas presas.
Abstract
Las presas de materiales sueltos funcionan como presas de gravedad, el peso propio de los materiales que la componen aportan la resistencia necesaria para soportar el empuje [...]
- 8
- read
A mixed Lagrangian and Eulerian FSI approach for simulation of overtopping on embankment dams
A. LARESE, R. Rossi

(2009). "III International Conference on Computational Methods for Coupled Problems in Science and Engineering". Ischia: CIMNE (Centro Internacional de Métodos Numéricos en Ingeniería), p. 1-4.

Abstract
The present work is concerned with some recent advances in FSI simulations with particular emphasis on the simulation of rockfill dams during extreme phenomena. The final objective is the study of the consequences of an overtopping on the stability of such structure and its possible failure mode. A mixed Lagrangian and Eulerian approach is used. The fluid behavior is described using a modified form of the Navier Stokes equations in order to consider the effect of a variable porosity. A non linear Darcy law is included in the momentum equation. A level set function is chosen to follow the movement of the free surface inside and outside the porous medium. The structure is described using a purely lagrangian PFEM formulation. The specific features of PFEM make it appropriate to treat the rockfill material and its large deformations and shape changes. A projection technique allows to perform the data transfer between the fluid and the structure non matching meshes
Abstract
The present work is concerned with some recent advances in FSI simulations with particular emphasis on the simulation of rockfill dams during extreme phenomena. The final [...]
- 8
- read
Particle finite element method and level set method for the simulation of the failure of rockfill dams due to overtopping phenomena
A. LARESE, R. Rossi, S. Idelsohn

(2010). 8th. World Congress on Computational Mechanics (WCCM8) 5th European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 2008)
- 8
- read
On the Validation of The Particle Finite Element Method (PFEM) for complex engineering fluid flow problems
A. LARESE, R. Rossi, S. Idelsohn

(2010). European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds)

Abstract
Several comparisons between experiments and computational models are presented in the following pages. The objective is to verify the ability of Particle Finite Elements Methods (PFEM) to reproduce hydraulic phenomena involving large deformation of the fluid domain.
Abstract
Several comparisons between experiments and computational models are presented in the following pages. The objective is to verify the ability of Particle Finite Elements Methods [...]
- 10
- read
Metodología de análisis del comportamiento de presas de escollera frente a un sobrevertido
A. LARESE, R. Rossi, M. Toledo

(2010). Actas de IX Jornadas Españolas de Presas

Abstract
Con el objetivo de llegar a un mejor conocimiento del principio del proceso de rotura de las presas de material suelto y obtener una herramienta que permita un análisis más riguroso de este fenómeno, se está desarrollando el proyecto XPRES financiado por el Ministerio de Educación y Ciencia [XPRES 2007-2010]. El proyecto consta de una parte experimental, llevada a cabo por CEDEX y UPM conjuntamente, y de una parte numérica, desarrollada en el Centro Internacional de Métodos Numéricos en Ingeniería (CIMNE). En el presente trabajo se quiere presentar lo que CIMNE está realizando para crear una herramienta de cálculo para el análisis del comportamiento de una presa de escollera durante un sobrevertido, que, por su incidencia en la seguridad de la presa y su entorno, merece un tratamiento más preciso. La nueva metodología permite analizar con detalle el estado de tensión y deformación de la presa y su posible rotura parcial, o catastrófica, bajo la acción del flujo de agua y combina métodos de elementos finitos avanzados con métodos de partículas.
Abstract
Con el objetivo de llegar a un mejor conocimiento del principio del proceso de rotura de las presas de material suelto y obtener una herramienta que permita un análisis [...]
- 5
- read
Physical and numerical modelization of the behaviour of rockfill dams during overtopping scenarios
A. Larese, R. Rossi, M. Toledo

"Dam maintenance and rehabilitation II: proceedings of the 2nd international congress on dam maintenance and rehabilitation, Zaragoza, Spain, 23-25 November 2010". Boca Raton, FL: CRC Press. Taylor & Francis Group, p. 479-487

Abstract
In recent years the technology on embankments dams has developed sensibly due to the advances in soil mechanic, and in all related sciences. Nevertheless their vulnerability to overtopping is still an open theme of discussion. Such an extreme phenomenon can lead to catastrophic consequences with loss of stability and failure of the entire structure. In spite of that, the consequences of an overtopping are still difficult to analyze and to prevent optimally. The possibility to define a numerical instrument to provide support for analyzing the failure of a dam is an important step ahead for organizing the intervention measures and for optimizing the economic plan. This is the principal reason why the Spanish ministry of Science and Innovation is funding the XPRES project, between the Polytechnic University of Madrid (UPM), the Centre for Hydrographical Studies of CEDEX and the International Centre for Numerical Methods in Engineering (CIMNE). The characterization of the failure process is dealt with by combining new numerical models, developed by CIMNE, with experimental tests on physical models of different sizes and scales carried out by UPM and CEDEX. In the present paper the work carried out by CIMNE is presented. The numerical model is based on combination of Lagrangian Particle Finite Element Method (PFEM) and Eulerian level set techniques in order to optimize the coupling between the dynamic effect of water inside and outside the dam, with the structural deformation and collapse induced on the down stream shoulder. The breaching formation and the consequent beginning of the failure are analyzed with particular interest. The specific features of the adopted method make it appropriate to treat the rockfill material and its large deformations and shape changes. A variable viscosity law and an energetic erosion criterion are used to take into account mass sliding and erosion acting in combined or alternative way
Abstract
In recent years the technology on embankments dams has developed sensibly due to the advances in soil mechanic, and in all related sciences. Nevertheless their vulnerability [...]
- 26
- read
Estudios sobre el comportamiento de presas de escollera ante vertidos por coronación
C. Lechuga, A. Lara, M. Berga, P. Viña, M. Toledo, R. Morán Moya, H. Campos, J. García-Palacios, A. Larese

(2010). Actas de IX Jornadas Españolas de Presas

Abstract
Como consecuencia del desarrollo de un Proyecto de Investigación, el CEDEX, la Universidad Politécnica de Madrid y el CIMNE han trabajado conjuntamente en el desarrollo de un método para el estudio del proceso de rotura de presas de escollera por sobrevertido combinando técnicas de elementos finitos y partículas. Para ello se han realizado y analizado 80 ensayos en modelo físico en las instalaciones del CEDEX y UPM, utilizando seis tamaños de material, cuatro taludes distintos (1,5H:1V; 2,2H:1V; 3,0H:1V y 3,5H:1V) y tres opciones de impermeabilización (sin pantalla, núcleo interno y pantalla externa). Estos ensayos están permitiendo conocer mejor el fenómeno y calibrar el modelo matemático. El modelo matemático basa su desarrollo en el acoplamiento de dos fenómenos: la filtración del agua en un medio granular (simulada mediante una malla euleriana) y la erosión del propio material granular (simulada mediante una malla lagrangiana). Durante 2010 se prevé finalizar la construcción, en el CEDEX, de un canal de grandes dimensiones que permitirá ensayar presas de hasta 4 m de altura en condiciones de laboratorio, que servirán como prototipo para el análisis del efecto de escala. Este prototipo aumentará considerablemente la fiabilidad de las conclusiones obtenidas.
Abstract
Como consecuencia del desarrollo de un Proyecto de Investigación, el CEDEX, la Universidad Politécnica de Madrid y el CIMNE han trabajado conjuntamente en el [...]
- 7
- read
PFEM application in fluid structure interaction problems
M. Celigueta, A. LARESE, S. Latorre

4th CONFERENCE ON ADVANTAGES AND APPLICATIONS OF GID. The Personal Pre and Post Processor

Abstract
In the current paper the Particle Finite Element Method (PFEM), an inno-vative numerical method for solving a wide spectrum of problems involving the interaction of ﬂuid and structures, is brieﬂy presented. Many examples of the use of the PFEM with GiD support are shown. GiD framework provides a useful pre and post processor for the speciﬁc features of the method. Its advantages and shortcomings are pointed out in the present work
Abstract
In the current paper the Particle Finite Element Method (PFEM), an inno-vative numerical method for solving a wide spectrum of problems involving the interaction of ﬂuid [...]
- 8
- read
Analysis of stability of earth dams in overtopping scenarios with the particle finite element method
A. LARESE, R. Rossi

(2009). International Conference on Particle-Based Methods Particles, E. Oñate and D.R.J. Owen (Eds), Barcelona

Abstract
The aim of this work is to study the consequences of an overtopping on the stability of earth dams and its possible failure mode. Rehabilitation and safety analysis of existing dams is nowadays an open field of research. The possibility to define a numerical instrument to provide support for analyzing the failure of a dam is a big step ahead for organizing the intervention measures and for optimizing the economic plan. In fact many existing dams have now to be modified due to the revision of previous design criteria in order to increase their safety. The objective of our work is to develop and validate a new computational method of general applicability that allows treating the above problems. The method will combine advanced finite element and particle techniques. Fluid-structure interaction effects and non-linear geometrical and mechanical effects in the dam material are considered. A mixed Lagrangian and Eulerian approach is used. The fluid behavior is described using a modified form of the Navier Stokes equations in order to consider the effect of a variable porosity. A non linear Darcy law is included in the momentum equation. A level set function is chosen to follow the movement of the free surface inside and outside the porous medium. The structure is described using a purely lagrangian PFEM formulation [1] and [2]. The specific features of PFEM make it appropriate to treat the rockfill material and its large deformations and shape changes. A projection technique allows to perform the data transfer between the fluid and the structure non matching meshes. Peer Reviewed Postprint (published version)
Abstract
The aim of this work is to study the consequences of an overtopping on the stability of earth dams and its possible failure mode. Rehabilitation and safety analysis of existing [...]
- 11
- read
Coupling eulerian and lagrangian models to simulate seepage and evolution of failure in prototype rockfill dams
A. LARESE, R. Rossi

(2011). ICOLD Benchmark Workshop on Numerical Analysis of Dams. "XI ICOLD Benchmark Workshop on Numerical Analysis of Dams". València: 2011, p. 1-20.

Abstract
This work presents a model which allows considering the interaction between rockfill dams and the water. The fluid behaviour is analyzed on an Eulerian fixed mesh. A level set technique is employed to tack the evolution of the free surface. An edge based approach is used to solve efficiently the modified form of the Navier‐Stokes equations. The dam response is evaluated on a Lagrangian moving mesh using PFEM. A visco‐rigid constitutive model is used to describe the structural behaviour and collapse of rockfill under fluid dynamic forces. Mohr‐Coulomb is the failure criteria adopted in the calculation. Finally a nodal projection algorithm allows the transferring of information through non matching meshes.
Abstract
This work presents a model which allows considering the interaction between rockfill dams and the water. The fluid behaviour is analyzed on an Eulerian fixed mesh. A level [...]
- 2
- read
Safety analysis of rockfill dams during overspill phenomena
A. LARESE, R. Rossi, S. Idelsohn

(2012). Annual Meeting of ICOLD. "76th Annual Meeting of ICOLD". Sofia: Bulgarian Committee on Large Dams.

Abstract
An increasing interest is rising on the study of rockfill embankment dams during overspilling. Design criteria of earth dams are being reviewed in many countries to guarantee an increasing safety level in front of an exceptional flooding. The possibility of studying the failure process is currently limited by the lack of a precise knowledge and by the absence of a suitable computational method which can predict accurately the onset and evolution of breaching. The analysis of the possible consequences of an accidental overspill is then impossible or very imprecise and the necessary economical measures for solving the problem are then inefficient. An appropriate computational method will help to reduce the economic impact of the investments in dam safety and in emergency plans for embankment dams. The dam failure is mainly due to the hydrodynamic forces acting on the rockfill material in the downstream face. Two failure mechanisms occur in an alternative or combined mode: dragging of rockfill particles with erosive effects, and loss of stability of a part of the downstream region due to the slippage of a bulk mass of the rockfill material. The loss of a relevant part of the rockfill material in the downstream side of the dam, can lead to mechanical failure of the impervious region, either if this is formed by a concrete (or asphalt) face, or by a clay core. In the latter case, erosion of the cohesive material can also influence the failure process, but the scale time of appearance of this erosion phenomenon is so different from the principal erosion of the rockfill part, that this second one is not taken into account. The objectives of the current work is to develop and validate a new computational method of general applicability that allows treating above problems. The method will combine advanced finite element and particle techniques. Fluid-structure interaction effects and non-linear geometrical and mechanical effects in the dam material will be considered. The specific features of the analysis method proposed make it appropriate to treat the rockfill material adjacent to the downstream surface where a transfer of momentum between the water running over the dam surface and that moving through the interior due to seepage occurs
Abstract
An increasing interest is rising on the study of rockfill embankment dams during overspilling. Design criteria of earth dams are being reviewed in many countries to guarantee [...]
- 15
- read
Theme B: simulation of the behavior of prototypes of rockfill dams during overtopping scenarios: seepage evolution and beginning of failure
A. LARESE, R. Rossi, M. Toledo, R. Morán, H. Campos, A. Lara

(2011). "XI ICOLD Benchmark Workshop on Numerical Analysis of Dams". València: 2011, p. 1-25.

Abstract
In recent years the technology on embankments dams has developed sensibly due to the advances in soil mechanic, and in all related sciences. Nevertheless their vulnerability to overtopping still remains their weakest point in comparison with concrete structures. The principal aim of this theme is the simulation of the initial stage of failure of the dam when an overtopping or an exceptional flood occurs. A fluid-structure coupled problem has to be considered. The sudden variation of the upstream conditions induces a quick evolution of a seepage line in the downstream shoulder. Non linear Darcy law has to be taken into account. On the other hand the water, emerging from the toe of the dam, induces dragging of particles and possible mass sliding, depending on geometrical and material conditions. UPM and CEDEX have carried out more than 70 experiments during the last two years. They analyzed the influence of a series of parameters on the failure mechanisms. These parameters are, for instance, the dimension of the rocks, the slope of the downstream part of the dam, the type of impermeable element used, and so on. We propose to reproduce numerically three experiments analyzing the evolution of seepage and following initial stage of failure in the case of a homogeneous dam, of a dam with impermeable upstream face and of a dam with internal core. Experimental data, bottom pressure distribution and topographic analysis of the geometry of the dam during the failure can be compared with numerical results.
Abstract
In recent years the technology on embankments dams has developed sensibly due to the advances in soil mechanic, and in all related sciences. Nevertheless their vulnerability [...]
- 22
- read
Simulation of the beginning of failure in rockfill dams caused by overtopping
A. LARESE, R. Rossi

(2015). Dam Protections against Overtopping and Accidental Leakage, Ed: CRC Press, 111 p.

Abstract
This paper presents the work carried out during the last years for the creation of a numerical model to simulate transient seepage phenomena in rockfill dams when overtopping occurs. The objective is to analyze the trigger of failure mechanism in the downstream slope due to hydrodynamic forces. The current work aims to give a contribution to this field proposing a new computational method of general applicability for simulating, with a unique formulation, the flow throughout and over the dam while failure occurs together with the dam structural response. The numerical results are validated with experiments on scale models rockfill dams.
Abstract
This paper presents the work carried out during the last years for the creation of a numerical model to simulate transient seepage phenomena in rockfill dams when overtopping [...]
- 5
- read
Simulation of two and three-dimensional viscoplastic flows using adaptive mesh refinement
J. Cotela, R. Rossi, A. LARESE

International Journal for Numerical Methods in Engineering (2017). Vol. 112 (11), pp. 1636-1658

Abstract
This paper presents a finite element solver for the simulation of steady non‐Newtonian flow problems, using a regularized Bingham model, with adaptive mesh refinement capabilities. The solver is based on a stabilized formulation derived from the variational multiscale framework. This choice allows the introduction of an a posteriori error indicator based on the small scale part of the solution, which is used to drive a mesh refinement procedure based on element subdivision. This approach applied to the solution of a series of benchmark examples, which allow us to validate the formulation and assess its capabilities to model 2D and 3D non‐Newtonian flows.
Abstract
This paper presents a finite element solver for the simulation of steady non‐Newtonian flow problems, using a regularized Bingham model, with adaptive mesh refinement capabilities. [...]
- 5
- read
Comparison of a Material Point Method and a Galerkin Meshfree Method for the Simulation of Cohesive-Frictional Materials
I. Iaconeta, A. LARESE, R. Rossi, Z. Guo

Materials (2017). Vol. 10 (10), pp. 1150

Abstract
The simulation of large deformation problems, involving complex history-dependent constitutive laws, is of paramount importance in several engineering fields. Particular attention has to be paid to the choice of a suitable numerical technique such that reliable results can be obtained. In this paper, a Material Point Method (MPM) and a Galerkin Meshfree Method (GMM) are presented and verified against classical benchmarks in solid mechanics. The aim is to demonstrate the good behavior of the methods in the simulation of cohesive-frictional materials, both in static and dynamic regimes and in problems dealing with large deformations. The vast majority of MPM techniques in the literatrue are based on some sort of explicit time integration. The techniques proposed in the current work, on the contrary, are based on implicit approaches, which can also be easily adapted to the simulation of static cases. The two methods are presented so as to highlight the similarities to rather than the differences from "standard" Updated Lagrangian (UL) approaches commonly employed by the Finite Elements (FE) community. Although both methods are able to give a good prediction, it is observed that, under very large deformation of the medium, GMM lacks robustness due to its meshfree natrue, which makes the definition of the meshless shape functions more difficult and expensive than in MPM. On the other hand, the mesh-based MPM is demonstrated to be more robust and reliable for extremely large deformation cases.
Abstract
The simulation of large deformation problems, involving complex history-dependent constitutive laws, is of paramount importance in several engineering fields. Particular attention [...]
- 14
- read
A Finite Element Implementation for the Three Variable Biot Model and an Introduction to its Application to Sugar Cane Milling Processes
R. Rossi, J. Amaya, A. Larese, M. Pohl

IEEE 38th Central America and Panama Convention (CONCAPAN XXXVIII)

Abstract
This document presents a step-by-step derivation for the formulation of a Biot model for saturated porous media, this formulation is apt to be implemented both in two and three dimensions. Afterwards, the spatial and temporal discretization for the model is also presented, focusing on a symbolic representation around the Gauss points for a particular element in an implementation. Lastly, the results for a couple of test cases are shown, and of special interest, the test case that constitutes the first approximation towards an implementation of this model for the simulation of sugar cane milling processes.
Abstract
This document presents a step-by-step derivation for the formulation of a Biot model for saturated porous media, this formulation is apt to be implemented both in two and [...]
- 11
- read
Soil-structure interaction simulation of landslides impacting a structure using an implicit material point method
B. Chandra, A. Larese, I. Iaconeta, R. Rossi, R. Wüchner

2nd International Conference on the Material Point Method for Modelling Soil-Water-Structure Interaction

Abstract
Global warming and climate changes have been one of the many causes that triggered numerous catastrophic landslide and mudflow disasters in the past twenty years. The increase of earth temperature has contributed to the increase of precipitation and undisputedly affected the soil slope stability, which by further may cause landslides in various scale and speed. This large soil deformation phenomenon carries along huge rocks and heavy materials that often results in extensive damage in civil infrastructures both directly or indirectly. Coming with this motivation, a soil-structure interaction simulation based on the implicit material point method (MPM) has been implemented within the Kratos Multiphysics framework for the objective of predicting structural deformation and, furthermore, structural failure caused by environmental flow problems such as landslides. In the current study, the soil is modeled using a non-associated Mohr-Coulomb-based elastoplastic law, while the structure is modeled as elastic and Neo-Hookean hyperelastic materials. In the numerical tests conducted, the equivalent stress and displacement measured on both rigid and flexible structures show a good qualitative agreement. In the future works, a more adequate consideration of the soil and structural model will be investigated before conducting a real-scale landslide simulation.
Abstract
Global warming and climate changes have been one of the many causes that triggered numerous catastrophic landslide and mudflow disasters in the past twenty years. The increase [...]
- 41
- read
Corrigendum to “Modelling of Bingham and Herschel-Bulkley flows with mixed P1/P1 finite elements stabilized with Orthogonal Subgrid Scale” [J. Nonnewt Fluid Mech, 228, 2016, p. 1–16]
E. Moreno, A. LARESE, M. Cervera

Journal of Non-Newtonian Fluid Mechanics (2019). Vol. 267, pp. 100

Abstract
This corrigendum provides some missing data and a reference for the paper Modelling of Bingham and Herschel-Bulkley flows with mixed P1/P1 finite elements stabilized with Orthogonal Subgrid Scale, as well as it corrects a mistaken unit of measure in the caption of a figure. The authors apologise for the inconveniences caused. Corrigendum – The caption of Figure 14 should read: Cylinder in an Herschel-Bulkley fluid. Drag force coefficient in terms of the relation L:R for Bn∗ = 10. – The flow index n in the sub-captions of Figure 17(a)-17(l) is nondimensional. – The following Reference should be added: [71] E. Moreno, M. Cervera "Elementos finitos mixtos estabilizados para flujos confinados de Bingham y de Herschel-Bulkley Parte II: soluciones numéricas" Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 2016; 32(3); 131-138. – The captions of Figures 2, 4-7, 22-25 should include the following text at the end: (Reproduced with permission from [71]).
Abstract
This corrigendum provides some missing data and a reference for the paper Modelling of Bingham and Herschel-Bulkley flows with mixed P1/P1 finite elements stabilized with [...]
- 11
- read
Coupled Soil-Structure Interaction Modeling and Simulation of Landslide Protective Structures
B. Chandra, A. LARESE, P. Bucher, R. Wüchner

(2019). VIII International Conferenceon Computational Methods for Coupled Problems in Science and Engineering

Abstract
Within the past two decades, mass movements hazards involving fast and large soil deformation have increased significantly in frequency and magnitude due to their strong relation to climate changes and global warming. These phenomena often bring along rocks, debris, and heavy materials that can extensively damage and destroy the landscape and infrastructures, causing devastating economic loss, and often, human casualties. The risk of future disasters continues to escalate with the increase of real estate development in suburban areas, including mountainous regions. Further assessment and prediction on such disasters and their countermeasures are, therefore, in high economic demands. One of the most intuitive ways is to install protective structures in mountain slopes and valleys that can hold the materials brought by the moving landslides. While the current state of the art of landslide prediction using numerical methods has been mainly dominated by the development of advanced geomechanical models suited for different types of soil materials, e.g. multi-phase unsaturated soil model, this study focuses more on the interaction of such phenomena with the installed protective structures. Here, an implicit formulation of material point method (MPM) is implemented to model the landslides considering finite strain assumption. Furthermore, a staggered coupling scheme with traditional Finite Element Method (FEM) is proposed to simulate accurately and robustly the dynamic force and displacement coupling of soil-structure interaction (SSI). All developments of the method are implemented within the Kratos-Multiphysics framework and available under the BSD license (https://github.com/KratosMultiphysics/Kratos/wiki). In the future works, more adequate consideration of coupling scheme and material models considering damage and fracture will be investigated before conducting a real-scale landslide simulation.
Abstract
Within the past two decades, mass movements hazards involving fast and large soil deformation have increased significantly in frequency and magnitude due to their strong relation [...]
- 19
- read
A procedure to design toe protections for rockfill dams against extreme through-flows
R. Morán Moya, M. Toledo, A. LARESE, R. Monteiro-Alves

Engineering Structures (2019). Vol. 195, pp. 400-412

Abstract
This article presents a procedure for designing rockfill toe protections with a highly permeable downstream shoulder to avoid mass sliding in dams during extremely high through-flow episodes. This accidental through-flow may be caused by reasons such as overtopping or leakage flow due to the loss of sealing at the impervious element of the dam or its foundation. The proposed protection is located at the downstream toe of the dam, and it is composed of highly permeable material, typically compacted rockfill. This material can be different from the rockfill forming the downstream shoulder of the main dam. The work is a result of a combination of numerical and experimental research carried out to analyze the influence of the geometry of the toe protection for given material properties of both dam shoulder and protection, when they can be considered highly isotropic. As a result of these studies, a design methodology is obtained. A series of validation tests are presented to support the reliability of the method.
Abstract
This article presents a procedure for designing rockfill toe protections with a highly permeable downstream shoulder to avoid mass sliding in dams during extremely high through-flow [...]
- 15
- read
Estudio de la capacidad de desagüe de aliviaderos con compuertas mediante modelación numérica. Aplicación a la presa de Oliana
F. Salazar, R. Rossi, R. Morán Moya, A. LARESE

(2011). Resúmenes y Comunicaciones JIA 2011

Abstract
En el presente trabajo se estudia la capacidad de desagüe de aliviaderos con compuertas radiales mediante la aplicación de un código de cálculo numérico al caso del aliviadero de la presa de Oliana.
Abstract
En el presente trabajo se estudia la capacidad de desagüe de aliviaderos con compuertas radiales mediante la aplicación de un código de cálculo numérico [...]
- 4
- read
Numerical modeling of the hydraulic performance of Oliana dam spillway using Kratos
F. Salazar, R. Morán Moya, R. Rossi, A. LARESE

(2011). "XI ICOLD Benchmark Workshop on Numerical Analysis of Dams". València, p. 1-6.

Abstract
Head-discharge curves in gated spillways are commonly obtained by means of the application of empirical formulations. The most popular one, recommended by SPANCOLD [1], was developed by the U.S. Army Corps of Engineers [2]. It is based in several experimental campaigns whose results were used to define a chart which provides the value of the discharge coefficient depending on the gate position. This formula does not take into account the effect of side contractions, which depends on the length-height ratio of the bays, as well as the head, the inflow conditions and spillway geometry. Thus, it is necessary to carry out specific experimental tests for every new spillway. The aim of the work is the analysis of the discharge capacity of radial gated spillways. The accuracy of empirical formulations has been studied, and Oliana Dam spillway has been investigated via numerical simulation.
Abstract
Head-discharge curves in gated spillways are commonly obtained by means of the application of empirical formulations. The most popular one, recommended by SPANCOLD [1], was [...]
- 7
- read
Applications of the particle finite element method in dam engineering
F. Salazar, D. Pozo, R. Morán Moya

"II Int. Conference on Particle-based Method - Fundamentals and Applications", pp. 1-11

Abstract
The paper presents the results of the application of the Particle Finite Element Method (PFEM) to the analysis of some of the more complex phenomena related to dam hydraulics: shock waves in spillways, aeration in bottom outlets, and erosion in the downstream river bed. Furthermore, the method has been applied to the study of the consequences of landslides in reservoirs: the wave generation, its propagation, and the affection to the dam.
Abstract
The paper presents the results of the application of the Particle Finite Element Method (PFEM) to the analysis of some of the more complex phenomena related to dam hydraulics: [...]
- 7
- read
Numerical modeling of air-water interaction in bottom outlets. Application to susqueda dam
F. Salazar, R. Morán Moya, M. Toledo

4th IAHR International Symposium on Hydraulic Structures

Abstract
The convenience of installing an aeration system in dam bottom outlets is well known nowadays. Otherwise, damages due to cavitation and vibration are frequently serious, as could be observed in several dams built in the beginning of the 20 th century. Several empirical formulas have been developed for calculating the maximum air demand, but their results are very sensitive to the variation of the geometrical characteristics of the installation. The Ministry of Science and Innovation of Spain has funded a project to address this problem using numerical modeling. The interaction of two fluids (air and water) with very different physical characteristics is the greatest difficulty to be overcome. To deal with this challenge, the Particle Finite Element Method (PFEM) has been used. The method has being applied to verify the performance of the new flat-seat circular-section valves, which design is being developed by INHISA as an alternative to the traditional Bureau-Type. Firstly, a case comprising the mentioned valve in a standard facility was analyzed for validation purposes, as well as for studying the influence of the different parameters in the air flow. Then, the method has been used for the analysis of the performance of the bottom outlets of Susqueda Dam.
Abstract
The convenience of installing an aeration system in dam bottom outlets is well known nowadays. Otherwise, damages due to cavitation and vibration are frequently serious, as [...]
- 7
- read
Protecciones de presas mediante repiés de escollera y bloques en forma de cuña
R. Morán Moya, R. Monteiro, M. Toledo, F. Salazar, A. LARESE, J. Mauro

(2013). Jornada Técnica. Avances en investigación aplicada en seguridad hidráulica de presas.

Abstract
Las conferencias presentadas durante la Jornada Técnica sobre Avances en investigación aplicada en seguridad hidráulica de presas organizada por el Centro de Estudios Hidrográficos del CEDEX y la Universidad Politécnica de Madrid el mes de junio de 2013 versan sobre la mejora de la seguridad hidráulica de las presas, una materia de especial importancia en nuestro país (habida cuenta del elevado número de infraestructuras de regulación en servicio y de su antigüedad media) y constituye una tarea en la que tanto el CEDEX como la Universidad Politécnica de Madrid han trabajado activamente en los últimos años. En este artículo se van a tratar dos tipos de protecciones de carácter muy diferente: las protecciones tipo repié y las protecciones mediante bloques prefabricados en forma de cuña (en adelante, BFC). En el caso de las protecciones tipo repié, aplicables a presas de materiales sueltos con espaldón de aguas abajo de escollera, se realiza una breve descripción del estado del arte de la tecnología y de sus aplicaciones así como de la innovación conseguida en este campo a partir de las investigaciones desarrolladas en la Universidad Politécnica de Madrid (UPM), en colaboración con el Centro de Estudios Hidrográficos del CEDEX y el Centro Internacional de Métodos Numéricos en la Ingeniería (CIMNE), en el marco del proyecto de investigación del Plan Nacional denominado EDAMS, y cuyo resultado directo ha sido la tesis doctoral titulada “Mejora de la seguridad de las presas de escollera frente a percolación accidental mediante protecciones tipo repié” (Morán, 2013). En la parte dedicada a la protección mediante BFC`s se describe la tecnología, sus posibilidades de utilización dentro de la ingeniería de presas y se presenta el caso de la presa de Barriga (Burgos), como primera aplicación mundial de esta tecnología para el diseño de un aliviadero de servicio sobre el propio cuerpo de una presa de materiales sueltos en explotación.
Abstract
Las conferencias presentadas durante la Jornada Técnica sobre Avances en investigación aplicada en seguridad hidráulica de presas organizada por el Centro [...]
- 6
- read
Interpretación de los datos de auscultación de presas por métodos no convencionales
F. Salazar, R. Morán Moya, L. Morera González, A. Roa, D. Sánchez, M. Toledo

(2013). Jornada Técnica. Avances en investigación aplicada en seguridad hidráulica de presas.

Abstract
Las conferencias presentadas durante la Jornada Técnica sobre Avances en investigación aplicada en seguridad hidráulica de presas organizada por el Centro de Estudios Hidrográficos del CEDEX y la Universidad Politécnica de Madrid el mes de junio de 2013 versan sobre la mejora de la seguridad hidráulica de las presas, una materia de especial importancia en nuestro país (habida cuenta del elevado número de infraestructuras de regulación en servicio y de su antigüedad media) y constituye una tarea en la que tanto el CEDEX como la Universidad Politécnica de Madrid han trabajado activamente en los últimos años. Todo profesional familiarizado con las presas conoce la dificultad de analizar su comportamiento, a veces sorprendente, con el grado de aproximación deseado. Esta dificultad emana de la complejidad de todos los elementos que afectan al modo en que la presa y el terreno de apoyo responden a las solicitaciones, y también a la complejidad de las propias solicitaciones. Muchas de las características del conjunto presa-terreno escapan a las posibilidades de modelación mediante los modelos disponibles de base teórica.
Abstract
Las conferencias presentadas durante la Jornada Técnica sobre Avances en investigación aplicada en seguridad hidráulica de presas organizada por el Centro [...]
- 6
- read
Nuevos métodos para el análisis de datos de auscultación de presas: la herramienta iComplex
F. Salazar, J. Piazzese, L. Morera González, M. Toledo, J. Mauro

VIII Congreso Argentino de Presas y Aprovechamientos Hidroeléctricos

Abstract
Los avances en los instrumentos de medida y en las técnicas de transmisión y almacenamiento de información han permitido aumentar el control de la seguridad de las presas, con medidas más fiables, precisas y frecuentes. En general, estos avances no han sido acompañados por una mejora en las técnicas de tratamiento y análisis de los datos de auscultación, que con frecuencia se limitan a gráficos y modelos estadísticos sencillos, como el modelo HST. Los autores están trabajando en la integración de tecnologías avanzadas de análisis de datos para extraer la mayor información posible acerca del estado de seguridad de la presa. Como resultado, se está desarrollando una herramienta que contiene: a) un módulo gráfico flexible y dinámico, que permite explorar los datos de forma interactiva, b) herramientas para el pretratamiento de datos para filtrar valores anómalos, completar series y calcular series modificadas, y c) una librería de herramientas estadísticas y de inteligencia artificial (HST, redes complejas y redes neuronales, entre otras) para seleccionar el modelo predictivo más adecuado en cada caso. Palabras clave: Seguridad de presas, auscultación, análisis de datos, modelos de predicción, bases de datos.
Abstract
Los avances en los instrumentos de medida y en las técnicas de transmisión y almacenamiento de información han permitido aumentar el control de la seguridad [...]
- 21
- read
Procedimiento de diseño de protecciones tipo repié para evitar la rotura por deslizamiento en masa de presas de escollera sometidas a percolaciones extremas
R. Morán Moya, R. Monteiro, M. Toledo, A. LARESE, F. Salazar, R. Rossi

(2013). "III Jornadas de Ingeniería del Agua. La protección contra los riesgos hídricos ( JIA 2013)", 23/10/2013 - 24/10/2013, Valencia (España). ISBN 978-84-267-2071-9. pp. 341-348.

Abstract
Síntesis del procedimiento de diseño formulado en la tesis doctoral de Rafael Morán. La demanda de seguridad en infraestructuras críticas ha aumentado en los últimos años, especialmente en los países con mayor nivel de desarrollo económico y social. En estos países, el control del estado de la seguridad de estas infraestructuras y las medidas para disminuir los riesgos asociados en caso de avería están adquiriendo una gran relevancia, lo que se ha plasmado en unas normativas más exigentes. Las presas están dentro de este tipo de infraestructuras tanto por su interés estratégico como por los daños que pueden llegar a ocasionar en caso de un mal funcionamiento o de su eventual rotura. La tecnología de protecciones de presas surge como respuesta a esta demanda de seguridad y, por este motivo, ha experimentado un gran avance en las últimas décadas. Dentro de este artículo se va a describir la protección tipo repié y el procedimiento de diseño propuesto para evitar el deslizamiento en masa de presas con espaldón de aguas abajo de escollera, cuando se produce una filtración anormalmente elevada debido a un fallo en la impermeabilidad en el conjunto presa-cimiento. El artículo describe brevemente el estado del arte de la tecnología y de sus aplicaciones así como de la innovación conseguida en este campo a partir de las investigaciones desarrolladas en la Universidad Politécnica de Madrid (UPM), en colaboración con el Centro de Estudios Hidrográficos del CEDEX y el Centro Internacional de Métodos Numéricos en la Ingeniería (CIMNE), en el marco del proyecto de investigación del Plan Nacional denominado EDAMS.
Abstract
Síntesis del procedimiento de diseño formulado en la tesis doctoral de Rafael Morán. La demanda de seguridad en infraestructuras críticas ha aumentado [...]
- 46
- read
Discussion on “Thermal displacements of concrete dams: Accounting for water temperature in statistical models”
F. Salazar, M. Toledo

Engineering Structures (2018). Vol. 171, pp. 1071-1072

Abstract
A discussion on the paper by Tatin et al. (2015) is presented. The paper described an innovative statistical model to interpret dam behaviour, which was validated with artificial data and then applied to seven dams in operation. This discussion provides several comments about the model performance evaluation, as well as suggestions for further analysis of the monitoring data.
Abstract
A discussion on the paper by Tatin et al. (2015) is presented. The paper described an innovative statistical model to interpret dam behaviour, which was validated with artificial [...]
- 5
- read
Corrigendum to “An empirical comparison of machine learning techniques for dam behaviour modelling” [Struct. Saf. 56 (2015) 9–17]
F. Salazar, M. Toledo, R. Morán Moya

Structural Safety (2016). Vol. 59, pp. 118

Abstract
The authors wish to report an error in Fig. 1 of the above-mentioned article. The correct position of the monitoring devices is shown below.
- 4
- read
A methodology for dam safety evaluation and anomaly detection based on boosted regression trees
F. Salazar, J. González, M. Toledo

(2016). Proceedings on 8th European Workshop On Structural Health Monitoring (EWSHM 2016)

Abstract
Many countries are implementing new dam safety regulations that often include more restrictive standards. This, together with the increasing average age of dams, results in a greater need for dam control and maintenance works. The advances in information and communications technologies improved the performance of dam monitoring systems, so a large amount of information on the dam behaviour can be collected. This has led to the use of more powerful tools for its analysis, many of which were first developed in the field of machine learning (e. g. neural networks). They offer some advantages over the conventional statistical methods. However, their capacity for early detection of anomalies has seldom been studied. As a result, they are far from being fully accepted by practitioners, whose analyses are often restricted to the interpretation of simple plots of time series data, together with basic statistical models. The present work describes a methodology for anomaly detection in dam behaviour, with the following features: a) The prediction model is based on boosted regression trees (BRTs). b) Causal and auto-regressive models are combined to detect different types of anomalies. c) It is checked whether the values of the external variables fall within the range of the training data. The performance of the proposed methodology was assessed through its application to a test case corresponding to an actual 100-m height arch dam, in operation since 1980. Artificial data were generated by means of a finite element model. Different anomalies were later added in order to test the anomaly detection capability. The method can be applied to other response variables and dam typologies, due to the great flexibility of BRTs, which automatically select the most relevant inputs.
Abstract
Many countries are implementing new dam safety regulations that often include more restrictive standards. This, together with the increasing average age of dams, results in [...]
- 8
- read
Modelación física y numérica de aliviaderos en laberinto con fondo poliédrico
J. Mauro, F. Salazar, M. Toledo, F. Caballero, C. Ponce-Farfán, T. Ramos

Ingeniería del agua (2016). Vol. 20 (3), pp. 127-138

Abstract
Un número importante de presas en España deben ampliar sus aliviaderos para cumplir criterios de seguridad distintos a los de su construcción. En este contexto, cobran interés soluciones de aliviaderos no convencionales como los aliviaderos en laberinto, que permiten aumentar la capacidad de desagüe de estructuras existentes con una inversión moderada. Para la optimización del diseño de este tipo de aliviaderos se desarrolla el proyecto POLILAB, en cuyo marco se llevaron a cabo las campañas de modelación física y numérica que se exponen en este artículo. Se describen los resultados más relevantes en lo relativo a la inclusión de fondos poliédricos en aliviaderos en laberinto como refuerzo estructural, contemplando además las modificaciones que implican con respecto a soluciones en laberinto con fondo plano, en lo relativo a capacidad de desagüe, patrón de flujo y acciones hidrodinámicas sobre los contornos del laberinto.
Abstract
Un número importante de presas en España deben ampliar sus aliviaderos para cumplir criterios de seguridad distintos a los de su construcción. En este [...]
- 6
- read
Hydraulic and stability analysis of the supporting layer of wedge-shaped blocks
F. Salazar, J. Mauro, A. Laresa, M. Toledo

(2016). Protections 2016: 2nd International Seminar on Dam Protection Against Overtopping". Fort Collins, CO: Colorado State University, p. 13-15

Abstract
Wedge shaped blocks (WSB) are attracting increasing attention as protection against overtopping for earth and rock-fill dams. However, there are limited examples of application and some aspects of the technology merit additional research and improvement. One key issue is the design of drainage and supporting layer for WSB protections. During overtopping, part of the overflow leaks through the joints between blocks, hence circulating through the granular material. The permeability and thickness of the supporting layer must be sufficient to prevent the flow from generating pressure on the bottom side of the blocks, which contributes to its destabilization. However, it must also be structurally stable to avoid undesirable deformations on the downstream face. Both the material permeability and the layer thickness determine the hydraulic behavior of this element. These, together with the weight of the blocks and the slope of the downstream face, directly influence mass and block stability. These aspects should be taken into account for the numerical modeling of seepage through the supporting layer. To this end, an application of the open source software Kratos Multi-physics was employed. A parametric study was conducted to quantify the influence of each design variable in the safety factor against mass sliding of the supporting and drainage layer.
Abstract
Wedge shaped blocks (WSB) are attracting increasing attention as protection against overtopping for earth and rock-fill dams. However, there are limited examples of application [...]
- 2
- read
Applications of numerical methods in design and evaluation of overtopping protection systems
F. Salazar, J. Mauro, A. Larese, M. Toledo

(2016). Protections 2016: 2nd International Seminar on Dam Protection Against Overtopping". Fort Collins, CO: Colorado State University, p. 1-10

Abstract
The development of overtopping protection systems often requires detailed analyses of complex physical phenomena. This hinders the comprehensive knowledge of their behavior, and therefore the development of suitable design criteria. In recent years, the authors have developed and validated different methods, combining continuous, particle and discrete numerical techniques, to obtain accurate and reliable solutions of different numerical problems involving fluid-soil-structure interaction. In this contribution, some applications of these methods to the study of dam protection against overtopping are presented. The main advantages of this approach include the ability to extract results of the governing variables (pressure, velocity) at any location of the domain, and the possibility to consider scenarios without the restrictions of the experimental facilities (flow rates, size, scale effects). In particular, the contribution gathers examples of application of numerical methods in a) analysis of rockfill dam stability against overtopping, including seepage evolution and deformation of the downstream shoulder, and b) stability analysis of wedge-shaped-blocks subjected to vandalism.
Abstract
The development of overtopping protection systems often requires detailed analyses of complex physical phenomena. This hinders the comprehensive knowledge of their behavior, [...]
- 2
- read