Effect of Time Step Size and Turbulence Model on the Open Water Hydrodynamic Performance Prediction of Contra-Rotating Propellers

A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R ＆ D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.

A Model of Steady-State Convective Combustion of Micropore Propellants

Aim To develop a physical and mathematical model related to micropore para- meters of steady-state convective combustion of micropore propellants(MPP). Methods The micropore parameters were measured by WXT-88 mage analysis apparatus and the convective combustion characteristic of MPP was measured by a large volume closed bomb, respectively. Rasults Statistical physical model of burning in the micropore and granular burning were developed. The burning rate equation of steady-state convective combustion of MPP was obtained. Conclusions This model correlates the convective burning rate with micropores para- meters for the first time,and the calculating values of convective burning rate are in agreement with test results.The model also can be used to estimate the effects of microporous parame- ters, basic mass burning rate, MPP density and pressure in combustion chamber on the convective combustion characteristics of MPP.

Scale Effects on Propeller Cavitating Hydrodynamic and Hydroacoustic Performances with Non-uniform Inflow

Considering the lack of theoretical models and ingredients necessary to explain the scaling of the results of propeller cavitation inception and cavitating hydroacoustics from model tests to full scale currently, and the insufficient reflection of the nuclei effects on cavitation in the numerical methods, the cavitating hydrodynamics and cavitation low frequency noise spectrum of three geometrically similar 7-bladed highly skewed propellers with non-uniform inflow are addressed. In this process, a numerical bridge from the multiphase viscous simulation of propeller cavitation hydrodynamics to its hydro-acoustics is built, and the scale effects on performances and the applicability of exist scaling law are analyzed. The effects of non-condensable gas（NCG） on cavitation inception are involved explicitly in the improved Sauer＇s cavitation model, and the cavity volume acceleration related to its characteristic length is used to produce the noise spectrum. Results show that, with the same cavitation number, the cavity extension on propeller blades increases with diameter associated with an earlier shift of the beginning point of thrust decline induced by cavitation, while the three decline slopes of thrust breakdown curves are found to be nearly the same. The power of the scaling law based on local Reynolds number around 0.9R section is determined as 0.11. As for the smallest propeller, the predominant tonal noise is located at blade passing frequency（BPF）, whereas 2BPF for the middle and both 2BPF and 3BPF for the largest, which shows the cavitating line spectrum is fully related to the interaction between non-uniform inflow and fluctuated cavity volume. The predicted spectrum level exceedance from the middle to the large propeller is 6.65 dB at BPF and 5.94 dB at 2BPF. Since it just differs less than 2 dB to the increment obtained by empirical scaling law, it is inferred that the scale effects on them are acceptable with a sufficient model scale, and so do the scaling law. The numerical implementation of cavitating hydrodynamics and hydro-acoustics prediction of propeller in big scale in wake has been completed.

Research on the influences of confining pressure and strain rate on NEPE propellant:Experimental assessment and constitutive model

In order to study the influences of confining pressure and strain rate on the mechanical properties of the Nitrate Ester Plasticized Polyether(NEPE)propellant,uniaxial tensile tests were conducted using the selfmade confining pressure system and material testing machine.The stress-strain responses of the NEPE propellant under different confining pressure conditions and strain rates were obtained and analyzed.The results show that confining pressure and strain rate have a remarkably influence on the mechanical responses of the NEPE propellant.As confining pressure increases(from 0 to 5.4 MPa),the maximum tensile stress and ultimate strain increase gradually.With the coupled effects of confining pressure and strain rate,the value of the maximum tensile stress and ultimate strain at 5.4 MPa and 0.0667 s^(-1)is 2.03 times and 2.19 times of their values under 0 MPa and 0.00333 s^(-1),respectively.Afterwards,the influence mechanism of confining pressure on the NEPE propellant was analyzed.Finally,based on the viscoelastic theory and continuous damage theory,a nonlinear constitutive model considering confining pressure and strain rate was developed.The damage was considered to be rate-dependent and pressuredependent.The constitutive model was validated by comparing experimental data with predictions of the constitutive model.The whole maximum stress errors of the model predictions are lower than 4%and the corresponding strain errors are lower than 7%.The results show that confining pressure can suppress the damage initiation and evolution of the NEPE propellant and the nonlinear constitutive model can describe the mechanical responses of the NEPE propellant under various confining pressure conditions and strain rates.This research can lay a theoretical foundation for analyzing the structural integrity of propellant grain accurately under working pressure loading.

