Spurious transients of projection methods in microflow simulations

Autor(es)

Sousa, Fabricio S. ^[1] ; Oishi, Cassio M. ^[2] ; Buscaglia, Gustavo C. ^[1]

Afiliação do(s) autor(es)

^[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

• Tipo de documento: Artigo Científico
• Fonte: COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING; v. 285, p. 659-693, MAR 1 2015.
• Citações Web of Science: 2

Resumo

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)

• Processo FAPESP: 11/00538-5 - Microfluídica computacional: uma abordagem lagrangeana-euleriana arbitrária
• Beneficiário: Fabrício Simeoni de Sousa
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 12/14481-8 - Aproximação numérica de interfaces microfluídicas com acoplamento mecânico-bioquímico
• Beneficiário: Gustavo Carlos Buscaglia
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 09/15892-9 - Estudo de métodos numéricos estáveis e acurados em escoamentos transientes: aperfeiçoamentos, implementações, aplicações com superfícies livres e modelos viscoelásticos
• Beneficiário: Cassio Machiaveli Oishi
• Modalidade de apoio: Auxílio à Pesquisa - Jovens Pesquisadores