Search

This work presents a partitioned method for landslide-generated wave events. The proposed strategy combines a Lagrangian Navier Stokes multi-fluid solver with an Eulerian method based on the Boussinesq shallow water equations. The Lagrangian solver uses the Particle Finite Element Method to model the landslide runout, its impact against the water body and the consequent wave generation. The results of this fully-resolved analysis are stored at selected interfaces and then used as input for the shallow water solver to model the far-field wave propagation. This one-way coupling scheme reduces drastically the computational cost of the analyses while maintaining high accuracy in reproducing the key phenomena of the cascading natural hazard. Several numerical examples are presented to show the accuracy and robustness of the proposed coupling strategy and its applicability to large-scale landslide-generated wave events. The validation of the partitioned method is performed versus available results of other numerical methods, analytical solutions and experimental measures.

The final publication is available at Springer via http://dx.doi.org/10.1007/s40571-021-00440-y A numerical technique based on a CFD-DEM method is presented for the analysis of particle impact drilling (PID) systems. The method is built from a preexisting finite element Navier–Stokes solver for the fluid phase and a discrete element method module for the steel particles which this drilling technology utilizes to enhance the penetration rate. We provide a detailed description of the most relevant implementation issues, including our choice of the hydrodynamic forces appropriate for power-law fluids. We also discuss several critical aspects related to the validity of the simplifying assumptions that will be helpful to simulation engineers. We apply our simple, one-way coupled approach on designs provided by an industrial partner to illustrate its potential as an analysis tool for this promising drilling technology. The goal of this work is, on the one hand, to provide evidence for the usefulness of the numerical approach as a design tool for PID systems, as well as a detailed discussion of the different aspects to be assessed for an effective simulation campaign. On the other hand, a series of modeling aspects that require further work are identified. The focus of the simulation campaign presented is on the particles and fluid flow characteristics within the drill bit. We acknowledge the financial support to CIMNE via the CERCA Programme/Generalitat de Catalunya, and also that provided by the Spanish Ministry of Economy and Competitiveness, through the “Severo Ochoa Programme for Centres of Excellence in R&D” (CEX2018-000797-S). We would also like to thank Particle Drilling Technologies and MindMesh Inc. for technical and financial support and for granting permission to publish this work.

A strategy to couple continuous Galerkin (CG) and hybridizable discontinuous Galerkin (HDG) discretizations based only on the HDG hybrid variable is presented for linear thermal and elastic problems. The hybrid CG-HDG coupling exploits the definition of the numerical flux and the trace of the solution on the mesh faces to impose the transmission conditions between the CG and HDG subdomains. The continuity of the solution is imposed in the CG problem via Nitsche's method, whereas the equilibrium of the flux at the interface is naturally enforced as a Neumann condition in the HDG global problem. The proposed strategy does not affect the core structure of CG and HDG discretizations. In fact, the resulting formulation leads to a minimally-intrusive coupling, suitable to be integrated in existing CG and HDG libraries. Numerical experiments in two and three dimensions show optimal global convergence of the stress and superconvergence of the displacement field, locking-free approximation, as well as the potential to treat structural problems of engineering interest featuring multiple materials with compressible and nearly incompressible behaviors., 28 pages, 14 figures

The final publication is available at Springer via http://dx.doi.org/10.1007/s40571-019-00229-0 A semi-explicit Lagrangian scheme for the simulation of thermally coupled incompressible flow problems is presented. The model relies on combining an explicit multi-step solver for the momentum equation with an implicit heat equation solver. Computational cost of the model is reduced via application of an efficient strategy adopted for the solution of momentum/continuity system by the authors in their previous work. The applicability of the method to solving thermo-mechanical problems is studied via various numerical examples.