Modeling and simulation of solvent behavior and temperature distribution within long stick propellants with large web thickness undergoing drying

Drying is a complicated physical process which involves simultaneous heat and mass transfer in the removal of solvents inside propellants.Inappropriate drying techniques may result in the formation of a hard skin layer near the surface to block the free access of most solvent through for long stick propellants with large web thickness,which lead to lower drying efficiency and worse drying quality.This study aims to gain a comprehensive understanding of drying process and clarify the mechanism of the blocked layer near the propellant surface.A new three-dimensional coupled heat and mass transfer(3D-CHMT)model was successfully developed under transient conditions.The drying experiment results show that the 3DCHMT model could be applied to describe the drying process well since the relative error of the content of solvent between simulation and experiment values is only 5.5%.The solvent behavior simulation demonstrates that the mass transfer process can be divided into super-fast(SF)and subsequent minorfast(MF)stages,and the SF stage is vital to the prevention of the blocked layer against the free access for solvent molecules inside propellant grains.The effective solvent diffusion coefficient(Deff)of the propellant surface initially increases from 3.4×10^(-6)to 5.3×10^(-6)m^(2)/s as the temperature increases,and then decreases to 4.1×10^(-8)m^(2)/s at 60-100 min.The value of Deffof surface between 0-1.4 mm has a unique trend of change compared with other regions,and it is much lower than that of the internal at100 min under simulation conditions.Meanwhile,the temperature of the propellant surface increases rapidly at the SF stage(0-100 min)and then very slowly thereafter.Both the evolution of Deffand temperature distribution demonstrate that the blocked layer near the propellant surface has been formed in the time period of approximately 0-100 min and its thickness is about 1.4 mm.To mitigate the formation of blocked layer and improve its drying quality of finial propellant products effectively,it should be initially dried at lower drying temperature(30-40℃)in 0-100 min and then dried at higher drying temperature(50-60℃)to reduce drying time for later drying process in double base gun propellants.The present results can provide theoretical guidance for drying process and optimization of drying parameters for long stick propellants with large web thickness.

Numerical Study on the Unsteady Characteristics of the Propeller Cavitation in Uniform and Nonuniform Wake Flows

Propeller cavitation is a problematic issue because of its negative effects, such as performances losses, noise,vibration and erosion. Numerical methodology is an effective and efficient technical tool for the study of propeller cavitation, however, it is hard to capture tip-vortex cavitation in the previous work by using common turbulence models based on turbulent-viscosity hypothesis. In this work, the Reynolds-Averaged Naiver-Stokes(RANS)approach, adopting the Reynolds stress turbulence model(RSM), is taken to study the unsteady characteristics of the cavitation on the four-bladed INSEAN E779 A model propeller. The numerical simulation was carried out using the commercial CFD software ANSYS Fluent 14.0. One kind of uniform wake flow and two kinds of nonuniform wake flows are considered here. The results in the uniform flow show a good agreement with previous experimental results on both the sheet cavitation and the tip vortex cavitation and prove the ability of the RSM on capturing the tip vortex cavitation. Two kinds of nonuniform wake flows are designed based on the previous experimental researches and the unsteady characteristics of the propeller cavitation are analyzed by comparing the results in the uniform and two nonuniform wake flows together.

Research on Propeller Dynamic Load Simulation System of Electric Propulsion Ship

A dynamic marine propeller simulation system was developed, which is utilized for meeting the experimental requirement of theory research and engineering design of marine electric propulsion system. By applying an actual ship parameter and its accurate propeller J＇ -KT＇ and J＇ - Kp＇ curve data, functional experiments based on the simulation system were carried out. The experiment results showed that the system can correctly emulate the propeller characteristics, produce the dynamic and steady performances of the propeller under different navigation modes, and present actual load torque for electric propulsion motor.

Application of transient burning rate model of solid propellant in electrothermal-chemical launch simulation

