International Center for Numerical Methods in Engineering

The paper presents a fully meshless procedure fo solving partial differential equations. The approach termed generically the ‘finite point method’ is based on a weighted least square interpolation of point data and point collocation for evaluating the approximation integrals. Some examples showing the accuracy of the method for solution of adjoint and non-self adjoint equations typical of convective-diffusive transport and also to the analysis of compressible fluid mechanics problem are presented.

Abstract Particle Methods are those in which the problem is represented by a discrete number of particles. Each particle moves accordingly with its own mass and the external/internal forces applied to it. Particle Methods may be used for both, discrete and continuous problems. In this paper, a Particle Method is used to solve the continuous fluid mechanics equations. To evaluate the external applied forces on each particle, the incompressible Navier–Stokes equations using a Lagrangian formulation are solved at each time step. The interpolation functions are those used in the Meshless Finite Element Method and the time integration is introduced by an implicit fractional‐step method. In this manner classical stabilization terms used in the momentum equations are unnecessary due to lack of convective terms in the Lagrangian formulation. Once the forces are evaluated, the particles move independently of the mesh. All the information is transmitted by the particles. Fluid–structure interaction problems including free‐fluid‐surfaces, breaking waves and fluid particle separation may be easily solved with this methodology. Copyright © 2004 John Wiley & Sons, Ltd.

We present a general formulation for the analysis of fluid-structure interaction problems using the particle finite element method (PFEM). The key feature of the PFEM is the use of a Lagrangian description to model the motion of nodes (particles) in both the fluid and the structure domains. Nodes are thus viewed as particles which can freely move and even separate from the main analysis domain representing, for instance, the effect of water drops. A mesh connects the nodes defining the discretized domain where the governing equations, expressed in an integral form, are solved as in the standard FEM. The necessary stabilization for dealing with the incompressibility condition in the fluid is introduced via the finite calculus (FIC) method. A fractional step scheme for the transient coupled fluid-structure solution is described. Examples of application of the PFEM method to solve a number of fluid-structure interaction problems involving large motions of the free surface and splashing of waves are presented.

Abstract The aim of this chapter is to provide an in‐depth survey of arbitrary Lagrangian–Eulerian (ALE) methods, including both conceptual aspects of the mixed kinematical description and numerical implementation details. Applications are discussed in fluid dynamics, nonlinear solid mechanics, and coupled problems describing fluid–structure interaction. The need for an adequate mesh‐update strategy is underlined, and various automatic mesh‐displacement prescription algorithms are reviewed. This includes mesh‐regularization methods essentially based on geometrical concepts, as well as mesh‐adaptation techniques aimed at optimizing the computational mesh according to some error indicator. Emphasis is then placed on particular issues related to the modeling of compressible and incompressible flow and nonlinear solid mechanics problems. This includes the treatment of convective terms in the conservation equations for mass, momentum, and energy, as well as a discussion of stress‐update procedures for materials with history‐dependent constitutive behavior.

Abstract The paper outlines the formulation of a novel algorithm which can be used for the solution of both compressible and incompressible Navier‐Stokes or Euler equations. Full incompressibility can be dealt with if the algorithm is used in its semi‐explicit form and its structure permits arbitrary interpolation functions to be used avoiding the Babuška‐Brezzi restriction. In a fully explicit version it introduces a rational form of balancing dissipation avoiding the use of arbitrary parameters and forms for this.

Abstract Blended cements, where Portland cement clinker is partially replaced by supplementary cementitious materials (SCMs), provide the most feasible route for reducing carbon dioxide emissions associated with concrete production. However, lowering the clinker content can lead to an increasing risk of neutralisation of the concrete pore solution and potential reinforcement corrosion due to carbonation. carbonation of concrete with SCMs differs from carbonation of concrete solely based on Portland cement (PC). This is a consequence of the differences in the hydrate phase assemblage and pore solution chemistry, as well as the pore structure and transport properties, when varying the binder composition, age and curing conditions of the concretes. The carbonation mechanism and kinetics also depend on the saturation degree of the concrete and CO 2 partial pressure which in turn depends on exposure conditions (e.g. relative humidity, volume, and duration of water in contact with the concrete surface and temperature conditions). This in turn influence the microstructural changes identified upon carbonation. This literature review, prepared by members of RILEM technical committee 281-CCC carbonation of concrete with supplementary cementitious materials, working groups 1 and 2, elucidates the effect of numerous SCM characteristics, exposure environments and curing conditions on the carbonation mechanism, kinetics and structural alterations in cementitious systems containing SCMs.

