| Full text | |
| Author(s): |
Total Authors: 3
|
| Affiliation: | [1] Univ Sao Paulo, ICMC, Dept Matemat Aplicada & Estat, BR-13560970 Sao Carlos, SP - Brazil
[2] Univ Estadual Paulista, Dept Matemat & Comp, Fac Ciencias & Tecnol, BR-19060900 Presidente Prudente, SP - Brazil
Total Affiliations: 2
|
| Document type: | Journal article |
| Source: | COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING; v. 285, p. 659-693, MAR 1 2015. |
| Web of Science Citations: | 2 |
| Abstract | |
The temporal behavior of projection methods for viscous incompressible low-Reynolds-number flows is addressed. The methods considered result from algebraically splitting the linear system corresponding to each time step, in such a way that the computation of velocity is segregated from that of pressure. Each method is characterized by two (possibly equal) approximate inverses (B-1 and B-2) of the momentum-equation velocity matrix, plus a parameter gamma which renders the method non-incremental (if gamma = 0) or incremental (if gamma = 1). The classical first-order projection method, together with more sophisticated methods (Perot's second-order method, Yosida method, pseudo-exact factorization method) and their incremental variants are put into the same algebraic form and their accuracy numerically tested. Splitting errors of first, second and third order in the time step size delta t are obtained, depending on the method. The methods are then discussed in terms of their ability and efficiency to compute steady states. Non-incremental methods are impractical because extremely small time steps are required for the steady state, which depends on delta t, to be reasonably accurate. Incremental methods, on the other hand, either become unstable as delta t is increased or develop a remarkable spurious transient which may last an extremely long time (much longer than any physical time scale involved). These transients have serious practical consequences on the simulation of steady (or slowly varying), low-inertia flows. From the physical viewpoint, the spurious transients may interfere with true slow processes of the system, such as heat transfer or species transport, without showing any obvious symptoms (wiggly behavior in space or time, for example, do not occur). From the computational viewpoint, the limitation in time step imposed by the spurious transient phenomenon weighs against choosing projection schemes for microflow applications, despite the low cost of each time step. (c) 2014 Elsevier B.V. All rights reserved. (AU) | |
| FAPESP's process: | 11/00538-5 - Computational microfluidics: an arbitrary lagrangian-eulerian approach |
| Grantee: | Fabrício Simeoni de Sousa |
| Support Opportunities: | Regular Research Grants |
| FAPESP's process: | 12/14481-8 - Numerical approximation of microfluidic interfaces with mechano-bio-chemical coupling |
| Grantee: | Gustavo Carlos Buscaglia |
| Support Opportunities: | Regular Research Grants |
| FAPESP's process: | 09/15892-9 - Study of stable and accurate numerical methods for transient flows: improvements, implementations, free surface flow problems and viscoelastic models |
| Grantee: | Cassio Machiaveli Oishi |
| Support Opportunities: | Research Grants - Young Investigators Grants |