No Arabic abstract
Several applications in the scientific simulation of physical systems can be formulated as control/optimization problems. The computational models for such systems generally contain hyperparameters, which control solution fidelity and computational expense. The tuning of these parameters is non-trivial and the general approach is to manually `spot-check for good combinations. This is because optimal hyperparameter configuration search becomes impractical when the parameter space is large and when they may vary dynamically. To address this issue, we present a framework based on deep reinforcement learning (RL) to train a deep neural network agent that controls a model solve by varying parameters dynamically. First, we validate our RL framework for the problem of controlling chaos in chaotic systems by dynamically changing the parameters of the system. Subsequently, we illustrate the capabilities of our framework for accelerating the convergence of a steady-state CFD solver by automatically adjusting the relaxation factors of discretized Navier-Stokes equations during run-time. The results indicate that the run-time control of the relaxation factors by the learned policy leads to a significant reduction in the number of iterations for convergence compared to the random selection of the relaxation factors. Our results point to potential benefits for learning adaptive hyperparameter learning strategies across different geometries and boundary conditions with implications for reduced computational campaign expenses. footnote{Data and codes available at url{https://github.com/Romit-Maulik/PAR-RL}}
Reduced Order Modeling (ROM) for engineering applications has been a major research focus in the past few decades due to the unprecedented physical insight into turbulence offered by high-fidelity CFD. The primary goal of a ROM is to model the key physics/features of a flow-field without computing the full Navier-Stokes (NS) equations. This is accomplished by projecting the high-dimensional dynamics to a low-dimensional subspace, typically utilizing dimensionality reduction techniques like Proper Orthogonal Decomposition (POD), coupled with Galerkin projection. In this work, we demonstrate a deep learning based approach to build a ROM using the POD basis of canonical DNS datasets, for turbulent flow control applications. We find that a type of Recurrent Neural Network, the Long Short Term Memory (LSTM) which has been primarily utilized for problems like speech modeling and language translation, shows attractive potential in modeling temporal dynamics of turbulence. Additionally, we introduce the Hurst Exponent as a tool to study LSTM behavior for non-stationary data, and uncover useful characteristics that may aid ROM development for a variety of applications.
In this paper, a coupling lattice Boltzmann (LB) model for simulating thermal flows on the standard D2Q9 lattice is developed in the framework of the double-distribution-function (DDF) approach in which the viscous heat dissipation and compression work are considered. In the model, a density distribution function is used to simulate the flow field, while a total energy distribution function is employed to simulate the temperature field. The discrete equilibrium density and total energy distribution functions are obtained from the Hermite expansions of the corresponding continuous equilibrium distribution functions. The pressure given by the equation of state of perfect gases is recovered in the macroscopic momentum and energy equations. The coupling between the momentum and energy transports makes the model applicable for general thermal flows such as non-Boussinesq flows, while the existing DDF LB models on standard lattices are usually limited to Boussinesq flows in which the temperature variation is small. Meanwhile, the simple structure and basic advantages of the DDF LB approach are retained. The model is tested by numerical simulations of thermal Couette flow, attenuation-driven acoustic streaming, and natural convection in a square cavity with small and large temperature differences. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.
We develop an operator splitting method to simulate flows of isothermal compressible natural gas over transmission pipelines. The method solves a system of nonlinear hyperbolic partial differential equations (PDEs) of hydrodynamic type for mass flow and pressure on a metric graph, where turbulent losses of momentum are modeled by phenomenological Darcy-Weisbach friction. Mass flow balance is maintained through the boundary conditions at the network nodes, where natural gas is injected or withdrawn from the system. Gas flow through the network is controlled by compressors boosting pressure at the inlet of the adjoint pipe. Our operator splitting numerical scheme is unconditionally stable and it is second order accurate in space and time. The scheme is explicit, and it is formulated to work with general networks with loops. We test the scheme over range of regimes and network configurations, also comparing its performance with performance of two other state of the art implicit schemes.
Soft solids in fluids find wide range of applications in science and engineering, especially in the study of biological tissues and membranes. In this study, an Eulerian finite volume approach has been developed to simulate fully resolved incompressible hyperelastic solids immersed in a fluid. We have adopted the recently developed reference map technique (RMT) by Valkov et. al (J. Appl. Mech., 82, 2015) and assessed multiple improvements for this approach.These modifications maintain the numerical robustness of the solver and allow the simulations without any artificial viscosity in the solid regions (to stabilize the solver). This has also resulted in eliminating the striations (wrinkles) of the fluid-solid interface that was seen before and hence obviates the need for any additional routines to achieve a smooth interface. An approximate projection method has been used to project the velocity field onto a divergence free field. Cost and accuracy improvements of the modifications on the method have also been discussed.
The general synthetic iteration scheme (GSIS) is extended to find the steady-state solution of nonlinear gas kinetic equation, removing the long-standing problems of slow convergence and requirement of ultra-fine grids in near-continuum flows. The key ingredients of GSIS are that the gas kinetic equation and macroscopic synthetic equations are tightly coupled, and the constitutive relations in macroscopic synthetic equations explicitly contain Newtons law of shear stress and Fouriers law of heat conduction. The higher-order constitutive relations describing rarefaction effects are calculated from the velocity distribution function, however, their constructions are simpler than our previous work (Su et al. Journal of Computational Physics 407 (2020) 109245) for linearized gas kinetic equations. On the other hand, solutions of macroscopic synthetic equations are used to inform the evolution of gas kinetic equation at the next iteration step. A rigorous linear Fourier stability analysis in periodic system shows that the error decay rate of GSIS can be smaller than 0.5, which means that the deviation to steady-state solution can be reduced by 3 orders of magnitude in 10 iterations. Other important advantages of the GSIS are (i) it does not rely on the specific form of Boltzmann collision operator and (ii) it can be solved by sophisticated techniques in computational fluid dynamics, making it amenable to large scale engineering applications. In this paper, the efficiency and accuracy of GSIS is demonstrated by a number of canonical test cases in rarefied gas dynamics.