A priori reduced order model of steady Navier‐Stokes equations for low Reynolds number

PAMM ◽  
2021 ◽  
Vol 20 (1) ◽  
Author(s):  
Christophe Hoareau ◽  
Lan Shang ◽  
Andreas Zilian
2019 ◽  
Vol 24 (2) ◽  
pp. 45 ◽  
Author(s):  
Nissrine Akkari ◽  
Fabien Casenave ◽  
Vincent Moureau

In the following paper, we consider the problem of constructing a time stable reduced order model of the 3D turbulent and incompressible Navier–Stokes equations. The lack of stability associated with the order reduction methods of the Navier–Stokes equations is a well-known problem and, in general, it is very difficult to account for different scales of a turbulent flow in the same reduced space. To remedy this problem, we propose a new stabilization technique based on an a priori enrichment of the classical proper orthogonal decomposition (POD) modes with dissipative modes associated with the gradient of the velocity fields. The main idea is to be able to do an a priori analysis of different modes in order to arrange a POD basis in a different way, which is defined by the enforcement of the energetic dissipative modes within the first orders of the reduced order basis. This enables us to model the production and the dissipation of the turbulent kinetic energy (TKE) in a separate fashion within the high ranked new velocity modes, hence to ensure good stability of the reduced order model. We show the importance of this a priori enrichment of the reduced basis, on a typical aeronautical injector with Reynolds number of 45,000. We demonstrate the capacity of this order reduction technique to recover large scale features for very long integration times (25 ms in our case). Moreover, the reduced order modeling (ROM) exhibits periodic fluctuations with a period of 2 . 2 ms corresponding to the time scale of the precessing vortex core (PVC) associated with this test case. We will end this paper by giving some prospects on the use of this stable reduced model in order to perform time extrapolation, that could be a strategy to study the limit cycle of the PVC.


1998 ◽  
Vol 120 (2) ◽  
pp. 257-262 ◽  
Author(s):  
Peter Gerlinger ◽  
Dieter Bru¨ggemann

A multigrid method for convergence acceleration is used for solving coupled fluid and turbulence transport equations. For turbulence closure a low-Reynolds-number q-ω turbulence model is employed, which requires very fine grids in the near wall regions. Due to the use of fine grids, convergence of most iterative solvers slows down, making the use of multigrid techniques especially attractive. However, special care has to be taken on the strong nonlinear turbulent source terms during restriction from fine to coarse grids. Due to the hyperbolic character of the governing equations in supersonic flows and the occurrence of shock waves, modifications to standard multigrid techniques are necessary. A simple and effective method is presented that enables the multigrid scheme to converge. A strong reduction in the required number of multigrid cycles and work units is achieved for different test cases, including a Mack 2 flow over a backward facing step.


Author(s):  
Imran Akhtar ◽  
Jeff Borggaard ◽  
John A. Burns ◽  
Lizette Zietsman

Functional gains are integral kernels of the standard feedback operator and are useful in control of partial differential equations (PDEs). These functional gains provide physical insight into how the control mechanism is operating. In some cases, these functional gains can provide information about the optimal placement of actuators and sensors. The study is motivated by fluid flow control and focuses on the computation of these functions. However, for practical purposes, one must be able to compute these functions for a wide variety of PDEs. For higher dimensional systems, computing these gains is at least as challenging as the original simulation problem. To reduce the complexity of the governing equations, reduced-order models are often developed by reducing the PDEs to ordinary-differential equations (ODEs). In this study, we use proper orthogonal decomposition (POD)-Galerkin based approach and develop a reduced-order model of a bluff body wake. We solve the incompressible Navier-Stokes equations, simulate the flow past a circular cylinder, and record the snapshots of the flow field. We compute the POD eigenfunctions and project the Navier-Stokes equations onto these few of these eigenfunctions to develop a reduced-order model. Later, we modify the model by introducing a control function simulating suction actuation on the cylinder surface. We linearize the model about the mean flow and apply feedback control to suppress vortex shedding. We then compute the functional gains for the applied control. We identify these gains at various stations in the wake region and suggest optimum locations for the sensors.


2021 ◽  
Vol 89 ◽  
pp. 1294-1315
Author(s):  
My Ha Dao ◽  
Hoang Huy Nguyen ◽  
Chin Chun Ooi ◽  
Quang Tuyen Le

2007 ◽  
Vol 23 (4) ◽  
pp. 423-432 ◽  
Author(s):  
S.-H. Lee ◽  
Tzuyin Wu

AbstractIn this paper, the drag force on a sphere moving constantly along the centerline of a circular pipe filled with viscous fluid (the falling-sphere problem) under low Reynolds number condition is investigated via numerical calculation. The incompressible Navier-Stokes equations are formulated in a pseudocompressibility form. The numerical scheme makes use of finite-volume method and the numerical flux terms are evaluated using the Total-Variation Diminishing (TVD) strategy commonly applied to the compressible flow. Steady solution is obtained by marching (iterating) in time until the artificial time derivative of pressure term in the continuity equation drops to zero.In the calculation, six different Reynolds number (Re) ranging from 0.1 to 1 and seven different pipe-to-sphere diameter ratios (D/d) ranging from 5 to 40 are selected to study the pipe-wall effect. In each case, the drag force on the sphere is evaluated and the results are compared with the existing approximate theoretical values derived from correcting the Stokes' formula. Both results agree in trend, but with noticeable deviation in values, particularly for cases with large pipe-to-sphere diameter ratios. The deviation is due to the fact that theoretical values were based on the solution to the linearized Navier-Stokes equations (Stokes' creeping-flow equations), while the fully nonlinear form of the Navier-Stokes equations are adopted in the present calculations. Finally, a least-square regression technique is applied to collapse the calculated results into a single expression exhibiting the functional relationship between the drag force, Reynolds number (Re), and the pipe-to-sphere diameter ratio (D/d).


Sign in / Sign up

Export Citation Format

Share Document