Para dar respuesta a los requerimientos en materia de aguas definidos en las directrices, reglamentos y recomendaciones existentes en la legislación española, los cuales están mayoritariamente basados en directivas europeas, se ha desarrollado una herramienta de modelización numérica del flujo de agua en lámina libre en 2 dimensiones. La herramienta, llamada Iber, combina un módulo hidrodinámico, un módulo de turbulencia y un módulo de transporte de sedimentos, y utiliza el método de volúmenes finitos para resolver las ecuaciones correspondientes. Al módulo de cálculo se le ha adaptado una interfaz que se basa en el software de preproceso y posproceso GiD, desarrollado por CIMNE. El resultado es una herramienta de modelización numérica del flujo de agua y sedimentos en ríos y estuarios, que utiliza esquemas numéricos avanzados especialmente estables y robustos en cualquier situación pero especialmente adecuados para flujos discontinuos y, en concreto, para cauces torrenciales y regímenes irregulares.

Abstract A meshless method is presented which has the advantages of the good meshless methods concerning the ease of introduction of node connectivity in a bounded time of order n , and the condition that the shape functions depend only on the node positions. Furthermore, the method proposed also shares several of the advantages of the finite element method such as: (a) the simplicity of the shape functions in a large part of the domain; (b) C 0 continuity between elements, which allows the treatment of material discontinuities, and (c) ease of introduction of the boundary conditions. Copyright © 2003 John Wiley & Sons, Ltd.

The paper explores the behaviour of compacted soils throughout the (dry density–water content) compaction plane by means of a conceptual framework that incorporates microstructural information. The engineering properties of compacted soils are described by an initial state in terms of a yielding stress, soil suction and a microstructural state variable. Microstructure is defined by the ratio of microvoid volume to total void volume. The pattern of variation of the microstructural parameter within the compaction plane has been determined, for some compacted soils, by analysing mercury intrusion porosimetry data. The microstructure of wet and dry compaction conditions can then be quantified. To ensure consistency, the framework is cast in the form of a constitutive model defined in terms of an effective suction and a constitutive stress that incorporate the microstructural variable. The model is shown to be consistent with a number of experimental observations and, in particular, it explains the intrinsic collapse potential of compacted soils. It predicts, for a common initial suction, a higher collapse potential for dry of optimum conditions than for wet compaction. It also predicts in a natural manner the observed evolution of soil compressibility during drained or undrained loading. Model capabilities are illustrated by application to a testing programme on statically compacted samples of low-plasticity silty clay. The compression behaviour of samples compacted wet and dry of optimum and the variation of collapse strains with confining stress have been successfully reproduced by the model.

Intestinal organoids capture essential features of the intestinal epithelium such as crypt folding, cellular compartmentalization and collective movements. Each of these processes and their coordination require patterned forces that are at present unknown. Here we map three-dimensional cellular forces in mouse intestinal organoids grown on soft hydrogels. We show that these organoids exhibit a non-monotonic stress distribution that defines mechanical and functional compartments. The stem cell compartment pushes the extracellular matrix and folds through apical constriction, whereas the transit amplifying zone pulls the extracellular matrix and elongates through basal constriction. The size of the stem cell compartment depends on the extracellular-matrix stiffness and endogenous cellular forces. Computational modelling reveals that crypt shape and force distribution rely on cell surface tensions following cortical actomyosin density. Finally, cells are pulled out of the crypt along a gradient of increasing tension. Our study unveils how patterned forces enable compartmentalization, folding and collective migration in the intestinal epithelium.

