Calculation of Hydrodynamic Forces for Unsteady Stokes Flows by Singularity Integral Equations Based on Fundamental Solutions

2013 ◽  
Vol 30 (2) ◽  
pp. 129-136 ◽  
Author(s):  
C. H. Hsiao ◽  
D. L. Young
Keyword(s):  
Time Domain ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Fundamental Solutions ◽  
Hydrodynamic Forces ◽  
Stokes Flows ◽  
Total Strength ◽  
Singularity Method ◽  
Unsteady Stokes Flow ◽  
Convolution Quadrature Method

ABSTRACTThe attractive feature of the singularity method for steady Stokes flows is that the hydrodynamic forces acting on the particle can be calculated by the total strength of distributed singularities. For unsteady Stokes flows, however we have to derive hydrodynamic forces acting on a solid body in terms of the strengths of both unsteady Stokeslets as well as unsteady potential dipoles if mass and force sources are both taken into consideration. Since the hydrodynamic force formulation results in a Volterra integral equation of the first kind, and the strengths are numerically approximated by means of the Lubich convolution quadrature method (CQM) in this study. As far as the numerical solutions of time-domain integral formulations of the unsteady Stokes equations are concerned, this paper requires only the Laplace-domain instead of the time- domain fundamental solutions of the governing equations. The stability and accuracy of the proposed method are verified through some well selected numerical examples. In total we include two examples presenting the accuracy of Lubich CQM, and another two examples for calculating general hydrody-namic forces of a sphere in oscillating or non-oscillating unsteady Stokes flows. It is concluded that this study is able to extend the unsteady Stokes flow theory to more general transient motions instead to limit to the oscillating flow assumption.

The Derivation and Application of Fundamental Solutions for Unsteady Stokes Equations

2015 ◽  
Vol 31 (6) ◽  
pp. 683-691 ◽  
Author(s):  
C-H. Hsiao ◽  
D.-L. Young
Keyword(s):  
Stokes Equations ◽  
Three Dimensional ◽  
Fundamental Solutions ◽  
Mass Source ◽  
Reynolds Number Flow ◽  
Singularity Method ◽  
Vector Calculus ◽  
Unsteady Stokes Flow ◽  
Number Flow ◽  
Unsteady Stokes Equations

AbstractIn this paper, two formulations in explicit form to derive the fundamental solutions for two and three dimensional unsteady unbounded Stokes flows due to a mass source and a point force are presented, based on the vector calculus method and also the Hörmander’s method. The mathematical derivation process for the fundamental solutions is detailed. The steady fundamental solutions of Stokes equations can be obtained from the unsteady fundamental solutions by the integral process. As an application, we adopt fundamental solutions: an unsteady Stokeslet and an unsteady potential dipole to validate a simple case that a sphere translates in Stokes or low-Reynolds-number flow by using the singularity method instead by the traditional method which in general limits to the assumption of oscillating flow. It is concluded that this study is able to extend the unsteady Stokes flow theory to more general transient motions by making use of the fundamental solutions of the linearly unsteady Stokes equations.

scholarly journals Hydrodynamics of self-propulsion near a boundary: predictions and accuracy of far-field approximations

2012 ◽  
Vol 700 ◽  
pp. 105-147 ◽  
Author(s):  
Saverio E. Spagnolie ◽  
Eric Lauga
Keyword(s):  
General Framework ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Fundamental Solutions ◽  
Reduced Model ◽  
Far Field ◽  
Free Surfaces ◽  
Source Dipole ◽  
Swimming Trajectories ◽  
Propulsive Mechanism

AbstractThe swimming trajectories of self-propelled organisms or synthetic devices in a viscous fluid can be altered by hydrodynamic interactions with nearby boundaries. We explore a multipole description of swimming bodies and provide a general framework for studying the fluid-mediated modifications to swimming trajectories. A general axisymmetric swimmer is described as a linear combination of fundamental solutions to the Stokes equations: a Stokeslet dipole, a source dipole, a Stokeslet quadrupole, and a rotlet dipole. The effects of nearby walls or stress-free surfaces on swimming trajectories are described through the contribution of each singularity, and we address the question of how accurately this multipole approach captures the wall effects observed in full numerical solutions of the Stokes equations. The reduced model is used to provide simple but accurate predictions of the wall-induced attraction and pitching dynamics for model Janus particles, ciliated organisms, and bacteria-like polar swimmers. Transitions in attraction and pitching behaviour as functions of body geometry and propulsive mechanism are described. The reduced model may help to explain a number of recent experimental results.

