Sea-based rocket launches encounter significant challenges stemming from dynamic marine environmental interactions.During the hot launch phase,characterized by low-velocity ascent,the departure of the rocket from the ...Sea-based rocket launches encounter significant challenges stemming from dynamic marine environmental interactions.During the hot launch phase,characterized by low-velocity ascent,the departure of the rocket from the oscillatory platform exhibits heightened sensitivity to external disturbances.In the development stage,assessing the launch dynamics and the clearance between the rocket and framed launcher are crucial for improving the reliability of sea-based rocket launches in rough sea conditions.This study presents a high-fidelity dynamic model of maritime hot launch system,demonstrating 3.21%prediction error through rigorous validation against experimental datasets from comprehensive modal analyses and the full-scale rocket flight test.To mitigate collision risks,we develop a computational method employing spatial vector analysis for dynamic measurement of rocket-launcher clearance during departure.Systematic investigations reveal that in rough sea conditions,optimal departure dynamics are achieved at θ_(thrust)=270°nozzle azimuth configuration,reducing failure probability compared to conventional orientations.The developed assessment framework not only resolves critical safety challenges in current sea launch systems but also establishes foundational principles for optimizing adapter axial configuration patterns in future designs.展开更多
The rocket sled system is not only a high-speed dynamic ground test system,but also one of the future aerospace horizontal launch schemes.The winged load,as a common type of payload,has greater vibration and noise int...The rocket sled system is not only a high-speed dynamic ground test system,but also one of the future aerospace horizontal launch schemes.The winged load,as a common type of payload,has greater vibration and noise intensity than the wingless load.Due to the severe aerodynamic instability prior to separation,the head-up or head-down phenomena are more evident and the test accuracy significantly decreases.The high-precision computer fluid dynamics and aeroacoustic analysis are employed to explore the multifield coupling mechanism of a rocket sled with the winged payload in the wide speed range(Ma=0.5–2).The results show that as the incoming velocity increases,the cone angle of the shock wave of the rocket sled decreases,the shock pressure increases quickly,and the vortex between the slippers splits and gradually shrinks in size.The velocity of the rocket sled exerts little influence on the modal resonance frequency.The wing has a significant impact on aerodynamic noise,and as the sound pressure level rises,the propagation direction gradually shifts towards the rear and upper regions of the wing.展开更多
In liquid rocket engines,regenerative cooling technology is essential for preserving structural integrity under extreme thermal loads.However,non-uniform coolant flow distribution within the cooling channels often lea...In liquid rocket engines,regenerative cooling technology is essential for preserving structural integrity under extreme thermal loads.However,non-uniform coolant flow distribution within the cooling channels often leads to localized overheating,posing serious risks to engine reliability and operational lifespan.This study employs a three-dimensional fluid–thermal coupled numerical model to systematically investigate the influence of geometric parameters-specifically the number of inlets,the number of channels,and inlet manifold configurations-on flow uniformity and thermal distribution in non-pyrolysis zones.Key findings reveal that increasing the number of inlets from one to three significantly enhances flow uniformity,reducing mass flow rate deviation from 1.2%to below 0.3%.However,further increasing the inlets to five yields only marginal improvements indicating diminishing(<0.1%),returns beyond three inlets.Additionally,temperature non-uniformity at the combustion chamber throat decreases by 37%-from 3050 K with 18 channels to 1915 K with 30 channels-highlighting the critical role of channel density in effective thermal regulation.Notably,while higher channel counts improve cooling efficiency,they also result in increased pressure losses of approximately 18%–22%,emphasizing the need to balance thermal performance against hydraulic resistance.An optimal configuration comprising 24 channels and three inlets was identified,providing minimal temperature gradients while maintaining acceptable pressure losses.The inlet manifold structure also plays a pivotal role in determining flow distribution.Configuration 3(Config-3),which features an enlarged manifold and reduced inlet velocity,achieves a 40%reduction in velocity fluctuations compared to Configuration 1(Config-1).This improvement leads to a more uniform mass flow distribution,with a relative standard deviation(RSD)of less than 0.15%.Furthermore,this design effectively mitigates localized hot spots near the nozzle-where temperature gradients are most severe-achieving a reduction of approximately 1135 K.展开更多
The reuse of liquid propellant rocket engines has increased the difficulty of their control and estimation.State and parameter Moving Horizon Estimation(MHE)is an optimization-based strategy that provides the necessar...The reuse of liquid propellant rocket engines has increased the difficulty of their control and estimation.State and parameter Moving Horizon Estimation(MHE)is an optimization-based strategy that provides the necessary information for model predictive control.Despite the many advantages of MHE,long computation time has limited its applications for system-level models of liquid propellant rocket engines.To address this issue,we propose an asynchronous MHE method called advanced-multi-step MHE with Noise Covariance Estimation(amsMHE-NCE).This method computes the MHE problem asynchronously to obtain the states and parameters and can be applied to multi-threaded computations.In the background,the state and covariance estimation optimization problems are computed using multiple sampling times.In real-time,sensitivity is used to quickly approximate state and parameter estimates.A covariance estimation method is developed using sensitivity to avoid redundant MHE problem calculations in case of sensor degradation during engine reuse.The amsMHE-NCE is validated through three cases based on the space shuttle main engine system-level model,and we demonstrate that it can provide more accurate real-time estimates of states and parameters compared to other commonly used estimation methods.展开更多
Taking a C1x motor with a backward-facing step which can generate a typical corner vortex as a reference,a numerical methodology using large eddy simulation was established in this study.Based on this methodology,the ...Taking a C1x motor with a backward-facing step which can generate a typical corner vortex as a reference,a numerical methodology using large eddy simulation was established in this study.Based on this methodology,the position of the backward-facing step of the motor was computed and analyzed to determine a basic configuration.Two key geometrical parameters,the head cavity angle and submerged nozzle cavity height,were subsequently introduced.Their effects on the corner vortex motion and their interactions with the acoustic pressure downstream of the backward-facing step were analyzed.The phenomena of vortex acoustic coupling and characteristics of pressure oscillations were further explored.The results show that the maximum error between the simulations and experimental data on the dominant frequency of pressure oscillations is 5.23%,which indicates that the numerical methodology built in this study is highly accurate.When the step is located at less than 5/8 of the total length of the combustion chamber,vortex acoustic coupling occurs,which can increase the pressure oscillations in the motor.Both the vorticity and the scale of vortices in the downstream step increase when the head cavity angle is greater than 24°,which increases the amplitude of the pressure oscillation by maximum 63.0%.The submerged nozzle cavity mainly affects the vortices in the cavity itself rather than those in the downstream step.When the height of the cavity increases from 10 to 20 mm,the pressure oscillation amplitude under the main frequency increases by 39.1%.As this height continues to increase,the amplitude of pressure oscillations increases but the primary frequency decreases.展开更多
A two-dimensional large eddy simulation numerical model is proposed to study the transient vortex flow and pressure oscillation of a large-aspect-ratio solid rocket motor.The numerical model is validated through exper...A two-dimensional large eddy simulation numerical model is proposed to study the transient vortex flow and pressure oscillation of a large-aspect-ratio solid rocket motor.The numerical model is validated through experimental data,finite element analysis and cumulative error analysis.The numerical simulations are executed to obtain the characteristics of the vortex-acoustic and pressure oscillation.The results show that the burning surface regression decreases the motor aspect ratio,increasing the corresponding natural frequency from 260 Hz to 293 Hz.The pressure oscillation phenomenon is formed due to the vortex-acoustic coupling.Decreasing the corner vortex shedding intensity shows negative effects on the dimensionless amplitude of the pressure oscillation.The head cavity without the injection can decrease the vortex-acoustic coupling level at the acoustic pressure antinode.The modified motor with head cavity can obtain a lower dimensionless oscillating pressure amplitude 0.00149 in comparison with 0.00895 of the original motor.The aspect ratio and volume of the head cavity without the injection have great effects on the pressure oscillation suppression,particularly at the low aspect ratio or large volume.The reason is that the mass in the region around the acoustic pressure antinode is extracted centrally,reducing the energy contribution to the acoustic system.With the volume increasing,the acoustic energy capacity increases.展开更多
As a critical component of pulse solid rocket motors(SRMs),the soft pulse separation device(PSD)is vital in enabling multi-pulse propulsion and has become a breakthrough in SRM engineering applications.To investigate ...As a critical component of pulse solid rocket motors(SRMs),the soft pulse separation device(PSD)is vital in enabling multi-pulse propulsion and has become a breakthrough in SRM engineering applications.To investigate the opening performance of the PSD,an axial PSD incorporating a star-shaped prefabricated defect was designed.The opening process was simulated using peridynamics,yielding the strain field distribution and the corresponding failure mode.A single-opening verification test was conducted.The simulation results showed good agreement with the experimental data,demonstrating the reliability of the peridynamic modeling approach.Furthermore,the effects of the prefabricated defect shape and depth on the opening performance of the PSD were analyzed through simulation.The research results indicate that the established constitutive model and failure criteria based on peridynamics can reasonably predict the failure location and the opening pressure of the soft PSD.Under the impact loading,the weak zone of the soft PSD firstly ruptures,and the damaged area gradually propagates along with the prefabricated defect,eventually leading to complete separation.A smaller prefabricated defect depth or a wider prefabricated defect distribution can cause a reduction in opening pressure.These research results provide valuable guidance for the preliminary design and optimization of PSDs in pulse solid rocket motors.展开更多
As the performance of the box-type multiple launch rocket system(BMLRS)improves,its mechanical structures,particularly the plane clearance design between the slider on the rocket and the guide inside the launch canist...As the performance of the box-type multiple launch rocket system(BMLRS)improves,its mechanical structures,particularly the plane clearance design between the slider on the rocket and the guide inside the launch canister,have grown increasingly complex.However,deficiencies still exist in the current launch modeling theory for BMLRS.In this study,a multi-rigid-flexible-body launch dynamics model coupling the launch platform and rocket was established using the multibody system transfer matrix method and the Newton-Euler formulation.Furthermore,considering the bending of the launch canister,a detection algorithm for slider-guide plane clearance contact was proposed.To quantify the contact force and friction effect between the slider and guide,the contact force model and modified Coulomb model were introduced.Both the modal and launch tests were conducted.Additionally,the modal convergence was verified.By comparing the modal experiments and simulation results,the maximum relative error of the eigenfrequency is 3.29%.thereby verifying the accuracy of the developed BMLRS dynamics model.Furthermore,the launch test validated the proposed plane clearance contact model.Moreover,the study investigated the influence of various model parameters on the dynamic characteristics of BMLRS,including launch canister bending stiffness,slider and guide material,slider-guide clearance,slider length and layout.This analysis of influencing factors provides a foundation for future optimization in BMLRS design.展开更多
This study investigates the potential of metal additives in acrylonitrile butadiene styrene(ABS)polymer fuel to enhance hybrid rocket motor(HRM)performance through computational analysis,Chemical Equilibrium with Appl...This study investigates the potential of metal additives in acrylonitrile butadiene styrene(ABS)polymer fuel to enhance hybrid rocket motor(HRM)performance through computational analysis,Chemical Equilibrium with Applications(CEA),software.ABS was selected as the base fuel due to its thermoplastic nature,which allows for the creation of complex fuel geometries through 3D printing,offering significant flexibility in fuel design.Hybrid rockets,which combine a solid fuel with a liquid oxidiser,offer advantages in terms of operational simplicity and safety.However,conventional polymer fuels often exhibit low regression rates and suboptimal combustion efficiencies.In this research,we evaluated a range of metal additives-aluminium(Al),boron(B),nickel(Ni),copper(Cu),and iron(Fe)-at chamber pressures ranging from 1 to 30 bar and oxidiser-to-fuel(O/F)ratios between 1.1 and 12,resulting in 1800 unique test conditions.The main performance parameters used to assess each formulation were characteristic velocity(C^(*))and adiabatic flame temperature.The results revealed that each test produced a different optimum O/F ratio,with most ratios falling between 4 and 6.The highest performance was achieved at a chamber pressure of 30 bar across all formulations.Among the additives,Al and B demonstrated significant potential for improved combustion performance with increasing metal loadings.In contrast,Fe,Cu,and Ni reached optimal performance at a minimum loading of 1%.Future work includes investigating B-Al metal composites as additives into the ABS base polymer fuel,and doing experimental validation tests where the metallised ABS polymer fuel is 3D printed.展开更多
Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision mod...Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx<sub>1</sub>, and slow-jet, vx<sub>2</sub>. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.展开更多
During the supersonic re-entry of multi-nozzle heavy rockets into the atmosphere,the basic flow state becomes increasingly complex due to the coupling effect between the retropropulsion plumes and the freestream.A num...During the supersonic re-entry of multi-nozzle heavy rockets into the atmosphere,the basic flow state becomes increasingly complex due to the coupling effect between the retropropulsion plumes and the freestream.A numerical method using the hybrid Reynolds-Averaged Navier-Stokes and Large Eddy Simulation(RES)method and discrete coordinate method is developed to accurately estimate the thermal environment.In addition,finite rate chemical kinetics is used to calculate the afterburning reactions.The numerical results agree well with wind tunnel data,which confirms the validity and accuracy of the numerical method.Computations are conducted for the heavy carrier rocket re-entry from 53.1 km to 39.5 km altitude with 180°angle of attack by using three different Supersonic Retro-Propulsion(SRP)modes.The numerical results reveal that these three SRP flow fields are all Short Penetration Models(SPM).As the re-entry altitudes decrease,both the plume-plume interaction and the plume-freestream interaction become weaker.The highest temperatures in the plume shear layers of the three SRP modes increase by 8.36%,7.33%and 6.92%respectively after considering afterburning reactions,and all occur at a reentry altitude of 39.5 km.As the rocket re-enters the atmosphere,the maximum heat flux on the rocket base plate of three SRP modes stabilizes at 290,170 and 200 kW/m^(2) respectively,but the maximum heat flux on the side wall increases significantly.When the altitude declines to 39.5 km,the extreme heat flux of the three modes increase by 84.16%,49.45%and 62.97%respectively compared to that at 53.1 km.展开更多
This study investigates the end-burning hybrid rocket motors with polyethylene fuel by the numerical simulation and experiment.Based on computational fluid dynamics,a numerical model is developed.The model is validate...This study investigates the end-burning hybrid rocket motors with polyethylene fuel by the numerical simulation and experiment.Based on computational fluid dynamics,a numerical model is developed.The model is validated by two firing tests in this hybrid rocket motor,which uses oxygen and polyethylene as propellants.The results show that the numerical and experimental data are in good agreement,and the error of the chamber pressure is less than 2.63%.Based on the simulation mode,the blowoff limit of the end-burning hybrid rocket motors is investigated.When the nozzle throat diameter and the inner diameter of grain are large,it is more difficult for the hybrid rocket motor to achieve end-burning mode,i.e.,the flame spreading is prevented in the narrow duct.The main reason is that when the nozzle throat and the grain port are large,chamber pressure and oxidizer flow velocity are low.Therefore,the friction velocity considering the pressure and flow velocity is proposed.The critical friction velocity is about 4.054–4.890 m/s in the hybrid rocket motors.When the friction velocity exceeds the critical friction velocity,the combustion mode in hybrid rocket motors changes from the flame spreading mode to the end-burning mode.Moreover,the regression rate formula is obtained by fitting,which shows that the regression rate has a good correlation with combustion chamber pressure.The critical friction velocity and regression rate formula can provide an important reference for end-burning hybrid rocket motors.展开更多
How to effectively evaluate the firing precision of weapon equipment at low cost is one of the core contents of improving the test level of weapon system.A new method to evaluate the firing precision of the MLRS consi...How to effectively evaluate the firing precision of weapon equipment at low cost is one of the core contents of improving the test level of weapon system.A new method to evaluate the firing precision of the MLRS considering the credibility of simulation system based on Bayesian theory is proposed in this paper.First of all,a comprehensive index system for the credibility of the simulation system of the firing precision of the MLRS is constructed combined with the group analytic hierarchy process.A modified method for determining the comprehensive weight of the index is established to improve the rationality of the index weight coefficients.The Bayesian posterior estimation formula of firing precision considering prior information is derived in the form of mixed prior distribution,and the rationality of prior information used in estimation model is discussed quantitatively.With the simulation tests,the different evaluation methods are compared to validate the effectiveness of the proposed method.Finally,the experimental results show that the effectiveness of estimation method for firing precision is improved by more than 25%.展开更多
Solid rocket motors have important applications in the propulsion of trans-media vehicles and underwater launched rockets.In this paper,the ignition start-up process of an underwater solid rocket motor across a wide d...Solid rocket motors have important applications in the propulsion of trans-media vehicles and underwater launched rockets.In this paper,the ignition start-up process of an underwater solid rocket motor across a wide depth range has been numerically studied.A novel multi-domain integrated model has been developed by combining the solid propellant ignition and combustion model with the volume of fluid multiphase model.This integrated model enables the coupled simulation of the propellant combustion and gas flow inside the motor,along with the gas jet evolution in the external water environment.The detailed flow field developments in the combustion chamber,nozzle,and wake field are carefully analyzed.The variation rules of the internal ballistics and thrust performance are also obtained.The effects of environmental medium and operating depth on the ignition start-up process are systematically discussed.The results show that the influence of the operating environment on the internal ballistic characteristics is primarily reflected in the initial period after the nozzle closure opens.The development of the gas jet in water lags significantly compared with that in air.As the water depth increases,the ignition delay time of the motor is shortened,and the morphology evolution of the gas jet is significantly compressed and accelerated.Furthermore,the necking and bulging of the jet boundary near the nozzle outlet and the consequent shock oscillations are intensified,resulting in stronger fluctuations in the wake pressure field and motor thrust.展开更多
The spatiotemporal distribution characteristics of the regression rate are crucial aspects of the research on Hybrid Rocket Motor(HRM). This study presents a pioneering effort in achieving a comprehensive numerical si...The spatiotemporal distribution characteristics of the regression rate are crucial aspects of the research on Hybrid Rocket Motor(HRM). This study presents a pioneering effort in achieving a comprehensive numerical simulation of fluid dynamics and heat transfer in both the fluid and solid regions throughout the entire operation of an HRM. To accomplish this, a dynamic grid technique that incorporates fluid–solid coupling is utilized. To validate the precision of the numerical simulations, a firing test is conducted, with embedded thermocouple probes being used to measure the inner temperature of the fuel grain. The temperature variations in the solid fuel obtained from both experiment and simulations show good agreement. The maximum combustion temperature and average thrust obtained from the simulations are found to deviate from the experimental results by only 3.3% and 2.4%, respectively. Thus, it can be demonstrated that transient numerical simulations accurately capture the fluid–solid coupling characteristics and transient regression rate. The dynamic simulation results of inner flow field and solid region throughout the entire working stage reveal that the presence of vortices enhances the blending of combustion gases and improves the regression rate at both the front and rear ends of the fuel grain. In addition, oscillations of the regression rate obtained in the simulation can also be well corresponded with the corrugated surface observed in the experiment. Furthermore, the zero-dimension regression rate formula and the formula describing the axial location dependence of the regression rate are fitted from the simulation results, with the corresponding coefficients of determination(R^(2)) of 0.9765 and 0.9298, respectively.This research serves as a reference for predicting the performance of HRM with gas oxygen and polyethylene, and presents a credible way for investigating the spatiotemporal distribution of the regression rate.展开更多
The thermal protection of rocket engines is a crucial aspect of rocket engine design.In this paper,the gas film/regenerative composite cooling of the liquid oxygen/liquid methane(LOX/LCH4)rocket engine thrust chamber ...The thermal protection of rocket engines is a crucial aspect of rocket engine design.In this paper,the gas film/regenerative composite cooling of the liquid oxygen/liquid methane(LOX/LCH4)rocket engine thrust chamber was investigated.A gas film/regenerative composite cooling model was developed based on the Grisson gas film cooling efficiency formula and the one-dimensional regenerative cooling model.The accuracy of the model was validated through experiments conducted on a 6 kg/s level gas film/regenerative composite cooling thrust chamber.Additionally,key parameters related to heat transfer performance were calculated.The results demonstrate that the model is sufficiently accurate to be used as a preliminary design tool.The temperature rise error of the coolant,when compared with the experimental results,was found to be less than 10%.Although the pressure drop error is relatively large,the calculated results still provide valuable guidance for heat transfer analysis.In addition,the performance of composite cooling is observed to be superior to regenerative cooling.Increasing the gas film flow rate results in higher cooling efficiency and a lower gas-side wall temperature.Furthermore,the position at which the gas film is introduced greatly impacts the cooling performance.The optimal introduction position for the gas film is determined when the film is introduced from a single row of holes.This optimal introduction position results in a more uniform wall temperature distribution and reduces the peak temperature.Lastly,it is observed that a double row of holes,when compared to a single row of holes,enhances the cooling effect in the superposition area of the gas film and further lowers the gas-side wall temperature.These results provide a basis for the design of gas film/regenerative composite cooling systems.展开更多
Bayesian estimation is applied to the analysis of backflow vortex instabilities in typical three-and four bladed liquid propellant rocket(LPR)engine inducers.The flow in the impeller eye is modeled as a set of equally...Bayesian estimation is applied to the analysis of backflow vortex instabilities in typical three-and four bladed liquid propellant rocket(LPR)engine inducers.The flow in the impeller eye is modeled as a set of equally intense and evenly spaced 2D axial vortices,located at the same radial distance from the axis and rotating at a fraction of the impeller speed.The circle theorem and the Bernoulli’s equation are used to predict the flow pressure in terms of the vortex number,intensity,rotational speed,and radial position.The theoretical spectra so obtained are frequency broadened to mimic the dispersion of the experimental data and parametrically fitted to the measured pressure spectra by maximum likelihood estimation with equal and independent Gaussian errors.The method is applied to three inducers,tested in water at room temperature and different loads and cavitation conditions.It successfully characterizes backflow instabilities using the signals of a single pressure transducer flush-mounted on the casing of the impeller eye,effectively by-passing the aliasing and data acquisition/reduction complexities of traditional multiple-sensor cross correlation methods.The identification returns the estimates of the model parameters and their standard errors,providing the information necessary for assessing the accuracy and statistical significance of the results.The flowrate is found to be the major factor affecting the backflow vortex instability,which,on the other hand,is rather insensitive to the occurrence of cavitation.The results are consistent with the data reported in the literature,as well as with those generated by the auxiliary models specifically developed for initializing the maximum likelihood searches and supporting the identification procedure.展开更多
The solid rocket motor driven system is one of the common ways for submarines to launch underwater missiles. It has significant advantages in improving the missile’s water exit speed, anti-interference capability, an...The solid rocket motor driven system is one of the common ways for submarines to launch underwater missiles. It has significant advantages in improving the missile’s water exit speed, anti-interference capability, and enemy striking power. The prediction of the underwater loading is a preliminary factor for the power system design of the underwater vehicle. This paper presents a rapid prediction method and validated by the experimental study for the underwater thrust of the solid rocket motor. Based on the potential flow assumption of the water field, a model of the bubble and a one-dimensional quasi-steady model of the nozzle are established to directly solve the flow status of the nozzle. The aerodynamic thrust and hydrodynamic thrust have been calculated and analyzed. The calculation results are within 5% error of the experimental results. Moreover, a design platform to predict the underwater thrust of the solid rocket motor has been developed based on Python and the PyQt library, which shows excellent system adaptability and computational efficiency.展开更多
基金the experimental technology support provided by the China Academy of Launch Vehicle Technology
文摘Sea-based rocket launches encounter significant challenges stemming from dynamic marine environmental interactions.During the hot launch phase,characterized by low-velocity ascent,the departure of the rocket from the oscillatory platform exhibits heightened sensitivity to external disturbances.In the development stage,assessing the launch dynamics and the clearance between the rocket and framed launcher are crucial for improving the reliability of sea-based rocket launches in rough sea conditions.This study presents a high-fidelity dynamic model of maritime hot launch system,demonstrating 3.21%prediction error through rigorous validation against experimental datasets from comprehensive modal analyses and the full-scale rocket flight test.To mitigate collision risks,we develop a computational method employing spatial vector analysis for dynamic measurement of rocket-launcher clearance during departure.Systematic investigations reveal that in rough sea conditions,optimal departure dynamics are achieved at θ_(thrust)=270°nozzle azimuth configuration,reducing failure probability compared to conventional orientations.The developed assessment framework not only resolves critical safety challenges in current sea launch systems but also establishes foundational principles for optimizing adapter axial configuration patterns in future designs.
基金supported by the National Natural Science Foundation of China(No.12104047)。
文摘The rocket sled system is not only a high-speed dynamic ground test system,but also one of the future aerospace horizontal launch schemes.The winged load,as a common type of payload,has greater vibration and noise intensity than the wingless load.Due to the severe aerodynamic instability prior to separation,the head-up or head-down phenomena are more evident and the test accuracy significantly decreases.The high-precision computer fluid dynamics and aeroacoustic analysis are employed to explore the multifield coupling mechanism of a rocket sled with the winged payload in the wide speed range(Ma=0.5–2).The results show that as the incoming velocity increases,the cone angle of the shock wave of the rocket sled decreases,the shock pressure increases quickly,and the vortex between the slippers splits and gradually shrinks in size.The velocity of the rocket sled exerts little influence on the modal resonance frequency.The wing has a significant impact on aerodynamic noise,and as the sound pressure level rises,the propagation direction gradually shifts towards the rear and upper regions of the wing.
基金supported by the Key project of Hunan Provincial Education Department(Grant Number:22A0485)The Natural Science Foundation of Hunan(Grant Number:2024JJ5293)The Key project of Hunan University of Arts and Science(Grant Number:23ZZ08).
文摘In liquid rocket engines,regenerative cooling technology is essential for preserving structural integrity under extreme thermal loads.However,non-uniform coolant flow distribution within the cooling channels often leads to localized overheating,posing serious risks to engine reliability and operational lifespan.This study employs a three-dimensional fluid–thermal coupled numerical model to systematically investigate the influence of geometric parameters-specifically the number of inlets,the number of channels,and inlet manifold configurations-on flow uniformity and thermal distribution in non-pyrolysis zones.Key findings reveal that increasing the number of inlets from one to three significantly enhances flow uniformity,reducing mass flow rate deviation from 1.2%to below 0.3%.However,further increasing the inlets to five yields only marginal improvements indicating diminishing(<0.1%),returns beyond three inlets.Additionally,temperature non-uniformity at the combustion chamber throat decreases by 37%-from 3050 K with 18 channels to 1915 K with 30 channels-highlighting the critical role of channel density in effective thermal regulation.Notably,while higher channel counts improve cooling efficiency,they also result in increased pressure losses of approximately 18%–22%,emphasizing the need to balance thermal performance against hydraulic resistance.An optimal configuration comprising 24 channels and three inlets was identified,providing minimal temperature gradients while maintaining acceptable pressure losses.The inlet manifold structure also plays a pivotal role in determining flow distribution.Configuration 3(Config-3),which features an enlarged manifold and reduced inlet velocity,achieves a 40%reduction in velocity fluctuations compared to Configuration 1(Config-1).This improvement leads to a more uniform mass flow distribution,with a relative standard deviation(RSD)of less than 0.15%.Furthermore,this design effectively mitigates localized hot spots near the nozzle-where temperature gradients are most severe-achieving a reduction of approximately 1135 K.
基金supported by the National Natural Science Foundation of China(Nos.62120106003 and 62173301)。
文摘The reuse of liquid propellant rocket engines has increased the difficulty of their control and estimation.State and parameter Moving Horizon Estimation(MHE)is an optimization-based strategy that provides the necessary information for model predictive control.Despite the many advantages of MHE,long computation time has limited its applications for system-level models of liquid propellant rocket engines.To address this issue,we propose an asynchronous MHE method called advanced-multi-step MHE with Noise Covariance Estimation(amsMHE-NCE).This method computes the MHE problem asynchronously to obtain the states and parameters and can be applied to multi-threaded computations.In the background,the state and covariance estimation optimization problems are computed using multiple sampling times.In real-time,sensitivity is used to quickly approximate state and parameter estimates.A covariance estimation method is developed using sensitivity to avoid redundant MHE problem calculations in case of sensor degradation during engine reuse.The amsMHE-NCE is validated through three cases based on the space shuttle main engine system-level model,and we demonstrate that it can provide more accurate real-time estimates of states and parameters compared to other commonly used estimation methods.
基金Sponsored by the Natural Science Foundation of Shaanxi Province (Grant No. S2025-JC-YB-0532)the Practice and Innovation Funds for Graduate Students of Northwestern Polytechnical University (PF2024044)
文摘Taking a C1x motor with a backward-facing step which can generate a typical corner vortex as a reference,a numerical methodology using large eddy simulation was established in this study.Based on this methodology,the position of the backward-facing step of the motor was computed and analyzed to determine a basic configuration.Two key geometrical parameters,the head cavity angle and submerged nozzle cavity height,were subsequently introduced.Their effects on the corner vortex motion and their interactions with the acoustic pressure downstream of the backward-facing step were analyzed.The phenomena of vortex acoustic coupling and characteristics of pressure oscillations were further explored.The results show that the maximum error between the simulations and experimental data on the dominant frequency of pressure oscillations is 5.23%,which indicates that the numerical methodology built in this study is highly accurate.When the step is located at less than 5/8 of the total length of the combustion chamber,vortex acoustic coupling occurs,which can increase the pressure oscillations in the motor.Both the vorticity and the scale of vortices in the downstream step increase when the head cavity angle is greater than 24°,which increases the amplitude of the pressure oscillation by maximum 63.0%.The submerged nozzle cavity mainly affects the vortices in the cavity itself rather than those in the downstream step.When the height of the cavity increases from 10 to 20 mm,the pressure oscillation amplitude under the main frequency increases by 39.1%.As this height continues to increase,the amplitude of pressure oscillations increases but the primary frequency decreases.
基金supported by the Natural Science Foundation of Hunan Province of China(No.2023JJ40672)the Innovation Science Fund Project of National University of Defense Technology,China(No.ZK2023-039)。
文摘A two-dimensional large eddy simulation numerical model is proposed to study the transient vortex flow and pressure oscillation of a large-aspect-ratio solid rocket motor.The numerical model is validated through experimental data,finite element analysis and cumulative error analysis.The numerical simulations are executed to obtain the characteristics of the vortex-acoustic and pressure oscillation.The results show that the burning surface regression decreases the motor aspect ratio,increasing the corresponding natural frequency from 260 Hz to 293 Hz.The pressure oscillation phenomenon is formed due to the vortex-acoustic coupling.Decreasing the corner vortex shedding intensity shows negative effects on the dimensionless amplitude of the pressure oscillation.The head cavity without the injection can decrease the vortex-acoustic coupling level at the acoustic pressure antinode.The modified motor with head cavity can obtain a lower dimensionless oscillating pressure amplitude 0.00149 in comparison with 0.00895 of the original motor.The aspect ratio and volume of the head cavity without the injection have great effects on the pressure oscillation suppression,particularly at the low aspect ratio or large volume.The reason is that the mass in the region around the acoustic pressure antinode is extracted centrally,reducing the energy contribution to the acoustic system.With the volume increasing,the acoustic energy capacity increases.
基金supported by the National Natural Science Foundation of China(No.12202011)the Youth Research fund of Shanghai Academy of Spaceflight Technology(KJW-KT-QNKYJJ-2022-25)China Postdoctoral Science Foundation(Nos.2024T170009,2022M710190).
文摘As a critical component of pulse solid rocket motors(SRMs),the soft pulse separation device(PSD)is vital in enabling multi-pulse propulsion and has become a breakthrough in SRM engineering applications.To investigate the opening performance of the PSD,an axial PSD incorporating a star-shaped prefabricated defect was designed.The opening process was simulated using peridynamics,yielding the strain field distribution and the corresponding failure mode.A single-opening verification test was conducted.The simulation results showed good agreement with the experimental data,demonstrating the reliability of the peridynamic modeling approach.Furthermore,the effects of the prefabricated defect shape and depth on the opening performance of the PSD were analyzed through simulation.The research results indicate that the established constitutive model and failure criteria based on peridynamics can reasonably predict the failure location and the opening pressure of the soft PSD.Under the impact loading,the weak zone of the soft PSD firstly ruptures,and the damaged area gradually propagates along with the prefabricated defect,eventually leading to complete separation.A smaller prefabricated defect depth or a wider prefabricated defect distribution can cause a reduction in opening pressure.These research results provide valuable guidance for the preliminary design and optimization of PSDs in pulse solid rocket motors.
基金supported by National Natural Science Foundation of China(Grant No.92266201).
文摘As the performance of the box-type multiple launch rocket system(BMLRS)improves,its mechanical structures,particularly the plane clearance design between the slider on the rocket and the guide inside the launch canister,have grown increasingly complex.However,deficiencies still exist in the current launch modeling theory for BMLRS.In this study,a multi-rigid-flexible-body launch dynamics model coupling the launch platform and rocket was established using the multibody system transfer matrix method and the Newton-Euler formulation.Furthermore,considering the bending of the launch canister,a detection algorithm for slider-guide plane clearance contact was proposed.To quantify the contact force and friction effect between the slider and guide,the contact force model and modified Coulomb model were introduced.Both the modal and launch tests were conducted.Additionally,the modal convergence was verified.By comparing the modal experiments and simulation results,the maximum relative error of the eigenfrequency is 3.29%.thereby verifying the accuracy of the developed BMLRS dynamics model.Furthermore,the launch test validated the proposed plane clearance contact model.Moreover,the study investigated the influence of various model parameters on the dynamic characteristics of BMLRS,including launch canister bending stiffness,slider and guide material,slider-guide clearance,slider length and layout.This analysis of influencing factors provides a foundation for future optimization in BMLRS design.
文摘This study investigates the potential of metal additives in acrylonitrile butadiene styrene(ABS)polymer fuel to enhance hybrid rocket motor(HRM)performance through computational analysis,Chemical Equilibrium with Applications(CEA),software.ABS was selected as the base fuel due to its thermoplastic nature,which allows for the creation of complex fuel geometries through 3D printing,offering significant flexibility in fuel design.Hybrid rockets,which combine a solid fuel with a liquid oxidiser,offer advantages in terms of operational simplicity and safety.However,conventional polymer fuels often exhibit low regression rates and suboptimal combustion efficiencies.In this research,we evaluated a range of metal additives-aluminium(Al),boron(B),nickel(Ni),copper(Cu),and iron(Fe)-at chamber pressures ranging from 1 to 30 bar and oxidiser-to-fuel(O/F)ratios between 1.1 and 12,resulting in 1800 unique test conditions.The main performance parameters used to assess each formulation were characteristic velocity(C^(*))and adiabatic flame temperature.The results revealed that each test produced a different optimum O/F ratio,with most ratios falling between 4 and 6.The highest performance was achieved at a chamber pressure of 30 bar across all formulations.Among the additives,Al and B demonstrated significant potential for improved combustion performance with increasing metal loadings.In contrast,Fe,Cu,and Ni reached optimal performance at a minimum loading of 1%.Future work includes investigating B-Al metal composites as additives into the ABS base polymer fuel,and doing experimental validation tests where the metallised ABS polymer fuel is 3D printed.
文摘Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx<sub>1</sub>, and slow-jet, vx<sub>2</sub>. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.
基金co-supported by the National Level Project,China and Shanghai Municipal Major Science and Technology Project,China.In addition,the authors gratefully acknowledge the guidance on English writing and numerical methods of Professor Junfeng ZHANG from the Faculty of Engineering at Laurentian University.
文摘During the supersonic re-entry of multi-nozzle heavy rockets into the atmosphere,the basic flow state becomes increasingly complex due to the coupling effect between the retropropulsion plumes and the freestream.A numerical method using the hybrid Reynolds-Averaged Navier-Stokes and Large Eddy Simulation(RES)method and discrete coordinate method is developed to accurately estimate the thermal environment.In addition,finite rate chemical kinetics is used to calculate the afterburning reactions.The numerical results agree well with wind tunnel data,which confirms the validity and accuracy of the numerical method.Computations are conducted for the heavy carrier rocket re-entry from 53.1 km to 39.5 km altitude with 180°angle of attack by using three different Supersonic Retro-Propulsion(SRP)modes.The numerical results reveal that these three SRP flow fields are all Short Penetration Models(SPM).As the re-entry altitudes decrease,both the plume-plume interaction and the plume-freestream interaction become weaker.The highest temperatures in the plume shear layers of the three SRP modes increase by 8.36%,7.33%and 6.92%respectively after considering afterburning reactions,and all occur at a reentry altitude of 39.5 km.As the rocket re-enters the atmosphere,the maximum heat flux on the rocket base plate of three SRP modes stabilizes at 290,170 and 200 kW/m^(2) respectively,but the maximum heat flux on the side wall increases significantly.When the altitude declines to 39.5 km,the extreme heat flux of the three modes increase by 84.16%,49.45%and 62.97%respectively compared to that at 53.1 km.
基金supported by the National Natural Science Foundation of China(No.U20B2034)the Academic Excellence Foundation of Beihang University for PhD Students,China.
文摘This study investigates the end-burning hybrid rocket motors with polyethylene fuel by the numerical simulation and experiment.Based on computational fluid dynamics,a numerical model is developed.The model is validated by two firing tests in this hybrid rocket motor,which uses oxygen and polyethylene as propellants.The results show that the numerical and experimental data are in good agreement,and the error of the chamber pressure is less than 2.63%.Based on the simulation mode,the blowoff limit of the end-burning hybrid rocket motors is investigated.When the nozzle throat diameter and the inner diameter of grain are large,it is more difficult for the hybrid rocket motor to achieve end-burning mode,i.e.,the flame spreading is prevented in the narrow duct.The main reason is that when the nozzle throat and the grain port are large,chamber pressure and oxidizer flow velocity are low.Therefore,the friction velocity considering the pressure and flow velocity is proposed.The critical friction velocity is about 4.054–4.890 m/s in the hybrid rocket motors.When the friction velocity exceeds the critical friction velocity,the combustion mode in hybrid rocket motors changes from the flame spreading mode to the end-burning mode.Moreover,the regression rate formula is obtained by fitting,which shows that the regression rate has a good correlation with combustion chamber pressure.The critical friction velocity and regression rate formula can provide an important reference for end-burning hybrid rocket motors.
基金National Natural Science Foundation of China(Grant Nos.11972193 and 92266201)。
文摘How to effectively evaluate the firing precision of weapon equipment at low cost is one of the core contents of improving the test level of weapon system.A new method to evaluate the firing precision of the MLRS considering the credibility of simulation system based on Bayesian theory is proposed in this paper.First of all,a comprehensive index system for the credibility of the simulation system of the firing precision of the MLRS is constructed combined with the group analytic hierarchy process.A modified method for determining the comprehensive weight of the index is established to improve the rationality of the index weight coefficients.The Bayesian posterior estimation formula of firing precision considering prior information is derived in the form of mixed prior distribution,and the rationality of prior information used in estimation model is discussed quantitatively.With the simulation tests,the different evaluation methods are compared to validate the effectiveness of the proposed method.Finally,the experimental results show that the effectiveness of estimation method for firing precision is improved by more than 25%.
基金supported by the National Level Project of China.
文摘Solid rocket motors have important applications in the propulsion of trans-media vehicles and underwater launched rockets.In this paper,the ignition start-up process of an underwater solid rocket motor across a wide depth range has been numerically studied.A novel multi-domain integrated model has been developed by combining the solid propellant ignition and combustion model with the volume of fluid multiphase model.This integrated model enables the coupled simulation of the propellant combustion and gas flow inside the motor,along with the gas jet evolution in the external water environment.The detailed flow field developments in the combustion chamber,nozzle,and wake field are carefully analyzed.The variation rules of the internal ballistics and thrust performance are also obtained.The effects of environmental medium and operating depth on the ignition start-up process are systematically discussed.The results show that the influence of the operating environment on the internal ballistic characteristics is primarily reflected in the initial period after the nozzle closure opens.The development of the gas jet in water lags significantly compared with that in air.As the water depth increases,the ignition delay time of the motor is shortened,and the morphology evolution of the gas jet is significantly compressed and accelerated.Furthermore,the necking and bulging of the jet boundary near the nozzle outlet and the consequent shock oscillations are intensified,resulting in stronger fluctuations in the wake pressure field and motor thrust.
基金supported by the National Natural Science Foundation of China (No.U20B2034).
文摘The spatiotemporal distribution characteristics of the regression rate are crucial aspects of the research on Hybrid Rocket Motor(HRM). This study presents a pioneering effort in achieving a comprehensive numerical simulation of fluid dynamics and heat transfer in both the fluid and solid regions throughout the entire operation of an HRM. To accomplish this, a dynamic grid technique that incorporates fluid–solid coupling is utilized. To validate the precision of the numerical simulations, a firing test is conducted, with embedded thermocouple probes being used to measure the inner temperature of the fuel grain. The temperature variations in the solid fuel obtained from both experiment and simulations show good agreement. The maximum combustion temperature and average thrust obtained from the simulations are found to deviate from the experimental results by only 3.3% and 2.4%, respectively. Thus, it can be demonstrated that transient numerical simulations accurately capture the fluid–solid coupling characteristics and transient regression rate. The dynamic simulation results of inner flow field and solid region throughout the entire working stage reveal that the presence of vortices enhances the blending of combustion gases and improves the regression rate at both the front and rear ends of the fuel grain. In addition, oscillations of the regression rate obtained in the simulation can also be well corresponded with the corrugated surface observed in the experiment. Furthermore, the zero-dimension regression rate formula and the formula describing the axial location dependence of the regression rate are fitted from the simulation results, with the corresponding coefficients of determination(R^(2)) of 0.9765 and 0.9298, respectively.This research serves as a reference for predicting the performance of HRM with gas oxygen and polyethylene, and presents a credible way for investigating the spatiotemporal distribution of the regression rate.
基金supported by the National Science Fund Project(No.2019-JCJQ-ZQ-019)the Innovative Research Group Project of National Natural Science Foundation of China(No.T2221002).
文摘The thermal protection of rocket engines is a crucial aspect of rocket engine design.In this paper,the gas film/regenerative composite cooling of the liquid oxygen/liquid methane(LOX/LCH4)rocket engine thrust chamber was investigated.A gas film/regenerative composite cooling model was developed based on the Grisson gas film cooling efficiency formula and the one-dimensional regenerative cooling model.The accuracy of the model was validated through experiments conducted on a 6 kg/s level gas film/regenerative composite cooling thrust chamber.Additionally,key parameters related to heat transfer performance were calculated.The results demonstrate that the model is sufficiently accurate to be used as a preliminary design tool.The temperature rise error of the coolant,when compared with the experimental results,was found to be less than 10%.Although the pressure drop error is relatively large,the calculated results still provide valuable guidance for heat transfer analysis.In addition,the performance of composite cooling is observed to be superior to regenerative cooling.Increasing the gas film flow rate results in higher cooling efficiency and a lower gas-side wall temperature.Furthermore,the position at which the gas film is introduced greatly impacts the cooling performance.The optimal introduction position for the gas film is determined when the film is introduced from a single row of holes.This optimal introduction position results in a more uniform wall temperature distribution and reduces the peak temperature.Lastly,it is observed that a double row of holes,when compared to a single row of holes,enhances the cooling effect in the superposition area of the gas film and further lowers the gas-side wall temperature.These results provide a basis for the design of gas film/regenerative composite cooling systems.
文摘Bayesian estimation is applied to the analysis of backflow vortex instabilities in typical three-and four bladed liquid propellant rocket(LPR)engine inducers.The flow in the impeller eye is modeled as a set of equally intense and evenly spaced 2D axial vortices,located at the same radial distance from the axis and rotating at a fraction of the impeller speed.The circle theorem and the Bernoulli’s equation are used to predict the flow pressure in terms of the vortex number,intensity,rotational speed,and radial position.The theoretical spectra so obtained are frequency broadened to mimic the dispersion of the experimental data and parametrically fitted to the measured pressure spectra by maximum likelihood estimation with equal and independent Gaussian errors.The method is applied to three inducers,tested in water at room temperature and different loads and cavitation conditions.It successfully characterizes backflow instabilities using the signals of a single pressure transducer flush-mounted on the casing of the impeller eye,effectively by-passing the aliasing and data acquisition/reduction complexities of traditional multiple-sensor cross correlation methods.The identification returns the estimates of the model parameters and their standard errors,providing the information necessary for assessing the accuracy and statistical significance of the results.The flowrate is found to be the major factor affecting the backflow vortex instability,which,on the other hand,is rather insensitive to the occurrence of cavitation.The results are consistent with the data reported in the literature,as well as with those generated by the auxiliary models specifically developed for initializing the maximum likelihood searches and supporting the identification procedure.
文摘The solid rocket motor driven system is one of the common ways for submarines to launch underwater missiles. It has significant advantages in improving the missile’s water exit speed, anti-interference capability, and enemy striking power. The prediction of the underwater loading is a preliminary factor for the power system design of the underwater vehicle. This paper presents a rapid prediction method and validated by the experimental study for the underwater thrust of the solid rocket motor. Based on the potential flow assumption of the water field, a model of the bubble and a one-dimensional quasi-steady model of the nozzle are established to directly solve the flow status of the nozzle. The aerodynamic thrust and hydrodynamic thrust have been calculated and analyzed. The calculation results are within 5% error of the experimental results. Moreover, a design platform to predict the underwater thrust of the solid rocket motor has been developed based on Python and the PyQt library, which shows excellent system adaptability and computational efficiency.