The final publication is available at Springer via http://dx.doi.org/10.1007/s00466-018-1624-3 Potential flow solvers represent an appealing alternative for the simulation of non-viscous subsonic flows. In order to deliver accurate results, such techniques require prescribing explicitly the so called Kutta condition, as well as adding a special treatment on the “wake” of the body. The wake is traditionally modelled by introducing a gap in the CFD mesh, which requires an often laborious meshing work. The novelty of the proposed work is to embed the wake within the CFD domain. The approach has obvious advantages in the context of aeroelastic optimization, where the position of the wake may change due to evolutionary steps of the geometry. This work presents a simple, yet effective, method for the imposition of the embedded wake boundary condition. The presented method preserves the possibility of employing iterative techniques in the solution of the linear problems which stem out of the discretization. Validation and verification of the solver are performed for a NACA 0012 airfoil.

We present a simulation tool for transient events in complex hydraulic networks. The code includes modelling of the transport of suspended cuttings in near-vertical wells. An unstructured finite volume formulation with implicit time integration has been chosen. The unconditional stability of the integrator makes the method suitable for the simulation of transient events over a wide range of characteristic timescales. It handles both very fast transients (e.g. fluid hammer events) and the long-term evolution of the well (e.g. hole cleaning operations). The software has been developed to address the need of the oil industry for a robust and efficient predictive tool allowing effective well control in managed pressure drilling operations. The physical modelling follows the standard practices accepted by the industry (e.g. mud rheology computations). The mathematical foundation of the algorithm is described followed by validation cases that illustrate its capabilities and accuracy. Finally, a practical industrial application example is provided to demonstrate the real-world performance of the software.

The paper presents a robust algorithm, which allows to implicitly describe and track immersed geometries within a background mesh. The background mesh is assumed to be unstructured and discretized by tetrahedrons. The contained geometry is assumed to be given as triangulated surface. Within the background mesh, the immersed geometry is described implicitly using a discontinuous distance function based on a level-set approach. This distance function allows to consider both, “double-sided” geometries like membrane or shell structures, and “single-sided” objects for which an enclosed volume is univocally defined. For the second case, the discontinuous distance function is complemented by a continuous signed distance function, whereas ray casting is applied to identify the closed volume regions. Furthermore, adaptive mesh refinement is employed to provide the necessary resolution of the background mesh. The proposed algorithm can handle arbitrarily complicated geometries, possibly containing modeling errors (i.e., gaps, overlaps or a non-unique orientation of surface normals). Another important advantage of the algorithm is the embarrassingly parallel nature of its operations. This characteristic allows for a straightforward parallelization using MPI. All developments were implemented within the open source framework “KratosMultiphysics” and are available under the BSD license. The capabilities of the implementation are demonstrated with various application examples involving practice-oriented geometries. The results finally show, that the algorithm is able to describe most complicated geometries within a background mesh, whereas the approximation quality may be directly controlled by mesh refinement.

This paper presents a finite element approach for modelling three-dimensional crack propagation in quasi-brittle materials, based on the strain injection and the crack-path field techniques. These numerical techniques were already tested and validated by static and dynamic simulations in 2D classical benchmarks [Dias et al., in: Monograph CIMNE No-134. International Center for Numerical Methods in Engineering, Barcelona, (2012); Oliver et al. in Comput Methods Appl Mech Eng 274:289–348, (2014); Lloberas-Valls et al. in Comput Methods Appl Mech Eng 308:499–534, (2016)] and, also, for modelling tensile crack propagation in real concrete structures, like concrete gravity dams [Dias et al. in Eng Fract Mech 154:288–310, (2016)]. The main advantages of the methodology are the low computational cost and the independence of the results on the size and orientation of the finite element mesh. These advantages were highlighted in previous works by the authors and motivate the present extension to 3D cases. The proposed methodology is implemented in the finite element framework using continuum constitutive models equipped with strain softening and consists, essentially, in injecting the elements candidate to capture the cracks with some goal oriented strain modes for improving the performance of the injected elements for simulating propagating displacement discontinuities. The goal-oriented strain modes are introduced by resorting to mixed formulations and to the Continuum Strong Discontinuity Approach (CSDA), while the crack position inside the finite elements is retrieved by resorting to the crack-path field technique. Representative numerical simulations in 3D benchmarks show that the advantages of the methodology already pointed out in 2D are kept in 3D scenarios

