A High Resolution Simulation of a Single Shock-Accelerated Particle

Author(s):  
W. Curtis Maxon ◽  
Tanner Nielsen ◽  
Nicholas Denissen ◽  
Johnathan D. Regele ◽  
Jacob McFarland

Abstract Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted using the hydrocode FLAG. FLAG's capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated -- a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG's Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

2021 ◽  
Vol 11 (11) ◽  
pp. 4748
Author(s):  
Monika Balázsová ◽  
Miloslav Feistauer ◽  
Jaromír Horáček ◽  
Adam Kosík

This study deals with the development of an accurate, efficient and robust method for the numerical solution of the interaction of compressible flow and nonlinear dynamic elasticity. This problem requires the reliable solution of flow in time-dependent domains and the solution of deformations of elastic bodies formed by several materials with complicated geometry depending on time. In this paper, the fluid–structure interaction (FSI) problem is solved numerically by the space-time discontinuous Galerkin method (STDGM). In the case of compressible flow, we use the compressible Navier–Stokes equations formulated by the arbitrary Lagrangian–Eulerian (ALE) method. The elasticity problem uses the non-stationary formulation of the dynamic system using the St. Venant–Kirchhoff and neo-Hookean models. The STDGM for the nonlinear elasticity is tested on the Hron–Turek benchmark. The main novelty of the study is the numerical simulation of the nonlinear vocal fold vibrations excited by the compressible airflow coming from the trachea to the simplified model of the vocal tract. The computations show that the nonlinear elasticity model of the vocal folds is needed in order to obtain substantially higher accuracy of the computed vocal folds deformation than for the linear elasticity model. Moreover, the numerical simulations showed that the differences between the two considered nonlinear material models are very small.


2021 ◽  
pp. 146808742110131
Author(s):  
Xiaohang Fang ◽  
Li Shen ◽  
Christopher Willman ◽  
Rachel Magnanon ◽  
Giuseppe Virelli ◽  
...  

In this article, different manifold reduction techniques are implemented for the post-processing of Particle Image Velocimetry (PIV) images from a Spark Ignition Direct Injection (SIDI) engine. The methods are proposed to help make a more objective comparison between Reynolds-averaged Navier-Stokes (RANS) simulations and PIV experiments when Cycle-to-Cycle Variations (CCV) are present in the flow field. The two different methods used here are based on Singular Value Decomposition (SVD) principles where Proper Orthogonal Decomposition (POD) and Kernel Principal Component Analysis (KPCA) are used for representing linear and non-linear manifold reduction techniques. To the authors’ best knowledge, this is the first time a non-linear manifold reduction technique, such as KPCA, has ever been used in the study of in-cylinder flow fields. Both qualitative and quantitative studies are given to show the capability of each method in validating the simulation and incorporating CCV for each engine cycle. Traditional Relevance Index (RI) and two other previously developed novel indexes: the Weighted Relevance Index (WRI) and the Weighted Magnitude Index (WMI), are used for the quantitative study. The results indicate that both POD and KPCA show improvements in capturing the main flow field features compared to ensemble-averaged PIV experimental data and single cycle experimental flow fields while capturing CCV. Both methods present similar quantitative accuracy when using the three indexes. However, challenges were highlighted in the POD method for the selection of the number of POD modes needed for a representative reconstruction. When the flow field region presents a Gaussian distribution, the KPCA method is seen to provide a more objective numerical process as the reconstructed flow field will see convergence with an increasing number of modes due to its usage of Gaussian properties. No additional criterion is needed to determine how to reconstruct the main flow field feature. Using KPCA can, therefore, reduce the amount of analysis needed in the process of extracting the main flow field while incorporating CCV.


1974 ◽  
Vol 96 (4) ◽  
pp. 394-400 ◽  
Author(s):  
V. A. Marple ◽  
B. Y. H. Liu ◽  
K. T. Whitby

The flow field in an inertial impactor was studied experimentally with a water model by means of a flow visualization technique. The influence of such parameters as Reynolds number and jet-to-plate distance on the flow field was determined. The Navier-Stokes equations describing the laminar flow field in the impactor were solved numerically by means of a finite difference relaxation method. The theoretical results were found to be in good agreement with the empirical observations made with the water model.


2013 ◽  
Vol 291-294 ◽  
pp. 1981-1984
Author(s):  
Zhang Xia Guo ◽  
Yu Tian Pan ◽  
Yong Cun Wang ◽  
Hai Yan Zhang

Gunpowder was released in an instant when the pill fly out of the shell during the firing, and then formed a complicated flow fields about the muzzle when the gas expanded sharply. Using the 2 d axisymmetric Navier-Stokes equation combined with single equation turbulent model to conduct the numerical simulation of the process of gunpowder gass evacuating out of the shell without muzzle regardless of the pill’s movement. The numerical simulation result was identical with the experimental. Then simulated the evacuating process of gunpowder gass of an artillery with muzzle brake. The result showed complicated wave structure of the flow fields with the muzzle brake and analysed the influence of muzzle brake to the gass flow field distribution.


