Compressible Fluid Research Articles

Analysis of natural vibration of truncated conical shells partially filled with fluid

AbstractIn this paper, we study the vibrational behavior of shells in the form of truncated cones containing an ideal compressible fluid. The sloshing effect on the free surface of the fluid is neglected. The dynamic behavior of the elastic structure is investigated based on the classical shell theory, the constitutive relations of which represent a system of ordinary differential equations written for new unknowns. Small fluid vibrations are described in terms of acoustic approximation using the wave equation for hydrodynamic pressure written in spherical coordinates. Its transformation into the system of ordinary differential equations is carried out by applying the generalized differential quadrature method. The formulated boundary value problem is solved by Godunov's orthogonal sweep method. Natural frequencies of shell vibrations are calculated using the stepwise procedure and the Muller method. The accuracy and reliability of the obtained results are estimated by making a comparison with the known numerical and analytical solutions. The dependencies of the lowest frequency on the fluid level and cone angle of shells under different combinations of boundary conditions (simply supported, rigidly clamped, and cantilevered shells) have been studied comprehensively. For conical straight and inverted shells, a numerical analysis has been performed to estimate the possibility of finding configurations at which the lowest natural frequencies exceed the corresponding values of the equivalent cylindrical shell.

A space marching method for sonic boom near field predictions

The conventional sonic boom propagation prediction method widely adopted in supersonic aircraft design involves a two-step procedure. In the first step, a compressible viscous or inviscid Computational Fluid Dynamics analysis is applied to the aircraft geometry at flight conditions to produce a flow field solution near the aircraft. Then in the second step, a one dimensional nonlinear Burgers’ equation model is used to propagate sonic boom signature traces from the near field solution at flight altitude to the ground along a ray path. For an accurate ground signature prediction the near field signature must be accurately modeled at a sufficient distance from the aircraft flight path in order to minimize errors in the one dimensional propagation model. This is a very challenging task for general purpose CFD tools in a design environment because the cost of maintaining highly accurate off-body solutions increases dramatically as the radial distance is enlarged in the computational domain. It is also particularly difficult to apply these tools for wave propagation because the algorithms are normally lower order and numerically dissipative and dispersive. In this work a space marching procedure based on an optimized higher order finite difference method is developed and applied in conjunction with a CFD solution concentrated in the close vicinity of the aircraft. This new approach is much more efficient, compared to previous methods, in providing highly accurate near field signatures for full carpet ground predictions. Results indicate that near field signatures retain more waveform shape information farther from the aircraft geometry while reducing the CFD cost significantly. The predicted ground signature is also shown to converge in shape as the radial distance of the near field signature grows, which is indicative of a more ideal initial condition being supplied to the one dimensional wave propagation to ground. This feature is very difficult to replicate in a tractable manner with common CFD approaches used in design.

MSD-based transient fluid model of liquid circular angular accelerometer for vibration-induced error analysis

This paper investigates the errors induced by linear vibration in the fluidic system of liquid circular angular accelerometers (LCAAs). A two-dimensional mass-spring-damper model is proposed to analyze the pressure difference between both sides of LCAA transducer subjected to linear vibration in both transient and steady states, where the compressibility of fluid is taken into account. The complex flow patterns of fluid are simplified to the positional changes of fluid particles. Springs and dampers are utilized to simulate the forces inside the fluid and between the fluid and the tube wall. Simulations are carried out based on the linear acceleration step excitation to explore the influence of model parameters, LCAA structure parameters and fluid parameters. The proposed model is validated by comparing the characteristics of the simulated pressure difference with their corresponding theoretical values, and is further verified by the comparison with the results from computational fluid dynamics method. Furthermore, the pressure difference under linear acceleration excitations in various forms is analyzed. The results indicate that the abrupt changes in linear acceleration are the primary cause of pressure difference oscillation, and the compressibility of fluid leads to increased pressure difference.

Newtonian gravitational waves from a continuum

