Automatic Hex-Dominant Mesh Generation for Complex Flow Configurations

2018 ◽  
Vol 11 (6) ◽  
pp. 615-624
Author(s):  
Nihar Sawant ◽  
Soji Yamakawa ◽  
Satbir Singh ◽  
Kenji Shimada
Author(s):  
Yoram Rubin

Spatial variability and the uncertainty in characterizing the flow domain play an important role in the transport of contaminants in porous media: they affect the pathlines followed by solute particles, the spread of solute bodies, the shape of breakthrough curves, the spatial variability of the concentration, and the ability to quantify any of these accurately. This chapter briefly reviews some basic concepts which we shall later employ for the analysis of solute transport in heterogeneous media, and also points out some issues we shall address in the subsequent chapters. Our exposition in chapters 8-10 on contaminant transport is built around the Lagrangian and the Eulerian approaches for analyzing transport. The Eulerian approach is a statement of mass conservation in control volumes of arbitrary dimensions, in the form of the advection-dispersion equation. As such, it is well suited for numerical modeling in complex flow configurations. Its main difficulties, however, are in the assignment of parameters, both hydrogeological and geochemical, to the numerical grid blocks such that the effects of subgrid-scale heterogeneity are accounted for, and in the numerical dispersion that occurs in advection-dominated flow situations. Another difficulty is in the disparity between the scale of the numerical elements and the scale of the samples collected in the field, which makes the interpretation of field data difficult. The Lagrangian approach focuses on the displacements and travel times of solute bodies of arbitrary dimensions, using the displacements of small solute particles along streamlines as its basic building block. Tracking such displacements requires that the solute particles do not transfer across streamlines. Since such mass transfer may only occur due to pore-scale dispersion, Lagrangian approaches are ideally suited for advection-dominated situations. Let us start by considering the displacement of a small solute body, a particle, as a function of time. “Small” here implies that the solute body is much smaller than the characteristic scale of heterogeneity. At the same time, to qualify for a description of its movement using Darcy’s law, the solute body also needs to be larger than a few pores. The small dimension of the solute body ensures that it moves along a single streamline and that it does not disintegrate due to velocity shear.


2015 ◽  
Vol 67 (3) ◽  
Author(s):  
R. J. Lingwood ◽  
P. Henrik Alfredsson

Research on the von Kármán boundary layer extends back almost 100 years but remains a topic of active study, which continues to reveal new results; it is only now that fully nonlinear direct numerical simulations (DNS) have been conducted of the flow to compare with theoretical and experimental results. The von Kármán boundary layer, or rotating-disk boundary layer, provides, in some senses, a simple three-dimensional boundary-layer model with which to compare other more complex flow configurations but we will show that in fact the rotating-disk boundary layer itself exhibits a wealth of complex instability behaviors that are not yet fully understood.


2004 ◽  
Vol 126 (3) ◽  
pp. 382-387 ◽  
Author(s):  
Gordan R. Stuhne and ◽  
David A. Steinman

Background. Computational fluid dynamics (CFD) simulations are used to analyze the wall shear stress distribution and flow streamlines near the throat of a stented basilar side-wall aneurysm. Previous studies of stented aneurysm flows used low mesh resolution, did not include mesh convergence analyses, and depended upon conformal meshing techniques that apply only to very artificial stent geometries. Method of Approach. We utilize general-purpose computer assisted design and unstructured mesh generation tools that apply in principle to stents and vasculature of arbitrary complexity. A mesh convergence analysis for stented steady flow is performed, varying node spacing near the stent. Physiologically realistic pulsatile simulations are then performed using the converged mesh. Results. Artifact-free resolution of the wall shear stress field on stent wires requires a node spacing of approximately 1/3 wire radius. Large-scale flow features tied to the velocity field are, however, captured at coarser resolution (nodes spaced by about one wire radius or more). Conclusions. Results are consistent with previous work, but our methods yield more detailed insights into the complex flow dynamics. However, routine applications of CFD to anatomically realistic cases still depend upon further development of dedicated algorithms, most crucially to handle geometry definition and mesh generation for complicated stent deployments.


Author(s):  
Haiwen Zhu ◽  
Zimo Lin ◽  
Jianlin Peng ◽  
Hong-Quan Zhang ◽  
Jianjun Zhu ◽  
...  

Abstract The performance of multi-stage Electrical Submersible Pumps (ESPs) under different flow conditions and its life span with sand production are commonly predicted by the Computational Fluid Dynamics (CFD) simulations. The mesh generation methodology and optimum grid number are usually validated by pump water catalog curves. Then, the validated mesh geometry is adopted in high viscosity, multiphase flow, and sand erosion simulations to study the effects including but not limited to: discrete phase bubble diameter, turbulence model, body forces, and erosion models. However, the mesh validation by pump water curves is not enough in complex flow conditions, especially in the erosion simulations. Different from the pump hydraulic performance simulation, the accuracy of the erosion simulation can be affected by mesh boundary and inner layer grid thickness, especially for small particles. In addition, the mesh-type (hexahedral and tetrahedral) and size of the inner domain can also significantly affect the particle trajectory. A comprehensive mesh independent study is conducted for water, oil, and gas-liquid conditions of a mixed type ESP in this paper. Then the near-wall inflation layer thickness and inner domain grid size effect to ESP erosion simulation are well analyzed. The mesh generation methodology can be applied to other turbomachinery simulations to improve accuracy.