A 30 mm electrothermal-chemical(ETC) gun experimental system is employed to research the burning rate characteristics of 4/7 high-nitrogen solid propellant. Enhanced gas generation rates(EGGR) of propellants during and after electrical discharges are verified in the experiments. A modified 0D internal ballistic model is established to simulate the ETC launch. According to the measured pressure and electrical parameters, a transient burning rate law including the influence of EGGR coefficient by electric power and pressure gradient(dp/dt) is added into the model. The EGGR coefficient of 4/7 high-nitrogen solid propellant is equal to 0.005 MW-1. Both simulated breech pressure and projectile muzzle velocity accord with the experimental results well. Compared with Woodley's modified burning rate law, the breech pressure curves acquired by the transient burning rate law are more consistent with test results. Based on the parameters calculated in the model, the relationship among propellant burning rate, pressure gradient(dp/dt) and electric power is analyzed. Depending on the transient burning rate law and experimental data, the burning of solid propellant under the condition of plasma is described more accurately.

Experimental Study of the Collective and Cyclic Pitch Propeller for an Underwater Vehicle

A series of experimental studies of the innovative propulsor named Collective and Cyclic Pitch Propeller(CCPP) applied to an underwater vehicle were designed and performed at the Australian Maritime College, University of Tasmania. The bollard pull and captive model tests were conducted to investigate the characteristics of CCPP and to examine the effect of different parameter settings to its performance. The results show that the CCPP is able to generate effective manoeuvring forces in various operational condition. In addition, the obtained results in the form of force coefficients provide a useful empirical model for the simulation and control of an underwater vehicle equipped with this propulsor.

Flight Transient Estimation of VTOL Aircraft with Propellers

Methodological issues associated with the determination of the vertical take-off and landing aerodynamic parameters equipped with two rotary propellers during take-off and hovering, descent and landing are studied in the proposed article. During the computer simulation process, kinematics parameters diagrams were made, aerodynamic coefficients and propellers thrust components at all stages of aircraft take-off were estimated. That numerical data can be used in a preliminary stage of aerodynamic design for the vertical take-off and landing aircraft and electric drones at the determination of control and equalization elements geometric and kinematic parameters.

Effect of neutral polymeric bonding agent on tensile mechanical properties and damage evolution of NEPE propellant

Introducing Neutral Polymeric bonding agents(NPBA) into the Nitrate Ester Plasticized Polyether(NEPE)propellant could improve the adhesion between filler/matrix interface, thereby contributing to the development of new generations of the NEPE propellant with better mechanical properties. Therefore,understanding the effects of NPBA on the deformation and damage evolution of the NEPE propellant is fundamental to material design and applications. This paper studies the uniaxial tensile and stress relaxation responses of the NEPE propellant with different amounts of NPBA. The damage evolution in terms of interface debonding is further investigated using a cohesive-zone model(CZM). Experimental results show that the initial modulus and strength of the NEPE propellant increase with the increasing amount of NPBA while the elongation decreases. Meanwhile, the relaxation rate slows down and a higher long-term equilibrium modulus is reached. Experimental and numerical analyses indicate that interface debonding and crack propagation along filler-matrix interface are the dominant damage mechanism for the samples with a low amount of NPBA, while damage localization and crack advancement through the matrix are predominant for the ones with a high amount of NPBA. Finally, crosslinking density tests and simulation results also show that the effect of the bonding agent is interfacial rather than due to the overall crosslinking density change of the binder.

Estimation of surface geometry on combustion characteristics of AP/HTPB propellant under rapid depressurization

The 2D sandwich model serves as a potent tool in exploring the influence of surface geometry on the combustion attributes of Ammonium perchlorate/Hydroxyl-terminated polybutadiene(AP/HTPB)propellant under rapid pressure decay.The thickness of the sandwich propellant is derived from slicing the 3D random particle packing,an approach that enables a more effective examination of the micro-flame structure.Comparative analysis of the predicted burning characteristics has been performed with experimental studies.The findings demonstrate a reasonable agreement,thereby validating the precision and soundness of the model.Based on the typical rapid depressurization environment of solid rocket motor(initial combustion pressure is 3 MPa and the maximum depressurization rate is 1000 MPa/s).A-type(a flatter surface),B-type(AP recesses from the combustion surface),and C-type(AP protrudes from the combustion surface)propellant combustion processes are numerically simulated.Upon comparison of the evolution of gas-phase flame between 0.1 and 1 ms,it is discerned that the flame strength and form created by the three sandwich models differ significantly at the beginning stage of depressurization,with the flame structures gradually becoming harmonized over time.Conclusions are drawn by comparison extinction times:the surface geometry plays a pivotal role in the combustion process,with AP protrusion favoring combustion the most.