2016 ◽  
Vol 792 ◽  
pp. 5-35 ◽  
Author(s):  
Giuseppe A. Zampogna ◽  
Alessandro Bottaro

The interaction between a fluid flow and a transversely isotropic porous medium is described. A homogenized model is used to treat the flow field in the porous region, and different interface conditions, needed to match solutions at the boundary between the pure fluid and the porous regions, are evaluated. Two problems in different flow regimes (laminar and turbulent) are considered to validate the system, which includes inertia in the leading-order equations for the permeability tensor through a Oseen approximation. The components of the permeability, which characterize microscopically the porous medium and determine the flow field at the macroscopic scale, are reasonably well estimated by the theory, both in the laminar and the turbulent case. This is demonstrated by comparing the model’s results to both experimental measurements and direct numerical simulations of the Navier–Stokes equations which resolve the flow also through the pores of the medium.


1997 ◽  
Vol 119 (1) ◽  
pp. 122-128 ◽  
Author(s):  
S. L. Puterbaugh ◽  
W. W. Copenhaver

An experimental investigation concerning tip flow field unsteadiness was performed for a high-performance, state-of-the-art transonic compressor rotor. Casing-mounted high frequency response pressure transducers were used to indicate both the ensemble averaged and time varying flow structure present in the tip region of the rotor at four different operating points at design speed. The ensemble averaged information revealed the shock structure as it evolved from a dual shock system at open throttle to an attached shock at peak efficiency to a detached orientation at near stall. Steady three-dimensional Navier Stokes analysis reveals the dominant flow structures in the tip region in support of the ensemble averaged measurements. A tip leakage vortex is evident at all operating points as regions of low static pressure and appears in the same location as the vortex found in the numerical solution. An unsteadiness parameter was calculated to quantify the unsteadiness in the tip cascade plane. In general, regions of peak unsteadiness appear near shocks and in the area interpreted as the shock-tip leakage vortex interaction. Local peaks of unsteadiness appear in mid-passage downstream of the shock-vortex interaction. Flow field features not evident in the ensemble averaged data are examined via a Navier-Stokes solution obtained at the near stall operating point.


2009 ◽  
Vol 132 (1) ◽  
Author(s):  
Alessandro Armellini ◽  
Filippo Coletti ◽  
Tony Arts ◽  
Christophe Scholtes

The present contribution addresses the aerothermal, experimental, and computational studies of a trapezoidal cross-sectional model simulating a trailing edge cooling cavity with one rib-roughened wall. The flow is fed through tilted slots on one side wall and exits through straight slots on the opposite side wall. The flow field aerodynamics is investigated in Part I of the paper. The reference Reynolds number is defined at the entrance of the test section and set at 67,500 for all the experiments. A qualitative flow model is deduced from surface-streamline flow visualizations. Two-dimensional particle image velocimetry measurements are performed in several planes around midspan of the channel and recombined to visualize and quantify three-dimensional flow features. The crossing-jets issued from the tilted slots are characterized and the jet-rib interaction is analyzed. Attention is drawn to the motion of the flow deflected by the rib-roughened wall and impinging on the opposite smooth wall. The experimental results are compared with the numerical predictions obtained from the finite volume Reynolds-averaged Navier–Stokes solver, CEDRE.


2021 ◽  
Vol 11 (04) ◽  
pp. 1-11
Author(s):  
Wanwan Li

In mechanical engineering educations, simulating fluid thermodynamics is rather helpful for students to understand the fluid’s natural behaviors. However, rendering both high-quality and realtime simulations for fluid dynamics are rather challenging tasks due to their intensive computations. So, in order to speed up the simulations, we have taken advantage of GPU acceleration techniques to simulate interactive fluid thermodynamics in real-time. In this paper, we present an elegant, basic, but practical OpenGL/SL framework for fluid simulation with a heat map rendering. By solving Navier-Stokes equations coupled with the heat diffusion equation, we validate our framework through some real-case studies of the smoke-like fluid rendering such as their interactions with moving obstacles and their heat diffusion effects. As shown in Fig. 1, a group of experimental results demonstrates that our GPU-accelerated solver of Navier-Stokes equations with heat transfer could give the observers impressive real-time and realistic rendering results.


Author(s):  
Elisabeth Longatte

This work is concerned with the modelling of the interaction of a fluid with a rigid or a flexible elastic cylinder in the presence of axial or cross-flow. A partitioned procedure is involved to perform the computation of the fully-coupled fluid solid system. The fluid flow is governed by the incompressible Navier-Stokes equations and modeled by using a fractional step scheme combined with a co-located finite volume method for space discretisation. The motion of the fluid domain is accounted for by a moving mesh strategy through an Arbitrary Lagrangian-Eulerian (ALE) formulation. Solid dyncamics is modeled by a finite element method in the linear elasticity framework and a fixed point method is used for the fluid solid system computation. In the present work two examples are presented to show the method robustness and efficiency.


Sign in / Sign up

Export Citation Format

Share Document