scholarly journals Improved theory for shock waves using the OBurnett equations

2021 ◽  
Vol 929 ◽  
Author(s):  
Ravi Sudam Jadhav ◽  
Abhimanyu Gavasane ◽  
Amit Agrawal

The main goal of the present study is to thoroughly test the recently derived OBurnett equations for the normal shock wave flow problem for a wide range of Mach number ( $3 \leq Ma \leq 9$ ). A dilute gas system composed of hard-sphere molecules is considered and the numerical results of the OBurnett equations are validated against in-house results from the direct simulation Monte Carlo method. The primary focus is to study the orbital structures in the phase space (velocity–temperature plane) and the variation of hydrodynamic fields across the shock. From the orbital structures, we observe that the heteroclinic trajectory exists for the OBurnett equations for all the Mach numbers considered, unlike the conventional Burnett equations. The thermodynamic consistency of the equations is also established by showing positive entropy generation across the shock. Further, the equations give smooth shock structures at all Mach numbers and significantly improve upon the results of the Navier–Stokes equations. With no tweaking of the equations in any way, the present work makes two important contributions by putting forward an improved theory of shock waves and establishing the validity of the OBurnett equations for solving complex flow problems.

2018 ◽  
Vol 856 ◽  
Author(s):  
M. Borgnino ◽  
G. Boffetta ◽  
F. De Lillo ◽  
M. Cencini

We study the dynamics and the statistics of dilute suspensions of gyrotactic swimmers, a model for many aquatic motile microorganisms. By means of extensive numerical simulations of the Navier–Stokes equations at different Reynolds numbers, we investigate preferential sampling and small-scale clustering as a function of the swimming (stability and speed) and shape parameters, considering in particular the limits of spherical and rod-like particles. While spherical swimmers preferentially sample local downwelling flow, for elongated swimmers we observe a transition from downwelling to upwelling regions at sufficiently high swimming speed. The spatial distribution of both spherical and elongated swimmers is found to be fractal at small scales in a wide range of swimming parameters. The direct comparison between the different shapes shows that spherical swimmers are more clusterized at small stability and speed numbers, while for large values of the parameters elongated cells concentrate more. The relevance of our results for phytoplankton swimming in the ocean is briefly discussed.


2020 ◽  
Vol 8 (6) ◽  
pp. 419 ◽  
Author(s):  
Yun-Ta Wu ◽  
Shih-Chun Hsiao

In this article, the interaction of solitary waves and a submerged slotted barrier is investigated in which the slotted barrier consists of three impermeable elements and its porosity can be determined by the distance between the two neighboring elements. A new experiment is conducted to measure free surface elevation, velocity, and turbulent kinetic energy. Numerical simulation is performed using a two-dimensional model based on the Reynolds-Averaged Navier-Stokes equations and the non-linear k-ɛ turbulence model. A detailed flow pattern is illustrated by a flow visualization technique. A laboratory observation indicates that flow separations occur at each element of the slotted barrier and the vortex shedding process is then triggered due to the complicated interaction of those induced vortices that further create a complex flow pattern. During the vortex shedding process, seeding particles that are initially accumulated near the seafloor are suspended by an upward jet formed by vortices interacting. Model-data comparisons are carried out to examine the accuracy of the model. Overall model-data comparisons are in satisfactory agreement, but modeled results sometimes fail to predict the positions of the induced vortices. Since the measured data is unique in terms of velocity and turbulence, the dataset can be used for further improvement of numerical modeling.


1983 ◽  
Vol 50 (2) ◽  
pp. 265-269
Author(s):  
D. Nixon

The perturbation theory for transonic flow is further developed for solutions of the Navier-Stokes equations in two dimensions or for experimental results. The strained coordinate technique is used to treat changes in location of any shock waves or large gradients.


1969 ◽  
Vol 37 (4) ◽  
pp. 727-750 ◽  
Author(s):  
Gareth P. Williams