Author(s):  
Mostafa Esmaeili ◽  
Asghar Afshari

In this study, an Eulerian-Lagrangian computational methodology is utilized for large eddy simulation (LES) of mixing phenomena in jet in cross-flows. A high-order multi-block algorithm is used to solve Eulerian equations in a generalized coordinate system. The composition is formulated based on the filtered mass density function (FMDF) and its equivalent stochastic Lagrangian equations, which is solved by Lagrangian Monte-Carlo method. Parameters influencing mixing enhancement including jet velocity profile, and jet pulsation are investigated. A good consistency between Eulerian and Lagrangian components of the numerical scheme is established. In jet in cross-flow (JICF) simulations, the vortical structures and flow features are predicted with the current numerical scheme. The results also show that the jet velocity profile affects both trajectory and mixing condition and the jet pulsation can enhance mixing depending on the Strouhal numbers. The obtained results including concentration distributions are in good agreement with available experimental data ensuring the performance and reliability of LES/FMDF methodology to study mixing in relatively complex flow configurations such as JICF.


Author(s):  
G. Freskos ◽  
O. Penanhoat

The demand for efficiency in today’s and in future civil aircraft is such that experimental studies alone do not suffice to optimize aircraft aerodynamics. In this context, much effort has been spent in the past decade to develop numerical methods capable of reproducing the phenomena that occur in the engine flow field. This paper presents some studies in Computational Fluid Dynamics related to supersonic inlets. Two approaches are considered. First, there is a need for code capable of calculating in a cost-efficient way the entire flow field around a 2D or 3D inlet, e.g. to perform parametric studies. To this effect, a computing method based on grid construction by mesh generator dedicated to inlet shapes and on the discretization of the unsteady Euler equations with an explicit upwind scheme was developed. The treatment of complex geometries led us to adopt a multiblock grid approach. Therefore particular attention was paid to the treatment of the boundary conditions between the different domains. Secondly, there is a need for code that can capture local phenomena in order to get a better understanding of inlet behaviour (shock/shock, shock/boundary layer interactions, etc.). To this effect a 2D turbulent Navier-Stokes code is used. The 2 equations k-ε turbulence model included in the program seems to be one of the most successful models for calculating flow realistically. Correction of the near-wall influence extends its capability to complex flow configurations, e.g. those with separated zones.


1994 ◽  
Vol 116 (1) ◽  
pp. 116-123 ◽  
Author(s):  
G. Freskos ◽  
O. Penanhoat

The demand for efficiency in today’s and in future civil aircraft is such that experimental studies alone do not suffice to optimize aircraft aerodynamics. In this context, much effort has been spent in the past decade to develop numerical methods capable of reproducing the phenomena that occur in the engine flow field. This paper presents some studies in Computational Fluid Dynamics related to supersonic inlets. Two approaches are considered. First, there is a need for a code capable of calculating in a cost-efficient way the entire flow field around a two-dimensional or three-dimensional inlet, e.g., to perform parametric studies. To this effect, a computing method based on grid construction by mesh generator dedicated to inlet shapes and on the discretization of the unsteady Euler equations with an explicit upwind scheme was developed. The treatment of complex geometries led us to adopt a multiblock grid approach. Therefore particular attention was paid to the treatment of the boundary conditions between the different domains. Second, there is a need for a code that can capture local phenomena in order to get a better understanding of inlet behavior (shock/shock, shock/boundary layer interactions, etc.). To this effect a two-dimensional turbulent Navier-Stokes code is used. The two-equation k-ε turbulence model included in the program seems to be one of the most successful models for calculating flow realistically. Correction of the near-wall influence extends its capability to complex flow configurations, e.g., those with separated zones.


1983 ◽  
Vol 105 (1) ◽  
pp. 105-112 ◽  
Author(s):  
C. Arcoumanis ◽  
A. F. Bicen ◽  
J. H. Whitelaw

Measurements of the three components of velocity and their corresponding fluctuations have been obtained by laser-Doppler anemometry mainly near TDC of compression in a model IC engine motored at 200 rpm with compression ratio of 6.7. The flow configurations comprised an axisymmetric cylinder head with and without upstream induced swirl and each of a flat piston and two centrally located, cylindrical and re-entrant, bowl-in-piston arrangements. In the absence of swirl and squish, the intake-generated mean motion and turbulence decayed considerably by the end of compression. The two piston-bowl configurations, however, resulted in a compression-induced squish motion with consequent formation of a toroidal vortex occupying the whole bowl space. Interacton of swirl, carried from intake and persisting through compression, with squish generated near TDC profoundly altered the axial flow structure. In the case of the cylindrical bowl, the sense of the vortex was reversed by swirl and, in the reentrant bowl, increased the number of vortices to two. The swirling motion inside the cylindrical bowl was close to solid body rotation while the re-entrant bowl gave rise to more complex flow patterns. Squish, in the presence or absence of swirl, did not augment the turbulent energy inside the cylindrical bowl contrary to the reentrant configuration where turbulence generation was observed.


Sign in / Sign up

Export Citation Format

Share Document