scholarly journals Hydromechanics of low-Reynolds-number flow. Part 2. Singularity method for Stokes flows

1975 ◽  
Vol 67 (4) ◽  
pp. 787-815 ◽  
Author(s):  
Allen T. Chwang ◽  
T. Yao-Tsu Wu
Keyword(s):  
Exact Solutions ◽  
Stokes Flow ◽  
Fundamental Solutions ◽  
The Body ◽  
Circular Cylinders ◽  
Geometrical Parameters ◽  
Flow Problems ◽  
Singularity Method ◽  
Wide Range ◽  
Body Shapes

The present study further explores the fundamental singular solutions for Stokes flow that can be useful for constructing solutions over a wide range of free-stream profiles and body shapes. The primary singularity is the Stokeslet, which is associated with a singular point force embedded in a Stokes flow. From its derivatives other fundamental singularities can be obtained, including rotlets, stresslets, potential doublets and higher-order poles derived from them. For treating interior Stokes-flow problems new fundamental solutions are introduced; they include the Stokeson and its derivatives, called the roton and stresson.These fundamental singularities are employed here to construct exact solutions to a number of exterior and interior Stokes-flow problems for several specific body shapes translating and rotating in a viscous fluid which may itself be providing a primary flow. The different primary flows considered here include the uniform stream, shear flows, parabolic profiles and extensional flows (hyper-bolic profiles), while the body shapes cover prolate spheroids, spheres and circular cylinders. The salient features of these exact solutions (all obtained in closed form) regarding the types of singularities required for the construction of a solution in each specific case, their distribution densities and the range of validity of the solution, which may depend on the characteristic Reynolds numbers and governing geometrical parameters, are discussed.

Machine learning–based reduced-order modeling of hydrodynamic forces using pressure mode decomposition

2021 ◽  
pp. 095441002199986
Author(s):  
Hassan F Ahmed ◽  
Hamayun Farooq ◽  
Imran Akhtar ◽  
Zafar Bangash
Keyword(s):  
Machine Learning ◽  
Short Term Memory ◽  
Stokes Equations ◽  
Decomposition Analysis ◽  
Hydrodynamic Forces ◽  
Reduced Order Modeling ◽  
Navier Stokes ◽  
Reduced Order ◽  
Mode Decomposition ◽  
Lift And Drag

In this article, we introduce a machine learning–based reduced-order modeling (ML-ROM) framework through the integration of proper orthogonal decomposition (POD) and deep neural networks (DNNs), in addition to long short-term memory (LSTM) networks. The DNN is utilized to upscale POD temporal coefficients and their respective spatial modes to account for the dynamics represented by the truncated modes. In the second part of the algorithm, temporal evolution of the POD coefficients is obtained by recursively predicting their future states using an LSTM network. The proposed model (ML-ROM) is tested for flow past a circular cylinder characterized by the Navier–Stokes equations. We perform pressure mode decomposition analysis on the flow data using both POD and ML-ROM to predict hydrodynamic forces and demonstrate the accuracy of the proposed strategy for modeling lift and drag coefficients.

scholarly journals Modeling Transient Flows in Heterogeneous Layered Porous Media Using the Space–Time Trefftz Method

2021 ◽  
Vol 11 (8) ◽  
pp. 3421
Author(s):  
Cheng-Yu Ku ◽  
Li-Dan Hong ◽  
Chih-Yu Liu ◽  
Jing-En Xiao ◽  
Wei-Po Huang
Keyword(s):  
Porous Media ◽  
Time Domain ◽  
Numerical Solutions ◽  
Space Time ◽  
Two Dimensional ◽  
Transient Flows ◽  
Layered Porous Media ◽  
Time Marching ◽  
Time Basis ◽  
Marching Scheme