A method of numerically integrating the Navier-Stokes equations for certain three-dimensional incompressible flows is described. The technique is presented through application to the particular problem of describing thermal convection in a rotating annulus. The equations, in cylindrical polar co-ordinate form, are integrated with respect to time by a marching process, together with the solving of a Poisson equation for the pressure. A suitable form of the finite difference equations gives a computationally-stable long-term integration with reasonably faithful representation of the spatial and temporal characteristics of the flow.Trigonometric interpolation techniques provide accurate (discretely exact) solutions to the Poisson equation. By using an auxiliary algorithm for rapid evaluation of trigonometric transforms, the proportion of computation needed to solve the Poisson equation can be reduced to less than 25% of the total time needed to’ advance one time step. Computing on a UNIVAC 1108 machine, the flow can be advanced one time-step in 2 sec for a 14 × 14 × 14 grid upward to 96 sec for a 60 × 34 × 34 grid.As an example of the method, some features of a solution for steady wave flow in annulus convection are presented. The resemblance of this flow to the classical Eady wave is noted.


2018 ◽  
Vol 855 ◽  
pp. 43-66 ◽  
Author(s):  
Ke Wu ◽  
Bruno D. Welfert ◽  
Juan M. Lopez

The dynamic response to shear of a fluid-filled square cavity with stable temperature stratification is investigated numerically. The shear is imposed by the constant translation of the top lid, and is quantified by the associated Reynolds number. The stratification, quantified by a Richardson number, is imposed by maintaining the temperature of the top lid at a higher constant temperature than that of the bottom, and the side walls are insulating. The Navier–Stokes equations under the Boussinesq approximation are solved, using a pseudospectral approximation, over a wide range of Reynolds and Richardson numbers. Particular attention is paid to the dynamical mechanisms associated with the onset of instability of steady state solutions, and to the complex and rich dynamics occurring beyond.


1998 ◽  
Vol 120 (1) ◽  
pp. 72-75 ◽  
Author(s):  
V. N. Kurdyumov ◽  
E. Ferna´ndez

A correlation formula, Nu = W0(Re)Pr1/3 + W1(Re), that is valid in a wide range of Reynolds and Prandtl numbers has been developed based on the asymptotic expansion for Pr → ∞ for the forced heat convection from a circular cylinder. For large Prandtl numbers, the boundary layer theory for the energy equation is applied and compared with the numerical solutions of the full Navier Stokes equations for the flow field and energy equation. It is shown that the two-terms asymptotic approximation can be used to calculate the Nusselt number even for Prandtl numbers of order unity to a high degree of accuracy. The formulas for coefficients W0 and W1, are provided.


2008 ◽  
Vol 136 (11) ◽  
pp. 4045-4062 ◽  
Author(s):  
Robert L. Walko ◽  
Roni Avissar

Abstract The dynamic core of the Ocean–Land–Atmosphere Model (OLAM), which is a new global model that is partly based on the Regional Atmospheric Modeling System (RAMS), is described and tested. OLAM adopts many features of its predecessor, but its dynamic core is new and incorporates a global geodesic grid with triangular mesh cells and a finite-volume discretization of the nonhydrostatic compressible Navier–Stokes equations. The spatial discretization of horizontal momentum is based on a C-staggered grid and uses a method that has not been previously applied in atmospheric modeling. The temporal discretization uses a unique form of time splitting that enforces consistency of advecting mass flux among all conservation equations. OLAM grid levels are horizontal, and topography is represented by the shaved-cell method. Aspects of the shaved-cell method that pertain to the OLAM discretization on the triangular mesh are described, and a method of conserving momentum in shaved cells on a C-staggered grid is presented. The dynamic core was tested in simulations with multiple vertical model levels and significant vertical motion. The tests include an idealized global circulation simulation, a cold density current, and mountain-wave flow over an orographic barrier, all of which are well-known standard benchmark experiments. OLAM gave acceptable results in all tests, demonstrating that its dynamic core produces accurate and robust solutions.


Sign in / Sign up

Export Citation Format

Share Document