scholarly journals An Explicit Meshless Point Collocation Solver for Incompressible Navier-Stokes Equations

Fluids ◽  
2019 ◽  
Vol 4 (3) ◽  
pp. 164 ◽  
Author(s):  
Bourantas ◽  
Zwick ◽  
Joldes ◽  
Loukopoulos ◽  
Tavner ◽  
...  
Keyword(s):  
Stokes Equations ◽  
Two Dimensions ◽  
Navier Stokes ◽  
Benchmark Problems ◽  
Implicit Methods ◽  
Time Step ◽  
Flow Equations ◽  
Point Collocation ◽  
Flow Past ◽  
Flow Past A Cylinder

We present a strong form, meshless point collocation explicit solver for the numerical solution of the transient, incompressible, viscous Navier-Stokes (N-S) equations in two dimensions. We numerically solve the governing flow equations in their stream function-vorticity formulation. We use a uniform Cartesian embedded grid to represent the flow domain. We discretize the governing equations using the Meshless Point Collocation (MPC) method. We compute the spatial derivatives that appear in the governing flow equations, using a novel interpolation meshless scheme, the Discretization Corrected Particle Strength Exchange (DC PSE). We verify the accuracy of the numerical scheme for commonly used benchmark problems including lid-driven cavity flow, flow over a backward-facing step and unbounded flow past a cylinder. We have examined the applicability of the proposed scheme by considering flow cases with complex geometries, such as flow in a duct with cylindrical obstacles, flow in a bifurcated geometry, and flow past complex-shaped obstacles. Our method offers high accuracy and excellent computational efficiency as demonstrated by the verification examples, while maintaining a stable time step comparable to that used in unconditionally stable implicit methods. We estimate the stable time step using the Gershgorin circle theorem. The stable time step can be increased through the increase of the support domain of the weight function used in the DC PSE method.

scholarly journals Computational Algorithms for Solving Spectral/$hp$ Stabilized Incompressible Flow Problems

2016 ◽  
Vol 8 (4) ◽  
pp. 21 ◽  
Author(s):  
Rakesh Ranjan ◽  
Anthony Theodore Chronopoulos ◽  
Yusheng Feng
Keyword(s):  
Stokes Equations ◽  
Three Dimensional ◽  
Higher Order ◽  
Navier Stokes ◽  
Krylov Methods ◽  
Higher Order Methods ◽  
The Past ◽  
Flow Past ◽  
Flow Past A Cylinder ◽  
Computationally Intensive

In this paper we implement the element-by-element preconditioner and inexact Newton-Krylov methods (developed in the past) for solving stabilized computational fluid dynamics (CFD) problems with spectral methods. Two different approaches are implemented for speeding up the process of solving both steady and unsteady incompressible Navier-Stokes equations. The first approach concerns the application of a scalable preconditioner namely the element by element LU preconditioner, while the second concerns the application of Newton-Krylov (NK) methods for solving non-linear problems. We obtain good agreement with benchmark results on standard CFD problems for various Reynolds numbers. We solve the Kovasznay flow and flow past a cylinder at Re-$100$ with this approach. We also utilize the Newton-Krylov algorithm to solve (in parallel) important model problems such as flow past a circular obstacle in a Newtonian flow field, three dimensional driven cavity, flow past a three dimensional cylinder with different immersion lengths. We explore the scalability and robustness of the formulations for both approaches and obtain very good speedup. Effective implementations of these procedures demonstrate for relatively coarse macro-meshes<br />the power of higher order methods in obtaining highly accurate results in CFD. While the procedures adopted in the paper have been explored in the past the novelty lies with applications with higher order methods which have been known to be computationally intensive.

scholarly journals Comparison of Higher-Order Numerical Schemes and Several Filtering Methods Applied to Navier-Stokes Equations with Applications to Computational Aeroacoustics

2009 ◽  
Vol 3 (3) ◽  
pp. 443-459
Author(s):  
L.S. Lai ◽  
G.S. Djambazov ◽  
C.-H. Lai ◽  
K.A. Pericleous
Keyword(s):  
Stokes Equations ◽  
Small Time ◽  
High Order ◽  
Navier Stokes ◽  
Acoustic Analogy ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Flow Equations ◽  
High Order Schemes ◽  
Fine Mesh

In computational acoustics, fluid-acoustic coupling methods for the computation of sound have been widely used by researchers for the last five decades. In the first part of the coupling procedure, the fully unsteady incompressible or compressible flow equations for the near-field of the unsteady flow are solved by using a Computational Fluid Dynamics (CFD) technique, such as Direct Numerical Simulation (DNS), Large Eddy Simulation (LES) or unsteady Reynolds averaged Navier-Stokes equations (RANS) the CFD predictions are then used to calculate sound sources using the acoustic analogy or solving a set of acoustic perturbation equations (APE) leading to the solution of the acoustic field. It is possible to use a 2-D reduced problem to provide a preliminary understanding of many acoustic problems. Unfortunately 2-D CFD simulations using a fine-mesh-small-time-step-LES-alike numerical method cannot be considered as LES, which applies to 3-D simulations only. Therefore it is necessary to understand the similarities and the effect between filters applied to unsteady compressible Navier-Stokes equations and the combined effect of high-order schemes and mesh size. The aim of this study is to provide suitable LES-alike methods for 2-D simulations. An efficient software implementation of high-order schemes is also proposed. Numerical examples are provided to illustrate these statistical similarities.

