Pressure Wave Speed Research Articles

Seasonal nitrate variation as a tracer of preferential flow in bedrock aquifers

Short-term nitrate variations in streams are common, but the extent of groundwater contributions to such variability has been unclear. Analysis of well, spring and stream nitrate concentrations plus well water levels and stream flow for a chalk aquifer in Dorset (UK) shows that nitrate variation in streams lags discharge peaks by days to weeks. Spring monitoring also shows a short lag time, showing that high nitrate concentrations are rapidly transmitted from aquifers to streams. The lag time reflects the difference between the celerity (pressure wave speed) and the advective velocity of groundwater flow through the fracture network. Results give celerities of ∼500 m/d, advective velocities of 100–200 m/d, and kinematic porosities of ∼0.001. In addition, long-term nitrate data gives time lags of several decades, reflecting the matrix porosity of ∼0.3. Consequently, nitrate can be used an environmental tracer for both preferential flow through fracture networks and for retarded matrix flow and storage in this dual-porosity bedrock aquifer. Dual porosity is common in bedrock aquifers, so it is possible that this methodology may be widely applicable to understand the hydraulic characteristics of bedrock aquifers and to explain short-term nitrate variation in streams.

Comparative Analysis of Water Hammer Performance in Different Pipe Parameters with FSI

Water hammer (WH) is a critical phenomenon in fluid-filled piping systems that can lead to severe pressure surges and structural damage. The characteristics of the pipe material, geometry, and support conditions play a crucial role in the fluid–structure interaction (FSI) during WH events. This study investigates the impact of various pipe parameters, including material, length, thickness, and diameter, on the WH behavior using an FSI-based numerical approach. A comprehensive computational model was developed based on the algorithm presented in Delft Hydraulics Benchmark Problem (A) to simulate the WH phenomenon in pipes made of different materials, such as steel, copper, ductile iron, PPR (polypropylene random copolymer), and GRP (glass-reinforced plastic). This study examines the influence of pipe parameters on WH performance in pipelines, utilizing FSI to analyze the phenomenon. The results show that the pipe material has a significant influence on the pressure wave speed, stress wave propagation, and the overall system response during WH. Pipes with lower modulus of elasticity, such as PPR and GRP, exhibit lower pressure wave speeds but higher stress wave speeds compared with steel pipes. Increasing the elastic modulus, pipe wall thickness, length, and diameter enhances the pipe’s stiffness and impacts the timing, magnitude of pressure surges, and the likelihood of cavitation. The findings of this study provide valuable insights into the design and mitigation of WH in piping systems.

Experimental Verification of 1D-Simulation Method of Water Hammer Induced in Two Series-Connected Pipes of Different Diameters: Determination of the Pressure Wave Speed

Experimental Verification of 1D-Simulation Method of Water Hammer Induced in Two Series-Connected Pipes of Different Diameters: Determination of the Pressure Wave Speed

A finite volume model for maintaining stationarity and reducing spurious oscillations in simulations of sewer system filling and emptying

This paper introduces a model for simulating the unsteady dynamics of sewer systems filling and emptying, offering greater accuracy and stability. This article presents two novel contributions: first, the HLLS scheme (Harten–Lax–van Leer + Source term) is adapted to ensure the preservation of stationary conditions not only in free surface flows, as originally conceived, but also in pressurized flows and mixed flow scenarios. Second, a new method is proposed for the treatment of open channel flow cells near pressurization or adjacent to pressurized cells to minimize spurious oscillations when utilizing the two-component pressure approach (TPA) model. To verify the new model's effectiveness, it was tested for various conditions against the outcomes of the Open Source Field Operation and Manipulation (OpenFOAM) computational fluid dynamics (CFD) model. Furthermore, to demonstrate the model's potential for simulating real systems, the model was applied to three sewer systems that closely resemble real-world conditions, each of which had been intentionally modified for confidentiality purposes. The results show that the improved model successfully maintains stationary conditions within a sloped pipe across various flow conditions, while also preventing spurious oscillations at mixed flow interfaces even when using a pressure wave speed of 1000 m s − 1 .

Rotating detonation combustor performance informed through a novel megahertz-rate stagnation pressure measurement

