Analyse du fonctionnement d'un propulseur a éjection modulée

ONERA developed, for studying the response of a propellant to a pressure or velocity fluctuation, an experimental rocket engine whose nozzle throat area can be modulated by a toothed disk.The paper presents a linearized theory of the functioning of this engine in the low frequency domain, i.e. when there is no wave propagation within the combuster.To describe the functioning of this motor, the Ryazantsev-Novozhilov method, which assumes that the gas response is instantaneous, is used. This analysis takes into account the erosion and radiation effects, the combustion efficiency and the thermal losses through the walls.Two particular cases are described, for two values of the Damköhler parameter $\text{D}{}_1\text{=}\text{t}{}_{\mathrm{c}}\text{t}{}_{\mathrm{th}}$, where t_c is the residence time in the combuster and t_th the characteristic thermal time of the heat penetration into the solid propellant. These two cases correspond, one to a classical propellant D₁ > 1, the other to a particular propellant of low burning rate ($\text{J}{}_{\mathrm{b}}\text{? 0.2}\text{to}\text{0.4}\text{mm s}{}^{\mathrm{?1}}$) D₁ < 1. The stability conditions are analysed as well as the pressure amplitute and phase as a function of the nozzle throat modulation frequency.Still in linearized theory, the complete solutions of the problem are presented, using a method of numerical resolution.     

A small parameter method in problems of maneuvering space vehicles with aerodynamic efficiency

The motion of a satellite with aerodynamic efficiency along a low near-circular orbit is considered in the paper. The controls of bank angle γ and lift coefficient C_y are used as control functions. The introduction of a small parameter ($\text{? = (}\text{?}{}_0\text{· S · g}{}_0\text{2G}\text{)}$) makes it possible to integrate an adjoint system of equations and to obtain an approximate solution to the complete problem in the class of piecewise-constant control functions. Maximum values for the coordinates of heading angle η and lateral derivation from the plane of a reference orbit ?, which are connected with orbit plane angle by the relation cos i = cos ? · cos h, are used as criteria of maneuvering capability for a satellite with aerodynamic efficiency. Optimal programs for bank angle and incidence variation are derived and the influence of lift-to-drag ratio on the vehicle maneuvering capabilities has been estimated.It is shown that the process of the optimal motion is a special kind of gravitational skipping similar to the Keplerian motion but with continuous descent.     

Order of magnitude analysis of convective effects in dendritic solidification

The different types of convective phenomena which may occur during the dendritic solidification of metallic alloys are discussed from an order of magnitude analysis. Bulk thermal convection and/or interdendritic solutal convection have to be considered according to the values of the experimental data. Scaling laws for the solute boundary layer resulting from bulk thermal convection have already been derived. It is shown here that the interdendritic flow depends on a solutal Grashof number Gr based on the horizontal density gradient and a characteristic length L_s which is of the order of the liquid channels width. For Gr < 1, which is generally verified in practical cases, the interdendritic flow velocity U_r is proportional to the Grashof number. This a priori law compares favorably with the results of horizontal solidification experiments where the mean interdendritic flow velocity has been estimated from the resulting measured macrosegregation. In these experiments, as well as for most horizontal dendritic solidifications of metallic alloys at 1 g, the ratio $\text{U}{}_{\mathrm{r}}\text{R}$ (R is the growth rate) is of order one. In order to cancel the interdendritic flow effects, this ratio has to be lowered by one order of magnitude. According to our analysis, this can be obtained by performing the experiments either at a slightly reduced g level (～10^?1 g), or at 1 g in a vertical stable configuration with a sufficiently low residual horizontal thermal gradient.     

Optimal aeroassisted return from high earth orbit with plane change

This paper gives a complete analysis of the problem of aeroassisted return from a high Earth orbit to a low Earth orbit with plane change. A discussion of pure propulsive maneuver leads to the necessary change for improvement of the fuel consumption by inserting in the middle of the trajectory an atmospheric phase to obtain all or part of the required plane change. The variational problem is reduced to a parametric optimization problem by using the known results in optimal impulsive transfer and solving the atmospheric turning problem for storage and use in the optimization process. The coupling effect between space maneuver and atmospheric maneuver is discussed. Depending on the values of the plane change i, the ratios of the radii, $\text{n =}\text{r}{}_1\text{r}{}_2$ between the orbits and $\text{a =}\text{r}{}_2\text{R}$ between the low orbit and the atmosphere, and the maximum lift-to-drag ratio E1 of the vehicle, the optimal maneuver can be pure propulsive or aeroassisted. For aeroassisted maneuver, the optimal mode can be parabolic, which requires only drag capability of the vehicle, or elliptic. In the elliptic mode, it can be by one-impulse for deorbit and one or two-impulse in postatmospheric flight, or by two-impulse for deorbit with only one impulse for final circularization. It is shown that whenever an impulse is applied, a plane change is made. The necessary conditions for the optimal split of the plane changes are derived and mechanized in a program routine for obtaining the solution.     