Computational Simulation of the Tethered Undersea Kites for Power Generation

2015 ◽  
Author(s):  
Amirmahdi Ghasemi ◽  
David J. Olinger ◽  
Gretar Tryggvason
Keyword(s):  
Power Output ◽  
High Speed ◽  
Stokes Equations ◽  
Computational Simulation ◽  
Hydrodynamic Forces ◽  
Two Dimensions ◽  
Ocean Current ◽  
Navier Stokes ◽  
Immersed Boundary ◽  
Flow Equations

The dynamic motion of tethered undersea kites (TUSK) is studied using numerical simulations. TUSK systems consist of a rigid-winged kite moving in an ocean current. The kite is connected by tethers to a platform on the ocean surface or anchored to the seabed. Hydrodynamic forces generated by the kite are transmitted through the tethers to a generator on the platform to produce electricity. TUSK systems are being considered as an alternative to marine turbines since the kite can move in high speed motions to increase power production compared to conventional marine turbines. The two-dimensional Navier-Stokes equations are solved on a regular structured grid that comprises the ocean current flow, and an immersed boundary method is used for the rigid kite. A two-step projection method along with Open Multi-Processing (OpenMP) is employed to solve the flow equations. The reel-out and reel-in velocities of the two tethers are adjusted to control the kite angle of attack and the resultant hydrodynamic forces. A baseline simulation was studied where a high net power output was achieved during successive kite power and retraction phases. System power output, vorticity flow fields, tether tensions, and hydrodynamic coefficients for the kite are determined. The power output results are in good agreement with established theoretical results for a kite moving in two dimensions.

scholarly journals The Sensitivity Analysis for the Flow Past Obstacles Problem with Respect to the Reynolds Number

2012 ◽  
Vol 4 (1) ◽  
pp. 21-35 ◽  
Author(s):  
Kazufumi Ito ◽  
Zhilin Li ◽  
Zhonghua Qiao
Keyword(s):  
Sensitivity Analysis ◽  
Reynolds Number ◽  
Critical Reynolds Number ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Irregular Domains ◽  
Time Step ◽  
Immersed Interface Method ◽  
Navier Stokes Equations ◽  
Flow Past

AbstractIn this paper, numerical sensitivity analysis with respect to the Reynolds number for the flow past obstacle problem is presented. To carry out such analysis, at each time step, we need to solve the incompressible Navier-Stokes equations on irregular domains twice, one for the primary variables; the other is for the sensitivity variables with homogeneous boundary conditions. The Navier-Stokes solver is the augmented immersed interface method for Navier-Stokes equations on irregular domains. One of the most important contribution of this paper is that our analysis can predict the critical Reynolds number at which the vortex shading begins to develop in the wake of the obstacle. Some interesting experiments are shown to illustrate how the critical Reynolds number varies with different geometric settings.

Flow Past a Forced Oscillating Cylinder: A Three-Dimensional Numerical Study

2019 ◽  
Author(s):  
Huan Ping ◽  
Yan Bao ◽  
Dai Zhou ◽  
Zhaolong Han
Keyword(s):  
Stokes Equations ◽  
Numerical Study ◽  
Three Dimensional ◽  
Flow Dynamics ◽  
Navier Stokes ◽  
Two Dimensional ◽  
Navier Stokes Equations ◽  
Flow Past ◽  
Flow Past A Cylinder ◽  
Two Dimensional Flow

Abstract In this paper, we conducted a three-dimensional investigation of flow past a cylinder undergoing forced oscillation. The flow configuration is similar to the work of Blackburn & Henderson (1999) [1], in which Reynolds number equals to 500 and a fixed motion amplitude of A/D = 0.25. The oscillation frequencies are varied in the range near to the natural shedding frequency of a stationary cylinder. The flow dynamics are governed by Navier-Stokes equations and the solutions are obtained by employing high-order spectral/hp element method. It is found that the flow dynamics are significantly distinguished from the study of two-dimensional flow by Blackburn & Henderson (1999) [1]. The values of hydrodynamic forces are smaller compared to that in the two-dimensional study. However, lock-in boundary we identified is broader. In addition, a different type of hysteresis loop of energy transfer coefficient is obtained.

An All-Speed Asymptotic-Preserving Method for the Isentropic Euler and Navier-Stokes Equations

2012 ◽  
Vol 12 (4) ◽  
pp. 955-980 ◽  
Author(s):  
Jeffrey Haack ◽  
Shi Jin ◽  
Jian‐Guo Liu
Keyword(s):  
Mach Number ◽  
Hyperbolic System ◽  
Stokes Equations ◽  
Two Dimensions ◽  
Shock Capturing ◽  
Navier Stokes ◽  
Time Step ◽  
Asymptotic Preserving ◽  
Navier Stokes Equations ◽  
Projection System

