Conservation equations of gasdynamics in curvilinear coordinate systems☆
Abstract
This paper presents a new method of writing the conservation equations of gasdynamics in curvilinear coordinates which eliminates undifferentiated terms. One canthus readily apply difference schemes derived for Cartesian coordinates which conserve mass, momentum, and energy in the total flow field. The method is derived for orthogonal coordinates, and then extended to cover the most general class of coordinate transformations, using general tensor analysis. Several special features of the equations are discussed.
References (11)
- J.L. Anderson et al.
J. Computational Phys.
(1968) - G. Moretti et al.
AIAA Paper No. 72-192
(June 1972) - P.D. Thomas et al.
AIAA J.
(1972) - P. Kutler et al.
J. Spacecraft Rockets
(1971) - R.D. Richtmyer et al.
Difference Methods for Initial-Value Problems
Cited by (355)
Artificial viscosity-based shock capturing scheme for the Spectral Difference method on simplicial elements
2024, Journal of Computational PhysicsAn artificial viscosity-based shock capturing scheme is extended to the context of high-order triangular Spectral Difference method for solving gas dynamics problems featuring discontinuities and shock waves. The equations are regularised thanks to an artificial diffusivity method combined with a shock sensor based on dilatation and vorticity fields valid for any polynomial order of approximation. The methodological novelty of this paper is an projection based iterative Restriction Prolongation filtering operation introduced for high-order triangular elements. This filtering operation applied to the artificial bulk viscosity and thermal conductivity stabilises the simulations and improves the quality of the shock capturing approach. The methodology is validated on various canonical relevant compressible test cases. Numerical results demonstrate the performance and the flexibility of the triangular high-order Spectral Difference method for solving simultaneously compressible flows and shock waves.
Strongly conservative discretization of governing equations in cylindrical coordinates
2024, Computers and FluidsCylindrical coordinates are often used in computational fluid dynamics (CFD), particularly when one considers axisymmetric flow. Although cylindrical coordinates have several advantages in accurately describing the rotational flow as well as the complex boundary, there are limits to how far they can be taken. These are (1) the governing equations are mostly in a non-conservative form in current standard models, which cannot preserve the conservation of numerical simulations, (2) the singularity along the axis at the coordinate origin (-axis), leading to difficulty reproducing the flow across the -axis, and (3) the time step is extremely shortened to satisfy the CFL condition near the -axis because the numerical cell thereof is narrow in the azimuthal direction for a given angular resolution. A strongly conservative model is proposed in this paper to overcome these difficulties. First, the strongly conservative form of the governing equations is derived by considering changes in the direction of the base vector within a cell, and a concept of direction error is defined to reveal the inherent error of the standard model. Second, inspired by the derivation of the first, a strongly conservative method for evaluating the flux across each cell surface is proposed, preserving the conservation of discretization and enlarging the time step. Finally, the proposed model is tested by three classical problems, i.e., uniform flow, Sod shock tube problem, and expanding outflows. The results imply that the proposed model is strongly conservative for the cylindrical coordinate system, and it can overcome the limits of the standard model and eliminate the inherent error in the standard model. Furthermore, the strongly conservative model and method are extended to arbitrary coordinates for the numerical conservation in general curvilinear coordinates of the problem with complex boundaries.
Advanced modelling of lubricated interfaces in general curvilinear grids
2023, Tribology InternationalTackling fluid-flow problems involving intricate surface geometries has been the catalyst for a plethora of numerical investigations aimed at accommodating curved complex boundaries. An example is the application of body-fitted curvilinear coordinate transformation, where the one-to-one correspondence of grid points from a non-orthogonal curvilinear grid in the physical domain to an orthogonal grid in the computational domain is achieved. In lubricated interfaces, such conversion is challenging due to the complexity of the governing equations in the mapped-grid, the numerical instabilities exhibited by their non-linearities and the severity of the operating conditions. The present contribution proposes a Reynolds-based, finite volume fluid–structure interaction (FSI) framework for solving thermal elastohydrodynamic lubrication (TEHL) problems mapped onto orthogonal grids in the computational domain. We demonstrate how the strong conservation form of the pertinent governing equations can be expressed in three-dimensional curvilinear grids and discretised using the finite volume method to ensure fluid-flow conservation and enforce mass-conserving cavitation conditions. Numerical and experimental benchmarks showcase the robustness and versatility of the proposed framework to simulate a diverse range of lubrication problems, hence achieving a predictive computational tool that would enable a shift towards tribology-aware design.
Research on printed circuit board external open-circuit failure gas discharge tube under short-circuit failure
2023, Electric Power Systems ResearchThe short-circuit failure of gas discharge tubes poses a potential danger to electrical equipment. It is also one of the causes of fire and combustion of electrical equipment. In this study, the above-mentioned short-circuit failure was investigated, and a printed circuit board external open-circuit failure gas discharge tube was developed. The structure and working mechanism of the printed circuit board external open-circuit failure gas discharge tube were primarily presented. The steady-state numerical model of the fault arc of the discharge tube was built according to magnetohydrodynamics. The arc temperature distribution and the relationship between the short-circuit current and the arc temperature were explored. The impact and power frequency follow current coupling experiments of the open-circuit failure gas discharge tube were performed using 8/20 µs impulse current generator coupled with an alternating current power supply. The open-circuit failure response time of the external open-circuit failure gas discharge tube was tested, and the effects of the short-circuit follow current and the external open-circuit failure location on the open-circuit failure time of the open-circuit failure gas discharge tube were examined through repetitive experiments. The high-speed camera was adopted to shoot the open-circuit failure action picture of the open-circuit failure gas discharge tube. The experimental results indicated that the size of the short-circuit follow current was positively correlated with the open-circuit response time of the open-circuit failure gas discharge tube, which was significantly affected by the power frequency follow current and the location of the trip point. The experimental results are consistent with the action picture and simulation results of the open-circuit failure gas discharge tube, such that the short-circuit follow current failure can be effectively suppressed. The conclusion of this study provides a certain reference for the application of the open-circuit failure gas discharge tube in real lightning.
Immersed boundary method for high-order flux reconstruction based on volume penalization
2022, Journal of Computational PhysicsIn the last decade, there has been a lot of interest in developing high-order methods as a viable option for unsteady scale-resolving-simulations which are increasingly important in the industrial design process. High-order methods offer the advantages of low numerical dissipation, high efficiency on modern architectures and quasi mesh-independence. Despite significant advance in high-order solution methods, the general CFD workflow (geometry, CAD preparation, meshing, solution, post-processing) has largely remained unchanged, with mesh generation being a significant bottleneck and often determining the overall quality of the solution. In this work, we aim to combine the numerical advantages of the high-order Flux Reconstruction (FR) method and the simplicity of the mesh generation (or lack thereof) of the Immersed Boundary Method (IBM) for steady and unsteady problems over moving geometries. The volume penalization (penalty-IBM) method is selected for its ease of implementation and robustness. Detailed discussions about numerical implementation, including the boundary representation, mask function, data reconstruction, and selection of the penalization parameter are given. Advantages of combining volume penalization in the high-order framework are shown by various numerical test cases. The approach is firstly demonstrated for the linear advection-diffusion equation by investigating the numerical convergence for the coupled FR-IBM approach. Thereafter, the accuracy of the approach is demonstrated for canonical (static) test cases in 2D and 3D when compared to a standard body-fitted unstructured simulation. Finally, the efficiency of the method to handle moving geometries is demonstrated for the flow around an airfoil with pitching and plunging motions.
Effects of blowing upon dynamic stability of blunt nosed re entry vehicles pitching in hypersonic flow
2023, Aeronautical Journal
- ☆
This research was supported by NASA-Ames Research Center under the Consortium Agreement, Collaboration Research Interchange No. SCC-301.