An experimental stagnation pressure measurement technique is presented for a rotating detonation combustor (RDC). Schlieren imaging enables rotating detonation wave passage to be correlated with oscillations observed in the under-expanded exhaust plume. By measuring the spatiotemporal variation in exhaust plume divergence angle, stagnation pressure measurements of the RDC were acquired at a rate of 1 MHz. Combustor mass flux was varied between 202 and 783 kg/m2s, producing equivalent available pressures (EAPs) in the range of 3.42–13.5 bar. Time-averaged stagnation pressure measurements gathered using this technique were in agreement with the measured EAP within ±1.5%. Time-resolved stagnation pressure measurements allow for the pressure ratio produced across detonation wave cycles to be determined. For the conditions tested, detonation pressure ratios and wave speeds decreased while increasing the mean operating pressure of the combustor. Numerical modeling of the conditions tested indicates that the decrease in pressure ratio and wave speed is a result of elevated levels of combustion prior to the detonation wave arrival (i.e., “preburning”). Simultaneous OH* chemiluminescence measurements within the combustion chamber show an increase in preburned heat release relative to detonative heat release for increasing operating pressures of the RDC, in agreement with the results of the numerical model. Modeled chemical kinetic timescales decrease by approximately the same magnitude by which the preburning mass fraction increased in the range of operating pressures tested, suggesting that the faster reaction rates associated with higher pressure combustion may be the reason for increased preburning within the combustor.

Numerical study of the transient behaviors of vapour-liquid equilibrium state CO2 decompression

The vapor–liquid equilibrium (VLE) state of CO2 commonly appears in pipelines for its transport, refrigeration systems, and large-scale trans-critical cycle systems. However, the behavior of sudden leaks in such systems remains unclear, which poses a challenge in assessing the risk of leakage. This work is focused on the decompression behavior of VLE state CO2 with different volume fractions of vapor phases in high-pressure pipeline leakage, specifically with the development of a computational fluid dynamics (CFD) model with a nonequilibrium phase transition and the real gas model. The results suggest that the Peng-Robinson Equation of State (PR EoS) is more conservative in predicting the degree of superheat than that of Span-Wagner (S-W) EoS and GERG-2008 EoS. Moreover, it is found that the initial volume fraction of the vapor phase in VLE state CO2 plays a crucial role in determining the characteristics of sudden leakage, such as pressure, temperature, and decompression wave speed both inside and outside the pipe. However, the position of the Mach disk remains unaffected by the initial volume fraction of the vapor phase. The initial vapor volume fraction of 0.2 has the most significant impact on the transient behaviors of the leakage, that is, the speed of the initial decompression wave is approximately 1.46 times slower than that of pure liquid CO2, the pressure and temperature start to decrease approximately 0.46 ms after those of pure liquid CO2, which is about 5.1 times later. Additionally, the peak pressure of the near-field jet is 6.3% higher and the maximum velocity is 11.4% higher than that of pure liquid CO2. It hopes that this work will contribute to the improvement of research models that assess the consequences of potential high-pressure pipeline rupture scenarios.

New physics of supersonic ruptures

AbstractUntil recently, it is believed that the rupture speed above the pressure wave is impossible since spontaneously propagating ruptures are driven by the energy released due to the rupture motion, which is transferred through the medium to the rupture tip region at the maximum speed equal to the pressure wave speed. However, the apparent violation of classic theories has been revealed by new experimental results demonstrating supersonic shear ruptures. This paper presents a detailed analysis of the recently discovered shear rupture mechanism (fan hinged), which suggests a new physics of energy supply to the tip of supersonic ruptures. The key element of this mechanism is the fan‐shaped structure of the head of extreme ruptures, which is formed as a result of an intense tensile cracking process with the creation of intercrack slabs that act as hinges between the shearing rupture faces. The fan structure is featured with the following extraordinary properties: extremely low friction approaching zero; amplification of shear stresses above the material strength at low applied shear stresses; creation of a self‐disbalancing stress state causing a spontaneous rupture growth; abnormally high energy release; generation of driving energy directly at the rupture tip which excludes the need to transfer energy through the medium. The fan mechanism operates in intact rocks at stress conditions corresponding to seismogenic depths and in pre‐existing extremely smooth interfaces due to identical tensile cracking processes at these conditions. This is Paper 1 (of two companion papers) which discusses the fan theory and extreme ruptures in experiments on extremely smooth interfaces. Paper 2 entitled “Fan‐hinged shear instead of frictional stick‐slip as the main and most dangerous mechanism of natural, induced and volcanic earthquakes in the earth's crust” considers extreme ruptures in intact rocks. Further study of this subject is a major challenge for deep underground science, earthquake and fracture mechanics, physics, and tribology.