Gravitational waves are being shown to derive directly from Newtonian dynamics for a continuous mass distribution, e.g. compressible fluids or equivalent. It is shown that the equations governing a continuous mass distribution, i.e. the inviscid Navier–Stokes equations for a general variable gravitational field g ( t , x ) , are equivalent to a form identical to the Maxwell equations from electromagnetism, subject to a specified condition. The consequence of this equivalence is the creation of gravity waves that propagate at a finite speed. The latter implies that Newtonian gravitation as presented in this paper is not ‘spooky action at a distance’ but rather is similar to electromagnetic waves propagating at finite speed, despite the apparent form appearing in the integrated field formula. In addition, this proves that, in analogy to the Maxwell equations, the Newtonian gravitation equations are Lorentz invariant for waves propagating at the speed of light. Since gravitational waves were so far derived only from Einstein’s general relativity theory, it becomes appealing to check if there is a connection between the Newtonian waves presented in this paper and the general relativity type of waves at least in a certain limit of overlapping validity, i.e. as a flat-space approximation. The latter is left for follow-up research.

Research on the prediction model of energy separation effect of vortex tube based on trajectory deflection characteristics and its numerical solution method

In this paper, in order to establish the energy separation mechanism of the vortex tube, the hydrodynamic behavior of the compressible fluid in the asymmetric cavity space is investigated, and a numerical model of the trajectory deflection behavior is deduced and established; in order to form the optimal design method of the structural parameters of the vortex tube, the force situation of the fluid microelements entering different regions of the vortex chamber of the vortex tube is analyzed, and the trajectory deflection equations are corrected by combining with the expansion behavior of the fluid and the characterizing equations of vortex strength, transportability, and vortex initiation characteristics are given. The characterization equations of vortex strength, transportability and vortex initiation characteristics are given, and the numerical simulation of their influence parameters is carried out; in order to realize the prediction of the vortex tube performance of a given structure, the multifactor Pearson thermodynamic map is used to correlate and analyze the experimental data of vortex tubes reported publicly in the past years, and the polynomial regression equations are designed and established for the prediction of the vortex tube's energy separation effect and the confidence level and the degree of coincidence of the prediction results are examined. The confidence level and degree of agreement of the prediction results were examined. It is found that: the trajectory deflection motion of the compressible fluid in the asymmetric cavity space is the result of the combined effect of structural air pressure bias and the expansion behavior of the incident fluid; in order to improve the vortex strength in the vortex tube, the vortex initiation chamber space should be as small as possible; the increase of the diameters of the hot-end pipe and the cold-end pipe is conducive to the enhancement of vortex strength, but at the same time, it weakens the vortex transport in the heat pipe; the vortex initiation chamber size has a negative correlation with the hot-end temperature rise, and the inlet fluid pressure has a The negative correlation between the size of the vortex chamber and the temperature rise at the hot end, the positive correlation between the increase of inlet fluid pressure and the resulting temperature rise, and the strong correlation between the inlet fluid pressure and the friction coefficient on the effect of energy separation; the predictive equations for the effect of energy separation obtained by the fitting are in good agreement with the real situation.

Hydrodynamic density functional theory of simple dissipative fluids

In this paper, a statistical physical derivation of thermodynamically consistent fluid mechanical equations is presented for non-isothermal viscous molecular fluids. The coarse-graining process is based on (i) the adiabatic expansion of the one-particle probability density function around local thermodynamic equilibrium, (ii) the assumption of decoupled particle positions and momenta, and (iii) the variational principle. It is shown that there exists a class of free energy functionals for which the conventional thermodynamic formalism can be naturally adopted for non-equilibrium scenarios, and describes entropy monotonic fluid flow in isolated systems. Furthermore, the analysis of the general continuum equations revealed the possibility of a non-local transport mode of energy in highly compressible dense fluids.

A new high-order shock-capturing TENO scheme combined with skew-symmetric-splitting method for compressible gas dynamics and turbulence simulation