Small cryogenic propulsion unit for upper stage application

The feasibility study was conducted to use the cryogenic propulsion system for the third stage of the future H-1 vehicle. While the $\text{LO}{}_2\text{LH}{}_2$ third-stage mass fraction is less than the current solid propellant third stage, the 50% higher Isp results in a significantly higher payload. Two basic configurations of the new propulsion system were proposed: one pressure-fed system and two pump-fed systems. The first is a pressure-fed system providing a 700 kg thrust at an Isp of 441 sec with restart capability. The second is a pump-fed system, operating on an expander cycle principle. A midget turbopump with a 90 000 rpm shaft speed feeds the thrust chamber which delivers 1 ton of thrust at an Isp of 471 sec. The third proposed system is also a pump-fed design using a unique expander bleed cycle, and delivers a 1 ton thrust at an Isp of 470 sec with a turbopump speed of 80 000 rpm. The results of engine testing predict the performance feasibility of respective propulsion system designs.     

Theoretical study about the Marangoni convection in a liquid column in zero-gravity

The thermal Marangoni effect on the surface of a liquid bridge induces a convection inside the liquid. For an imposed arbitrary periodic axial circumferential temperature distribution on the liquid surface the velocity distributions in radial-, angular- and axial direction are determined theoretically by solving the linearized Navier-Stokes equations. Of particular interest is the effect of the viscosity parameter $\text{v}\text{a}{}^2$ and axial wave length to diameter ratio $\text{l}\text{a}$. It was found that the increase of viscosity decreases the magnitude of the velocity distributions and that for small axial wave length to diameter ratios the radial- and axial velocities exhibit peak values close to the free surface of the liquid. This is in a less pronounced way also true for the angular velocity, which shows for increasing moderate values $\text{l}\text{a}\text{(0 ≤}\text{l}\text{a}\text{≤ 2}$) a strong increase in magnitude and for larger axial wavelength a decrease again. For increasing axial wavelength the peak value of the radial- and axial velocity shifts towards the center of the liquid bridge, of which for a further increase a decrease of the magnitude appears.     

Effects of stand-off bursts on rubble-pile targets: Evaluation of a hazardous asteroid mitigation strategy

We explore the aftereffects of stand-off burst mitigation on kilometer-scale rubble pile asteroids. We use a simple model of X-ray energy deposition to calculate the impulse transferred to the target, in particular to burst-facing blocks on the target surface. The impulse allows us to estimate an initial velocity field for the blocks on the outer side of the target facing the burst. We model the dynamics using an N-body polyhedron program built on the Open Dynamics Engine, a “physics engine” that integrates the dynamical equations for objects of general shapes and includes collision detection, friction, and dissipation.We tested several different models for target objects: rubble piles with different mass distributions, a “brick-pile” made of closely fitting blocks and zero void space, and a non-spherical “contact binary” rubble pile. Objects were bound together by self-gravity and friction/inelastic restitution with no other cohesive forces. Our fiducial cases involved objects of m=3.5×10¹² kg (corresponding to a radius of 0.7 km for the bulk object), an X-ray yield of 1 megaton, and stand-off burst distances of R=0.8–2.5 km from the target center of mass.Kilometer-scale rubble piles are robust to stand-off bursts of a yield ($Y\sim 1$ megaton) that would be sufficient to provide an effective velocity change ($\mathrm{\Delta }v\sim 0.05\mathrm{m}{\mathrm{s}}^{-1}$). Disaggregation involving some tens of percent of the target mass happens immediately after the impulse; the bulk of the object re-accretes on a few gravitational timescales, and the final deflected target contains over 95% (typically, 98–99%) of the original mass. Off-center components of the mitigation impulse and the target mass distribution cause a small amount of induced spin and off-axis components of velocity change. The off-axis velocity component amounts to an angular deviation of $\sim$ 0.05–0.1 radians from the nominal impulse vector, which may be important for mitigation planning.