Nodally exact Ritz discretizations of 1D diffusion–absorption and Helmholtz equations by variational FIC and modified equation methods

The final publication is available at Springer via http://dx.doi.org/10.1007/s00466-005-0011-zThis article presents the first application of the Finite Calculus (FIC) in a Ritz-FEM variational framework. FIC provides a steplength parametrization of mesh dimensions, which is used to modify the shape functions. This approach is applied to the FEM discretization of the steady-state, one-dimensional, diffusion–absorption and Helmholtz equations. Parametrized linear shape functions are directly inserted into a FIC functional. The resulting Ritz-FIC equations are symmetric and carry a element-level free parameter coming from the function modification process. Both constant- and variable-coefficient cases are studied. It is shown that the parameter can be used to produce nodally exact solutions for the constant coefficient case. The optimal value is found by matching the finite-order modified differential equation (FOMoDE) of the Ritz-FIC equations with the original field equation. The inclusion of the Ritz-FIC models in the context of templates is examined. This inclusion shows that there is an infinite number of nodally exact models for the constant coefficient case. The ingredients of these methods (FIC, Ritz, MoDE and templates) can be extended to multiple dimensions.Peer ReviewedPostprint (author's final draft

Posibilidades de los métodos numéricos en el mundo industrial

Postprint (published version

A finite point method to solve shallow water equations

Preprin

The CIMNE model for generating knowledge on computational engineering and its transfer to society

We present an overview of the model implemented by the International Centre for Numerical Methods in Engineering (CIMNE, www.cimne.com ) for the generation of scientific and technical knowledge on computational engineering, understood in the broad sense, the subsequent generation of product resulting from the research activities and the transfer of these products to society for their exploitation in the industrial sector. We present examples of application of the CIMNE model to a number of academic, scientific and industry activities of CIMNE .Preprin

The cycle of ideas in research development and technology transfer

When people talks about research, development and technology transfer (RTD in short) and commercial industrial activities, it is quite usual to mix up objectives, resources and responsibilities. There is a growing opinion, spread by some public administrators and the media, that research groups at universities and RTD centers must be "profitable". In few words, many persons would wish that research is funded with competitive funds or even loans (either public or private) and that the outcomes of the research get into the market rapidly, so that the profits from marketing the corresponding products will allow research groups to be financially self-sustaining. It is also usually argued that companies should increase their research activity and invest more resources in finding new discoveries. The humoristic limit of this situation is pretending that companies will produce Nobel laureates and that a research group would be among the 10 first positions in the ranking of economic organizations in a country. Certainly, both these situations are not impossible, although they are highly improbable. The opposite case is, unfortunately, more frequent, ie. that researchers compelled by their needs to guarantee the financial survival of their groups, abandon or considerably reduce their fundamental research activities, and that companies, confusing what is a Technical Department with an Innovation Department, consider themselves self-sufficient and underestimate or ignore the contact with the real RTD world.Preprin

A FIC-based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems

This is the peer reviewed version of the following article: [de-Pouplana, I., and Oñate, E. (2017) A FIC-based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems. Int. J. Numer. Anal. Meth. Geomech., 41: 110–134. doi: 10.1002/nag.2550], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nag.2550/abstract. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving."A new mixed displacement-pressure element for solving solid–pore fluid interaction problems is presented. In the resulting coupled system of equations, the balance of momentum equation remains unaltered, while the mass balance equation for the pore fluid is stabilized with the inclusion of higher-order terms multiplied by arbitrary dimensions in space, following the finite calculus (FIC) procedure. The stabilized FIC-FEM formulation can be applied to any kind of interpolation for the displacements and the pressure, but in this work, we have used linear elements of equal order interpolation for both set of unknowns. Examples in 2D and 3D are presented to illustrate the accuracy of the stabilized formulation for solid–pore fluid interaction problems.Peer ReviewedPostprint (author's final draft

Discrete/finite element modelling of rock cutting with a TBM disc cutter

The final publication is available at Springer via http://dx.doi.org/10.1007/s00603-016-1133-7This paper presents advanced computer simulation of rock cutting process typical for excavation works in civil engineering. Theoretical formulation of the hybrid discrete/finite element model has been presented. The discrete and finite element methods have been used in different subdomains of a rock sample according to expected material behaviour, the part which is fractured and damaged during cutting is discretized with the discrete elements while the other part is treated as a continuous body and it is modelled using the finite element method. In this way, an optimum model is created, enabling a proper representation of the physical phenomena during cutting and efficient numerical computation. The model has been applied to simulation of the laboratory test of rock cutting with a single TBM (tunnel boring machine) disc cutter. The micromechanical parameters have been determined using the dimensionless relationships between micro- and macroscopic parameters. A number of numerical simulations of the LCM test in the unrelieved and relieved cutting modes have been performed. Numerical results have been compared with available data from in-situ measurements in a real TBM as well as with the theoretical predictions showing quite a good agreement. The numerical model has provided a new insight into the cutting mechanism enabling us to investigate the stress and pressure distribution at the tool–rock interaction. Sensitivity analysis of rock cutting performed for different parameters including disc geometry, cutting velocity, disc penetration and spacing has shown that the presented numerical model is a suitable tool for the design and optimization of rock cutting process.Peer ReviewedPostprint (published version

Simulation of flows with violent free surface motion and moving objects using unstructured grids

This is the peer reviewed version of the following article: [Löhner, R. , Yang, C. and Oñate, E. (2007), Simulation of flows with violent free surface motion and moving objects using unstructured grids. Int. J. Numer. Meth. Fluids, 53: 1315-1338. doi:10.1002/fld.1244], which has been published in final form at https://doi.org/10.1002/fld.1244. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.A volume of fluid (VOF) technique has been developed and coupled with an incompressible Euler/Navier–Stokes solver operating on adaptive, unstructured grids to simulate the interactions of extreme waves and three-dimensional structures. The present implementation follows the classic VOF implementation for the liquid–gas system, considering only the liquid phase. Extrapolation algorithms are used to obtain velocities and pressure in the gas region near the free surface. The VOF technique is validated against the classic dam-break problem, as well as series of 2D sloshing experiments and results from SPH calculations. These and a series of other examples demonstrate that the ability of the present approach to simulate violent free surface flows with strong nonlinear behaviour.Peer ReviewedPostprint (author's final draft

Variational formulation of the finite calculus equations in solid mechanics and diffusion-reaction problems

We present a variational formulation of the finite calculus (FIC) equations for problems in mechanics governed by differential equations with symmetric operators. Applications considered include solid mechanics, diffusion-transport and diffusion-reaction problems. The key of the variational formulation is the identification of the FIC governing equations with the classical differential equations of mechanics written in terms of modified non-local variables. A total potential energy (TPE) functional is found in terms of the modified variables. The FIC equations in the domain and the boundary are recovered as the Euler-Lagrange equations and the natural boundary condition of the TPE functional, respectively. Symmetric finite element equations are obtained after discretization of the TPE functional, therefore preserving the symmetry of the governing infinitesimal equations. The variational FIC expression is reinterpreted as a Petrov Galerkin weighted residual form of the original FIC equations with non-local weighting functions. The analogy of the variational FIC-FEM formulation with a discontinuous Galerkin method is recognized. Extensions to multidimensional linear elastostatics and diffusion-reaction problems are presented.Preprin

Membrane structures formed by low pressure inflatable tubes: new analysis methods and recent constructions

This paper shows applications of a recently developed thin shell element adequate for the analysis of membrane and inflatable structures. The element is a three node triangle with only translational degrees of freedom that uses the configuration of the three adjacent elements to evaluate the strains in terms of the nodal displacements only. This allows to compute (constant) bending strains and (linear) membrane strains using a total Lagrangian formulation. Several examples, including inflation and deflation of membranes and some practical applications to the analysis, design and construction of membrane structures formed by low pressure inflatable tubes are presented.Preprin