scholarly journals ISWFoam: a numerical model for internal solitary wave simulation in continuously stratified fluids

2022 ◽  
Vol 15 (1) ◽  
pp. 105-127
Author(s):  
Jingyuan Li ◽  
Qinghe Zhang ◽  
Tongqing Chen

Abstract. A numerical model, ISWFoam, for simulating internal solitary waves (ISWs) in continuously stratified, incompressible, viscous fluids is developed based on a fully three-dimensional (3D) Navier–Stokes equation using the open-source code OpenFOAM®. This model combines the density transport equation with the Reynolds-averaged Navier–Stokes equation with the Coriolis force, and the model discrete equation adopts the finite-volume method. The k–ω SST turbulence model has also been modified according to the variable density field. ISWFoam provides two initial wave generation methods to generate an ISW in continuously stratified fluids, including solving the weakly nonlinear models of the extended Korteweg–de Vries (eKdV) equation and the fully nonlinear models of the Dubreil–Jacotin–Long (DJL) equation. Grid independence tests for ISWFoam are performed, and considering the accuracy and computing efficiency, the appropriate grid size of the ISW simulation is recommended to be 1/150th of the characteristic length and 1/25th of the ISW amplitude. Model verifications are conducted through comparisons between the simulated and experimental data for ISW propagation examples over a flat bottom section, including laboratory scale and actual ocean scale, a submerged triangular ridge, a Gaussian ridge, and slope. The laboratory test results, including the ISW profile, wave breaking location, ISW arrival time, and the spatial and temporal changes in the mixture region, are well reproduced by ISWFoam. The ISWFoam model with unstructured grids and local mesh refinement can effectively simulate the evolution of ISWs, the ISW breaking phenomenon, waveform inversion of ISWs, and the interaction between ISWs and complex topography.

2021 ◽  
Author(s):  
Jingyuan Li ◽  
Qinghe Zhang ◽  
Tongqing Chen

Abstract. A numerical model, ISWFoam, for simulating internal solitary waves (ISWs) in continuously stratified, incompressible, viscous fluids is developed based on a fully three-dimensional (3D) Navier-Stokes equation using the open source code OpenFOAM. This model combines the density transport equation with the Reynolds-averaged Navier-Stokes equation with the Coriolis force, and the model discrete equation adopts the finite volume method. The k-ω SST turbulence model has also been modified accordingly to the variable density field. ISWFoam provides two initial wave generation methods to generate an ISW in continuously stratified fluids, including solving the weakly nonlinear models of the extended Korteweg–de Vries (eKdV) equation and the fully nonlinear models of the Dubreil-Jacotin-Long (DJL) equation. Grid independence tests for ISWFoam are performed, considering the accuracy and computing efficiency, the appropriate grid size of the ISW simulation is recommended to be one-one hundred and fiftieth of the characteristic length and one-twenty fifth of the ISW amplitude. Model verifications are conducted through comparisons between the simulated and experimental data for ISW propagation examples over a flat bottom section, including laboratory scale and actual ocean scale, a submerged triangular ridge, a Gaussian ridge and slope. The laboratory test results, including the ISW profile, wave breaking location, ISW arrival time, and the spatial and temporal changes in the mixture region, are well reproduced by ISWFoam. The ISWFoam model with unstructured grids and local mesh refinement can accurately simulate the generation and evolution of ISWs, the ISW breaking phenomenon and the interaction between ISWs and complex structures and topography.


2021 ◽  
Vol 9 (12) ◽  
pp. 1374
Author(s):  
Jingyuan Li ◽  
Qinghe Zhang ◽  
Tongqing Chen

A numerical model of internal solitary waves in continuously stratified fluids is developed by introducing a density transport equation to the three-dimensional Navier–Stokes equation and adopting the fully nonlinear models of the Dubreil-Jacotin-Long equation to obtain the initial field of the ISW. The corresponding turbulence model has also been modified to ensure that it considers the variable density field. Comparisons between numerical simulation results and experimental results show that the total resistance, the nondimensional pressure coefficient, and the nondimensional friction coefficient for the standard submarine model proposed by the Defense Advanced Research Projects Agency under different flow field conditions are highly consistent with the experimental results. The model established is used to numerically analyse the forces and moments of the standard submarine model encountering ISWs at different submergence depths. The influence of the rotation centre position on the moment is discussed, and the position range of the optimal rotation centre is proposed.


1998 ◽  
Vol 115 (1) ◽  
pp. 18-24 ◽  
Author(s):  
G.W. Wei ◽  
D.S. Zhang ◽  
S.C. Althorpe ◽  
D.J. Kouri ◽  
D.K. Hoffman

Symmetry ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 288
Author(s):  
Alexei Kushner ◽  
Valentin Lychagin

The first analysis of media with internal structure were done by the Cosserat brothers. Birkhoff noted that the classical Navier–Stokes equation does not fully describe the motion of water. In this article, we propose an approach to the dynamics of media formed by chiral, planar and rigid molecules and propose some kind of Navier–Stokes equations for their description. Examples of such media are water, ozone, carbon dioxide and hydrogen cyanide.


1973 ◽  
Vol 59 (2) ◽  
pp. 391-396 ◽  
Author(s):  
N. C. Freeman ◽  
S. Kumar

It is shown that, for a spherically symmetric expansion of a gas into a low pressure, the shock wave with area change region discussed earlier (Freeman & Kumar 1972) can be further divided into two parts. For the Navier–Stokes equation, these are a region in which the asymptotic zero-pressure behaviour predicted by Ladyzhenskii is achieved followed further downstream by a transition to subsonic-type flow. The distance of this final region downstream is of order (pressure)−2/3 × (Reynolds number)−1/3.


Sign in / Sign up

Export Citation Format

Share Document