The high-order shock-capturing scheme is one of the main building blocks for the simulation of the compressible fluid characterized by strong shockwaves and broadband length scales. However, the classical shock-capturing scheme fails to perform long-time stable and non-dissipative simulations since the quadratic invariants associated with the conservation equations cannot be conserved as a result of the inherent numerical dissipation. Additionally, the overall computational cost for classical shock-capturing schemes is quite expensive as a result of the time-consuming local characteristic decomposition and the nonlinear-weights computing process. In this work, based on a new efficient discontinuity indicator, which distinguishes the non-smooth high-wavenumber fluctuations and discontinuities from smooth scales in the wavenumber space, a paradigm of high-order shock-capturing scheme by recasting the non-dissipative skew-symmetric-splitting method with newly optimized dispersion property for smooth flow scales and invoking the nonlinear targeted ENO (TENO) schemes for non-smooth discontinuities is proposed. The resulting TENO-S scheme not only successfully performs long-time stable computations for smooth flows without numerical dissipation, but also recovers the robust shock-capturing capabilities with adaptive numerical dissipation. Without the necessity of parameter tuning case by case, extensive benchmark simulations involving a wide range of flow length scales and strong discontinuities demonstrate that the proposed TENO-S scheme performs significantly better than the straightforward deployment of WENO/TENO-family schemes with better spectral property and higher computational efficiency.

On the modeling of compressible viscous fluids via Burgers and Oldroyd equations

AbstractThe paper develops some generalizations of the Burgers and the Oldroyd equations for the dynamics of fluids. To account for compressibility, in addition to viscosity, first the Oldroyd derivative is replaced with the Truesdell derivative. Consequently, the two equations can be given a linear form within the Lagrangian formulation. Furthermore, possible anisotropies are modeled by replacing some (scalar) coefficients with tensors. To emphasize the compressibility property, generalizations are established to allow for nonzero longitudinal viscosity. Next, the thermodynamic consistency is investigated by regarding both types of equations as rate equations, of second order and first order. The requirements on the parameters entering the two equations are derived while the linearity of the two equations allow the free energy potential be quadratic. The Oldroyd equation is found to be compatible via appropriate restrictions of the tensor coefficients, through different pairs of free energy and entropy production.

A critical analysis and improvements of empirical models for predicting the performance of Electrical Submersible Pumps under viscous flow

One of the main challenges associated with Electrical Submersible Pumps (ESPs) system is the correct sizing of the pump. High viscosity fluids, emulsions, or compressible fluids significantly impact the pump’s ability to transfer energy to the fluid. Obtaining ESP characteristic curves in conditions resembling real applications is complex but essential for the correct design, operation, optimization, and safeguarding of ESP systems. The objective of this study is to critical assess the traditional methods for predicting the ESP performance under viscous flow. This study also improves these methods driven by ESP performance data. The results show that the original models have limitations because their databases are primarily based on data from conventional centrifugal pumps. Consequently, using these methods can result in significant errors in the sizing of the production system. The accuracy of the models improved significantly after modifying the methods. Reductions in errors were observed in all correction factors, and particularly in the prediction of global performance. The improvements were effective in mitigating the effects of increased deviations with the reduction in Reynolds number observed when applying the original models to predict flow rate and efficiency.

Numerical Investigations on the Jet Dynamics during Cavitation Bubble Collapsing between Dual Particles

The jet dynamics during cavitation bubble collapsing between unequal-sized dual particles are investigated utilizing a numerical model that combines the finite volume approach alongside the volume of fluid approach. The model incorporates the compressibility of the two-phase fluid and accounts for mass and heat transfer between two phases. The computational model utilizes an axisymmetric model, where the axis of symmetry is defined as the line that connects the centers of the particles and the bubble. A comprehensive analysis is presented on the influence of the particle radius and bubble–particle distance on the jet behavior. Furthermore, the variations of surface pressure on the particles induced by jet impingement are quantitatively analyzed. Four distinct jet behaviors are categorized, depending on the formation mechanism, as well as the number and the direction of the jets. For case 1, the bubble produces a single jet directed toward a small particle; for case 2, the bubble fragments produces double jets receding from each other; for case 3, the bubble produces double jets approaching each other; and for case 4, the bubble produces a single jet directed toward a large particle. The pressure perturbations induced by jet impingement upon the particles exceed those caused by shock wave impacts. The larger the bubble volume at the moment of jet formation, the longer the duration of the pressure variation caused by the jet impinging on the particles.