Advances in the production of optical fibers made possible the recent development of innovative sensing systems for the health monitoring of civil structures. The main reasons for this development are the reduced weight and dimensions of fiber optic sensors, the strong immunity to electromagnetic interference, the improved environmental resistance and the scale flexibility for small-gauge and long-gauge measurements. These systems can provide high-resolution and measurement capabilities that are not feasible with conventional technologies. In addition, they can be manufactured at a low cost and they offer a number of key advantages, including the ability to multiplex an appreciable number of sensors along a single fiber and interrogate such systems over large distances. For these reasons, it is evident that fiber optic sensors will change the instrumentation industry in the same way fiber optics has revolutionized communications. This paper provides an overview of the intensity modulated and spectrometric fiber optic sensors and techniques to assess the condition of existing structures in order to enhance the durability of the new bridges, increasing lifetime and reliability and decreasing maintenance activities. Application of these sensors to monitoring strain, temperature, inclination, acceleration, load measurements, ice detection, vehicles speeds and weights, and corrosion and cracking of reinforced and prestressed concrete structures will be described.

We have conducted a direct numerical simulation study of dilute turbulent particulate flow in a vertical plane channel by considering thousands of finite-size rigid particles with resolved phase interfaces. The particle diameter corresponds to approximately 11 wall units and their terminal Reynolds number is set to 136. The fluid flow with a bulk Reynolds number of 2700 is directed upward, which maintains the particles suspended on average. Two density ratios were simulated, differing by a factor of 4.5. The corresponding Stokes numbers of the two flow cases were O(10) in the near-wall region and O(1) in the outer flow. We have observed the formation of large-scale elongated streaklike structures with streamwise dimensions of the order of eight channel half-widths and cross-stream dimensions of the order of one half-width. At the same time, we have found no evidence of significant formation of particle clusters, which suggests that the large structures are due to an intrinsic instability of the flow, which is triggered by the presence of the particles. It was found that the mean fluid velocity profile tends toward a concave shape, and the turbulence intensity and the normal stress anisotropy are strongly increased. The effect of varying the Stokes number while maintaining the buoyancy, particle size, and volume fraction constant was relatively weak.

In this article, we discuss connections between sequential system identification and control for linear time-invariant systems, often termed indirect data-driven control, as well as a contemporary direct data-driven control approach seeking an optimal decision compatible with recorded data assembled in a Hankel matrix and robustified through suitable regularizations. We formulate these two problems in the language of behavioral systems theory and parametric mathematical programs, and we bridge them through a multicriteria formulation trading off system identification and control objectives. We illustrate our results with two methods from subspace identification and control: namely, subspace predictive control and low-rank approximation, which constrain trajectories to be consistent with a nonparametric predictor derived from (respectively, the column span of) a data Hankel matrix. In both cases, we conclude that direct and regularized data-driven control can be derived as convex relaxation of the indirect approach, and the regularizations account for an implicit identification step. Our analysis further reveals a novel regularizer and a plausible hypothesis explaining the remarkable empirical performance of direct methods on nonlinear systems.

Abstract The particle finite element method (PFEM) is a powerful and robust numerical tool for the simulation of multi-physics problems in evolving domains. The PFEM exploits the Lagrangian framework to automatically identify and follow interfaces between different materials (e.g. fluid–fluid, fluid–solid or free surfaces). The method solves the governing equations with the standard finite element method and overcomes mesh distortion issues using a fast and efficient remeshing procedure. The flexibility and robustness of the method together with its capability for dealing with large topological variations of the computational domains, explain its success for solving a wide range of industrial and engineering problems. This paper provides an extended overview of the theory and applications of the method, giving the tools required to understand the PFEM from its basic ideas to the more advanced applications. Moreover, this work aims to confirm the flexibility and robustness of the PFEM for a broad range of engineering applications. Furthermore, presenting the advantages and disadvantages of the method, this overview can be the starting point for improvements of PFEM technology and for widening its application fields.