HYBRID APPROACH FOR ESTIMATING COUPLING EFFECT BETWEEN PROPELLANT SLOSHING DYNAMIC AND SPACECRAFT GNC

An approach based on equivalent mechanics theory and computational fluid dynamics （CFD） technology is proposed to estimate dynamical influence of propellant sloshing on the spacecraft. A mechanical model is estab- lished by using CFD technique and packed as a ＂sloshing＂ block used in spacecraft guidance navigation and control （GNC） simulation loop. The block takes motion characteristics of the spacecraft as inputs and outputs of pertur- bative force and torques induced by propellant sloshing, thus it is more convenient for analyzing coupling effect between propellant sloshing dynamic and spacecraft GNC than using CFD packages. An example demonstrates the accuracy and the superiority of the approach. Then, the deducing process is applied to practical cases, and simulation results validate that the proposed approach is efficient for identifying the problems induced by sloshing and evaluating effectiveness of several typical designs of sloshing suppression.

Numerical Simulation of Particle/Matrix Interface Failure in Composite Propellant

Interface debonding between particle and matrix in composite propellant influences its macroscopic mechanical properties greatly. For this, the laws of interface cohesive damage and failure were analyzed. Then, its microscopic computational model was established. The interface mechanical response was modeled by the bilinear cohesive zone model. The effects of interface properties and particle sizes on the macroscopic mechanical behavior were investigated. Numerical simulation of debonding damage evolution of composite propellant under finite deformation was carried out. The debonding damage nucleation, propagation mechanism and non-uniform distribution of microscopic stress-strain fields were discussed. The results show that the finite element simulation method based on microstructure model can effectively predict the trend of macroscopic mechanical behavior and particle/matrix debonding evolution process. It can be used for damage simulation and failure assessment of composite propellants.

Multi-scale simulation of diffusion behavior of deterrent in propellant

Concentration distribution of the deterrent in single-base propellant during the process of firing plays an important role in the ballistic properties of gun propellant in weapons. However, the diffusion coefficient calculated by molecular dynamics(MD) simulation is 6 orders of magnitude larger than the experimental values. Meanwhile, few simple and comprehensive theoretical models can explain the phenomenon and accurately predict the concentration distribution of the propellant. Herein, an onion model combining with MD simulation and finite element method of diffusion in propellants is introduced to bridge the gap between the experiments and simulations, and correctly predict the concentration distribution of deterrent. Furthermore, a new time scale is found to characterize the diffusion process. Finally, the time-and position-depended concentration distributions of dibutyl phthalate in nitrocellulose are measured by Raman spectroscopy to verify the correctness of the onion model. This work not only provides guidance for the design of the deterrent, but could be also extended to the diffusion of small molecules in polymer with different crystallinity.

Modelling Study to Compare the Flow and Heat Transfer Characteristics of Low-Power Hydrogen,Nitrogen and Argon Arc-Heated Thrusters

A modelling study is performed to compare the plasma flow and heat transfer characteristics of low-power arc-heated thrusters （arcjets） for three different propellants： hydrogen, nitrogen and argon. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. The temperature, velocity and Mach number distributions calculated within the thruster nozzle obtained with different propellant gases are compared for the same thruster structure, dimensions, inlet-gas stagnant pressure and arc currents. The temperature distributions in the solid region of the anode-nozzle wall are also given. It is found that the flow and energy conversion processes in the thruster nozzle show many similar features for all three propellants. For example, the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip; the flow transition from the subsonic to supersonic regime occurs within the constrictor region; the highest axial velocity appears inside the nozzle; and most of the input propellant flows towards the thruster exit through the cooler gas region near the anode-nozzle wall. However, since the properties of hydrogen, nitrogen and argon, especially their molecular weights, specific enthMpies and thermal conductivities, are different, there are appreciable differences in arcjet performance. For example, compared to the other two propellants, the hydrogen arcjet thruster shows a higher plasma temperature in the arc region, and higher axial velocity but lower temperature at the thruster exit. Correspondingly, the hydrogen arcjet thruster has the highest specific impulse and arc voltage for the same inlet stagnant pressure and arc current. The predictions of the modelling are compared favourably with available experimental results.