On Tsunami Waves Induced by Atmospheric Pressure Shock Waves After the 2022 Hunga Tonga‐Hunga Ha'apai Volcano Eruption

AbstractEmploying a linear shallow water equation (LSWE) model in the spherical coordinates, this paper investigates the tsunami waves generated by the atmospheric pressure shock waves due to the explosion of the submarine volcano Hunga Tonga–Hunga Ha'apai on 15 January 2022. Using the selected 59 atmospheric pressure records in the Pacific Ocean, an empirical atmospheric pressure model is first constructed. Applying the atmospheric pressure model and realistic bathymetric data in the LSWE model, tsunami generation and propagation are simulated in the Pacific Ocean. The numerical results show clearly the co‐existence of the leading locked waves, propagating with the speed of the atmospheric pressure waves (∼1,100 km/hr), and the trailing free waves, propagating with long gravity ocean wave celerity (∼750 km/hr). During the event, tsunamis were reported by 41 Deep‐ocean Assessment and Reporting of Tsunamis (DART) buoys in the Pacific Ocean, which require corrections because of the occurrence of atmospheric pressure waves. The numerically simulated tsunami arrival time and the amplitudes of the wave crest and trough of the leading locked waves compare reasonably well with the corrected DART measurements. The comparisons for the trailing waves are less satisfactory, since free waves could also have been generated by other tsunami generation mechanisms, which have not been considered in the present model, and by the scattering of locked waves over changing bathymetry. In this regard, the numerical results show clearly that the deep Tonga trench (∼10 km) amplifies the trailing waves in the Southeast part of the Pacific Ocean via the Proudman resonance condition.

Ultrasonic Evaluation of Laser Scanning Speed Effect on the Spectral Properties of Three-Dimensional-Printed Metal Phononic Crystal Artifacts.

Additive manufacturing/three-dimensional printing (AM/3DP) processes promise a flexible production modality to fabricate a complex build directly from its digital design file with minimal postprocessing. However, some critical shortcomings of AM/3DP processes related to the build quality and process repeatability are frequently experienced and reported in the literature. In this study, an in situ real-time nondestructive monitoring framework based on the dispersive properties of phononic crystal artifacts (PCAs) to address such quality challenges is described. Similar to a witness coupon, a PCA is printed alongside a build while it is interrogated and monitored with ultrasound. A PCA is substantially smaller than the actual build. Due to its periodic internal structures, a PCA creates pass and stop bands in its spectral response, which are sensitive to the variations in its process and material parameters. These periodic structures, representing the geometric complexities of an actual build, are designed for a specific monitoring objective(s) in AM/3DP. As a model application, in this demonstration study, the effect of the laser scanning speed of a slective laser melting (SLM) printer on the spectral properties of metal PCAs (mPCAs) is ultrasonically evaluated offline. The dependency of the pressure and shear wave speeds, the apparent Young's and shear moduli, and Poisson's ratio on the scanning speed are quantified, and it is found that they are highly sensitive to the laser scanning speed of an SLM printer. The sensitivity of the peaks of the pressure and shear spectral waveforms acquired for the identical mPCA designs printed on the same build plate with the same process parameters is also quantified. For powder-based AM/3DP technologies, where scanning speed is among the crucial process parameters such as laser power and bed temperature, the reported correlations between scanning speeds and the mechanical and spectral features of the mPCAs are expected to be instrumental in developing in situ real-time monitoring systems.

Investigation of pressure wave propagation and attenuation characteristics in managed pressure drilling by fast switching throttle valve

Pressure waves possess many significant applications in the oil and gas drilling engineering field, such as mud pulse telemetry (MPT) and measurement while drilling (MWD). The focus of this research is to study the pressure wave propagation and attenuation characteristics of wellbore liquid-phase flow in managed pressure drilling (MPD) by fast switching throttle valve (FSTV). First, a mathematical model of transient pressure wave propagation along the wellbore in both upstream and downstream directions is proposed in MPD by FSTV based on the one-dimensional transient flow theory. The model considering the frictional shear effect between the pipe walls is solved by utilizing the method of characteristics. Meanwhile, boundary conditions at the drill string inlet and annulus outlet, at the throttle valve, at the junction of drill bit, and at the reducer are adequately taken into account according to the actual situation of fluid flow. Second, a laboratory experiment of excited pressure waves in a vertical wellbore is conducted to measure the variation of the pressure fluctuation with different pump rates by FSTV. Comparing with the measurement result, the calculation result is discovered that the overall change is consistent, where the maximum absolute relative error at the peak of the pressure wave is only 4.5%. Finally, it further analyzes influential factors affecting the propagation and attenuation behaviors of wellbore pressure waves in liquid-phase flow based on the model. The results indicate that pressure waveforms present sinusoidal wave propagation, and pump rates, pressure wave speed, excitation time, fluid type, mud density, fluid viscosity, and borehole size exert varying degrees of influence on downhole pressure fluctuations. The proposed model achieves accurate quantitative interpretation and analysis of downhole pressure wave in MPD by FSTV, which has important significance for the realization of safe and efficient drilling.