AbstractThe computation of compressible flows becomes more challenging when the Mach number has different orders of magnitude. When the Mach number is of order one, modern shock capturing methods are able to capture shocks and other complex structures with high numerical resolutions. However, if the Mach number is small, the acoustic waves lead to stiffness in time and excessively large numerical viscosity, thus demanding much smaller time step and mesh size than normally needed for incompressible flow simulation. In this paper, we develop an all-speed asymptotic preserving (AP) numerical scheme for the compressible isentropic Euler and Navier-Stokes equations that is uniformly stable and accurate for all Mach numbers. Our idea is to split the system into two parts: one involves a slow, nonlinear and conservative hyperbolic system adequate for the use of modern shock capturing methods and the other a linear hyperbolic system which contains the stiff acoustic dynamics, to be solved implicitly. This implicit part is reformulated into a standard pressure Poisson projection system and thus possesses sufficient structure for efficient fast Fourier transform solution techniques. In the zero Mach number limit, the scheme automatically becomes a projection method-like incompressible solver. We present numerical results in one and two dimensions in both compressible and incompressible regimes.

scholarly journals Revisiting Surface-Subsurface Exchange at Intertidal Zone with a Coupled 2D Hydrodynamic and 3D Variably-Saturated Groundwater Model

Water ◽  
2021 ◽  
Vol 13 (7) ◽  
pp. 902
Author(s):  
Zhi Li ◽  
Ben R. Hodges
Keyword(s):  
Intertidal Zone ◽  
Coastal Wetlands ◽  
High Performance ◽  
Stokes Equations ◽  
Wave Model ◽  
Surface Flow ◽  
Richards Equation ◽  
Navier Stokes ◽  
Dynamic Surface ◽  
Flow Equations

A new high-performance numerical model (Frehg) is developed to simulate water flow in shallow coastal wetlands. Frehg solves the 2D depth-integrated, hydrostatic, Navier–Stokes equations (i.e., shallow-water equations) in the surface domain and the 3D variably-saturated Richards equation in the subsurface domain. The two domains are asynchronously coupled to model surface-subsurface exchange. The Frehg model is applied to evaluate model sensitivity to a variety of simplifications that are commonly adopted for shallow wetland models, especially the use of the diffusive wave approximation in place of the traditional Saint-Venant equations for surface flow. The results suggest that a dynamic model for momentum is preferred over diffusive wave model for shallow coastal wetlands and marshes because the latter fails to capture flow unsteadiness. Under the combined effects of evaporation and wetting/drying, using diffusive wave model leads to discrepancies in modeled surface-subsurface exchange flux in the intertidal zone where strong exchange processes occur. It indicates shallow wetland models should be built with (i) dynamic surface flow equations that capture the timing of inundation, (ii) complex topographic features that render accurate spatial extent of inundation, and (iii) variably-saturated subsurface flow solver that is capable of modeling moisture change in the subsurface due to evaporation and infiltration.

A finite element method for the Navier-Stokes equations in moving domain with application to hemodynamics of the left ventricle

2017 ◽  
Vol 32 (4) ◽  
Author(s):  
Alexander Danilov ◽  
Alexander Lozovskiy ◽  
Maxim Olshanskii ◽  
Yuri Vassilevski
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Left Ventricle ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Time Step ◽  
Navier Stokes Equations ◽  
Computed Tomography Images ◽  
Time Dependent Domain ◽  
Element Method

AbstractThe paper introduces a finite element method for the Navier-Stokes equations of incompressible viscous fluid in a time-dependent domain. The method is based on a quasi-Lagrangian formulation of the problem and handling the geometry in a time-explicit way. We prove that numerical solution satisfies a discrete analogue of the fundamental energy estimate. This stability estimate does not require a CFL time-step restriction. The method is further applied to simulation of a flow in a model of the left ventricle of a human heart, where the ventricle wall dynamics is reconstructed from a sequence of contrast enhanced Computed Tomography images.

Essential Ingredients for the Computation of Steady and Unsteady Blade Boundary Layers

1991 ◽  
Vol 113 (4) ◽  
pp. 608-616 ◽  
Author(s):  
H. M. Jang ◽  
J. A. Ekaterinaris ◽  
M. F. Platzer ◽  
T. Cebeci
Keyword(s):  
Boundary Layer ◽  
Boundary Layers ◽  
Stokes Equations ◽  
Inviscid Flow ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Transitional Region ◽  
Low Reynolds Numbers ◽  
Flow Equations ◽  
Low Reynolds

Two methods are described for calculating pressure distributions and boundary layers on blades subjected to low Reynolds numbers and ramp-type motion. The first is based on an interactive scheme in which the inviscid flow is computed by a panel method and the boundary layer flow by an inverse method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier–Stokes equations with an embedded grid technique that permits accurate calculation of boundary layer flows. Studies for the Eppler-387 and NACA-0012 airfoils indicate that both methods can be used to calculate the behavior of unsteady blade boundary layers at low Reynolds numbers provided that the location of transition is computed with the en method and the transitional region is modeled properly.

Sign in / Sign up

Export Citation Format

Share Document