Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The full numerical simulation of turbulent flows in the context of industrial applications remains a very challenging problem. One of the difficulties is the presence of a large number of interacting scales of different orders of magnitude, ranging from the macroscopic scale to the Kolmogorov dissipation scale. In order to better understand and simulate these interactions, new algorithms of incremental type have recently been introduced, stemming from the theory of infinite-dimensional dynamical systems (see, e.g., the algorithms of Foias, Jolly et al. (1988), Foias, Manley and Temam (1988), Laminie et al. (1993, 1994), Marion and Temam (1989, 1990), and Temam (1990)). These algorithms are based on decompositions of the unknown functions into a large- and a small-scale component, one of the underlying ideas being to approximate the corresponding attractor. In the context of finite element methods, the decomposition of a solution into small- and large-scale components appears when we consider hierarchical bases (Yserentant, 1986; Zienkiewicz et al., 1982). In the present article we derive several estimates of the size of the structures for the linear and nonlinear terms which correspond to interactions of different hierarchical components of the velocity field, and also their time variation. The one-step time variation of these terms can be smaller than the expected accuracy of the computation. Using this remark, we implement an adaptive spatial and temporal multilevel algorithm which treats differently the small- and large-scale components of the flow. We derive several a priori estimates in order to study the perturbation introduced into the approximated equations. All the interaction terms between small- and large-scale structures are frozen during several time steps. Finally, we implement the multilevel method in order to simulate a bidimensional flow described by the Burgers equations. We perform a parametric study of the procedure and its efficiency. The gain in CPU time is significant and the trajectories computed by our multi-scale method remain close to the trajectories obtained with the classical Galerkin method.
ISSN: 0168-9274

hal-00730710
https://hal.archives-ouvertes.fr/hal-00730710

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Algorithms issued from the NonLinear Galerkin method have been used in many situations and with different discretizations for the resolution of evolutionary nonlinear equations. The main idea of these methods is to use a splitting of the solution in order to model the equation. According to the splitting of the solution, a splitting of the equation is obtained. The modeling principle is to freeze terms which have a small time variation. In this work we use wavelet discretizations of the 2-D Burgers equations and compare the results with the hierarchical finite elements method. The numerical tests indicate that wavelets give better results than finite elements
ISSN: 0885-7474

hal-00383423
https://hal.archives-ouvertes.fr/hal-00383423
DOI : 10.1007/s10915-004-4806-4

Éditeur(s) : HAL CCSD
Résumé : International audience
In this paper, via the algebra of generalized functions, we investigate the generalized Riemann's problem associated to conservation laws with analytical coefficients. This allows us to transform the problem into a system of ordinary differential equations. In some particular cases, such that Burgers' and conservative Richard's equation, approximated solutions are obtained by the truncation of the so called Hugoniot-Maslov's chain and numerical simulations are also presented in the case of equations with polynomial coefficients.
Novi Sad Journal of Mathematics

hal-00771405
https://hal.archives-ouvertes.fr/hal-00771405

Éditeur(s) : HAL CCSD Society for Industrial and Applied Mathematics
Résumé : International audience
In this paper, we developed a theoretical study for nonconservative sytems in one dimension in order to construct numerical schemes for solving the Riemann problem. The nonconservative form of our model system required the use of a well-adapted theory in order to give us a sense of our problem. We chose a framework of generalized functions for solving a scalar hyperbolic equation with a discontinuous coefficient $sigma_t +usigma_x approx 0$, where u is the velocity solution of a Burgers's equation. After an explicit solution of the Riemann problem, we derived Godunov split schemes for computing an approximate solution of the Cauchy problem. We applied our study to a system modeling elasticity and a system modeling gas dynamics. Some stability properties of a scheme and its convergence to a generalized solution are proved for the first model. Numerical experiments confirmed this convergence result. For the second model, calculations of flows containing weak-to-moderate shocks showed that conservation errors are reduced when the mesh is refined but were not entirely eliminated.
ISSN: 0036-1429

hal-00770266
https://hal.univ-antilles.fr/hal-00770266
DOI : 10.1137/S0036142900378637

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
This article is intended as a preliminary report on the implementation of a finite volume multilevel scheme for the discretization of the incompressible Navier-Stokes equations. As is well known, the use of staggered grids (e.g., MAC grids [M. Perić, R. Kessler and G. Scheuerer, Comput. & Fluids 16 (1988), no. 4, 389-403; Zbl 0672.76018]) is a serious impediment for the implementation of multilevel schemes in the context of finite differences. This difficulty is circumvented here by the use of a colocated finite volume discretization [S. Faure, J Laminie and R. Temam, "Colocated finite volume schemes for fluid flows'', submitted; M. Perić, R. Kessler and G. Scheuerer, op. cit.], for which the algebra of multilevel methods is much simpler than in the context of MAC type finite differences. The general ideas and the numerical simulations are presented in this article in the simplified context of a two-dimensional Burgers equation; the two- and three-dimensional Navier-Stokes equations, which introduce new difficulties related to the incompressibility condition and the time discretization, will be considered elsewhere [see S. Faure, J Laminie and R. Temam, op. cit.; "Finite volume discretization and multilevel methods for the Navier-Stokes equations'', in preparation].
ISSN: 0885-7474

hal-00730193
https://hal.archives-ouvertes.fr/hal-00730193

Éditeur(s) : HAL CCSD Wiley
Résumé : International audience
We study if the multilevel algorithm introduced by Debussche, T. Dubois and R. Temam [Theoret. Comput. Fluid Dynam. 7 (1995), no. 4, 279-315; Zbl 838.76060] and Dubois, F. Jauberteau and Temam [J. Sci. Comput. 8 (1993), no. 2, 167-194; MR1242960 (94f:65098)] for the 2D Navier-Stokes equations with periodic boundary conditions and spectral discretization can be generalized to more general boundary conditions and to finite elements. We first show that a direct generalization, as in [C. Calgaro, J. Laminie and R. Temam, Appl. Numer. Math. 23 (1997), no. 4, 403-442; MR1453424 (98d:76124)], for the Burgers equation, would not be very efficient. We then propose a new approach where the domain of integration is decomposed into subdomains. This enables us to define localized small-scale components and we show that, in this context, there is a good separation of scales. We conclude that all the ingredients necessary for the implementation of the multilevel algorithm are present.
ISSN: 0271-2091

hal-00730709
https://hal.archives-ouvertes.fr/hal-00730709