Published in Comput. Meth. Appl. Mech. Engng., Vol. 192 (28-30), pp. 3043-3059, 2003
doi: 10.1016/S0045-7825(03)00340-2
Abstract
The paper introduces a general procedure for computational analysis of a wide class of multiscale problems in mechanics using a finite calculus (FIC) formulation. The FIC approach is based in expressing the governing equations in mechanics accepting that the domain where the standard balance laws are established has a finite size. This introduces naturally additional terms into the classical equations of infinitesimal theory in mechanics which are useful for the numerical solution of problems involving different scales in the physical parameters. The discrete nodal values obtained with the FIC formulation and the finite element method (FEM) can be effectively used as the starting point for obtaining a more refined solution in zones where high gradients of the relevant variables occur using hierarchical or enriched FEM. Typical multiscale problems in mechanics which can be solved with the FIC method include convection-diffusion-reaction problems with high localized gradients, incompressible problems in solid and fluid mechanics, localization problems such as prediction of shear bands in solids and shock waves in compressible fluids, turbulence, etc. The paper presents an introduction of the treatment of multiscale problems using the FIC approach in conjunction with the finite element method (FEM). Examples of application of the FIC/FEM formulation to the solution of simple multiscale convection-diffusion problems are given.
1 INTRODUCTION
Multiscale problems in mechanics are characterized by the existence of numerical values of the geometrical and/or physical parameters differing in several orders of magnitude. Typical examples are convection-diffusion-absorption problems where high values of the convection or the absorption terms lead to sharp boundary and/or internal layers along which the numerical solution can change in many orders of magnitude. A similar problem is found in the analysis of thin discontinuity layers in solids, such as in the case of fracture in concrete, rock or geomaterials, or in shock waves in compressible fluids. A traditional multiscale problem is that of a turbulence fluid, where high gradients of the velocity field occur along random directions of the flow. The transition from a compressible to an incompressible solid or fluid can also be considered as a multiscale problem, in the sense that the propagation speed of sound changes from a finite value (for the compressible case) to an infinite value (in the case of incompressibility). Indeed, problems where many different scales of the material properties co-exist in a solid (such as in composites) are also classical examples of multiscale situations in solid mechanics. For the sake of conciseness, however, the case of composites will not be dealt with in the paper, although it can be also tackled with the general FIC formulation here described.
Recent attempts to develop a unified relativity theory in physics adequate for dealing with the whole spectrum of scales in the universe, from the Plank atomistic scale (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 10^{-33}}
cm) to the cosmological scale (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle 10^{+28}}
cm) indicate that a suitable mathematical model should incorporate the effect of the different scales within the governing equations of the model [1]. Failing to include these scale terms may lead to unstable behaviour of the equations when solved numerically for very different values of the physical, geometrical and/or time parameters of the problem.
The Finite Calculus (FIC) method developed by Oñate and co-workers [2,3,4,5,6,7] is a consistent procedure to re-formulate the governing equation in mechanics introducing new terms involving characteristic space and time dimensions into the equations. The modified equations are derived by invoking the balance laws in mechanics in a space-time domain of finite size. The new terms introduced by the FIC approach are essential to obtain physical (stable) numerical solutions for all ranges of the parameters governing the physical problem.
This paper presents an extension of the FIC method to deal with multiscale problems in mechanics using the finite element method (FEM) [10]. The discrete (nodal) values provided by the numerical FEM solution of the FIC equations are taken as the “correct macroscale” solution for a given discretization. A more refined numerical solution in zones where high gradients of the relevant variables occur can be simply obtained by using a hierarchical FE approximations or by introducing enriched interpolation functions accounting for the “correct” shape of the solution sought. The paper describes the basic ideas of the multiscale FIC procedure proposed and two examples of application to simple 1D and 2D convection-diffusion problems showing the applicability of the method. The possibilities of the FIC method for strain localization problems are also briefly highlighted.
2 THE FINITE CALCULUS METHOD IN MECHANICS
Let us consider a 1D space-time domain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar \Omega : \Omega \times t} , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega }
is the domain length and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle t}
is the time.
We can now invoke any of the classical balance laws in mechanics, to be satisfied within a space-time slab of finite size Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle [x-d,x]\times [t-T ,t] \subset \bar \Omega }
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d}
is the length of the space balance domain and and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle T}
is a time increment defining the size of the balance domain in the time axis. Using the standard procedure of expanding in Taylor series the changes of the chosen governing variables and retaining the linear terms in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle T}
gives the FIC balance equations as
|
(1) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r}
denotes the governing differential equations of the infinitessimal theory. In Eq.(1) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
are respectively characteristic length and characteristic time parameters satisfying
|
(2) |
For an advection-diffusion thermal problem, Eq.(1) is the heat balance equation with
|
(3) |
where the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
is the temperature, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \rho }
is the density, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v}
is the velocity, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle c}
are the thermal diffusivity and the specific heat parameters, respectively and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle Q}
is the heat source term.
Let us consider now a 1D solid mechanics problem. In this case Eq.(1) is the force equilibrium equation with
|
(4) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u}
is the displacement, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
is the axial force and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle b}
is the axial load.
The derivation of Eq.(1) for a simple steady state problem is given in the Appendix. Details of the derivation of the FIC equations for the transient case can be found in [6,7].
Note that as the dimensions of the balance slab tend to zero, also the values of the characteristic parameters Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
tend to zero. In the limit case Eq.(1) tends to
|
(5) |
which is the standard form of the governing differential equations of the infinitessimal theory.
Eq.(1) can be generalized for a multidimensional problem as
|
(6) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n_b}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n_d}
are the number of balance equations and the number of space dimensions, respectively.
A typical example of Eq.(6) are the FIC equations for an incompressible fluid written for the 3D case as
|
In above Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r_{m_i}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r_d}
are the ith momentum equation and the mass balance equation given by
|
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v_i}
are the velocities, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \sigma _{ij}=\tau _{ij}-p\delta _{ij}}
are the total stresses, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p}
is the pressure, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tau _{ij}}
are the viscous stresses and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle b_i}
are the body forces in the fluid [3,4].
Note that Eq.(8) is the result of applying the FIC concepts to the equations expressing the balance of mass over a space-time domain of finite size [2,4].
For a general solid mechanics problem the governing FIC equations stating the equilibrium of forces will be identical to Eq.(7) where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r_{m_i}}
is given by Eq.(9), dropping now the convective terms [8,9].
2.1 Boundary and initial conditions
The Dirichlet and initial conditions are the same as in the standard infinitessimal theory. The Neumann boundary condition expressing balance of fluxes across the boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma _q}
is however modified in the FIC formulation introducing, for consistency, the concept of balance over a domain of finite size adjacent to the Neumann boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma _q}
[2]. The new boundary condition is written generally as
|
(11) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_{n_i}}
is the normal “flow” across the boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma _q}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar q_{n_i}}
is the prescribed flow value and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n_i}
are the components of the unit normal to the boundary.
For example, in a 3D advection-diffusion problem, the boundary condition for the “diffusive” flux is
|
(12) |
For 3D fluid or solid mechanics problems, Eq.(11) would be
|
(13) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar t_i}
are the prescribed boundary tractions and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r_{m_i}}
is given by Eq.(9).
Applications of the FIC method to problems in fluid and solid mechanics using the finite element method and the meshless finite point method can be found in [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19].
3 WHAT DO THE CHARACTERISTIC PARAMETERS MEAN?
The characteristic parameters in the space and time dimension (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h_i}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
) can be interpreted as free intrinsic parameters giving the “exact” expression of the balance equations (up to first order terms) in a space-time domain of finite size.
Let us consider now the discretized solution of the modified governing equations. For simplicity we will focuss here in the simplest scalar 1D problem. The variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
is approximated as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi \simeq \hat \phi }
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat \phi }
denotes the approximated solution. The values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat \phi }
are now expressed in terms of a finite set of parameters using any discretization procedure (finite elements, finite diferences, etc.). The discretized FIC governing equation (1) would read now
|
(14) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat r:=r (\hat \phi )}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r_{\bar \Omega }}
is the residual of the discretized FIC equations.
The meaning of the characteristic parameters Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
in the discretized equation (14) changes (although the same simbols than in Eq.(1) have been kept). The Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
parameters become now of the order of magnitude of the discrete domain where the balance laws are satisfied. In practice, this means that
|
(15) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle l}
is the grid dimension in the space discretization (i.e. the element size in the FEM or the cell size in the FDM) and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Delta t}
is the time step used to solve the transient problem.
The values of the characteristic parameters can be found now in order to obtain an a correct numerical solution. The meaning of “correct solution” must be obviously properly defined. Ideally, this would be a solution giving “exact” values at a discrete number of points (i.e. the nodes in a FE mesh). This is infortunatelly impossible for practical problems (with the exception of very simple 1D and 2D cases) and, in practice, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
are computed making use of some chosen optimality rule, such as ensuring that the error of the numerical solution diminishes for appropriate values of the characteristic parameters [2,12,13]. This suffices in practice to obtain physically sound (stable) numerical results for any range of the physical parameters of the problem, always within the limitation of the discretization method chosen.
The concept of discretization of the FIC governing equations is explained in more detail with an example in the next section.
4 FINITE ELEMENT DISCRETIZATION OF THE FIC EQUATIONS. AN EXAMPLE: THE ADVECTIVE-DIFFUSIVE PROBLEM
Let us consider the FIC governing equations for the steady-state advective-diffusive problem defined in vector form for the multidimensional case as
|
(16) |
with
|
(17a) |
|
(17b) |
with
|
(18) |
In above equations Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h}
is the characteristic length vector, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol \nabla }}
is the gradient vector, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{D}}
is the diffusivity matrix, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{n}}
is the normal vector and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{v}}
is the velocity vector. As usual Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma _\phi }
is the Dirichlet boundary where the unknown is prescribed to a fixed value Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \bar \phi }
.
A finite element interpolation of the unknown Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
can be written as
|
(19) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N_i}
are the shape functions and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat \phi _i}
are the nodal values of the approximate function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat \phi }
[10].
Application of the Galerkin FE method to Equations (12)-(13b) gives, after integrating by parts the term Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol \nabla }r}
(and neglecting the space derivatives of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}}
) [2,10]
|
(20) |
The last integral in Equation (20) has been expressed as the sum of the element contributions to allow for interelement discontinuities in the term Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol \nabla }\hat r} , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat r = r (\hat \phi )}
is the residual of the FE solution of the infinitesimal equations.
Note that the residual terms have disappeared from the Neumann boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma _q} . This is due to the consistency of the FIC terms in Equation (17b).
Integrating by parts the diffusive terms in the first integral of (20) leads to
|
(21) |
In matrix form
|
(22) |
Matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{K}}
and vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{f}}
are assembled from the element contributions given by
|
(23) |
|
(24) |
Note that the method introduces in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol K}
an additional diffusivity matrix given by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {1\over 2} \mathbf{h}^T \mathbf{v}}
. Also the second integral of Equation (23) vanishes for linear approximations. The same happens with the second term of the first integral of Equation (24) when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N_i}
is linear and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle Q}
is constant. The evaluation of these integrals is mandatory in any other case.
4.1 The discrete 1D advective-diffusive problem
The role of the stabilization parameter can be better clarified with a simple 1D steady-state advection-diffusion problem governed by the FIC equation
|
(25) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
is the length of the analysis domain and
|
(26) |
The form of the element stiffness matrix and the element flux vector for the 1D case are obtained by simplification of Eqs.(23) and (24) as
|
(27a) |
|
(27b) |
For simplicity we will assume Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle Q=0}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi=0}
at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=0}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi = \bar \phi }
at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=L}
.
We will solve the problem with the simplest mesh of two linear finite element of dimensions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle l^e=L/2} . The stencil for this problem is
|
(28) |
The solution for the central node is obtained making Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _1=0}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _3=\bar \phi }
in Eq.(28) as
|
(29) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \gamma =\displaystyle{vl^e\over 2k}}
is the element Peclet number and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha =\displaystyle{h\over l^e}}
is a dimensionless characteristic parameter.
We note that for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha =0}
(i.e. the discrete solution of the standard infinitessimal equations) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _2 = {1\over 2} (1-\gamma ) \bar \phi }
which gives non physical (unstable) negative values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _2}
for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \gamma > 1}
. A critical value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha }
can be deduced from Eq.(29) ensuring that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _2 \ge 0}
for all range of the parameters Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u,k}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle l}
. This gives
|
(30) |
which is the well known expression of the critical stabilization parameter for this problem [10,11].
The value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha }
given by Eq.(30) allow us to obtain the best possible physically meaningful solution for the problem. Indeed for this simple case we could aim to obtain the exact analytical solution at node Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _2}
. This is simply granted is Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha }
is taken as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha =\coth \gamma - 1/\gamma }
which can be again deduced from Eq.(28) [6,10]. Note that the values of this expression for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha }
and those of Eq.(30) are quite similar for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \gamma > 2}
[10].
The merit of the expression of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha }
of Eq.(30) is that it has been deduced with no previous knowledge of the exact solution. Indeed, the approximated solution within the element adjacent to the boundary at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=L}
will differ substantially from the exact one: the numerical solution is linear whereas the analytical solution gives a boundary layer adjacent to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=L}
. The stabilized FIC numerical solution can however be taken now as the starting point to obtain a refined solution giving a more accurate description of the physical problem within the element where the boundary layer is located. This can be obtained by mesh refinement or either by enriching the approximation within that element in an adaptive manner. Here a number of possibilities can be explored using mesh refinement, hierarchical approximations, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle p} -refinement or enrichment techniques, etc. Some of these options are discussed in a next section.
4.2 The characteristic length vector
From the FIC expressions in Section 4 we note that in multidimensional problems the space stabilization parameters can be grouped in a vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h} , where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h}
has the space dimensions of the problem (i.e. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h} =[h_1,h_2,h_3]^T}
for 3D situations). Computation of vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h}
is a crucial step as the quality of the stabilized solution depends on the choice of the module and direction of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h}
.
We could, for instance, assume that the direction of vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h}
is parallel to the velocity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol v}
, i.e. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}= h{\mathbf{v}\over \vert \mathbf{v}\vert }}
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
is a characteristic length. Under these conditions, Equation (20) reads
|
(31) |
Equation (31) coincides precisely with the so called Streamline-Upwind-Petrov-Galerkin (SUPG) method [2,10,20]. The ratio Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \displaystyle{h\over 2\vert \mathbf{v}\vert }}
has dimensions of time and it is termed element intrinsic time parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \tau }
.
It is important to note that the SUPG expression is a particular case of the more general FIC formulation. This explains the limitations of the SUPG method to provide stabilized numerical results in the vicinity of sharp gradients of the solution transverse to the flow direction. In general, the adequate direction of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{ h}}
is not coincident with that of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{v}}
and the components of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}}
introduce the necessary stabilization along the streamlines and the transverse directions to the flow. In this manner, the FIC method reproduces the best-features of the so-called stabilized discontinuity-capturing schemes [5,12,13].
The computation of the characteristic vector Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h}
is a topioc which fails outside the scope of this paper. We will just mention that in practice it is convenient to split Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol h}
as
|
(32) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_s}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_t}
are called streamline and transverse characteristic lenght vectors. A useful option is to take
|
(33) |
The characteristic lengths Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {h}_s}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {h}_t}
are computed at element level by imposing that the residual of the numerical solution diminishes over each element [2,5,12,13]. In practice, these lengths can be computed for each element as the maximum value of the scalar product between the element sides and the velocity vector (for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h_s}
) and the gradient vector (for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h_t} ) [14,15,16,17].
Note that the gradient of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
in the expression of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_t}
introduces a non linearity in the computation process which is typical of discontinuity capturing techniques. A procedure to linealize the expression of the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}_t}
in element next to boundary layers by approximating Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\boldsymbol \nabla }\phi }
by the normal vector to the boundary is presented in [5].
5 INTERPRETATION OF THE DISCRETE SOLUTION OF THE FIC EQUATIONS
Let us consider the solution of a physical problem in a space domain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega } , governed by a differential equation Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r(\phi ) =0}
in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega }
with the corresponding boundary conditions. The “exact” (analytical) solution of the problem will be a function giving the sought distribution of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
for any value of the geometrical and physical parameters of the problem. Obviously, since the analytical solution is difficult to find (practically impossible for real situations), an approximate numerical solution is found Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi \simeq \hat \phi }
by solving the problem Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat r =0}
, with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat r = r(\hat \phi )} , using a particular discretization method (such as the FEM). The distribution of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega }
is now obtained for specific values of the geometrical and physical parameters. The accuracy of the numerical solution depends on the discretization parameters, such as the number of elements and the approximating functions chosen in the FEM. Figure 1 shows a schematic representation of the distribution of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat \phi }
along a line for different discretizations Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_1,M_2,\cdots , M_n}
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_n}
are the coarser and finer meshes, respectively. Obviously for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle n}
sufficiently large a good approximation of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
will be obtained and for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_\infty }
the numerical solution Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat \phi }
will coincide with the “exact” (and probably unreachable) analytical solution Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
at all points. Indeed in some problems the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_\infty }
solution can be found by a “clever” choice of the discretization parameters.
An unstable solution will occur when for some (typically coarse) discretizations, the numerical solution provides non physical or very unaccurate values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat \phi } . A situation of this kind is represented by curves Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_2}
of the left hand side of Figure 1. These unstabilities would disappear by an appropriate mesh refinement (curves Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_3,M_4 \cdots }
in Figure 1) at the obvious increase of the computational cost.
In the FIC formulation the starting point are the modified differential equations of the problem and the corresponding modified boundary conditions as previously described. The modified equations are not useful to find an analytical solution, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi (x)} , for the physical problem. However, the numerical solution of the FIC equation can be readily found. Moreover, by adequately choosing the values of the characteristic length parameter Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h} , the numerical solution of the FIC equations will be always stable (physically sound) for any discretization level chosen.
This process is schematically represented in Figure 1 where it is shown that the numerical oscillations for the coarser discretizations Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_2}
dissapear when using the FIC procedure.
We can therefore conclude the FIC approach allows us to obtain a better numerical solution for a given discretization. Indeed, as in the standard infinitesimal case, the choice of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle M_\infty }
will yield the (unreachable) exact analytical solution and this ensures the consistency of the method.
|
| Figure 1: Schematic representation of the numerical solution of a physical problem using standard infinitesimal calculus and finite calculus. |
6 MULTISCALE ADAPTIVE FINITE ELEMENT ANALYSIS USING THE FIC FORMULATION
6.1 Adaptive solution by updating the stability parameters
The FIC formulation can be applied in conjunction with any numerical method (say the finite element method) in order to obtain a “correct” solution within the limit of the discretization chosen. The accuracy of the finite element solution depends greatly of the values chosen for the stabilization parameters in the FIC formulation. Indeed if these values are taken equal to zero, an unstable non physical solution is obtained in many cases. The space and time stabilization parameters Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h_i}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
in the FIC formulation can be given adequate values such that the desired accuracy in the numerical solution is achieved. In [2,5,12,13] a procedure is presented to adaptively compute the values of the stabilization parameters so that the element error, measured in terms of the average residual of the governing equations over the element, diminishes until it reaches a prescribed value. Other procedures to compute the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h_i}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta }
parameters can be devised, all aiming to obtain the “best” possible nodal solution for a given mesh.
Once this solution has been obtained, the problem of accurately reproducing sharp gradients in selected zones remains. Some methods aiming to solve this problem are proposed next.
The nodal finite element values computed with the stabilized FIC formulation can be taken as the starting point for obtaining an enhanced solution in zones where high gradients of the sought variables occur. Three options will be considered here: mesh adaptivity, hierarchical shape functions and enriched interpolations incorporating the shape of the analytical solution.
6.2 Mesh refinement
This is the more standard procedure, although it is still unpractical for many practical 3D problems involving complex geometries. The mesh adaption process can be based on mesh regeneration or mesh enrichment procedures. The adaptivity strategy can use global or local error estimation based on adequate error norms. Details of mesh adaptivity procedures based on different error measures can be found elsewhere and will not be discussed here [10,21,22]. Given the intrinsic difficulties for modifying the meshes for 3D problems we will favour the two procedures described next versus the mesh adaptivity process.
6.3 Hierarchical interpolation
The nodal values provided by the FIC method can be considered “correct” (in the sense that they are physically sound) whithin the limitations of the mesh discretization chosen. A procedure to obtain an enhanced approximation within a particular element, or a patch of elements, is to introduce hierarchical interpolations in the selected elements. The hierarchical interpolation can be simply written for an individual element
|
(34) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N_i}
denote the standard shape functions, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _i}
are the nodal values, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N_j^h}
are the hierarchical shape functions and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle a_j^h}
are the hierarchical values. As usual, the functions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N_j^h}
should vanish along the element edges. The simplest hierarchical interpolations (satisfying orthogonality conditions) can be based on Legendre polynomials or on standard trigonometric functions [10]. More advanced hierarchical finite element interpolations based on the concept of the partition of unity can be effectively used [23].
A good approximation of the hierarchical variables can be obtained taking the nodal values obtained with the FIC numerical solution as fixed. A simple set of equations needs to be solved at each element level only for the hierarchical variables Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle a_j^h} . Indeed, this computation can be performed within a more general iterative procedure, where the nodal variables and the hierarchical values are updated in an staggered manner.
6.4 Enrichment interpolation using an approximation to the analytical shape of the solution
An effective procedure to obtain an enhanced solution starting with the FIC nodal values is to use a local interpolation over each element injecting into the interpolating functions a good approximation to the shape of the sought solution. This can be effectively carried out by using as the new approximation the free-space solution of the homogeneous and constant-coefficient version of the standard governing partial differential equations. A procedure of this kind has been used by Farhat et al. [24] who added the enriched field to the initial polynomial field. The resulting field is not continuous in this case and continuity is enforced by means of Lagrange multipliers.
In this work it is proposed to apply the “analytical” enrichment in order to obtain an enhanced solution within the selected elements where high gradients occur, while keeping the nodal values given by the FIC solution fixed.
Indeed the three adaptive procedures above suggested can be combined and the efficiency of the different alternatives should be investigated and tested in practice.
The concepts explained are clarified next with two simple examples of application.
7 EXAMPLES
7.1 1D advective-diffusive problem
We will consider the simplest 1D advective-diffusive problem governed by Eqs.(25) and (26) with boundary conditions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi =0}
at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=0}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi =\bar \phi }
at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=L}
. Again for simplicity we will consider a mesh of two linear 1D elements.
Figure 2 shows the unstable numerical solution obtained for the case of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha =0}
(Galerkin) and the “correct” FIC solution for a critical value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha =(1- {1\over \gamma }) (1+\varepsilon )}
with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \varepsilon =10^{-10}}
. Note that the numerical solution, although nodally correct, is far from the “exact” analytical distribution giving a boundary layer next to the end node at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=L} .
|
| Figure 2: 1D advection-diffusion problem solved with two linear elements (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \gamma =50
). Unstable solution (- - -) obtained for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \alpha =0 (Galerkin). Stabilized solution (–-) obtained with the FIC method for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \alpha = (1+ 10^{-10})(1- {1\over \gamma })
|
We will consider next a hierarchical approximation of the solution over the element 2–3 using Eq.(34) with hierarchical interpolation functions defined as
|
(35) |
Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N_j^h}
vanishes for the element ends located at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \xi =\pm 1}
. Figure 3 shows the approximation of the numerical solution over the element 2–3 obtained for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle m=3}
keeping the nodal values given by the original FIC solution fixed. The improvement in the solution is noticeable. An almost identical solution is obtained if the value for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _2}
is recomputed until convergence is achieved. A similar improvement is obtained using Legendre polynomials up to a quartic order.
A considerable better solution over the element 2–3 can be obtained by using for approximating the unknown within this element the following interpolation
|
(36) |
with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle P_1 = \displaystyle{e^{\gamma \xi }-e^{\gamma }\over e^{-\gamma } - e^{\gamma }}\quad ,\quad P_2=\displaystyle{e^{-\gamma } - e^{\gamma \xi }\over e^{-\gamma } - e^{\gamma }}} .
The new shape functions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle P_1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle P_2}
are obtained from the analytical solution of the homogeneous differential equation (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle r=0}
) for constant values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle u}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
. Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle P_i (\xi _j ) =1}
for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle i=j}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle =0}
for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle i\not = j}
.
Application of the approximation (36) over the element 2–3 with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _2^{(e)} =\bar \phi }
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _1^{(e)} = \phi _2}
, where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi _2}
is the nodal value at node 2 obtained with the FIC formulation, reproduces accurately the boundary layer within this element, as expected (see Figure 3). Recall that Eq.(36) has been applied “a posteriori”, once the nodal value at node 2 has been found using the original linear interpolation and the FIC method with the critical value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha }
above given.
|
| Figure 3: 1D advection-diffusion problem solved with two linear elements. Initial stabilized FIC solution (curve A). Enhanced solutions over element 2–3 using a hierarchical trigonometric interpolation (curve B) and an exponential interpolation (curve C) |
This simple example shows that the nodal values provided by the FIC method are a good starting point for obtaining an enhanced numerical solution in the vecinity of sharp gradients. The results from this simple example indicate that it suffices with a straightforward “a posteriori” analysis to obtain a better solution in these regions using either a hierarchical interpolation or (ideally) an approximation enriched with the shape of the solution of the homogeneous equation for constant coefficients.
7.2 2D stationary convection-diffusion problem with diagonal velocity, zero source and uniform Dirichlet conditions
|
|
|
|
Figure 4: Convection-diffusion problem with uniform Dirichlet conditions. One step solution. Distribution of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \phi
along the lines Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): AA'
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): BB'
| |
The steady-state convection-diffusion equation is solved in a domain of unit size (Figure 4) with
|
(37) |
The following Dirichlet conditions are assumed
|
|
The expected solution is a uniform distribution of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi =0}
over the domain except in the vecinity of the boundaries Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y=1}
where a boundary layer is formed.
The domain is discretized with a uniform mesh of 400 four node quadrilaterals (Figure 4).
The solution has been obtained in one single step using the value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathbf{h}}
given by Eqs.(32) and (33) approximating the gradient of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
in the elements next to the outflow boundaries by the normal vector. For details see [5]. Figure 4 shows that the sharp gradients expected are perfectly captured over the boundary elements without any oscillation.
An enhanced solution over the elements adjacent to the outflow boundaries can now be obtained using hierarchical interpolations or exponential interpolations similarly as described for the 1D example in the previous section. The efficiency of these procedures is similar to that shown in Figure 4 for the 1D case giving in both cases an enhanced reproduction of the boundary layer within the elements adjacent to the lines Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle x=1}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y=1}
.
8 POSSIBILITIES OF FINITE CALCULUS FOR STRAIN LOCALIZATION
The FIC method introduces naturally higher order derivative terms of the displacements in the equilibrium equations. These terms resemble those introduced by the Cosserat model [25] and the so called “non local” constitutive models [25,26,27,28,29]. These models are typically used in order to preserve the ellipticity of the solid mechanics equations in the presence of localized high displacement gradient zones, such as shear bands and fracture lines.
For instance, the FIC governing equations for a simple 1D bar can be written (in absence of external body forces) as [5]
|
(38) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
is the axial force.
The relationship between the axial force and the elongation Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \varepsilon =\displaystyle{du\over dx}}
in the softening branch can be written in incremental form as
|
(39) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle E_T}
is the tangent modulus of the material and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \delta u}
is the displacement increment.
Substituting Equation (39) into (38) gives
|
(40) |
Equation (40) resembles the governing equations derived by Lasry and Belytschko [26] using localization limiters in the stress-strain relationship and by Schreier and Chen [27] using gradient-dependent plasticity constitutive models.
Many similarities can be found between the governing equations derived by the FIC method and those resulting from the non local constitutive models [30]. Although in some situations the resulting governing differential equations are the same, the basic difference between the two approaches is that non local constitutive models aim to enhance the constitutive equation by adding new terms (typically strain gradient terms), whereas the FIC method defines a new equilibrium equation applicable to all the discrete scales while preserving the features of the constitutive equation. This opens a world of possibilities for the FIC equations in solid mechanics to solve strain localization problems using standard constitutive models.
9 CONCLUSIONS
The paper has presented some basic ideas to combine the finite calculus (FIC) method with adaptive finite element strategies in order to solve problems where sharp gradients of the solution occur in localized zones. The key argument of the procedures proposed is that the discretized FE solution of the FIC equations leads to nodal values which are “as good as possible” by virtue of the adequate selection of the stabilization parameters. The FIC solution can be used as the starting point to compute an “a posteriori” enhanced solution over selected elements using hierarchical or enriched approximations. Indeed, mesh adativity can be included here as a third alternative and the combination of the different procedures opens new areas of research.
The same ideas can be used to solve a wide class of multiscale problems in mechanics where the solution changes in several orders of magnitude in specific regions.
REFERENCES
[1] Nottale, L. Relativity in all its states. (in French), Hachette, 1998.
[2] Oñate E. Derivation of stabilized equations for advective-diffusive transport and fluid flow problems. Comput. Methods Appl. Mech. Engrg. 1998; 151:1-2:233–267.
[3] Oñate E, García J. A finite element method for fluid-structure interaction with surface waves using a finite calculus formulation. in Comput. Methods Appl. Mech. Engrg. 2001; 191:6-7:635-660.
[4] Oñate E. A stabilized finite element method for incompressible viscous flows using a finite increment calculus formulation. Comput. Methods Appl. Mech. Engrg. 2000; 182:1–2:355–370.
[5] Oñate E. Possibilities of finite calculus in computational mechanics. To be published in Int. J. Num. Meth. Engng., 2003.
[6] Oñate E, Manzan M. Stabilization techniques for finite element analysis of convection diffusion problems. In Comput. Anal. of Heat Transfer, WIT Press, G. Comini and B. Sunden (Eds.), 2000.
[7] Oñate E, Manzan M. A general procedure for deriving stabilized space-time finite element methods for advective-diffusive problems. Int. J. Num. Meth. Fluids 1999; 31:203–221.
[8] Oñate E, Rojek J, Taylor RL, Zienkiewicz OC. Finite calculus formulation for analysis of incompressible solids using linear triangles and tetrahedra. Submitted to Int. J. Num. Meth. Engng., May 2002.
[9] Oñate E, Rojek J, Taylor RL, Zienkiewicz OC. Linear triangles and tetrahedra for incompressible problems using a finite calculus formulation. In European Conference on Computational Mechanics (ECCM2001), Cracow, Poland, 26–29 June, 2001.
[10] Zienkiewicz OC, Taylor RL. The finite element method. 5th Edition, 3 Volumes, Butterworth–Heinemann, 2000.
[11] Hirsch C. Numerical computation of internal and external flow, J. Wiley, Vol. 1 1988, Vol. 2, 1990.
[12] Oñate E, García J, Idelsohn SR. Computation of the stabilization parameter for the finite element solution of advective-diffusive problems. Int. J. Num. Meth. Fluids 1997;25:1385–1407.
[13] Oñate E, García J, Idelsohn SR. An alpha-adaptive approach for stabilized finite element solution of advective-diffusive problems with sharp gradients. New Advances in Adaptive Comput. Methods in Mech., Elsevier, P. Ladeveze, J.T. Oden (Eds.), 1998.
[14] García J, Oñate E. An unstructured finite element solver for ship hydrodynamics. J. Appl. Mechanics 2002.
[15] Oñate E, García J, Bugeda G, Idelsohn SR. A general stabilized formulation for incompressible fluid flow using finite calculus and the FEM. In Towards a New Fluid Dynamics with its Challenges in Aeronautics, J. Periaux, D. Champion, O. Pironneau and Ph. Thomas (Eds.), CIMNE, Barcelona, 2002.
[16] Oñate E, García J, Idelsohn SR. Ship Hydrodynamics. Encyclopedia of Computational Mechanics, T. Hughes, R. de Borst and E. Stein (Eds.), J. Wiley, 2004 (to be published).
[17] Oñate E, Taylor RL, Zienkiewicz OC, Rojek J. A residual correction method based on finite calculus. FEMClass of 42 Meeting, 30–31 May, Ibiza, 2002, www.congress.cimne.upc.es/femclass42. To be published in Engineering Computation.
[18] Oñate E, Idelsohn SR, Zienkiewicz OC, Taylor RL. A finite point method in computational mechanics. Applications to convective transport and fluid flow. Int. J. Num. Meth. Engng. 1996;39:3839–3866.
[19] Oñate E, Sacco C, Idelsohn S. A finite point method for incompressible flow problems. Computing and Visualization in Science 2000; 2:67–75.
[20] Brooks A, Hughes TJR. Streamline upwind/Petrov-Galerkin formulation for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations. Comput. Methods Appl. Mech. Engrg. 1982; 32:199–259.
[21] Ladeveze P., Pelle J. La maïtrise du calcul en mecanique linearise et non lineaire. Hermes, 2001.
[22] Oñate E, Bugeda G. A study of mesh adaptivity criteria in adaptive finite element computations. Engng. Computations 1993; 10:307–321.
[23] Taylor RL, Zienkiewicz OC, Oñate E. A hierarchical finite element method based on the partition of unity. Comput. Methods Appl. Mech. Engrg. 1998; 152:73–84
[24] Farhat C, Harari I, Franca L. The discontinuous enrichment method. Comput. Methods Appl. Mech. Engrg. 2001; 190:6455–6479.
[25] de Borst R. Simulation of strain localization: A reappraisal of the Cosserat continuum. Continuum Computation 1991; 8:317–332.
[26] Larsry D, Belytschko T. Localization limiters in transient problems. Int. J. Solids and Structures 1988; 24:581–597.
[27] Schreier HL, Chen Z. One dimensional softening with localization. ASME J. Appl. Mech. 1980; 53:791–797.
[28] de Borst R. Gradient-dependent plasticity: Formulation and algorithmic aspects. Int. J. Num. Meth. Engng 1992; 35:521–539.
[29] Peric D, Yu J, Owen DRJ. On error estimates and adaptivity in elastoplastic solids: Application to the numerical simulation of strain localization in classical and Cosserat continua. Int. J. Num. Meth. Engng. 1994; 37:1351–1379.
[30] Borja R. Possibilities of the FIC method for strain localization problems. Private communication, 2001.
APPENDIX
We will consider a convection-diffusion problem in a 1D domain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega }
of length Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle L}
. The equation of balance of fluxes in a subdomain of size Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d}
belonging to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega }
(Figure A1) is written as
|
(41) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_A}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_B}
are the incoming and outgoing fluxes at points Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle B}
, respectively. The flux Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q}
includes both convective and diffusive terms; i.e. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q=v\phi - k{d\phi \over dx}}
, where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
is the transported variable (i.e. the temperature in a thermal problem), Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle v}
is the velocity and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
is the diffusitivity of the material.
|
| Figure 5: Equilibrium of fluxes in a space balance domain of finite size |
Let us express now the fluxes Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_A}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_B}
in terms of the flux at an arbitrary point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C}
within the balance domain (Figure A.1). Expanding Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_A}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_B}
in Taylor series around point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C}
up to second order terms gives
|
(42) |
Substituting Equation (A.2) into Equation (A.1) gives after simplification
|
(43) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h=d_1-d_2}
and all the derivatives are computed at point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C}
.
Standard calculus theory assumes that the domain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d}
is of infinitesimal size and the resulting balance equation is simply Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {dq\over dx}=0}
. We will relax this assumption and allow the space balance domain to have a finite size. The new balance equation (A.3) incorporates now the underlined term which introduces the characteristic length Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h} . Obviously, accounting for higher order terms in Equation (A2) would lead to new terms in Equation (A.3) involving higher powers of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h} .
Distance Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
in Equation (A.3) can be interpreted as a free parameter depending on the location of point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C}
within the balance domain. Note that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle -d\le h\le d}
and, hence, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
can take a negative value. At the discrete solution level the domain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d}
should be replaced by the balance domain around a node. This gives for an equal size discretization Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle -l^e\le h \le l^e}
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle l^e}
is the element or cell dimension. The fact that Equation (A.3) is the exact balance equation (up to second order terms) for any 1D domain of finite size and that the position of point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle C}
is arbitrary, can be used to derive numerical schemes with enhanced properties simply by computing the characteristic length parameter from an adequate “optimality” rule leading to an smaller error in the numerical solution.
Consider, for instance, Equation (A.3) applied to the 1D convection-diffusion problem. Neglecting third order derivatives of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi } , Equation (A.3) can be rewritten in terms of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
as
|
(44) |
We see clearly that the FIC method introduces naturally an additional diffusion term in the standard convection-diffusion equation. This is the basis of the popular “artificial diffusion” procedure [10,11] where the characteristic length Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle h}
is typically expressed as a function of the cell or element dimension.
Equation (A.3) can be extended to account for source terms. The modified governing equation can then be written in compact form as
|
(45) |
with
|
(46) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle Q}
is the external source. For consistency a “finite” form of the Neumann boundary condition should be used. This can be readily obtained by invoking balance of fluxes in a domain of finite size next to the boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma _q}
where the external flux is prescribed to a value Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle q_p}
. The modified Neumann boundary condition is [2]
|
(47) |
The governing equations of the problem are completed with the standard Dirichlet condition prescribing the value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \phi }
at the boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma _\phi }
.
The underlined terms in Equations (A.5) and (A.7) introduce the necessary stabilization in the discrete solution using whatever numerical scheme.
The time dimension can be simply accounted for the FIC method by considering the balance equation in a space-time slab domain as described in [2,6,7].
Document information
Published on 01/01/2003
DOI: 10.1016/S0045-7825(03)00340-2
Licence: CC BY-NC-SA license
