SECOND ORDER NUMERICAL METHODS FOR NAVIER-STOKES AND DARCY-BRINKMAN EQUATIONS

Download
2022-6-28
DEMİR, Medine
In this thesis, second-order, efficient and reliable numerical stabilization methods are considered for approximating solutions to the incompressible, viscous fluid flow driven by the Navier-Stokes equations and for the Darcy-Brinkman equations driven by double-diffusive convection. The standard Galerkin finite element method remains insufficient for accurately solving these complex nonlinear equations that creates some problems such as numerical instabilities and unphysical oscillations in the solution. A good numerical algorithm should resolve all the scales in the solution to avoid these problems which requires too much computational effort. Thus, developing proper and efficient numerical algorithm that exhibits correct physical behaviour of the flow and accurately approximates solutions over a finite time interval remains a great challenge in computational fluid dynamics. First, this thesis proposes a numerical scheme which tests and analyzes a subgrid artificial viscosity method to model the incompressible Navier-Stokes equations along a linearly extrapolated BDF2 time discretization method. The method considers the viscous term as a combination of the vorticity and the grad-div stabilization term. The method introduces global stabilization by adding a term, then antidiffuses through the extra mixed variables. A detailed analysis of conservation laws, including both energy and helicity balance of the method is presented. It is shown that the approximate solutions of the method are unconditionally stable and optimally convergent. Several numerical tests are presented for validating the support of the derived theoretical results. Second, this thesis considers the backward Euler based linear time filtering method for the developed energy-momentum-angular momentum conserving formulation of the time dependent, incompressible Navier-Stokes equations in the case of weakly enforced divergence constraint. The method adds time filtering as a post-processing step to the energy-momentum-angular momentum conserving formulation to enhance the accuracy and to improve the approximate solutions. It is shown that in comparison with the backward-Euler based energy-momentum-angular momentum conserving formulation without any filter, the proposed method not only leads to a 2-step, unconditionally stable and second order accurate method but also increases numerical accuracy of solutions. Numerical studies verify the theoretical findings and demonstrate preeminence of the proposed method over the unfiltered case. Third, this thesis studies an efficient, accurate, effective and unconditionally stable time stepping scheme for the Darcy-Brinkman equations in double-diffusive convection. The stabilization within the proposed method uses the idea of stabilizing the curvature for velocity, temperature and concentration equations. Accuracy in time is proven and the convergence results for the fully discrete solution of problem variables is given. Several numerical examples including a convergence study are provided that support the derived theoretical results and demonstrate the efficiency and the accuracy of the method.

Suggestions

SECOND ORDER OSCILLATION OF MIXED NONLINEAR DYNAMIC EQUATIONS WITH SEVERAL POSITIVE AND NEGATIVE COEFFICIENTS
ÖZBEKLER, ABDULLAH; Zafer, Ağacık (2011-09-01)
New oscillation criteria are obtained for superlinear and sublinear forced dynamic equations having positive and negative coefficients by means of nonprincipal solutions.
Linear-linear basis functions for MLFMA solutions of magnetic-field and combined-field integral equations
Ergül, Özgür Salih (2007-04-01)
We present the linear-linear (LL) basis functions to improve the accuracy of the magnetic-field integral equation (MFIE) and the combined-field integral equation (CFIE) for three-dimensional electromagnetic scattering problems involving closed conductors. We consider the solutions of relatively large scattering problems by employing the multilevel fast multipole algorithm. Accuracy problems of MFIE and CFIE arising from their implementations with the conventional Rao-Wilton-Glisson (RWG) basis functions can...
Hybrid Surface Integral Equations for Optimal Analysis of Perfectly Conducting Bodies
Karaosmanoglu, Bariscan; Ergül, Özgür Salih (2015-07-24)
We consider hybrid formulations involving simultaneous applications of the electric-field integral equation (EFIE), the magnetic-field integral equation (MFIE), and the combined-field integral equation (CFIE) for the electromagnetic analysis of three-dimensional conductors with arbitrary geometries. By selecting EFIE, MFIE, and CFIE regions on a given object, and optimizing these regions in accordance with the simulation requirements, one can construct an optimal hybrid-field integral equation (HFIE) that p...
Generalized Hybrid Surface Integral Equations for Finite Periodic Perfectly Conducting Objects
Karaosmanoglu, Bariscan; Ergül, Özgür Salih (2017-01-01)
Hybrid formulations that are based on simultaneous applications of diversely weighted electric-field integral equation (EFIE) and magnetic-field integral equation (MFIE) on periodic but finite structures involving perfectly conducting surfaces are presented. Formulations are particularly designed for closed conductors by considering the unit cells of periodic structures as sample problems for optimizing EFIE and MFIE weights in selected regions. Three-region hybrid formulations, which are designed by geneti...
Oscillatory behavior of integro-dynamic and integral equations on time scales
Grace, S. R.; Zafer, Ağacık (2014-02-01)
By making use of asymptotic properties of nonoscillatory solutions, the oscillation behavior of solutions for the integro-dynamic equation
Citation Formats
M. DEMİR, “SECOND ORDER NUMERICAL METHODS FOR NAVIER-STOKES AND DARCY-BRINKMAN EQUATIONS,” Ph.D. - Doctoral Program, Middle East Technical University, 2022.