The cover image is based on the Research Article Shape-preserving finite elements in cylindrical and spherical geometries: The double Jacobian approach by Jason P. Morgan, https://doi.org/10.1002/fld.4799 .
We have recently proposed a data-driven correction reduced-order model (DDC-ROM) framework for the numerical simulation of fluid flows, which can be formally written as follows.
The cover image is based on the Research Article Pressure boundary condition in a multi-phase lattice Boltzmann method and its applications on simulations of two-phase flows by Lei Wang and Ze-Rui Peng, https://doi.org/10.1002/fld.4800 .
Studies on the unphysical increase of turbulent quantities for RANS simulation induced by shock waves in hypersonic flows are carried out. Numerical experiments on the hypersonic flow over a blunt body reveal that the phenomenon of unphysical increase of turbulent quantities across the detached shock wave is induced by the strain-rate-based production terms of the k- and k- SST turbulence models, which leads to the over-prediction of aerothermal prediction. While this phenomenon does not occur for Spalart–Allmaras (S–A) turbulence model because of its vorticity-based production term. In order to eliminate this unphysical phenomenon, and to maintain the accuracy of the original models for boundary layer and separation flows, a new correction method for the k- and k- SST models is proposed: by comparing the orders of magnitude between the strain-rate-based and vorticity-based production terms, the vorticity-based production term is used near the shock waves, while the original strain-rate-based production term is still used in other regions. Finally, the correction method is applied to turbulence and transition flows over blunt bodies, and the numerical results show that the correction method effectively eliminates the unphysical increase of turbulent quantities across shock waves and improves the accuracy of aerothermal and transition onset location prediction. 相似文献
This paper presents a new spectral model for solving the fully nonlinear potential flow problem for water waves in a single horizontal dimension. At the heart of the numerical method is the solution to the Laplace equation which is solved using a variant of the ‐transform. The method discretizes the spatial part of the governing equations using the Galerkin method and the temporal part using the classical fourth‐order Runge‐Kutta method. A careful investigation of the numerical method's stability properties is carried out, and it is shown that the method is stable up to a certain threshold steepness when applied to nonlinear monochromatic waves in deep water. Above this threshold artificial damping may be employed to obtain stable solutions. The accuracy of the model is tested for: (i) highly nonlinear progressive wave trains, (ii) solitary wave reflection, and (iii) deep water wave focusing events. In all cases it is demonstrated that the model is capable of obtaining excellent results, essentially up to very near breaking. 相似文献
A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the ‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained. 相似文献
A 3D axisymmetric Galerkin boundary integral formulation for potential flow is employed to model two fluids of different densities, one fluid enclosed inside the other. The interface variables are the velocity potential and the normal velocity, and they can be solved for separately, the second linear system being symmetric. The algorithm is validated by comparing with the analytic solutions for a static interior spherical drop over a range of values for the relative densities of exterior and interior fluids and various boundary conditions. For time‐dependent simulations utilizing a level set method for the interface tracking, the accuracy has been checked by comparing against the known oscillation frequency of the sphere. Pinch‐off profiles corresponding to an initial two‐lobe geometry drop and D = 6 are also presented. Published in 2011 by John Wiley & Sons, Ltd. 相似文献
In this paper, the equivalent equations of DLR k−ε turbulent model in the boundary-fitted curvilinear coordinate are employed. Using the upwind idea that the contribution of the difference coefficients to the main node is positive contribution, and the other nodes are negative contribution or no contribution, new five-point difference schemes with a full diagonally dominant coefficient matrix (5-point-DD difference scheme) are constructed. Finally, taking the equation in the DLR k−ε turbulent model as an example, the mathematical characteristics of the 5-point-DD difference scheme are analyzed, and the uniform boundedness and convergence theorems of the Gauss-Seidel iterative sequence are given. Numerical simulations show that the five-point schemes are strictly diagonally dominant, and the calculated results are in good agreement with the experimental results. 相似文献