In this study, we developed a novel boundary-type meshless approach for dealing with two-dimensional transient flows in heterogeneous layered porous media. The novelty of the proposed method is that we derived the Trefftz space–time basis function for the two-dimensional diffusion equation in layered porous media in the space–time domain. The continuity conditions at the interface of the subdomains were satisfied in terms of the domain decomposition method. Numerical solutions were approximated based on the superposition principle utilizing the space–time basis functions of the governing equation. Using the space–time collocation scheme, the numerical solutions of the problem were solved with boundary and initial data assigned on the space–time boundaries, which combined spatial and temporal discretizations in the space–time manifold. Accordingly, the transient flows through the heterogeneous layered porous media in the space–time domain could be solved without using a time-marching scheme. Numerical examples and a convergence analysis were carried out to validate the accuracy and the stability of the method. The results illustrate that an excellent agreement with the analytical solution was obtained. Additionally, the proposed method was relatively simple because we only needed to deal with the boundary data, even for the problems in the heterogeneous layered porous media. Finally, when compared with the conventional time-marching scheme, highly accurate solutions were obtained and the error accumulation from the time-marching scheme was avoided.

A numerical study of the interaction between unsteay free-stream disturbances and localized variations in surface geometry

1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar
Keyword(s):  
Numerical Solutions ◽  
Free Stream ◽  
Stokes Equations ◽  
Numerical Study ◽  
Navier Stokes ◽  
Surface Geometry ◽  
Navier Stokes Equations ◽  
S Waves ◽  
Tollmien Schlichting

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.

Theory of a Time Domain Boundary Element Development for the Dynamic Analysis of Coupled Multiphase Porous Media

2017 ◽  
Vol 08 (03n04) ◽  
pp. 1750007
Author(s):  
Pooneh Maghoul ◽  
Behrouz Gatmiri
Keyword(s):  
Boundary Element ◽  
Time Domain ◽  
Unsaturated Soils ◽  
Water Pressure ◽  
Domain Boundary ◽  
Fundamental Solutions ◽  
Element Formulation ◽  
Inertia Effects ◽  
Coupled Equations ◽  
The Time Domain

This paper presents an advanced formulation of the time-domain two-dimensional (2D) boundary element method (BEM) for an elastic, homogeneous unsaturated soil subjected to dynamic loadings. Unlike the usual time-domain BEM, the present formulation applies a convolution quadrature which requires only the Laplace-domain instead of the time-domain fundamental solutions. The coupled equations governing the dynamic behavior of unsaturated soils ignoring contributions of the inertia effects of the fluids (water and air) are derived based on the poromechanics theory within the framework of a suction-based mathematical model. In this formulation, the solid skeleton displacements [Formula: see text], water pressure [Formula: see text] and air pressure [Formula: see text] are presumed to be independent variables. The fundamental solutions in Laplace transformed-domain for such a dynamic [Formula: see text] theory have been obtained previously by authors. Then, the BE formulation in time is derived after regularization by partial integrations and time and spatial discretizations. Thereafter, the BE formulation is implemented in a 2D boundary element code (PORO-BEM) for the numerical solution. To verify the accuracy of this implementation, the displacement response obtained by the boundary element formulation is verified by comparison with the elastodynamics problem.

Spin-up from rest of a compressible fluid in a rapidly rotating cylinder

1992 ◽  
Vol 237 ◽  
pp. 413-434 ◽  
Author(s):  
Jae Min Hyun ◽  
Jun Sang Park
Keyword(s):  
Mach Number ◽  
Compressible Fluid ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Axial Direction ◽  
Ekman Layer ◽  
Infinite Cylinder ◽  
Full Time ◽  
Viscous Effects ◽  
Initial State

Spin-up flows of a compressible gas in a finite, closed cylinder from an initial state of rest are studied, The flow is characterized by small reference Ekman numbers, and the peripheral Mach number is O(1). Comprehensive numerical solutions have been obtained for the full, time-dependent compressible Navier-Stokes equations. The details of the flow, temperature, and density evolution are described. In the early phase of spin-up, owing to the thermoacoustic disturbances caused by the compressible Rayleigh effect, the flows are oscillatory, and this oscillatory behaviour is pronounced at higher Mach numbers. The principal dynamical role of the Ekman layer is dominant over moderate times of orders of the homogeneous spin-up timescales. Owing to the density stratification in the radial direction, the Ekman layer is thicker in the central region of the interior. The interior azimuthal flows are mainly uniform in the axial direction. As the Mach number increases, the rate of spin-up in the interior becomes slower, and the propagating shear front is more diffusive. Explicit comparisons with the results for an infinite cylinder are made to ascertain the contributions of the endwall disks. In contrast to the usual incompressible spin-up from rest, the viscous effects are relatively more important for the case of a compressible fluid.

Sign in / Sign up

Export Citation Format

Share Document