The final publication is available at Springer via http://dx.doi.org/10.1007/s40571-016-0154-z A methodology to integrate geographical information system (GIS) data with large-scale pedestrian simulations has been developed. Advances in automatic data acquisition and archiving from GIS databases, automatic input for pedestrian simulations, as well as scalable pedestrian simulation tools have made it possible to simulate pedestrians at the individual level for complete cities in real time. An example that simulates the evacuation of the city of Barcelona demonstrates that this is now possible. This is the first step towards a fully integrated crowd prediction and management tool that takes into account not only data gathered in real time from cameras, cell phones or other sensors, but also merges these with advanced simulation tools to predict the future state of the crowd.

The application of the semi-Lagrangian particle finite element method (SL–PFEM) for the seakeeping simulation of the wave adaptive modular vehicle under spray generating conditions is presented. The time integration of the Lagrangian advection is done using the explicit integration of the velocity and acceleration along the streamlines (X-IVAS). Despite the suitability of the SL–PFEM for the considered seakeeping application, small time steps were needed in the X-IVAS scheme to control the solution accuracy. A preliminary proposal to overcome this limitation of the X-IVAS scheme for seakeeping simulations is presented.

The final publication is available at Springer via http://dx.doi.org/10.1007/s40571-015-0089-9 A modular discrete element framework is presented for large-scale simulations of industrial grain-handling systems. Our framework enables us to simulate a markedly larger number of particles than previous studies, thereby allowing for efficient and more realistic process simulations. This is achieved by partitioning the particle dynamics into distinct regimes based on their contact interactions, and integrating them using different time-steps, while exchanging phase-space data between them. The framework is illustrated using numerical experiments based on fertilizer spreader applications. The model predictions show very good qualitative and quantitative agreement with available experimental data. Valuable insights are developed regarding the role of lift vs drag forces on the particle trajectories in-flight, and on the role of geometric discretization errors for surface meshing in governing the emergent behavior of a system of particles.

The final publication is available at Springer via http://dx.doi.org/10.1007/s40571-016-0117-4 Dam bottom outlets play a vital role in dam operation and safety, as they allow controlling the water surface elevation below the spillway level. For partial openings, water flows under the gate lip at high velocity and drags the air downstream of the gate, which may cause damages due to cavitation and vibration. The convenience of installing air vents in dam bottom outlets is well known by practitioners. The design of this element depends basically on the maximum air flow through the air vent, which in turn is a function of the specific geometry and the boundary conditions. The intrinsic features of this phenomenon makes it hard to analyse either on site or in full scaled experimental facilities. As a consequence, empirical formulas are frequently employed, which offer a conservative estimate of the maximum air flow. In this work, the particle finite element method was used to model the air–water interaction in Susqueda Dam bottom outlet, with different gate openings. Specific enhancements of the formulation were developed to consider air–water interaction. The results were analysed as compared to the conventional design criteria and to information gathered on site during the gate operation tests. This analysis suggests that numerical modelling with the PFEM can be helpful for the design of this kind of hydraulic works.

Discretization processes leading to numerical schemes sometimes produce undesirable effects. One potentially serious problem is that a discretization may produce the loss of validity of some of the physical principles or mathematical properties originally present in the continuous equation. Such loss may lead to uncertain results such as numerical instabilities or unexpected non-physical solutions. As a consequence, the compatibility of a discrete formulation with respect to intrinsic physical principles might be essential for the success of a numerical scheme. This paper addresses such type of issue. Its main objective is to demonstrate that standard Finite Element discretizations of the heat conduction equation violate Clausius’s postulate of the second law of thermodynamics, at nodal level. The problem occurs because non-physical, reversed nodal heat-fluxes arise in such discretizations. Conditions for compatibility of discrete nodal heat-fluxes with respect to Clausius’s postulate are derived here and named discrete thermodynamic compatibility conditions (DTCC). Simple numerical examples are presented to show the undesirable consequences of such failure. It must be pointed out that such DTCCs have previously appeared in the context of the study of the conditions that make discrete solutions to satisfy the discrete maximum principle (DMP). However, the present article does not put attention on such mathematical principle but on the satisfaction of a fundamental physical one: the second law of thermodynamics. Of course, from the presented point of view, it is clear that the violation of such fundamental law will cause, among different problems, the violation of the DMP. Fil: Limache, Alejandro Cesar. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina Fil: Idelsohn, Sergio Rodolfo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina

To explore the capability of unstructured mesh to simulate detonation wave propagation phenomena, numerical simulation of H2/air detonation using unstructured mesh was conducted. The unstructured mesh has several adv- antages such as easy mesh adaptation and flexibility to the complicated configurations. To examine the resolution dependency of the unstructured mesh, several simulations varying the mesh size were conducted and compared with a computed result using a structured mesh. The results show that the unstructured mesh solution captures the detailed structure of detonation wave, as well as the structured mesh solution. To capture the detailed detonation cell structure, the unstructured mesh simulations required at least twice, ideally 5times the resolution of structured mesh solution.

Model order reduction (MOR) is one of the most appealing choices for real-time simulation of non-linear solids. In this work a method is presented in which real time performance is achieved by means of the o -line solution of a (high dimensional) parametric problem that provides a sort of response surface or computational vademecum. This solution is then evaluated in real-time at feedback rates compatible with haptic devices, for instance (i.e., more than 1kHz). This high dimensional problem can be solved without the limitations imposed by the curse of dimensionality by employing Proper Generalized Decomposition (PGD) methods. Essentially, PGD assumes a separated representation for the essential eld of the problem. Here, an error estimator is proposed for this type of solutions that takes into account the non-linear character of the studied problems. This error estimator allows to compute the necessary number of modes employed to obtain an approximation to the solution within a prescribed error tolerance in a given quantity of interest.

The aim of the present work is to present an overview of some numerical procedures for the simulation of free surface flows within a porous structure. A particular algorithm developed by the authors for solving this type of problems is presented. A modified form of the classical Navier–Stokes equations is proposed, with the principal aim of simulating in a unified way the seepage flow inside rockfill-like porous material and the free surface flow in the clear fluid region. The problem is solved using a semi-explicit stabilized fractional step algorithm where velocity is calculated using a 4th order Runge–Kutta scheme. The numerical formulation is developed in an Eulerian framework using a level set technique to track the evolution of the free surface. An edge-based data structure is employed to allow an easy OpenMP parallelization of the resulting finite element code. The numerical model is validated against laboratory experiments on small scale rockfill dams and is compared with other existing methods for solving similar problems.

The Proper Generalized Decomposition (PGD) requires separability of the input data (e.g. physical properties, source term, boundary conditions, initial state). In many cases the input data is not expressed in a separated form and it has to be replaced by some separable approximation. These approximations constitute a new error source that, in some cases, may dominate the standard ones (discretization, truncation…) and control the nal accuracy of the PGD solution. In this work the relation between errors in the separated input data and the errors induced in the PGD solution is discussed. Error estimators proposed for homogenized problems and oscillation terms are adapted to asses the behaviour of the PGD errors resulting from approximated input data. The PGD is stable with respect to error in the separated data, with no critical amplification of the perturbations. Interestingly, we identified a high sensitiveness of the resulting accuracy on the selection of the sampling grid used to compute the separated data. The separation has to be performed on the basis of values sampled at integration points: sampling at the nodes defining the functional interpolation results in an important loss of accuracy. For the case of a Poisson problem separated in the spatial coordinates (a complex diffusivity function requires a separable approximation), the final PGD error is linear with the truncation error of the separated data. This relation is used to estimate the number of terms required in the separated data, that has to be in good agreement with the truncation error accepted in the PGD truncation (tolerance for the stoping criteria in the enrichment procedure). A sensible choice or the prescribed accuracy of the PGD solution has to be kept within the limits set by the errors in the separated input data.