In 1995 the two senior authors of the present paper introduced a new algorithm designed to replace the Taylor–Galerkin (or Lax–Wendroff) methods, used by them so far in the solution of compressible flow problems. The new algorithm was applicable to a wide variety of situations, including fully incompressible flows and shallow water equations, as well as supersonic and hypersonic situations, and has proved to be always at least as accurate as other algorithms currently used. The algorithm is based on the solution of conservation equations of fluid mechanics to avoid any possibility of spurious solutions that may otherwise result. The main aspect of the procedure is to split the equations into two parts, (1) a part that is a set of simple scalar equations of convective–diffusion type for which it is well known that the characteristic Galerkin procedure yields an optimal solution; and (2) the part where the equations are self-adjoint and therefore discretized optimally by the Galerkin procedure. It is possible to solve both the first and second parts of the system explicitly, retaining there the time step limitations of the Taylor–Galerkin procedure. But it is also possible to use semi-implicit processes where in the first part we use a much bigger time step generally governed by the Peclet number of the system while the second part is solved implicitly and is unconditionally stable. It turns out that the characteristic-based-split (CBS) process allows equal interpolation to be used for all system variables without difficulties when the incompressible or nearly incompressible stage is reached. It is hoped that the paper will help to make the algorithm more widely available and understood by the profession and that its advantages can be widely realised. Copyright © 1999 John Wiley & Sons, Ltd.

We have performed direct numerical simulations of turbulent flows in a square duct considering a range of Reynolds numbers spanning from a marginal state up to fully developed turbulent states at low Reynolds numbers. The main motivation stems from the relatively poor knowledge about the basic physical mechanisms that are responsible for one of the most outstanding features of this class of turbulent flows: Prandtl's secondary motion of the second kind. In particular, the focus is upon the role of flow structures in its generation and characterization when increasing the Reynolds number. We present a two-fold scenario. On the one hand, buffer layer structures determine the distribution of mean streamwise vorticity. On the other hand, the shape and the quantitative character of the mean secondary flow, defined through the mean cross-stream function, are influenced by motions taking place at larger scales. It is shown that high velocity streaks are preferentially located in the corner region (e.g. less than 50 wall units apart from a sidewall), flanked by low velocity ones. These locations are determined by the positioning of quasi-streamwise vortices with a preferential sign of rotation in agreement with the above described velocity streaks' positions. This preferential arrangement of the classical buffer layer structures determines the pattern of the mean streamwise vorticity that approaches the corners with increasing Reynolds number. On the other hand, the centre of the mean secondary flow, defined as the position of the extrema of the mean cross-stream function (computed using the mean streamwise vorticity), remains at a constant location departing from the mean streamwise vorticity field for larger Reynolds numbers, i.e. it scales in outer units. This paper also presents a detailed validation of the numerical technique including a comparison of the numerical results with data obtained from a companion experiment.

We present a general framework for the dimensional reduction, in terms of number of degrees of freedom as well as number of integration points (“hyper-reduction”), of nonlinear parameterized finite element (FE) models. The reduction process is divided into two sequential stages. The first stage consists in a common Galerkin projection onto a reduced-order space, as well as in the condensation of boundary conditions and external forces. For the second stage (reduction in number of integration points), we present a novel cubature scheme that efficiently determines optimal points and associated positive weights so that the error in integrating reduced internal forces is minimized. The distinguishing features of the proposed method are: (1) The minimization problem is posed in terms of orthogonal basis vector (obtained via a partitioned Singular Value Decomposition) rather that in terms of snapshots of the integrand. (2) The volume of the domain is exactly integrated. (3) The selection algorithm need not solve in all iterations a nonnegative least-squares problem to force the positiveness of the weights. Furthermore, we show that the proposed method converges to the absolute minimum (zero integration error) when the number of selected points is equal to the number of internal force modes included in the objective function. We illustrate this model reduction methodology by two nonlinear, structural examples (quasi-static bending and resonant vibration of elastoplastic composite plates). In both examples, the number of integration points is reduced three order of magnitudes (with respect to FE analyses) without significantly sacrificing accuracy.

