You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
<!-- metadata commented in wiki content
2
3
4
==The Particle Finite Element Method (PFEM). An effective numerical technique for solving marine, naval and harbour engineering problems==
5
6
E. Oñate<sup>1</sup>, S.R. Idelsohn<sup>*</sup>, M.A. Celigueta<sup>1</sup> and B. Suárez<sup>2</sup>
7
8
''<sup>1</sup>'' ''International Center for Numerical Methods in Engineering (CIMNE), Spain''
9
10
''<sup>*</sup>'' ''ICREA Research Professor at CIMNE''
11
12
''<sup>2</sup>'' ''Universitat Politècnica de Catalunya (UPC), Spain''
13
-->
14
15
==Abstract==
16
17
We present some developments in the Particle Finite Element Method (PFEM) for the solution of complex coupled problems in marine, naval and harbour engineering involving fluid-soil-structure interaction (FSSI). The PFEM uses an updated Lagrangian description to model the motion of nodes (particles) in a continuum domain containing fluid, soil/rock and structures subdomains. A mesh connects the nodes defining the discretized domain where the governing equations for each of the constituent materials are solved with the FEM. The stabilization for dealing with an incompressibility material is introduced via the finite calculus (FIC) method. An incremental iterative scheme for solving the non linear transient FSSI problem is described. The procedure to model frictional contact conditions and material erosion at fluid-solid and solid-solid interfaces is described. We present examples of application of the PFEM to solve FSSI problems in marine, naval and harbour engineering such as the motion of rocks by water streams, the stability of breakwaters and constructions under sea waves, the sinking of ships and the collision of a ship with ice blocks.
18
19
==1 Introduction==
20
21
The analysis of problems involving the interaction of fluids, soil/rocks and structures is relevant in many areas of marine, naval and harbour engineering. Some representative examples are the study of off-shore and harbour structures under large waves, constructions hit by tsunamis, erosion of sea bed and landscape adjacent to sea shore, motion of ships under severe sea conditions, the simulation of the sinking of a ship, the collision of a ship with ice blocks, excavation and drilling problems in the sea in petroleum and gas engineering, etc.
22
23
The authors have developed in previous works a particular class of Lagrangian formulation for solving problems involving complex interactions between (free surface) fluids and solids. The so-called the particle finite element method ([[http://www.cimne.com/pfem PFEM]]), treats the mesh nodes in the fluid and solid domains as particles which can freely move and even separate from the main fluid domain representing, for instance, the effect of water drops. A mesh connects the nodes discretizing the domain where the governing equations are solved using a stabilized FEM.
24
25
An advantage of the Lagrangian formulation is that the non-linear and non symmetric convective terms disappear from the fluid equations <span id="citeF-1"></span>[[#cite-1|[1]],<span id="citeF-2"></span>[[#cite-2|2]]]. The difficulty is however transferred to the problem of adequately (and efficiently) moving the mesh nodes. The theory and applications of the PFEM are reported in <span id="citeF-1"></span>[[#cite-1|[1]]-<span id="citeF-18"></span>[[#cite-18|18]]].
26
27
The solution of a FSSI problem requires solving the momentum and mass balance equations for the underlying continuum that may include incompressible domains (as it is the case for a fluid and some soils). In our work we use a stabilized mixed FEM based on Finite Calculus (FIC) which allows us to solve incompressible continua with a linear approximation for the velocity and pressure variables <span id="citeF-6"></span>[[#cite-6|[6]],<span id="citeF-9"></span>[[#cite-9|9]],<span id="citeF-18"></span>[[#cite-18|18]],<span id="citeF-20"></span>[[#cite-20|20]]-<span id="citeF-23"></span>[[#cite-23|23]]].
28
29
The layout of the chapter is as follows. In the next section the key ideas of the PFEM are outlined. Next the basic equations for a general continuum using a Lagrangian description and the FIC formulation are schematically presented. Then an algorithm for the transient solution is briefly described. The treatment of the coupled FSSI problem and the method for mesh generation and for identification of the free surface nodes are outlined. The procedure for treating the frictional contact interaction between fluid, soil and structure interfaces is explained. We present several examples of application of the PFEM to solve FSSI problems in marine, naval and harbour engineering such as the motion of rocks by water streams, the stability of breakwaters and constructions hit by sea waves, the study of the sinking of ships and the collision of a ship with ice blocks.
30
31
==2 The basis of the particle finite element method==
32
33
Let us consider a domain containing both fluid and solid subdomains (the solid subdomain may include soil/rock materials and/or structural elements). The moving fluid particles interact with the solid boundaries, thereby inducing the deformation of the solid which in turn affects the flow motion and, therefore, the problem is fully coupled.
34
35
In the PFEM both the fluid and the solid domains are modelled using an ''updated'' ''Lagrangian formulation'' <span id="citeF-24"></span>[[#cite-24|[24]]]. That is, all variables are assumed to be known in the ''current configuration'' at time  <math>t</math> . The new set of variables in both domains is sought for in the ''next or updated configuration'' at time <math>t+\Delta t</math> . The finite element method (FEM) is used to solve the equations of continuum mechanics for each of the subdomains. Hence a mesh discretizing these domains must be generated in order to solve the governing equations for each subdomain in the standard FEM fashion <span id="citeF-24"></span> [[#cite-24|[24]]-<span id="citeF-27"></span>[[#cite-27|27]]].
36
37
The quality of the numerical solution depends on the discretization chosen as in the standard FEM. Adaptive mesh refinement techniques can be used to improve the solution.
38
39
===2.1 Basic steps of the PFEM===
40
41
For clarity purposes we will define the '' collection or cloud of nodes (<math display="inline">C</math>)'' pertaining to the fluid and solid domains, the ''volume (<math display="inline">V</math>)'' defining the analysis domain for the fluid and the solid and the ''mesh (<math display="inline">M</math>)'' discretizing both domains.
42
43
A typical solution with the PFEM involves the following steps.
44
<ol>
45
<span id="step-1></span>
46
<li> The starting point at each time step is the cloud of points in the fluid and solid domains. For instance  <math>{}^nC</math> denotes the cloud at time  <math>t=t_n</math> (Figure [[#img-1|1]]).</li>
47
<span id="step-2></span>
48
<li> Identify the boundaries for both the fluid and solid domains defining the analysis domain  <math>{}^nV</math> in the fluid and the solid. This is an essential step as some boundaries (such as the free surface in fluids) may be severely distorted during the solution, including separation and re-entering of nodes. The Alpha Shape method <span id="citeF-19"></span>[[#cite-19|[19]]] is used for the boundary definition.</li>
49
<span id="step-3></span>
50
<li> Discretize the fluid and solid domains with a finite element mesh <math>{}^nM.</math> We use an effect mesh generation scheme based on the extended Delaunay tesselation <span id="citeF-3"></span>[[#cite-3|[3]],<span id="citeF-4"></span>[[#cite-4|4]]].</li>
51
<span id="step-4></span>
52
<li> Solve the coupled Lagrangian equations of motion for the overall continuum. Compute the state variables in at the next (updated) configuration for <math>t+\Delta t</math> : velocities, pressure and viscous stresses in the fluid and displacements, stresses and strains in the solid.</li>
53
<span id="step-5></span>
54
<li> Move the mesh nodes to a new position  <math>{}^{n+1}C</math> where <math displya="inline">n+1</math> denotes the time <math>t_n+\Delta t</math>, in terms of the time increment size. This step is typically a consequence of the solution process of step [[#step-4|4]].</li>
55
<span id="step-6></span>
56
<li> Go back to step [[#step-1|1]] and repeat the solution for the next time step to obtain <math>{}^{n+2}C</math> (Figure [[#img-1|1]]).</li>
57
</ol>
58
59
<div id='img-x'></div>
60
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
61
|-
62
| [[Image:Draft_Samper_913569722-image5.png|400px]]
63
|- style="text-align: center; font-size: 75%;"
64
| colspan="x" | '''Figure 1:'''Sequence of steps to update a ``cloud´´ of nodes representing a domain containing a fluid and a solid part from time  <math>n</math> ( <math>t=t_n</math> ) to time <math displya="inline">n+2</math> ( <math>t=t_n+2\Delta t</math> )
65
|}
66
67
==3  FIC/FEM formulation for a Lagrangian continuum==
68
69
===3.1  Governing equations===
70
71
The equations to be solved are the standard ones in Lagrangian continuum mechanics <span id="citeF-24"></span>[[#cite-24|[24]]].
72
73
''Momentum''
74
<span id="eq-1"></span>
75
{| class="formulaSCP" style="width: 100%; text-align: center;" 
76
|-
77
| 
78
{| style="text-align: left; margin:auto;width: 100%;" 
79
|-
80
| style="text-align: center;" |<math>\rho \frac{\partial v_i}{\partial t}=\frac{\partial {\sigma }_{ij}}{\partial x_j}+</math><math>b_i\mbox{ }\mbox{ }\mbox{ }in\mbox{ }\mbox{ }V</math>
81
|}
82
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
83
|}
84
85
''Pressure-velocity relationship''
86
<span id="eq-2"></span>
87
{| class="formulaSCP" style="width: 100%; text-align: center;" 
88
|-
89
| 
90
{| style="text-align: left; margin:auto;width: 100%;" 
91
|-
92
| style="text-align: center;" |<math>\frac{1}{K}\frac{\partial p}{\partial t}-\frac{\partial v_i}{\partial x_i}=0\mbox{ }\mbox{ }\mbox{ }in\mbox{ }\mbox{ }V</math>
93
|}
94
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
95
|}
96
97
In above equations  <math>v_i</math> is the velocity along the <math>i</math>th global (cartesian) axis, <math>p</math> is the pressure (assumed to be positive in tension)  <math>\rho </math> and <math>K</math> are the density and bulk modulus of the material, respectively,  <math>b_i</math> and  <math>{\sigma }_{ij}</math> are the body forces and the (Cauchy) stresses. Eqns ([[#eq-1|1]]) and ([[#eq-2|2]]) are completed with the constitutive relationships <span id="citeF-12"></span>[[#cite-12|[12]]]:
98
99
'' Incompressible continuum''
100
101
<span id="eq-3"></span>
102
{| class="formulaSCP" style="width: 100%; text-align: center;" 
103
|-
104
| 
105
{| style="text-align: left; margin:auto;width: 100%;" 
106
|-
107
| style="text-align: center;" |<math>{}^{t+1}{\sigma }_{ij}=2\mu {\dot{\epsilon }}_{ij}+</math><math>^{t+1}p{\delta }_{ij}</math>
108
|}
109
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
110
|}
111
112
''Compressible/quasi-incompressible continuum''
113
<span id="eq-4"></span>
114
<span id="eq-4a"></span>
115
{| class="formulaSCP" style="width: 100%; text-align: center;" 
116
|-
117
| 
118
{| style="text-align: left; margin:auto;width: 100%;" 
119
|-
120
| style="text-align: center;" |<math>{}^{t+1}{\sigma }_{ij}=^t{\overset{\mbox{ˆ}}{\sigma }}_{ij}+</math><math>2\mu {\dot{\epsilon }}_{ij}+\lambda {\dot{\epsilon }}_{ii}{\delta }_{ij}</math>
121
|}
122
| style="width: 5px;text-align: right;white-space: nowrap;" | (4a)
123
|}
124
125
where  <math>{\overset{\mbox{ˆ}}{\sigma }}_{ij}</math> are the component of the Cauchy stress tensor  <math>\left[\overset{\mbox{ˆ}}{\sigma }\right]</math>
126
127
<span id="eq-4b"></span>
128
{| class="formulaSCP" style="width: 100%; text-align: center;" 
129
|-
130
| 
131
{| style="text-align: left; margin:auto;width: 100%;" 
132
|-
133
| style="text-align: center;" |<math>[\overset{\mbox{ˆ}}{\sigma }]=\frac{1}{J}\boldsymbol{F}^T\boldsymbol{SF}</math>
134
|}
135
| style="width: 5px;text-align: right;white-space: nowrap;" | (4b)
136
|}
137
138
where <math display="inline">S</math> is the second Piola-Kirchhoff stress tensor, <math display="inline">F</math> is the deformation gradient tensor and  <math display="inline">F</math> <span id="citeF-24"></span>[[#cite-24|[24]]]. Parameters  <math>\mu </math> and  <math>\lambda </math> take the following values for a fluid or solid material:
139
140
''Fluid    '' <math>\mu </math> : viscosity;  <math>\lambda =\Delta tK-\frac{2\mu }{3}</math>
141
142
''Solid    '' <math>\mu =\frac{\Delta tG}{J}</math> ;  <math>\lambda =\frac{2G\nu \Delta t}{J(1-2\nu )}</math> , where  <math>\nu </math> is the Poisson ratio, <math display="inline>G</math> is the shear modulus and  <math>\Delta t</math> the time increment.
143
144
In eqns ([[#eq-3|3]]) and ([[#eq-4|4]]),  <math>{\dot{\epsilon }}_{ij}</math> is the rate of deformation and  <math>{\delta }_{ij}</math> is the Kronecker delta.  <math>{}^t(\cdot )</math> denotes values at time <math display="inline">t</math>.
145
146
Note that <math>{}^t{\overset{\frown}{\sigma }}_{ij}=0</math> for a fluid in eqn.([[#eq-4a|4a]]), as the stresses depend on the rates of  deformation only.
147
148
Indexes in eqns ([[#eq-1|1]])-([[#eq-4|4]]) range from  <math>i,j=1,n_d</math> , where  <math>n_d</math> is the number of space dimensions. These equations are completed with the standard boundary conditions of prescribed velocities and surface tractions in the mechanical problem <span id="citeF-8"></span>[[#cite-8|[8]]-<span id="citeF-13"></span>[[#cite-13|13]]].
149
150
===3.2  Discretization of the equations===
151
152
A key problem in the numerical solution of eqns ([[#eq-1|1]])-([[#eq-4|4]]) is the satisfaction of the mass balance condition for the fully incompressible case (i.e.  <math>K=\infty </math> in eqn [[#eq-2|(2)]]). A number of procedures to solve this problem exists in the finite element literature <span id="citeF-25"></span>[[#cite-25|[25]],<span id="citeF-27"></span>[[#cite-27|27]]]. In our approach we use a stabilized formulation based in the so-called finite calculus (FIC) procedure <span id="citeF-9"></span>[[#cite-9|[9]],<span id="citeF-14"></span>[[#cite-14|14]],<span id="citeF-20"></span>[[#cite-20|20]]-<span id="citeF-23"></span>[[#cite-23|23]]]. The essence of this method is the solution of a ''modified mass balance'' equation which is written as
153
154
<span id="eq-5"></span>
155
{| class="formulaSCP" style="width: 100%; text-align: center;" 
156
|-
157
| 
158
{| style="text-align: left; margin:auto;width: 100%;" 
159
|-
160
| style="text-align: center;" |<math>\frac{1}{K}\frac{\partial p}{\partial t}-\frac{\partial v_i}{\partial x_i}+</math><math>\sum_{i=1}^3\tau \frac{\partial q}{\partial x_i}\left[\frac{\partial p}{\partial x_i}+\right. </math><math>\left. {\pi }_i\right]=0</math>
161
|}
162
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
163
|}
164
165
where <math display="inline">q</math> are weighting functions,  <math>\tau </math> is a stabilization parameter given by <span id="citeF-9"></span>[[#cite-9|[9]],<span id="citeF-14"></span>[[#cite-14|14]],<span id="citeF-21"></span>[[#cite-21|21]]-<span id="citeF-23"></span>[[#cite-23|23]]]
166
167
<span id="eq-6"></span>
168
{| class="formulaSCP" style="width: 100%; text-align: center;" 
169
|-
170
| 
171
{| style="text-align: left; margin:auto;width: 100%;" 
172
|-
173
| style="text-align: center;" |<math>\tau ={\left(\frac{2\rho \vert v\vert }{h}+\frac{8\mu }{3h^2}\right)}^{-1}</math>
174
|}
175
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
176
|}
177
178
In the above, <math display="inline">h</math> is a characteristic length of each finite element and  <math>\vert v\vert </math> is the modulus of the velocity vector. In eqn ([[#eq-5|5]])  <math>{\pi }_i</math> are auxiliary pressure projection variables chosen so as to ensure that the second term in eqn ([[#eq-5|5]]) can be interpreted as weighted sum of the residuals of the momentum equations and therefore it vanishes for the exact solution. The set of governing equations is completed by adding the following constraint equation <span id="citeF-9"></span>[[#cite-9|[9]]].
179
180
<span id="eq-7"></span>
181
{| class="formulaSCP" style="width: 100%; text-align: center;" 
182
|-
183
| 
184
{| style="text-align: left; margin:auto;width: 100%;" 
185
|-
186
| style="text-align: center;" |<math>{\int }_V\tau w_i\left(\frac{\partial p}{\partial x_i}+\right. </math><math>\left. {\pi }_i\right)dV=0\mbox{ }i=1,n_d\mbox{ }\mbox{ }(\mbox{no}\mbox{ }\mbox{sum}\mbox{ }\mbox{in}\mbox{ }\mbox{ }i)</math>
187
|}
188
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
189
|}
190
191
where  <math>w_i</math> are arbitrary weighting functions.
192
193
The rest of the integral equations are obtained by applying the weighted residual technique to the governing eqns ([[#eq-1|1]]), ([[#eq-2|2]]), ([[#eq-3|3]]) and ([[#eq-5|5]]) and the boundary conditions <span id="citeF-9"></span>[[#cite-9|[9]],<span id="citeF-12"></span>[[#cite-12|12]],<span id="citeF-14"></span>[[#cite-14|14]],<span id="citeF-18"></span>[[#cite-18|18]],<span id="citeF-22"></span>[[#cite-22|22]],<span id="citeF-26"></span>[[#cite-26|26]]].
194
195
We interpolate next in the finite element fashion the set of problem variables. For 3D problems these are the three velocities  <math>v_i</math> , the pressure <math>p</math> and the three pressure gradient projections  <math>{\pi }_i.</math> In our work we use equal order ''linear interpolation'' for all variables over meshes of 3-noded triangles (in 2D) and 4-noded tetrahedra (in 3D). The resulting set of discretized equations using the standard Galerkin technique has the following form
196
197
''Momentum''
198
199
<span id="eq-8"></span>
200
{| class="formulaSCP" style="width: 100%; text-align: center;" 
201
|-
202
| 
203
{| style="text-align: left; margin:auto;width: 100%;" 
204
|-
205
| style="text-align: center;" |<math>\boldsymbol{M\dot{\overline{v}}}+\boldsymbol{K\overline{v}}+</math><math>\boldsymbol{G\overline{p}}=\boldsymbol{f}</math>
206
|}
207
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
208
|}
209
210
''Pressure-velocity relationship''
211
212
<span id="eq-9"></span>
213
{| class="formulaSCP" style="width: 100%; text-align: center;" 
214
|-
215
| 
216
{| style="text-align: left; margin:auto;width: 100%;" 
217
|-
218
| style="text-align: center;" |<math>\boldsymbol{\overline{M}\dot{\overline{p}}}+\boldsymbol{L\overline{p}}-</math><math>\boldsymbol{G\overline{v}}+\boldsymbol{Q}\mbox{ }\boldsymbol{\overline{\pi }}\mbox{ }=0</math>
219
|}
220
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
221
|}
222
223
''Pressure gradient projection  ''
224
225
<span id="eq-10"></span>
226
{| class="formulaSCP" style="width: 100%; text-align: center;" 
227
|-
228
| 
229
{| style="text-align: left; margin:auto;width: 100%;" 
230
|-
231
| style="text-align: center;" |<math>\boldsymbol{\overset{\frown}{M}\overline{\pi }}\mbox{ }+</math><math>\boldsymbol{Q}^T\boldsymbol{\overline{p}}=\boldsymbol{0}</math>
232
|}
233
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
234
|}
235
236
In eqns ([[#eq-8|8]])-([[#eq-10|10]]) <math>\overline{\left(\cdot \right)}</math> denotes nodal variables,  <math display="inline">\overset{\cdot}{\overline{\left( \cdot \right) }}=</math><math>\frac{\partial }{\partial t}\overline{(\cdot )}</math>. The matrices and vectors are given in <span id="citeF-14"></span>[[#cite-14|[14]],<span id="citeF-18"></span>[[#cite-18|18]]].
237
238
The solution in time of eqns ([[#eq-8|8]])-([[#eq-10|10]]) can be performed using any time integration scheme typical of the updated Lagrangian FEM <span id="citeF-24"></span>[[#cite-24|[24]]]. A basic algorithm following the conceptual process described in Section [[#2 The basis of the particle finite element method|2]] is presented in Box [[#Box-1|I]].
239
240
<span id="Box-1"></span>
241
'''Box I'''. Basic PFEM algorithm for a Lagrangian continuum
242
243
{|  class="floating_tableSCP wikitable" style="text-align: right; margin: 1em auto;min-width:50%;"
244
245
|- style="border-top: 2px solid;"
246
| style="text-align: left;border-left: 2px solid;border-right: 2px solid;" | 1. LOOP OVER TIME STEPS, <math display="inline">t=1</math>, NTIME
247
248
Known values
249
250
<math display="inline">^{t} \bar{\boldsymbol x},{}^{t} \bar{\boldsymbol v},{}^{t} \bar{\boldsymbol p},{}^{t} \bar{\boldsymbol \pi },{}^{t} \bar{T},{}^{t} \mu ,{}^{t}{\boldsymbol f},{}^{t}\mathbf{ q},{}^{t} C,{}^{t} V,{}^{t} M</math> 
251
252
<span id='section-2'></span>
253
2. LOOP OVER NUMBER OF ITERATIONS, <math display="inline">i=1</math>, NITER 
254
255
<math display="inline">\bullet </math>  Compute  nodal velocities by solving Eq.([[#eq-8|8]])
256
257
<math display="inline">\displaystyle \left[\frac{1}{\Delta t} \mathbf{M}+\mathbf{K}\right]{}^{t+1} \bar{\boldsymbol v}^{i+1} ={}^{t+1} \mathbf{f} - {\boldsymbol G}^{t+1} {\boldsymbol p}^{i} +\frac{1}{\Delta t} \mathbf{M} {}^{t} \bar{\boldsymbol v}</math>
258
259
<math display="inline">\bullet </math>   Compute nodal pressures from Eq.([[#eq-9|9]])
260
261
<math display="inline">\displaystyle \left[\frac{1}{\Delta t}\bar {\boldsymbol M}-{\boldsymbol L}\right]{}^{t+1} \bar{\boldsymbol p}^{i+1} =\mathbf{G}{}^{T} {}^{t+1} \bar{\boldsymbol v}^{i+1} + \mathbf{Q}{}^{t+1} \bar{\boldsymbol \pi }^{i} +\frac{1}{\Delta t}\bar {\boldsymbol M} {}^t\bar{\boldsymbol p} </math> 
262
263
<math display="inline">\bullet </math> Compute nodal pressure gradient projections from Eq.([[#eq-10|10]])
264
265
<math display="inline">{}^{n+1} \bar{\boldsymbol \pi }^{i+1} =-\hat{\boldsymbol M}^{-1}_{D} \left[\mathbf{Q}^{T} \right]{}^{t+1} \bar{\boldsymbol p}^{i+1} \begin{array}{ccc} {} & {,} & {\hat{\boldsymbol M}_{D} =diag\left[\hat{\boldsymbol M}_{D} \right]} \end{array}</math> 
266
267
<math display="inline">\bullet </math>   Update position of analysis domain nodes:
268
269
<math display="inline">{}^{t+\Delta t} \bar{\boldsymbol x}^{i+1} ={}^{t} \mathbf{x}^{i} +{}^{t+\Delta t} \mathbf{v}^{i+1} \Delta t</math>
270
271
Define new “cloud” of nodes <math display="inline">{}^{t+1} C^{i+1} </math> 
272
273
<math display="inline">\bullet </math>   Update strain rate and strain values
274
275
<math display="inline">\bullet </math>   Update stress values 
276
277
Check convergence <math display="inline">\rightarrow </math> NO <math display="inline">\rightarrow </math> Next iteration <math display="inline">i\to i+1</math>
278
279
<math display="inline">\downarrow </math> YES
280
281
Next time step <math display="inline">t\to t+1</math> 
282
283
<math display="inline">\bullet </math> Identify new analysis domain boundary: <math display="inline">{}^{t+1} V</math> 
284
285
<math display="inline">\bullet </math> Generate mesh:<math display="inline">{}^{t+1} M</math> 
286
287
Go to 1
288
|- style="border-bottom: 2px solid;"
289
| style="text-align: left;border-left: 2px solid;border-right: 2px solid;"|
290
|}
291
292
==4  Generation of a new mesh==
293
294
A key point for the success of the PFEM is the fast regeneration of a mesh at every time step on the basis of the position of the nodes in the space domain. In our work the mesh is generated  using the so called extended Delaunay tesselation (EDT) presented in <span id="citeF-4"></span>[[#cite-4|[4]]].
295
296
The CPU time required for meshing grows linearly with the number of nodes. The CPU time for solving the equations exceeds that required for meshing as the number of nodes increases. As a general rule for large 3D problems meshing consumes around <math display="inline">15%</math> of the total CPU time per time step, while the solution of the equations (with typically 3 iterations per time step) and the system assembly consume approximately <math display="inline">70%</math> and <math display="inline">15%</math> of the CPU time per time step, respectively. These figures refer to analyses in a single processor Pentium IV PC and prove that the generation of the mesh has an acceptable cost in the PFEM. Indeed considerable speed can be gained using parallel computing techniques.
297
298
==5  Identification of boundary surfaces==
299
300
One of the main tasks in the PFEM is the correct definition of the boundary domain. Boundary nodes are sometimes explicitly identified. In other cases, the total set of nodes is the only information available and the algorithm must recognize the boundary nodes (Figure [[#img-2|2]]).
301
302
In our work we use an extended Delaunay partition for recognizing boundary nodes <span id="citeF-4"></span>[[#cite-4|[4]]][4]. Considering that the nodes follow a variable <math displya="inline">h(x)</math> distribution, where <math displya="inline">h(x)</math> is typically the minimum distance between two nodes. '' All nodes on an empty sphere with a radius greater than  <math>\alpha h</math>, are considered as boundary nodes''. In practice  <math>\alpha </math> is a parameter close to, but greater than one. Values of <math>\alpha </math> ranging around 1.3 have been found to be optimal in all examples analyzed. This criterion is coincident with the Alpha Shape concept <span id="citeF-19"></span>[[#cite-19|[19]]].
303
304
Once a decision has been made concerning which nodes are on the boundaries, the boundary surface is defined by all the polyhedral surfaces (or polygons in 2D) having all their nodes on the boundary and belonging to just one polyhedron.
305
306
The method also allows one to identify isolated fluid particles outside the main fluid domain. These particles are treated as part of the external boundary where the pressure is fixed to the atmospheric value. We recall that each particle is a material point characterized by the density of the solid or fluid domain to which it belongs. The mass lost when a boundary element is eliminated due to departure of a node from the analysis domain is regained when the  node falls down and a new boundary element is created by the Alpha Shape algorithm.
307
308
The boundary recognition method is useful for detecting contact conditions between the fluid domain and a boundary, as well as between different solids as detailed in the next section.
309
310
We emphasize that the key differences between the PFEM and the classical FEM are the remeshing technique and the identification of the domain boundary at each time step.
311
312
<div id='img-2'></div>
313
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
314
|-
315
|[[Image:Draft_Samper_913569722-image67.png|414px]]
316
|- style="text-align: center; font-size: 75%;"
317
| colspan="1" | '''Figure 2:'''Modelling of contact conditions at a solid-solid interface with the PFEM.
318
|}
319
320
==6  Treatment of contact conditions in the PFEM==
321
322
===6.1  Contact between the fluid and a fixed boundary===
323
324
Known velocities at boundaries in the PFEM are prescribed in strong form to the boundary nodes. These nodes might belong to fixed external boundaries or to moving boundaries linked to the interacting solids. Contact between fluid particles and fixed boundaries is accounted for by the incompressibility condition which  ''naturally prevents fluid nodes to penetrate into the solid boundaries'' <span id="citeF-6"></span>[[#cite-6|[6]],<span id="citeF-9"></span>[[#cite-9|9]],<span id="citeF-14"></span>[[#cite-14|14]]].
325
326
===6.2  Contact between solid-solid interfaces===
327
328
The contact between two solid interfaces is treated by introducing a layer of '' contact elements'' between the two interacting solid interfaces. This layer is ''automatically created during the mesh generation step'' by prescribing a minimum distance <math>\left(h_c\right)</math> between two solid boundaries. If the distance exceeds the minimum value <math>\left(h_c\right)</math> then the generated elements are treated as fluid elements. Otherwise the elements are treated as contact elements where a relationship between the tangential and normal forces and the corresponding displacement is introduced (Figure [[#img-2|2]]) <span id="citeF-6"></span>[[#cite-x|[6]],<span id="citeF-10"></span>[[#cite-10|10]]].
329
330
This algorithm allows us to identify and model complex frictional contact conditions between two or more interacting bodies moving in water in an extremely simple manner. The algorithm can also be used effectively to model frictional contact conditions between rigid or elastic solids in structural mechanics applications <span id="citeF-16"></span>[[#cite-16|[16]]].
331
332
==7  Modeling of bed erosion==
333
334
Prediction of bed erosion and sediment transport in open channel flows are important tasks in river and environmental engineering. Bed erosion can lead to instabilities of the river basin slopes. It can also undermine the foundation of bridge piles thereby favouring structural failure. Modeling of bed erosion is also relevant for predicting the evolution of surface material dragged in earth dams in overspill situations. Bed erosion is one of the main causes of environmental damage in floods.
335
336
In recent works we have proposed an extension of the PFEM to model bed erosion <span id="citeF-10"></span>[[#cite-10|[10]],<span id="citeF-14"></span>[[#cite-14|14]]]. The erosion model is based on the frictional work at the bed surface originated by the shear stresses in the fluid. The resulting erosion model resembles Archard law typically used for modeling abrasive wear in surfaces under frictional contact conditions <span id="citeF-28"></span>[[#cite-28|[28]]].
337
338
The algorithm for modeling bed erosion  is the following:
339
340
# Compute at  the bed surface the resultant tangential stress induced by the fluid motion.
341
# Compute the frictional work <math>{}^nW_f</math> originated by the tangential stresses at the bed surface.
342
# The onset of erosion at a bed point occurs when  <math>{}^nW_f</math> exceeds a critical threshold value  <math>W_c</math> .
343
# If  <math>{}^nW_f>W_c</math> at a bed node, then the node is detached from the bed region and it is allowed to move with the fluid. Also, the mass of the patch of bed elements surrounding the bed node is transferred to the new fluid node. This mass is subsequently transported with the fluid.
344
345
Figure [[#img-3|3]] shows an schematic view of the bed erosion algorithm described.
346
347
<div id='img-3'></div>
348
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
349
|-
350
| [[Image:Draft_Samper_913569722-image73.png|390px]]
351
|- style="text-align: center; font-size: 75%;"
352
| colspan="x" | '''Figure 3:'''Modeling of bed erosion with the PFEM by dragging of bed material.
353
|}
354
355
Sediment deposition can be modeled by an inverse process. Hence, a suspended node adjacent to the bed surface with a velocity below a threshold value is attached to the bed surface.
356
357
Examples of the bed erosion algorithm for modeling excavation and rock cutting problems are presented in <span id="citeF-16"></span>[[#cite-16|[16]]].
358
359
==8  Examples==
360
361
===8.1  Dragging of rocks by a water stream===
362
363
Predicting the critical speed at which a rock will be dragged by a water stream is of great importance in many problems in hydraulic, harbour, civil and environmental engineering.
364
365
The PFEM has been successfully applied to the study of the motion of a 1Tn quasi-spherical rock due to a water stream. The rock lays on a collection of rocks that are kept rigid.
366
367
Frictional conditions between the analyzed rock and the rest of the rocks have been assumed.
368
369
Figure [[#img-4|4a]] shows that a water stream of 1m/s is not able to displace the individual rock. An increase of the water speed to 2m/s induces the motion of the rock as shown in Figure [[#img-4|4b]].
370
371
<div id='img-4'></div>
372
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
373
|-
374
|[[Image:Draft_Samper_913569722-image74.png|258px]]
375
|[[Image:Draft_Samper_913569722-image75.png|258px]]
376
|- style="text-align: center; font-size: 75%;"
377
|colspan="1"|(a) Water speed of 1m/s. The individual rock can not be dragged by the stream
378
|colspan="1"|(b) Water speed of 2m/s. The individual rock is dragged by the stream
379
|- style="text-align: center; font-size: 75%;"
380
| colspan="2" | '''Figure 4''': Drag of a 1 Tn rock under a water stream at speeds of a) 1m/s and b) 2m/s.
381
|}
382
383
===8.2  Impact of sea waves on piers and breakwaters===
384
385
Figure [[#img-5|5]] shows the analysis of the effect of breaking waves on two different sites of a breakwater containing reinforced concrete blocks (each one of <math display="inline">4x4x4</math> mts). The figures correspond to the study of Langosteira harbour in A Coruña, Spain using PFEM.
386
387
<div id='img-5'></div>
388
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
389
|-
390
|[[Image:Draft_Samper_913569722-image76-c.png|300px]]
391
|[[File:Draft_Samper_913569722_6528_image77-c.png]]
392
|- style="text-align: center; font-size: 75%;"
393
| colspan="2" | '''Figure 5:''' Breaking waves on breakwater slopes containing reinforced concrete blocks.
394
|}
395
396
===8.3  Soil erosion problems===
397
398
Figure [[#img-6|6a]] shows the capacity of the PFEM for modelling soil erosion, sediment transport and material deposition in a river bed. The soil particles are first detached from the bed surface under the action of the jet stream. Then they are transported by the flow and eventually fall down due to gravity forces into the bed surface at a downstream point.
399
400
Figure [[#img-6|6b]] shows the progressive erosion of the unprotected part of a breakwater slope in the Langosteira harbour in A Coruña, Spain. The non protected upper shoulder zone is progressively eroded as it is hit by the sea waves.
401
402
403
<div id='img-6'></div>
404
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
405
|-
406
|[[Image:Draft_Samper_913569722-image78-c.png|382px]]
407
|[[Image:Draft_Samper_913569722-image79-c.png|294px]]
408
|- style="text-align: center; font-size: 75%;"
409
| colspan="1" | (a)
410
| colspan="1" | (b)
411
|- style="text-align: center; font-size: 75%;"
412
| colspan="2" | '''Figure 6''': (a) Erosion, transport and deposition of soil particles at a river bed due to an impacting jet stream (b) Erosion of an unprotected shoulder of a breakwater due to sea waves
413
|}
414
415
===8.4  Falling of a lorry into the sea by sea wave erosion of the road slope===
416
417
Figure [[#img-7|7]] shows a representative example of the progressive erosion of a soil mass adjacent to the shore due to sea waves and the subsequent falling into the sea of a 2D object representing the section of a lorry. The object has been modeled as a rigid solid.
418
419
This example and the previous ones, although still quite simple and schematic, show the possibilities of the PFEM for modeling complex FSSI problems involving soil erosion, free surface waves and rigid/deformable structures.
420
421
<div id='img-7'></div>
422
{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
423
|-
424
| [[Image:Draft_Samper_913569722-image80-c.png|300px]]
425
| [[Image:Draft_Samper_913569722-image81-c.png|300px]]
426
| [[Image:Draft_Samper_913569722-image82-c.png|300pc]]
427
|- style="text-align: center; font-size: 75%;"
428
| colspan="3" | '''Figure 7:''' Erosion of a soil mass due to sea waves and the subsequent falling into the sea operating in a road adjacent to the sea lorry
429
|}
430
431
432
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
433
 
434
{|
435
|-
436
 </div>
437
438
Figure 7:
439
440
===8.5  Impact of waves over transport vehicles in harbour===
441
442
Figures 8 and 9 show two examples of the study of the impact of large waves on lorries operating in roads adjacent to a breakwater.
443
444
Figure 8 shows the effect of a wave of 20mts amplitude on a 20 Tns lorry adjacent to a breakwater shoulder.
445
446
Figure 9 displays the impact of an overtopping wave on the same lorry placed on a road behing a breakwater.
447
448
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
449
 
450
{|
451
|-
452
| [[Image:Draft_Samper_913569722-image83.png|240px]]
453
| [[Image:Draft_Samper_913569722-image84.png|center|240px]]
454
|}
455
</div>
456
457
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
458
 
459
{|
460
|-
461
| [[Image:Draft_Samper_913569722-image85.png|240px]]
462
| [[Image:Draft_Samper_913569722-image86-c.png|center|222px]]
463
|}
464
</div>
465
466
Figure 8: 20mts amplitude wave hitting a 20Tn lorry placed close to a breakwater slope.
467
468
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
469
 
470
{|
471
|-
472
| [[Image:Draft_Samper_913569722-image87-c.png|228px]]
473
| [[Image:Draft_Samper_913569722-image88-c.png|center|222px]]
474
|}
475
</div>
476
477
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
478
 
479
{|
480
|-
481
| [[Image:Draft_Samper_913569722-image89-c.png|222px]]
482
| [[Image:Draft_Samper_913569722-image90-c.png|center|252px]]
483
|}
484
</div>
485
486
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
487
<span style="text-align: center; font-size: 75%;">Figure 9: Effect of an overtopping wave on a 20Tn lorry in a road adjacent to a breakwater.</span></div>
488
489
===8.6 Simulation of sinking of ships===
490
491
The PFEM can be effectively applied for simulating the sinking of ships under a variety of scenarios.
492
493
Figure 10 shows images of  the 2D simulation of the sinking of a cargo vessel induced by a breach in the bow region.
494
495
Figure 11 displays a 3D simulation of the skinking of a simple fisherman boat induced by a hole in the side of the hull.
496
497
These examples evidence the potential of PFEM for the study of the sinking of ships.
498
499
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
500
<span style="text-align: center; font-size: 75%;">(a)</span></div>
501
502
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
503
<span style="text-align: center; font-size: 75%;"> 
504
{|
505
|-
506
| [[Image:Draft_Samper_913569722-image91.png|246px]]
507
| [[Image:Draft_Samper_913569722-image92.png|center|246px]]
508
|}
509
</span></div>
510
511
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
512
<span style="text-align: center; font-size: 75%;"> 
513
{|
514
|-
515
| [[Image:Draft_Samper_913569722-image93.png|246px]]
516
| [[Image:Draft_Samper_913569722-image94.png|center|252px]]
517
|}
518
</span></div>
519
520
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
521
<span style="text-align: center; font-size: 75%;">(b)</span></div>
522
523
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
524
<span style="text-align: center; font-size: 75%;"> 
525
{|
526
|-
527
| [[Image:Draft_Samper_913569722-image95.png|246px]]
528
| [[Image:Draft_Samper_913569722-image96.png|center|252px]]
529
|}
530
</span></div>
531
532
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
533
<span style="text-align: center; font-size: 75%;"> 
534
{|
535
|-
536
| [[Image:Draft_Samper_913569722-image97.png|246px]]
537
| [[Image:Draft_Samper_913569722-image98.png|center|246px]]
538
|}
539
</span></div>
540
541
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
542
<span style="text-align: center; font-size: 75%;">Figure 10: 2D simulation of the sinking of a cargo vessel due to a breach in the bow region. (a) Water streamline at different times. (b) Water velocity pattern at different times during sinking.</span></div>
543
544
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
545
 
546
{|
547
|-
548
| [[Image:Draft_Samper_913569722-image99.png|192px]]
549
| [[Image:Draft_Samper_913569722-image100.png|center|192px]]
550
|}
551
 [[Image:Draft_Samper_913569722-image101.png|192px]] </div>
552
553
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
554
 
555
{|
556
|-
557
| [[Image:Draft_Samper_913569722-image102-c.png|186px]]
558
| [[Image:Draft_Samper_913569722-image103-c.png|center|186px]]
559
|}
560
 [[Image:Draft_Samper_913569722-image104-c.png|186px]] </div>
561
562
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
563
<span style="text-align: center; font-size: 75%;">Figure 11: 3D simulation of the sinking of a boat induced by a hole in the side of the hull.  Figures show different views of the water flow inside the boat during sinking.</span></div>
564
565
===8.7 Colision of boat with ice blocks===
566
567
Figures 12 and 13 show two examples of the application of PFEM to the study of the collision of a ship with ice blocks.
568
569
Figure 12 shows snapshots of the motion of a boat which collides with several ice blocks.
570
571
Figure 13 displays the interaction between a boat and two ice slabs that trap the boat in their motion.
572
573
We note that the boat and the ice blocks have been modelled as rigid bodies in these examples. Indeed, the deformation of the ship strucutre due to the ice-ship interaction forces cand be accounted for in the analysis.
574
575
576
{|
577
|-
578
| [[Image:Draft_Samper_913569722-image105-c.png|156px]]
579
| [[Image:Draft_Samper_913569722-image106-c.png|center|150px]]
580
|}
581
 
582
{|
583
|-
584
| [[Image:Draft_Samper_913569722-image107-c.png|150px]]
585
| [[Image:Draft_Samper_913569722-image108-c.png|center|144px]]
586
|}
587
588
589
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
590
<span style="text-align: center; font-size: 75%;">Figure 12: 3D simulation of a boat colliding with five ice blocks.</span></div>
591
592
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
593
 
594
{|
595
|-
596
| [[Image:Draft_Samper_913569722-image109.png|216px]]
597
| [[Image:Draft_Samper_913569722-image110.png|center|216px]]
598
|}
599
</div>
600
601
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
602
 
603
{|
604
|-
605
| [[Image:Draft_Samper_913569722-image111.png|216px]]
606
| [[Image:Draft_Samper_913569722-image112.png|center|222px]]
607
|}
608
</div>
609
610
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
611
<span style="text-align: center; font-size: 75%;">Figure 13: Simulation of the interaction of two adjacent ice slabs and boat. The ice slabs move towards the boat that ends up out of the water and over the slabs.</span></div>
612
613
<span id='GrindEQpgref5049b1f220'></span>
614
==9  Conclusions==
615
616
The particle finite element method (PFEM) is a promising numerical technique for solving fluid-soil-structure interaction (FSSI) problems in naval, marine and harbour engineering involving large motion of fluid and solid particles, surface waves, water splashing, frictional contact situations between fluid-solid and solid-solid interfaces and bed erosion, among other complex phenomena. The success of the PFEM lies in the accurate and efficient solution of the equations of an incompressible continuum using an updated Lagrangian formulation and a stabilized finite element method allowing the use of low order elements with equal order interpolation for all the variables. Other essential solution ingredients are the efficient regeneration of the finite element mesh, the identification of the boundary nodes using the Alpha-Shape technique and the simple algorithm to treat frictional contact conditions and erosion/wear at fluid-solid and solid-solid interfaces via mesh generation. The examples presented have shown the potential of the PFEM for solving a wide class of practical FSSI problems in naval, marine and harbour engineering.
617
618
==10 Acknowledgements==
619
620
This research was partially supported by projects SAFECON and REALTIME of the European Research Council of the European Commission (EC).
621
622
==11 References==
623
624
<div id="cite-1"></div>
625
[[#citeF-1|[1]]]   Idelsohn, S.R., Oñate, E., Del Pin, F. & Calvo, N., ''Lagrangian formulation: the only way to solve some free-surface fluid mechanics problems'', 5th World Congress on Comput. Mechanics, H.A. Mang, F.G. Rammerstorfer & J. Eberhardsteiner (Eds), July 7-12, Viena, Austria, 2002.
626
627
<div id="cite-2"></div>
628
[[#citeF-2|[2]]]   Idelsohn, S.R., Oñate, E. & Del Pin, F., A lagrangian meshless finite element method applied to fluid-structure interaction problems, ''Comput. and Struct''., 81, pp. 655-671, 2003.
629
630
<div id="cite-3"></div>
631
[[#citeF-3|[3]]] Idelsohn, S.R., Oñate, E., Calvo, N. & Del Pin, F., The meshless finite element method, ''Int. J. Num. Meth. Engng''., '''58(6)''', pp. 893-912, 2003.
632
633
<div id="cite-4"></div>
634
[[#citeF-4|[4]]]  Idelsohn, S.R., Calvo, N. & Oñate, E., Polyhedrization of an arbitrary point set, ''Comput. Method Appl. Mech. Engng''., '''192(22-24)''', pp. 2649-2668, 2003.
635
636
<div id="cite-5"></div>
637
[[#citeF-5|[5]]]  Idelsohn,  S.R., Oñate, E. & Del Pin, F., The particle finite element method: a powerful tool to solve incompressible flows with free-surfaces and breaking waves, ''Int. J. Num. Meth. Engng'',. '''61''', pp. 964-989, 2004.
638
639
<div id="cite-6"></div>
640
[[#citeF-6|[6]]]  Oñate, E., Idelsohn, S.R., Del Pin, F. & Aubry, R., The particle finite element method. An overview, ''Int. J. Comput. Methods'', '''1(2)''', pp. 267-307, 2004b.
641
642
<div id="cite-7"></div>
643
[[#citeF-7|[7]]]   Aubry, R., Idelsohn, S.R. & Oñate, E., Particle finite element method in fluid mechanics including thermal convection-diffusion, ''Computer & Structures'', '''83(17-18)''', pp. 1459-1475, 2005.
644
645
<div id="cite-8"></div>
646
[[#citeF-8|[8]]]  Idelsohn, S.R., Oñate, E., Del Pin, F. & Calvo, N., Fluid-structure interaction using the particle finite element method, ''Comput. Meth. Appl. Mech. Engng''., '''195''', pp. 2100-2113, 2006.
647
648
<div id="cite-9"></div>
649
[[#citeF-9|[9]]]   Oñate, E., García, J., Idelsohn, S.R. & Del Pin, F., FIC formulations for finite element analysis of incompressible flows. Eulerian, ALE and Lagrangian approaches, ''Comput. Meth. Appl. Mech. Engng''., '''195(23-24)''', pp. 3001-3037, 2006.
650
651
<div id="cite-10"></div>
652
[[#citeF-10|[10]]]  Oñate, E., Celigueta, M.A. & Idelsohn, S.R., Modeling bed erosion in free surface flows by the Particle Finite Element Method, ''Acta Geotechnia'', '''1(4)''', 237-252, 2006.
653
654
<div id="cite-11"></div>
655
[[#citeF-11|[11]]]  Del Pin, F., Idelsohn, S.R., Oñate, E. & Aubry, R., The ALE/Lagrangian particle finite element method: A new approach to computation of free-surface flows and fluid-object interactions, ''Computers & Fluids'', '''36''', pp. 27-38, 2007.
656
657
<div id="cite-12"></div>
658
[[#citeF-11|[11]]]  Idelsohn, S.R., Marti, J., Limache, A. & Oñate, E. Unified Lagrangian formulation for elastic solids and incompressible fluids: Application to fluid-structure interaction problems via the PFEM, ''Comput Methods Appl Mech Engrg''., '''197''', pp. 1762-1776, 2008.
659
660
<div id="cite-13"></div>
661
[[#citeF-13|[13]]] Larese, A., Rossi, R., Oñate, E. & Idelsohn, S.R., Validation of the Particle Finite Element Method (PFEM) for  free surface flows,  ''Engng. Computations'',  '''25(4)''', pp. 385-425, 2008.
662
663
<div id="cite-14"></div>
664
[[#citeF-14|[14]]]  Oñate, E., Idelsohn, S.R., Celigueta, M.A. & Rossi, R., Advances in the particle finite element method for the analysis of fluid-multibody interaction and bed erosion in free surface flows, ''Comput. Meth. Appl. Mech. Engng''., '''197(19-20)''', 1777-1800, 2008.
665
666
<div id="cite-15"></div>
667
[[#citeF-15|[15]]]  Idelsohn, S.R., Mier-Torrecilla, M. & Oñate, E., Multi-fluid flows with the Particle Finite Element Method, ''Comput Methods Appl Mech Engrg''., '''198''', pp. 2750-2767, 2009.
668
669
<div id="cite-16"></div>
670
[[#citeF-16|[16]]]  Carbonell, J.M,. Oñate, E. & Suárez, B., Modeling of ground excavation with the Particle Finite Element Method,  ''J. of Engineering Mechanics (ASCE)'',  '''136(4)''', pp. 455- 463, 2010.
671
672
<div id="cite-17"></div>
673
[[#citeF-17|[17]]]  Oñate, E., Rossi, R., Idelsohn, S.R. & Butler, K., Melting and spread of polymers in fire with the particle finite element method, ''Int. J. Numerical Methods in Engng''., '''81(8)''', 1046-1072, 2010.
674
675
<div id="cite-18"></div>
676
[[#citeF-18|[18]]] Oñate, E., Celigueta, M.A., Idelsohn, S.R., Salazar, F. & Suárez B., Possibilities of the particle finite element method for fluid–soil–structure interaction problems, ''Comput. Mech''., '''48''', pp. 307-318, 2011.
677
678
<div id="cite-19"></div>
679
[[#citeF-19|[19]]]  Edelsbrunner, H. & Mucke, E.P., Three dimensional alpha shapes, ''ACM Trans. Graphics'' '''13''', pp. 43-72, 1999.
680
681
<div id="cite-20"></div>
682
[[#citeF-20|[20]]]  Oñate, E., Derivation of stabilized equations for advective-diffusive transport and fluid flow problems, ''Comput. Meth. Appl. Mech. Engng''., '''151''', pp. 233-267, 1998.
683
<div id="cite-21"></div>
684
[[#citeF-21|[21]]]  Oñate, E., Possibilities of finite calculus in computational mechanics, ''Int. J. Num. Meth. Engng''. '''60(1)''', pp. 255-281, 2004.
685
686
<div id="cite-22"></div>
687
[[#citeF-22|[22]]]  Oñate, E. & García, J., A finite element method for fluid-structure interaction with surface waves using a finite calculus formulation, ''Comp. Meth. Appl. Mech. Eng''., '''191''', pp. 635-660, 2001.
688
689
<div id="cite-23"></div>
690
[[#citeF-23|[23]]]  Oñate, E., Valls, A. & García, J., FIC/FEM formulation with matrix stabilizing terms for incompressible flows at low and high Reynold's numbers, ''Comput. Mech''., '''38 (4-5)''', pp. 440-455, 2006.
691
692
<div id="cite-24"></div>
693
[[#citeF-24|[24]]]  Zienkiewicz, O.C. & Taylor, R.L.,  ''The finite element method for solid and structural mechanics'', Elsevier, 2005.
694
695
<div id="cite-25"></div>
696
[[#citeF-25|[25]]]  Donea, J. & Huerta, A., ''Finite element method for flow problems'', J. Wiley, 2003.
697
698
<div id="cite-26"></div>
699
[[#citeF-26|[26]]]  Zienkiewicz, O.C., Taylor, R.L. & Zhu, J.Z., ''The finite element method. Its basis and fundamentals'', Elsevier, 2005.
700
701
<div id="cite-27"></div>
702
[[#citeF-27|[27]]]  Zienkiewicz, O.C., Taylor, R.L. & Nithiarasu, P., ''The finite element method for fluid dynamics'', Elsevier, 2006.
703
704
<div id="cite-28"></div>
705
[[#citeF-28|[28]]]  Archard, J.F., Contact and rubbing of flat surfaces,  ''J. Appl. Phys.'' '''24(8)''', pp.981-988, 1953.
706

Return to Onate et al 2013a.

Back to Top