Effect of shape on the physical properties of pharmaceutical tablets

Despite a well-established process understanding, quality issues for compressed oral solid dosage forms are frequently encountered during various drug product development and production stages. In the current work, a non-destructive contact ultrasonic experimental rig integrated with a collaborative robot arm and an advanced vision system is presented and employed to quantify the effect of the shape of a compressed tablet on its mechanical properties. It is observed that these properties are affected by the tablet geometric shapes and found to be linearly sensitive to the compaction pressures. It is demonstrated that the presented approach significantly improves the repeatability of the experimental waveform acquisition. In addition, with the increased confidence levels in waveform acquisition accuracy and corresponding pressure and shear wave speeds due to improved measurement repeatability, we conclude that pharmaceutical compact materials can indeed have a negative Poisson's ratio, therefore can be auxetic. The presented technique and instrument could find critical applications in continuous tablet manufacturing, and its real-time quality monitoring as measurement repeatability has been significantly improved, minimizing product quality variations.

Illegal connection detection in a viscoelastic pipeline using inverse transient analysis in the time domain

Abstract Illegal connection (IC) in water supply systems and water network wastes water as well as energy and reduces water quality, which has negative technical and economic effects on water management. Transient-based defect detection is a powerful method applied in water pipelines. This paper investigates the efficiency of the transient-based inverse transient analysis (ITA) method in estimating the characteristics of ICs in the viscoelastic water supply system in the time domain. To better evaluate this method, an experimental transient model was developed using polyethylene pipelines (with a length of 158 meters and a nominal diameter of 2 inches). In the first step, the hydraulic transient solver was calibrated in the calibration approach of dynamic parameters of pressure wave speed and pipe wall viscoelasticity. Then, the sensitivity of the ITA method to the spatial step of the method of characteristics and signal sample size was assessed. Finally, the efficiency of the ITA method was evaluated for several experiments with different transient intensities and a noisy signal. The results indicated that the accuracy for locating the IC was higher than the accuracy for the IC's length.

Experimental visualization and characteristics of bubble nucleation during rapid decompression of supercritical and subcooled carbon dioxide

The rapid decompression of sub-cooled and supercritical carbon dioxide was investigated experimentally using schlieren and regular high-speed photography. An optically accessible expansion tube was designed and constructed. A test campaign with initial pressures ranging from 6.45–9.02 MPa and initial temperatures from 297–308 K was performed. High-speed schlieren imagery was used to visualize the decompression waves, evaporation waves, and an oscillatory phase change region present in the blowdown process. The time history of pressures at different locations along the expansion tube were recorded using fast response pressure transducers. Image and pressure signal analyses were performed to extract decompression wave speeds and evaporation wave speeds. Results were compared to analytical models and other experiments reported in the literature. The relationship between pressure undershoot and wave speed to the initial conditions was analysed.

Comparison of three testing methods of rock materials equations of state

For rock materials, it is very important to use effective and reliable testing methods to obtain accurate Hugoniot parameters and further determine the equation of state (EOS) reflecting the intrinsic properties of materials. In this paper, taking marble as the experimental sample, three test methods of the Hugoniot EOS (i.e., symmetric, asymmetric, and reverse impacts) are proposed based on the first-class light gas gun experimental platform. The Displacement Interferometer System for Any Reflector is used to measure the speed history of the target plate and the arrival time of the shock wave. The electric probe is used to measure the speed of the flyer. Under the same launch conditions, at least three repeated experiments were conducted for each scheme. The experimental results showed that when the relative standard deviation (RSD) of the launch speed was controlled within 2%, the RSD of post-wave pressure and shock wave speed in the reverse impact method was the smallest, making it the best choice for the measurement of the Hugoniot EOS of rock materials. However, the symmetric impact can evaluate the impact flatness of the flying plate and the target plate by the time difference between the shock wave to reach different positions on the same plane of the target.