Tsunamis are unpredictable and infrequent but potentially large impact natural disasters. To prepare, mitigate and prevent losses from tsunamis, probabilistic hazard and risk analysis methods have been developed and have proved useful. However, large gaps and uncertainties still exist and many steps in the assessment methods lack information, theoretical foundation, or commonly accepted methods. Moreover, applied methods have very different levels of maturity, from already advanced probabilistic tsunami hazard analysis for earthquake sources, to less mature probabilistic risk analysis. In this review we give an overview of the current state of probabilistic tsunami hazard and risk analysis. Identifying research gaps, we offer suggestions for future research directions. An extensive literature list allows for branching into diverse aspects of this scientific approach.

UNLABELLED: We developed a reporting guideline to provide authors with guidance about what should be reported when writing a paper for publication in a scientific journal using a particular type of research design: the single-case experimental design. This report describes the methods used to develop the Single-Case Reporting guideline In BEhavioural interventions (SCRIBE) 2016. As a result of 2 online surveys and a 2-day meeting of experts, the SCRIBE 2016 checklist was developed, which is a set of 26 items that authors need to address when writing about single-case research. This article complements the more detailed SCRIBE 2016 Explanation and Elaboration article (Tate et al., 2016) that provides a rationale for each of the items and examples of adequate reporting from the literature. Both these resources will assist authors to prepare reports of single-case research with clarity, completeness, accuracy, and transparency. They will also provide journal reviewers and editors with a practical checklist against which such reports may be critically evaluated. We recommend that the SCRIBE 2016 is used by authors preparing manuscripts describing single-case research for publication, as well as journal reviewers and editors who are evaluating such manuscripts. SCIENTIFIC ABSTRACT: Reporting guidelines, such as the Consolidated Standards of Reporting Trials (CONSORT) Statement, improve the reporting of research in the medical literature (Turner et al., 2012). Many such guidelines exist and the CONSORT Extension to Nonpharmacological Trials (Boutron et al., 2008) provides suitable guidance for reporting between-groups intervention studies in the behavioral sciences. The CONSORT Extension for N-of-1 Trials (CENT 2015) was developed for multiple crossover trials with single individuals in the medical sciences (Shamseer et al., 2015; Vohra et al., 2015), but there is no reporting guideline in the CONSORT tradition for single-case research used in the behavioral sciences. We developed the Single-Case Reporting guideline In BEhavioural interventions (SCRIBE) 2016 to meet this need. This Statement article describes the methodology of the development of the SCRIBE 2016, along with the outcome of 2 Delphi surveys and a consensus meeting of experts. We present the resulting 26-item SCRIBE 2016 checklist. The article complements the more detailed SCRIBE 2016 Explanation and Elaboration article (Tate et al., 2016) that provides a rationale for each of the items and examples of adequate reporting from the literature. Both these resources will assist authors to prepare reports of single-case research with clarity, completeness, accuracy, and transparency. They will also provide journal reviewers and editors with a practical checklist against which such reports may be critically evaluated.

The European parliament has declared a global “climate and environmental emergency” on 28 November 2019. Given that, climate change is a clear strategic issue all around the world. Then, greenhouse gas emissions are reported by each country to the United Nations Framework Convention on Climate Change (UNFCCC) every year. In addition, The Intergovernmental Panel on Climate Change (IPCC) in the “2006 IPCC Guidelines for National Greenhouse Gas Inventories” give the procedure to calculate and manage the national greenhouse gases (GHG) emissions. However, these guidelines do not provide any method to consider the net carbon dioxide emissions to the atmosphere (released in clinker fabrication minus those due to concrete carbonation) by the Portland cement clinker industry. This topic should be implemented in the climatic models of the next IPCC assessment report. This paper provides an easy procedure of estimating net CO2 emissions proposed in the “recarbonation project” (simplified method); that is to say, carbon dioxide uptake during the service-life stage is considered as the 20% of the CO2 released by the calcination (process emissions), whereas the end-of-life and secondary usage is only the 3% of the CO2 released by calcination. The outcome of this study reveals that 31,290.753 tons of carbon dioxide will be absorbed by the cement-based materials produced in Spain with the cements manufactured from 2005 to 2015.