Dynamic Simulation of the Stability of Annular Flows in Gas Wells

This paper provides a holistic view of the development and progression of multiphase gas-liquid annular flow in gas wells. With the help of a momentum balance, the stability of a fully developed annular flow is defined. The Euler-Euler model was used to model the multiphase flow system, while the Reynolds-Averaged Navier–Stokes (RANS) model was used to depict a time-dependent solution for a vertical fluid flow system, considering acceleration due to gravity and fluid compressibility. The resultant flow regime was tested for stability by gradually reducing the field velocity of flow and observing the effects using a computational fluid dynamics (CFD) software. The principal parameters considered are the superficial gas velocity and liquid film flow rate. The properties of the region where the flow becomes unstable are documented and related to the incidence of liquid loading in a typical gas well. The effect are critically observed and compared to flow regimes at lower superficial gas velocity. The resultant flow regime configurations are analyzed and compared to the annular flow system. This study will help researchers understand the behaviour of gas wells that have potential in developing liquid loading problems and the movement of liquid film on the wall of the tubing when in annular flow regime.

Method for assessing well interference at under-gas cap zone using CRM material balance model

Optimization of injection operation conditions is a primary task when designing the development of mature oil fields. To select optimal injectivities, solutions to the optimization problem based on a CRM (capacitance resistance model) analytical model are used. CRM models based on analytical solutions of the material balance equations of weakly compressible fluids due to their speed can be used as an alternative to flow simulation models in solving a number of problems to support oil field development. The main task of CRM models is to determine the well interference factor, i.e. the share of fluid produced due to a particular injection well. These factors can be used to analyze waterflooding and develop solutions for waterflood optimization. At the same time, in some cases reservoir fluids may have high gas content, which is a classical limitation of a CRM model and leads to underestimation of the interference factors when estimating the injection performance. In this paper, the authors propose an approach to solve the problems of distortion of such factors, and for the first time in CRM application practice a model has been adapted for use in the conditions of high GORs and has a potential for being applied for under-gas cap zones. This is achieved by including in the equations the physical properties of gas and their behavior under reservoir conditions. The complication of algorithms did not significantly affect the model run time, but allowed to match the model both separately and jointly for liquid and gas phases. Thus, the improved classical CRM model has significantly expanded the scope of application of operational waterflood analysis tools for assessing the current situation. No less complicated and urgent is the problem of forecasting the production of gas liberated in the conditions of the previously mentioned under-gas-cap zones, because it implies the need to estimate the future reservoir pressure. The authors intend to devote the next paper to finding a solution to this issue, taking into account the proposed new method of reservoir pressure estimation using a CRM model.

Incompressible versus compressible fluid flow models: A case study on furnace tap-hole lancing

Pyrometallurgical furnaces, essential for metal extraction, operate at temperatures exceeding 1600°C and represent complex multiphase systems that challenge direct industrial research. Multiphysics models play a key role in shedding light on their intricate behaviours, supporting the refinement of design and operational strategies. Integral to the operation are the tap-holes, which facilitate the removal of molten products and are routinely opened by lancing, a process comparable to the use of a cutting torch, where high temperatures result from oxygen reacting with an iron lance. When the lance pierces the clay, oxygen gas enters the furnace, which could influence the behaviour of the molten material inside. In this work, a multiphase fluid flow model was used to investigate bulk flow dynamics, with a focus on the effects of the lancing process on the inside of the furnace, immediately behind the tap-hole. Incompressible and compressible multiphase fluid solvers were used and compared with respect their performance - the intention was to assess whether using a compressible solver would yield a different solution to the incompressible one. It was concluded that there are negligible disparities in bulk fluid flow behaviour between the solvers for the case studies examined, indicating that solver selection might be less consequential for certain aspects of oxygen lancing.

Low Mach number limit on perforated domains for the evolutionary Navier–Stokes–Fourier system

We consider the Navier–Stokes–Fourier system describing the motion of a compressible, viscous and heat-conducting fluid on a domain perforated by tiny holes. First, we identify a class of dissipative solutions to the Oberbeck–Boussinesq approximation as a low Mach number limit of the primitive system. Secondly, by proving the weak–strong uniqueness principle, we obtain strong convergence to the target system on the lifespan of the strong solution.