In the present work an implementation of the Back and Forth Error Compensation and Correction (BFECC) algorithm specially suited for running on General-Purpose Graphics Processing Units (GPGPUs) through Nvidia`s Compute Unified Device Architecture (CUDA) is analyzed in order to solve transient pure advection equations. The objective is to compare it to a previous explicit version used in a Navier-Stokes solver fully written in CUDA. It turns out that BFECC could be implemented with unconditional stable stability using Semi-Lagrangian time integration allowing larger time steps than Eulerian ones. Fil: Costarelli, Santiago Daniel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico. Centro de Investigación de Métodos Computacionales; Argentina Fil: Storti, Mario Alberto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico. Centro de Investigación de Métodos Computacionales; Argentina Fil: Paz, Rodrigo Rafael. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico. Centro de Investigación de Métodos Computacionales; Argentina Fil: Dalcin, Lisandro Daniel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico. Centro de Investigación de Métodos Computacionales; Argentina Fil: Idelsohn, Sergio Rodolfo. Institució Catalana de Recerca i Estudis Avançats; España. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico. Centro de Investigación de Métodos Computacionales; Argentina

Trabajo presentado en el 4th International RILEM PhD Workshop, celebrado en Madrid (España), el 19 de noviembre de 2010, We have used Ab initio calculations to study the structural changes produced by the inclusion of light impurities into pure bcc iron. In order to obtain a clear picture of the mechanics of the phase changes Bain's pathway was studied in detail for pure iron. The position that hydrogen atoms tend to occupy at high densities favours octahedral sites inside the bcc matrix, producing an internal stress field that suggests a deformation that matches the prediction of martensitic transformation predicted by Bain's pathway. We have used Density Functional Theory in order to optimize the structures studied, obtaining the enthalpy of the configuration as a function of c/a, allowing a better understanding of the dynamics of the process of phase changes. © RILEM 2012.

The paper focuses on the Continuum Strong Discontinuity Approach (CSDA) to fracture mechanics, and the traction-separation cohesive laws induced from continuum dissipative models as their projections onto the failure interface. They are compared with the cohesive laws commonly used for the fracture simulation in quasi-brittle materials, typically concrete. Emphasis is placed in the analysis of the mechanical stress-strain states induced by the CSDA into the fracture process zone: first when the damage mechanism is initiated and, after, when the cohesive model determines the crack response. The influence of the material parameters, particularly the fracture energy and the initial continuum softening modulus, in the obtained phenomenological responses is also analyzed. Representative numerical solutions of fracture problems are finally presented. Fil: Huespe, Alfredo Edmundo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Desarrollo Tecnológico para la Industria Química. Universidad Nacional del Litoral. Instituto de Desarrollo Tecnológico para la Industria Química; Argentina Fil: Oliver, J.. Technical University of Catalonia; España Fil: Pulido, M. D. G.. Technical University of Catalonia; España Fil: Blanco, S.. Technical University of Catalonia; España Fil: Linero, D.L.. Technical University of Catalonia; España

This recommendation covers the description of non-destructive electrochemical test methods for the estimation in large size concrete structures of the instantaneous corrosion current densiti, icor, expressed in gA/cm2, by means of the so called Polarization Resistance technique, Rp, in order to asses the condition of embedded Steel reinforcement related to its corrosion. The values of icor, can be used to assess the rate of degradation of concrete structures affected by reinforcement corrosión. However, the cannot give information on the actual loss in steel cross section which, at present, only can be assessed by means of direct visual observation. Values of the free corrosión potential or half-cell potential, Ecorr [V], of the embedded reinforcing Steel and of the electrical concrete resistance, Re [f)], are obtained as preliminary steps of the Rp measurements. Values of the concrete resistivity, P [~m], can be calculated from Re values providing the geometrical arrangement of the electrodes enables this calculation. Both parameters, Ecorr and Re (or P) may be used to complement the reliability of the ico~r measurements.

The final publication is available at Springer via http://dx.doi.org/10.1007/BF00372173. Ket qualitative features of solutions exhibiting strong discontinuities in rate-independent inelastic solids are identified and exploited in the design of a new class of finite element approximations. The analysis shows that the softening law must be re-interpreted in a distributional sense for the continuum solutions to make mathematical sense and provides a precise physical interpretation to the softening modulus. These results are verified by numerical simulations employing a regularized discontinuous finite element method which circumvent the strong mesh-dependence exhibited by conventional methods, without resorting to viscosity or introducing additional ad-hoc parameters. The analysis is extended to a new class of anisotropic rate-independent damage models for brittle materials.