Study of Static and Dynamic Properties of Sand under Low Stress Compression

The aim of this work was to investigate a wide range of grain sizes of sand in the pre-yield regime during compression through the combined study of ultrasound (US) wave speed and acoustic emission (AE). The specific study was performed using modified oedometer and uniaxial compression experimental set-ups. The studied samples were natural dune sand (poorly graded on the poorly graded sand (SP) index) as well as its three extracted fractions as follows: 2.36–0.6 mm, 0.6–0.3 mm and 0.3–0.075 mm. The maximum compression stress during the modified oedometer experiments was <150 kPa, while during the modified uniaxial compression experiments, it was <400 kPa. Each sample was loaded while measuring the US pressure (P) wave speed and AE at each loading stage. The results show that the stiffer the soil is, the higher the value of the P wave speed measured, resulting in similar P wave velocity values achieved at a much lower applied stress during the oedometer experiments in comparison with the uniaxial compression tests. Regarding the AE results, it is seen that the higher the stress level is, causing more friction between the sand particles, the more AE events there are during their movement. The following parameters of AE were shown to be the most sensitive to the stress increase: the number of AE hits and the signals’ energy.

Natural shear wave propagation speed is influenced by both changes in myocardial structural properties as well as loading conditions

Abstract Funding Acknowledgements Type of funding sources: Public grant(s) – National budget only. Main funding source(s): Research Foundation - Flanders (FWO) Background Shear wave elastography (SWE) is a promising tool for the non-invasive assessment of myocardial stiffness. It is based on the evaluation of the propagation speed of shear waves by high frame rate echocardiography. These waves can be induced by for instance mitral valve closure (MVC) and the speed at which they travel is related to the instantaneous stiffness of the myocardium. Myocardial stiffness is defined by the local slope of the stress-strain relation and can therefore be altered by both changes in structural properties of the myocardium as well as loading conditions. Purpose The aim of this study was to investigate how changes in myocardial structural properties as well as loading conditions affect shear wave speed after MVC. Methods Until now, 8 pigs (weight: 33.6 ± 5.4 kg) were included. The following interventions were performed: 1) preload was reduced by balloon occlusion of the vena cava inferior, 2) afterload was increased by balloon occlusion of descending aorta, 3) preload was increased by intravenous administration of 500 ml of saline and 4) ischemia/reperfusion injury (I/R injury) was induced in the septal wall by balloon occlusion of the LAD for 90 min. with subsequent reperfusion for 40 min. Echocardiographic and left ventricular pressure recordings were simultaneously obtained during each intervention. Left ventricular parasternal long-axis views were acquired with an experimental high frame rate ultrasound scanner (average frame rate: 1279 ± 148 Hz). Shear waves were visualized on tissue acceleration maps by drawing an M-mode line along the interventricular septum. Shear wave propagation speed after MVC was calculated by assessing the slope of the wave pattern on the tissue acceleration map (Figure A). Results The change in left ventricular end-diastolic pressure (LVEDP) and shear wave speed after MVC between baseline and each intervention are shown in Figure B and C, respectively. Preload reduction resulted in significant lower LVEDP compared to baseline (p < 0.01), while the other loading changes did not have a significant effect. Shear wave speed after MVC significantly increased by afterload and preload increase (p < 0.01). I/R injury resulted in increased shear wave speed (p < 0.01) without significantly altering LVEDP. There was a good positive correlation between the change in LVEDP and the change in shear wave speed induced by loading changes (r = 0.76; p < 0.001) (Figure D). However, the correlation became less strong if data of I/R injury was taken into account as well (r = 0.63; p < 0.001). Conclusion Our results suggest that SWE is capable to characterize myocardial tissue properties and besides has the potential as a novel method for the estimation of left ventricular filling pressures. However, in the presence of structural changes of the myocardium, care should be taken when estimating filling pressures based on shear wave propagation speed. Abstract Figure.

A CFD computational model for high-pressure liquid CO2 decompression

A CFD decompression model with considering external flow field is established based on the non-equilibrium phase transition. The results of pressure and decompression wave speed calculated by CFD are compared with the results of “shock tube” test, showing good agreement. On this basis, the effects of external flow field on the transient behaviour of liquid CO2 decompression in the pipe is further investigated. It is found that the external flow field basically has no effect on the variation of pressure and decompression wave speed in the pipe. Moreover, a Mach disk is formed at about 7.35R (R is the radius of rupture opening, 19.05 mm) away from the rupture opening in the process of jet expansion. Although the external flow field has a significant influence on the static pressure at the rupture opening, the mass flow in the case without considering external flow field is larger than that with considering external flow field, and the overall difference between the two cases is within 6%.