A tensor basis neural network-based turbulence model for transonic axial compressor flows

Traditional turbulence models encounter limitations when simulating intricate flows within transonic axial compressors. In contrast, recent advancements in machine learning turbulence models have demonstrated enhanced potential in refining the precision of turbulence modeling. Notably, the tensor basis neural network (TBNN) methodology has successfully developed non-linear eddy viscosity turbulence models. These models possess the capability to capture the anisotropic characteristic of Reynolds stress. However, applying machine learning non-linear eddy viscosity models to aircraft turbomachinery remains relatively infrequent. Current research mainly focuses on linear eddy viscosity turbulence models based on the Boussinesq hypothesis that presupposes isotropy in Reynolds stress. In this work, we introduce a non-linear eddy viscosity turbulence model, denoted as the k-ω-SST-TBNN model based on the TBNN framework. This model has been employed to simulate the transonic axial compressor NASA Rotor 37. The TBNN is trained using large eddy simulation (LES) results datasets. The importance of input scalar features is analyzed using the random forest method, from which significant and low-degree variables are selected as inputs for the TBNN. Furthermore, this study proposes including the turbulent Mach number as one of the extra features, representing fluid compressibility, thereby extending the computational mass flow rate range of compressor. Additionally, this paper proposes a method involving the weighted average of Reynolds stress, combining the high-precision but less robust TBNN predictions with results based on the Boussinesq assumption to enhance the turbulence model's robustness. The trained k-ω-SST-TBNN model undergoes validation on Rotor 37, where it exhibits a marked improvement in the prediction of overall performance, tip-gap vortex, and radial distribution of flow parameters. The model also displays a commendable capacity for generalization.

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly.

We introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

Compressible FSI of elastic spikes for drag reduction under hypersonic flow

In this study, a new fluid structure interaction (FSI) algorithm is developed for the aerodynamic-elasticity problem, aiming to address the coupling problem between compressible fluids and deformable solid structures. To verify the feasibility of the proposed model, a two-dimensional cantilever panel case under shock waves was conducted. Our numerical results present good consistency with the data from literatures. Furthermore, the validated FSI algorithm is applied on an elastic spike under compressible hypersonic flow. The complicated interaction between the fluid and solid structure under various spike geometry sizes and materials were revealed in respect of the spike deformation, pressure, temperature, and flow fields distribution. The main findings show that the spike deformation largely affects pressure and temperature on the blunt, the aerodynamic drag and lift force, indicating the necessity of FSI. For the spike diameters and lengths, it is found that larger spike diameters and shorter spikes benefits for deformation reduction, thereby minimizing fluctuations in aerodynamic drag force. For spike materials, spikes with larger elastic moduli and lower densities are preferred for the drag reduction. Thus, the proposed FSI model can accurately predict flow fields around and deformations for elastic spikes and offer guidance for drag reduction design in aerospace engineering.

On the motion of a body with a cavity filled with magnetohydrodynamic fluid

We study the dynamics of a coupled system, formed by a rigid body with a cavity entirely filled with magnetohydrodynamic compressible fluid. Our aim is to derive the global existence of the unique classical solutions and weak solutions to this system. Moreover, we show the weak-strong uniqueness principle which means that a weak solution coincides with a strong solution on the time existence of a strong solution, provided they emanate from the same initial data.

On the existence and stability of stationary solutions for the compressible micropolar fluids

We consider the existence, uniqueness and nonlinear stability of stationary solutions for the three-dimensional compressible micropolar fluids with the external force of general form. More precisely, with the weighted L2 and L∞ estimates on solutions to the linearized problem, we show the existence and uniqueness of stationary solutions in some suitable function space by the contraction mapping principle. The stability of the steady flow is studied by an elementary energy method.

A New Blowup Criterion to the Cauchy Problem for the Three-Dimensional Compressible Viscous Micropolar Fluids with Vacuum

A New Blowup Criterion to the Cauchy Problem for the Three-Dimensional Compressible Viscous Micropolar Fluids with Vacuum