Dynamic Analysis of the Seafloor Pilot Miner Based on Single-Body Vehicle Model and Discretized Track-Terrain Interaction Model

In order to achieve the complex dynamic analysis of the self-propelled seafloor pilot miner moving on the seafloor of extremely cohesive soft soil and further to make it possible to integrate the miner system with some subsystems to form the complete integrated deep ocean mining pilot system and perform dynamic analysis, a new method for the dynamic modeling and analysis of the miner is proposed and developed in this paper, resulting in a simplified 3D single-body vehicle model with three translational and three rotational degrees of freedom, while the track-terrain interaction model is built by partitioning the track-terrain interface into discrete elements with parameterized force dements built on the theory of terramechanics acting on each discrete dement. To evaluate and verify the correctness and effectiveness of this new modeling and analysis method, typical comparative studies with regard to computational efficiency and solution accuracy are carried out between the traditional modeling method of building the tracked vehicle as a multi-body model and the new modeling method. In full consideration of the particMar structure design of the pilot miner, the special characteristics of the seafioor soil and the hydrodynamic force of near-seafloor currnt, the dynamic simulation analysis of the miner is performed and discussed, which can provide useful guidance and reference for the practical miner system in design and operation. This new method can not only realize the rapid dynamic simulation analysis of the miner but also make possible the integration and rapid dynamic analysis of the complete integrated deep ocean mining pilot system in further researches.

基于Modelica语言的某气垫船机桨匹配建模及仿真

为满足某型气垫船推进系统的性能指标,开展基于国产化软件平台的机桨匹配算法研究。以某型气垫船的推进系统为对象,建立数学模型,并以基于Modelica语言的国产仿真平台MWorks为基础对该系统进行建模仿真。通过研究变螺距角、变转速时气垫船推进系统的性能,对机桨匹配算法进行优化。仿真结果表明,优化算法可满足螺距角变距时间要求、螺距角稳态误差要求以及燃机机动性能曲线要求,所提出的机桨匹配控制算法满足性能要求,研究结果可为类似气垫船推进系统设计提供借鉴。

The Application of Erosive Burning to Propellant Charge Interior Ballistics

Erosive burning is a common burning phenomenon of the gunpowder with inner holes.The actual combustion law of the gunpowder with inner holes can be changed by erosive burning.Pressure difference between the inner and the outer of hole caused by loading density variation of the propellant charge makes erosive burning occur at inner holes during in-bore burning.The effect of erosive burning on burning speed of the propellant is studied by using the effects of flow rate,heat transfer and erosion of the combustion gas in inner holes on burning rate.The mathematic model of erosive burning of the propellant is established.The effects of the factors such as loading density,inner hole size and grain length on erosive burning and interior ballistic performance are analyzed.The method to improve the bore pressure for small charge mass and small firing range by erosive burning is proposed.

Traction Potential Analysis of Self-Propelled Seafloor Trencher Based on Mechanical Measurements in Bentonite Soil

Shear stress-displacement relationship model of soil is very important to predict the tractive performance of tracked vehicle. Most shear stress-displacement models were proposed for terrestrial field.However,they are not suit for soft seafloor with flow surface and high water content. Based on comprehensive analysis of seafloor soil shear deformation and track segment shear tests,a new empirical model of shear stressdisplacement relationship for saturated soft-plastic soil(SSP model) was proposed. To validate the SSP model and evaluate potential tractive force of self-propelled seafloor trencher,a test platform,where track segment shear test and drawbar pull test can be performed in seafloor soil substitute(bentonite water mixture),was built. Series of shear tests were carried out. Test results show that the SSP model can describe the mechanical behavior of track segment in seafloor soil substitute with good approximation. Through analyzing the main external forces,including environmental loads from seafloor soil and sea current applied to seafloor tracked trencher during the trenching process in a straight line,drawbar pull analysis model was deduced with the SSP model. A scale test model of seafloor tracked trencher was built,and the verification tests for drawbar pull analysis model were designed and carried out. Results of verification tests indicate that the drawbar pull analysis model was feasible and effective. The drawbar pull tests also indicated that the SSP model is valid from another side.