Characteristics of Very Low Frequency Sound Propagation in Full Waveguides of Shallow Water.

This work is concerned with the characteristics of very low frequency sound propagation (VLF, ≤100 Hz) in the shallow marine environment. Under these conditions, the classical hypothesis of considering the sea bottom as a fluid environment is no longer appropriate, and the sound propagation characteristics at the sea bottom should be also considered. Hence, based on the finite element method (FEM), and setting the sea bottom as an elastic medium, a proposed model which unifies the sea water and sea bottom is established, and the propagation characteristics in full waveguides of shallow water can be synchronously discussed. Using this model, the effects of the sea bottom topography and the various geoacoustic parameters on VLF sound propagation and its corresponding mechanisms are investigated through numerical examples and acoustic theory. The simulation results demonstrate the adaptability of the proposed model to complex shallow water waveguides and the accuracy of the calculated acoustic field. For the sea bottom topography, the greater the inclination angle of an up-sloping sea bottom, the stronger the leak of acoustic energy to the sea bottom, and the more rapid the attenuation of the acoustic energy in sea water. The effect of a down-sloping sea bottom on acoustic energy is the opposite. Moreover, the greater the pressure wave (P-wave) speed in the sea bottom, the more acoustic energy remains in the water rather than leaking into the bottom; the influence laws of the density and the shear wave (S-wave) speed in the sea bottom are opposite.

Effect of hydrogen enrichment on the dynamics of a lean technically premixed elevated pressure flame

Abstract The heat release distribution, combustion instability characteristics and flow dynamics of lean swirl-stabilized flames of hydrogen enriched natural gas were studied using time-resolved (10 kHz) stereo PIV, OH* chemiluminescence imaging and acoustic pressure measurements. The technically premixed gas turbine model combustor was operated at elevated pressure up to 5 bar with preheated air. The H2 vol. fraction in the fuel was varied up to 50%. An M-shaped aerodynamically stabilized flame persisted at 1 bar up to low H2 enrichment ratios. Increasing H2 enrichment caused the M-flame to first transition to a bistable flame then to a shear layer stabilized V-flame. Increasing pressure and H2 enrichment both decreased the flame length. The combustion instability frequency and thermoacoustic amplitude varied and were mapped across operating conditions. The effect of H2 enrichment on the thermoacoustic amplitude was explained by examining the convective coupling between the Helmholtz acoustics and heat release. H2 enrichment was shown to increase the phase delay between the pressure and heat release by decreasing the pressure wave travel times in the plenum and swirler sections of the burner and by decreasing the flame length while simultaneously increasing the pressure wave speed within the burner. Thus H2 enrichment increased the thermoacoustic amplitude when the phase delay at the initial condition was negative, and decreased it when the initial phase delay was positive. Flow dynamics were studied by characterizing the precessing vortex core frequency and energy. The PVC structure existed across all conditions, and its frequency reduced to a different Strouhal number for the M-flame and V-flame. H2 enrichment consistently slightly decreased the PVC energy.

Study on pressure wave propagation through cryogenic condensing two-phase flow in liquid rocket propellant feedline

This paper presents a numerical study on the pressure wave traveling in the condensing flow of the liquid and vapor oxygen in the tank-to-engine propellant feedline of liquid rocket. A computational fluid dynamics (CFD) model is developed based on Eulerian multiphase model, and the coupling of pressure wave and flow condensation is achieved through incorporating user-defined functions (UDF) with interphase mass, heat and momentum exchange models. Experimental validation is undertaken and the predicted pressure profiles by CFD simulation prove to be in good agreement with experimental data. The simulation results indicate that the distance required to condense the oxygen mixture to pure liquid fluctuates with propagation of the pressure wave which experiences attenuation and deceleration with a speed of the order of 1–14 m/s. As the local vapor content reduces, the pressure wave speed increases in an exponential way with an abrupt increase greater than an order of magnitude taking place on transition to the single-phase liquid oxygen, while the amplitude attenuation decreases linearly. The dependence of pressure wave speed and attenuation on frequency satisfy exponential and logarithmic relations, respectively. The wave speed is to a large extent independent of the high frequencies and exhibits a dispersed characteristic at low frequencies.