{"id":1753,"date":"2023-05-29T06:03:57","date_gmt":"2023-05-29T03:03:57","guid":{"rendered":"https:\/\/fti.dp.ua\/conf\/?p=1753"},"modified":"2023-05-29T06:03:59","modified_gmt":"2023-05-29T03:03:59","slug":"05291-0602","status":"publish","type":"post","link":"https:\/\/fti.dp.ua\/conf\/2023\/05291-0602\/","title":{"rendered":"Statement of the problem of designing a liquid rocket engine dual bell nozzle of the maximum thrust using the direct method of the calculus of variations"},"content":{"rendered":"\n<h1 class=\"wp-block-heading citation_title\">Statement of the problem of designing a liquid rocket engine dual bell nozzle of the maximum thrust using the direct method of the calculus of variations<\/h1>\n\n\n\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<h5 class=\"wp-block-heading citation_author\">Ivan Dubrovskyi<\/h5>\n\n\n\n<p class=\"citation_author_url\"><em>ORCID: <a href=\"https:\/\/orcid.org\/0000-0002-0707-0074\" target=\"_blank\" rel=\"noopener\" title=\"\">https:\/\/orcid.org\/0000-0002-0707-0074<\/a><\/em><\/p>\n\n\n\n<p><em><em>Oles Honchar Dnipro National University, Dnipro<\/em><\/em><\/p>\n<\/div><\/div>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<h5 class=\"wp-block-heading citation_author\">Valerii Bucharskyi<\/h5>\n\n\n\n<p class=\"citation_author_url\"><em>ORCID: <a href=\"https:\/\/orcid.org\/0000-0002-8245-5652\" target=\"_blank\" rel=\"noopener\" title=\"\">https:\/\/orcid.org\/0000-0002-8245-5652<\/a><\/em><\/p>\n\n\n\n<p><em><em><em>Oles Honchar Dnipro National University, Dnipro<\/em><\/em><\/em><\/p>\n<\/div><\/div>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">Introduction<\/h3>\n\n\n\n<p>A rocket engine with a standard single-section profiled nozzle provides maximum thrust only at a certain value of atmospheric pressure [1]. To expand the range of atmospheric pressure values at which the optimum engine operation mode is achieved, a dual bell nozzle can be used [2]. A typical diagram of such a nozzle is shown in Fig. 1. Its main elements are three cross-sections <strong>a<\/strong>, <strong>b<\/strong>, <strong>c<\/strong> and two profiled section elements <strong>a-b<\/strong>, <strong>b-c<\/strong> bounded by them. The fig. 1 shows the general case of the contour of such a nozzle, for which a corner in the generatrix is located at points <strong>a<\/strong> and <strong>b<\/strong>. In addition, the section <strong>a<\/strong> coincides with the critical section of the engine chamber.<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-dual-bell-nozzle-contour.png\" alt=\"\" class=\"wp-image-1754\" width=\"396\" height=\"203\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-dual-bell-nozzle-contour.png 792w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-dual-bell-nozzle-contour-300x154.png 300w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-dual-bell-nozzle-contour-768x394.png 768w\" sizes=\"auto, (max-width: 396px) 100vw, 396px\" \/><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 1 \u2013 Dual bell nozzle contour<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Let\u2019s consider the principle of operation of such a nozzle. When operating at low altitude, with a high value of atmospheric pressure <em>p<sub>e1<\/sub><\/em>, the thrust is generated only by the first profiled section <strong>a<\/strong>&#8211;<strong>b<\/strong>. With an increase in flight altitude and a corresponding drop in atmospheric pressure to <em>p<sub>e2<\/sub> <\/em>(<em>p<sub>e1 <\/sub><\/em>&gt;<em> p<sub>e2<\/sub><\/em>), the flow of combustion products expands, begins to flow around the second profiled section and create additional thrust. Thus, an engine with a dual bell nozzle will generate more average flight time thrust than an engine with a single nozzle designed to operate at one of the atmospheric pressures <em>p<sub>e1 <\/sub><\/em>or <em>p<sub>e2<\/sub><\/em>.<\/p>\n\n\n\n<p>Usually, the method of characteristics [3] is used for profiling a dual bell nozzle. However, it has the following disadvantages:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>the inability to profile the maximum thrust nozzle with explicit restrictions on its dimensions, weight, etc.;<\/li>\n\n\n\n<li>the requirement for the absence of shock waves in the flow of combustion products inside the nozzle.<\/li>\n<\/ul>\n\n\n\n<p>As an alternative to the method of characteristics, this paper proposes to use the direct method of calculus of variations [4] for profiling a dual bell nozzle.<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">The aim and tasks of the study<\/h3>\n\n\n\n<p>The aim of the work is to obtain the thrust functional for its further maximization using the direct method of the calculus of variations. To achieve this goal, it is necessary to solve the following tasks:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>determine the initial data necessary to derive the thrust functional and to implement its further maximization;<\/li>\n\n\n\n<li>derive the expression of the thrust functional.<\/li>\n<\/ul>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">Initial data<\/h3>\n\n\n\n<p>The initial data for the problem of profiling the axisymmetric dual bell nozzle by the direct method of the calculus of variations are:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>known contour of the subsonic part of the rocket engine chamber;<\/li>\n\n\n\n<li>mathematical model of combustion products used to describe their flow inside the nozzle;<\/li>\n\n\n\n<li>parameters of combustion products in the critical section: pressure, density and velocity vector;<\/li>\n\n\n\n<li>atmospheric pressure values at both engine operating modes <em>p<sub>e1<\/sub><\/em> and <em>p<sub>e2<\/sub><\/em>,(<em>p<sub>e1 <\/sub><\/em>&gt;<em> p<sub>e2<\/sub><\/em>).<\/li>\n<\/ul>\n\n\n\n<p>As the main mathematical model describing the flow of combustion products in the liquid rocket engine chamber, a model of an inviscid compressible ideal gas of constant chemical composition was chosen, consisting of a system of stationary Euler equations (1), which was closed by the Mendeleev-Clapeyron equation of state [5]. System (1) is written in differential form in a cylindrical coordinate system under the assumption of axial symmetry of the flow of combustion products inside the LRE chamber.<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"83\" height=\"33\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image003.png\" alt=\"\" class=\"wp-image-1755\"\/><\/figure><\/div>\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"274\" height=\"17\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image004-1.png\" alt=\"\" class=\"wp-image-1756\"\/><\/figure><\/div>\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"281\" height=\"17\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image005-1.png\" alt=\"\" class=\"wp-image-1757\"\/><\/figure><\/div>\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"102\" height=\"17\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image006.png\" alt=\"\" class=\"wp-image-1758\"\/><\/figure><\/div>\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"296\" height=\"35\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image007-1.png\" alt=\"\" class=\"wp-image-1759\"\/><\/figure><\/div>\n\n\n<p>where <strong><em>F<\/em> \u2013 <\/strong>axial flux vector;<\/p>\n\n\n\n<p><strong><em>G<\/em> \u2013 <\/strong>radial flux vector;<\/p>\n\n\n\n<p><strong><em>S<\/em><\/strong><strong> \u2013 <\/strong>source term;<\/p>\n\n\n\n<p><em>R \u2013 <\/em>radius, m;<\/p>\n\n\n\n<p><em>\u03c1<\/em> \u2013 density, kg\/m<sup>3<\/sup>;<\/p>\n\n\n\n<p><em>p \u2013<\/em> pressure, Pa;<\/p>\n\n\n\n<p><em>u<\/em> \u2013 axial velocity, m\/s;<\/p>\n\n\n\n<p><em>v<\/em> \u2013 radial velocity, m\/s;<\/p>\n\n\n\n<p><em>E<\/em> \u2013 total energy, J\/kg;<\/p>\n\n\n\n<p><em>k<\/em> \u2013 heat capacity ratio.<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">Thrust functional derivation<\/h3>\n\n\n\n<p>The thrust of an arbitrary axisymmetric nozzle can be defined as the resultant of the pressure forces applied to the side surface of the nozzle [1]:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"266\" height=\"54\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image008.png\" alt=\"\" class=\"wp-image-1760\"\/><\/figure><\/div>\n\n\n<p>where <em>P<\/em> \u2013 thrust, H;<\/p>\n\n\n\n<p><em>F \u2013 <\/em>side surface area of a nozzle, m<sup>2<\/sup>;<\/p>\n\n\n\n<p><em>p<sub>in<\/sub> \u2013<\/em> combustion products pressure inside the nozzle, Pa;<\/p>\n\n\n\n<p><em>p<sub>e<\/sub> \u2013<\/em> atmospheric pressure at the given altitude, Pa.<\/p>\n\n\n\n<p>Let us select on the nozzle contour an annular element with a width <em>dx<\/em>, a height <em>dR<\/em> and located at a distance <em>R<\/em> from the axis of symmetry of the nozzle (fig. 2).<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-2-\u2013-annual-element-of-a-nozzle-contour.png\" alt=\"\" class=\"wp-image-1761\" width=\"330\" height=\"296\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-2-\u2013-annual-element-of-a-nozzle-contour.png 660w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-2-\u2013-annual-element-of-a-nozzle-contour-300x269.png 300w\" sizes=\"auto, (max-width: 330px) 100vw, 330px\" \/><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 2 \u2013 Annual element of a nozzle contour<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Its side surface area equals <em>dF<\/em> = <em>2\u03c0Rds<\/em>. Using this equality, the integral (2) can be represented as:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"427\" height=\"54\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image010.png\" alt=\"\" class=\"wp-image-1762\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image010.png 427w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image010-300x38.png 300w\" sizes=\"auto, (max-width: 427px) 100vw, 427px\" \/><\/figure><\/div>\n\n\n<p>Due to the axial symmetry of the nozzle, all the radial components of the thrust will be mutually balanced, therefore, in what follows, we will take into account only its axial component. In view of this, expression (3) is transformed to the following form:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"325\" height=\"54\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image011.png\" alt=\"\" class=\"wp-image-1763\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image011.png 325w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image011-300x50.png 300w\" sizes=\"auto, (max-width: 325px) 100vw, 325px\" \/><\/figure><\/div>\n\n\n<p>From the fig. 2 one can deduce that <em>cos\u03c6 = cos(90-\u03b8) = sin\u03b8 = dR\/ds<\/em>, and since <em>R`<\/em> = <em>dR\/dx<\/em>, then rewrite (4) in this way:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"433\" height=\"109\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image012.png\" alt=\"\" class=\"wp-image-1764\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image012.png 433w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image012-300x76.png 300w\" sizes=\"auto, (max-width: 433px) 100vw, 433px\" \/><\/figure><\/div>\n\n\n<p>where x<sub>1<\/sub> \u2013 the coordinate of the beginning of the nozzle at the horizontal axis, m;<\/p>\n\n\n\n<p>x<sub>2 <\/sub>\u2013 the coordinate of the end of the nozzle at the horizontal axis, m;<\/p>\n\n\n\n<p><em>R <\/em>=<em> R(x)<\/em> \u2013 nozzle contour that is being designed.<\/p>\n\n\n\n<p>Let\u2019s transform the integrand from (5):<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"406\" height=\"55\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image013.png\" alt=\"\" class=\"wp-image-1765\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image013.png 406w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image013-300x41.png 300w\" sizes=\"auto, (max-width: 406px) 100vw, 406px\" \/><\/figure><\/div>\n\n\n<p>Since the atmospheric pressure does not depend on the integration variable, we take it out of the brackets and integrate the second term into (6):<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"413\" height=\"56\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image014.png\" alt=\"\" class=\"wp-image-1766\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image014.png 413w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image014-300x41.png 300w\" sizes=\"auto, (max-width: 413px) 100vw, 413px\" \/><\/figure><\/div>\n\n\n<p>Let\u2019s substitute (6) and (7) into (5). As a result, we get the final expression for determining the thrust of an arbitrary nozzle:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"385\" height=\"55\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image015.png\" alt=\"\" class=\"wp-image-1767\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image015.png 385w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image015-300x43.png 300w\" sizes=\"auto, (max-width: 385px) 100vw, 385px\" \/><\/figure><\/div>\n\n\n<p>In (8), there are clearly two thrust components: an internal one, which depends on the pressure of the combustion products, and an external one, which depends only on the pressure of the atmosphere.<\/p>\n\n\n\n<p>Now use (8) to calculate the thrust of the dual bell nozzle shown in<br>fig. 1. In the first operation mode with high atmospheric pressure <em>p<sub>e1<\/sub><\/em>, only the first section of the nozzle <strong>a<\/strong>&#8211;<strong>b<\/strong> will create thrust, in the second mode &#8211; both sections <strong>a<\/strong>&#8211;<strong>b<\/strong> and <strong>b<\/strong>&#8211;<strong>c<\/strong>. Therefore, (8) for the second and the first mode will take the form (9):<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"402\" height=\"55\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image016.png\" alt=\"\" class=\"wp-image-1768\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image016.png 402w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image016-300x41.png 300w\" sizes=\"auto, (max-width: 402px) 100vw, 402px\" \/><\/figure><\/div>\n\n\n<p>Here the number <em>m<\/em> can get values 1 and 2 depending on the mode, <em>R<sub>1<\/sub> = R<sub>b<\/sub>, R<sub>2<\/sub> = R<sub>c<\/sub><\/em>.<\/p>\n\n\n\n<p>Let&#8217;s combine the separate thrust equations taken from (9) to obtain the total thrust equation (10) for the both modes simultaneously. To do this, we introduce the coefficient <em>\u03b1<\/em> equal to 0 in the first mode and 1 in the second. Then the expression for the total thrust of the dual bell nozzle takes the form:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"336\" height=\"18\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image017.png\" alt=\"\" class=\"wp-image-1769\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image017.png 336w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image017-300x16.png 300w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image017-330x18.png 330w\" sizes=\"auto, (max-width: 336px) 100vw, 336px\" \/><\/figure><\/div>\n\n\n<p>where <em>\u03b1<\/em> \u2013 the coefficient, that equals to 0 in the first operation mode and 1 \u2013 in the second;<\/p>\n\n\n\n<p><em>P<sub>\u03a3<\/sub><\/em> \u2013 dual bell nozzle total thrust, H.<\/p>\n\n\n\n<p>We will obtain the final form of the thrust functional by adding to (10) the restrictions imposed by the chosen mathematical model of combustion products (1). To do this, we use the Lagrange multipliers [6] &#8211; we introduce four unknown functions <em>\u03bb<sub>1<\/sub><\/em>, <em>\u03bb<sub>2<\/sub><\/em>, <em>\u03bb<sub>3<\/sub><\/em>, <em>\u03bb<sub>4<\/sub><\/em>, that depend on <em>x<\/em> and <em>R<\/em>, multiply them by the corresponding equations from the system (1) and sum them up:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"303\" height=\"49\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image018.png\" alt=\"\" class=\"wp-image-1770\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image018.png 303w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image018-300x49.png 300w\" sizes=\"auto, (max-width: 303px) 100vw, 303px\" \/><\/figure><\/div>\n\n\n<p>To add (11) to (10), we integrate (11) over the region <strong><em>K<\/em><\/strong> bounded from below by the axis of symmetry of the nozzle and from above by its contour. After that, the desired thrust functional will take the final form:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"363\" height=\"54\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image019.png\" alt=\"\" class=\"wp-image-1771\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image019.png 363w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/image019-300x45.png 300w\" sizes=\"auto, (max-width: 363px) 100vw, 363px\" \/><\/figure><\/div>\n\n\n<p>The values of functional (12) can only be obtained as a result of numerical simulation of the flow of combustion products in the chamber of a liquid rocket engine. To implement this, various numerical methods can be used [7, 8].<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">Conclusions<\/h3>\n\n\n\n<p>In the work, the initial data for the problem of profiling a dual bell nozzle were determined. An expression for the thrust functional was also obtained for such a nozzle. This will allow in the future to solve the problem of its maximization using the direct method of calculus of variations, the result of which will be the dual bell contour of the, producing the maximum thrust.<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center citation_references\">References<\/h3>\n\n\n\n<ol type=\"1\" class=\"citation_references wp-block-list\">\n<li>Biblarz&nbsp;O., Sutton&nbsp;G.&nbsp;P. Rocket propulsion elements. Wiley &amp; Sons, Incorporated, John, 2016. 792&nbsp;p.<\/li>\n\n\n\n<li>Future space-transport-system components under high thermal and mechanical loads\u00a0\/ ed. by N.\u00a0A.\u00a0Adams et al. Cham\u00a0: Springer International Publishing, 2021. URL:\u00a0https:\/\/doi.org\/10.1007\/978-3-030-53847-7\u00a0(date of access: 25.05.2023).<\/li>\n\n\n\n<li>Zucrow&nbsp;M.&nbsp;J., Hoffman&nbsp;J.&nbsp;D. Gas dynamics, vol. 2: multidimensional flow. John Wiley &amp; Sons Inc, 1977. 490&nbsp;p.<\/li>\n\n\n\n<li>Giusti&nbsp;E. Direct methods in the calculus of variations. River Edge, NJ&nbsp;: World Scientific, 2003. 403&nbsp;p.<\/li>\n\n\n\n<li>Chorin&nbsp;A.&nbsp;J., Marsden&nbsp;J.&nbsp;E. Mathematical introduction to fluid mechanics. Springer London, Limited, 2013. 172 p.<\/li>\n\n\n\n<li>Kot&nbsp;M. A first course in the calculus of variations. Providence, Rhode Island&nbsp;: American Mathematical Society, 2014. 298&nbsp;p.<\/li>\n\n\n\n<li>Dubrovskyi I., Bucharskyi V. Development of a method of extended cells for the formulation of boundary conditions in numerical integration of gas dynamics equations in the domains of a curvilinear shape. Eastern-European Journal of Enterprise Technologies. 2020. Vol. 5, No 7. P. 72-84.<\/li>\n\n\n\n<li>Ferziger J. H., Milovan P., Street R. L. Computational Methods for Fluid Dynamics. Germany: Springer-Verlag. 606 p<\/li>\n<\/ol>\n\n\n\n<div style=\"height:3em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity is-style-default\"\/>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<div class=\"wp-block-group is-vertical is-content-justification-right is-layout-flex wp-container-core-group-is-layout-b6c475e2 wp-block-group-is-layout-flex\">\n<div class=\"wp-block-group is-content-justification-right is-nowrap is-layout-flex wp-container-core-group-is-layout-fd526d70 wp-block-group-is-layout-flex\"><div class=\"taxonomy-post_tag wp-block-post-terms\"><a href=\"https:\/\/fti.dp.ua\/conf\/tag\/cims-2023\/\" rel=\"tag\">CIMS 2023<\/a><\/div>\n\n<div class=\"wp-block-post-date\"><time datetime=\"2023-05-29T06:03:57+03:00\">May 29, 2023<\/time><\/div><\/div>\n\n\n<div class=\"taxonomy-category wp-block-post-terms\"><a href=\"https:\/\/fti.dp.ua\/conf\/session\/aerospace\/\" rel=\"tag\">Aerospace vehicles<\/a><\/div>\n\n\n<p><em>URI: <a href=\"\/conf\/2023\/05291-0602\/\" title=\"\">https:\/\/fti.dp.ua\/conf\/2023\/05291-0602\/<\/a><\/em><\/p>\n<\/div>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity is-style-default\"\/>\n","protected":false},"excerpt":{"rendered":"<p>Statement of the problem of designing a liquid rocket engine dual bell nozzle of the maximum thrust using the direct method of the calculus of variations Ivan Dubrovskyi ORCID: https:\/\/orcid.org\/0000-0002-0707-0074 Oles Honchar Dnipro National University, Dnipro Valerii Bucharskyi ORCID: https:\/\/orcid.org\/0000-0002-8245-5652 Oles Honchar Dnipro National University, Dnipro Introduction A rocket engine with a standard single-section profiled nozzle provides maximum thrust only at a certain value of atmospheric pressure [1]. To expand the range of atmospheric pressure values at which the optimum &hellip; <\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[15],"tags":[23],"class_list":["post-1753","post","type-post","status-publish","format-standard","hentry","category-aerospace","tag-cims-2023"],"_links":{"self":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts\/1753","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/comments?post=1753"}],"version-history":[{"count":1,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts\/1753\/revisions"}],"predecessor-version":[{"id":1772,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts\/1753\/revisions\/1772"}],"wp:attachment":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/media?parent=1753"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/categories?post=1753"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/tags?post=1753"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}