Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID INJECTION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

Destabilization of binary mixing layer in supercritical conditions

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Geochemical Changes in Response to CO 2 Injection in a CO 2 -EOR Complex in Northern Michigan (Volume II.1)

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the geochemical monitoring program under MRCSP to use stable and radiogenic isotope geochemistry in concert with analysis of general geochemical parameters of fluids and gases and analysis of core samples to determine geochemical processes occurring in the reef structure because of CO 2 injection. Specifically, brine and gas samples were collected and analyzed to determine changes occurring between reefs prior to and following CO 2 injection. The analytical results for general geochemical parameters were modeled with chemical equilibrium models to determine if the injection of CO 2 resulted in the mineral dissolution or precipitation. Finally, core samples were collected and analyzed to determine if there was evidence of dissolution features or mineral precipitation.

01 COAL, LIGNITE, AND PEAT↗

Experimental study of crossover from capillary to viscous fingering for supercritical CO 2 - water displacement in a homogeneous pore network

Carbon sequestration in saline aquifers involves displacing resident brine from the pore space by supercritical CO 2 (scCO 2 ). The displacement process is considered unstable due to the unfavorable viscosity ratio (logM < 0). The unstable mechanisms that affect scCO 2 - water displacement under reservoir conditions (i.e., 41 °C, 9 MPa) were investigated in a homogeneous micromodel. A wide range of injection rates (logCa = -7.61~-4.73) was studied in two sets of experiments: discontinuous-rate injection, where the micromodel was first cleaned and saturated with water before each injection rate was imposed, and continuous-rate injection, where the rate was increased after quasi-steady conditions were reached for a certain rate. For the discontinuous-rate experiments, capillary fingering and viscous fingering are the dominant mechanisms for low (logCa <= -6.61) and high injection rates (logCa >= -5.21), respectively. Crossover from capillary to viscous fingering was observed for logCa = -5.91~-5.21, resulting in a large decrease in scCO 2 saturation. The discontinuous-rate experimental results confirmed the decrease in nonwetting fluid saturation during crossover from capillary to viscous fingering predicted by numerical simulations by Lenormand et al. (1988).1 Capillary fingering was the only mechanism that dominates all injection rates in the continuous-rate experiment, and resulted in monotonic increase in scCO 2 saturation.

Environmental Molecular Sciences Laboratory↗

Numerical estimations of storage efficiency for the prospective CO 2 storage resource of shales

Here, hydrocarbon-bearing shale formations might be an attractive geologic reservoir for permanent carbon dioxide (CO 2 ) storage. Shale formations applicable for storage require previous hydrocarbon production to deplete a reservoir, hydraulic fracturing to provide a highly-permeable stimulated zone, sufficient depths (generally > 800 m or 2600 feet) to maintain CO 2 in a supercritical state, and an overlying seal to prevent CO 2 migration to underground sources of drinking water and into the atmosphere. CO 2 is stored in shale as a free phase within fractures and matrix pores, and as a sorbed phase on organic matter and clay minerals. Recently in our previous work, we presented a screening-level assessment methodology for CO 2 storage in shales using a volumetric approach. The approach deals with reduction of CO 2 storage through estimations of efficiency factors based on petrophysical properties (i.e., bulk volume, porosity, sorption, etc.) and their limitations on fluid transport. Here, we conducted numerical simulations using the FRACGEN/NFFLOW simulator to study the CO 2 injection into a depleted hydro-fractured shale reservoir to estimate storage efficiencies using a range of reservoir parameters and injection scenarios. Specifically, the ranges for two efficiency factors, E Φ and E S , measure the effectiveness of free and sorbed CO 2 storage. These efficiency factors were estimated to have P 10 to P 90 probability ranges of 0.15 to 0.36 for and 0.11 to 0.24 for E S , reported after 60 years of CO 2 injection.

58 GEOSCIENCES↗

PARTICLE ACCELERATION AT A FLARE TERMINATION SHOCK: EFFECT OF LARGE-SCALE MAGNETIC TURBULENCE

We investigate the acceleration of charged particles (both electrons and protons) at collisionless shocks predicted to exist in the vicinity of solar flares. The existence of standing termination shocks has been examined by flare models and numerical simulations. We study electron energization by numerically integrating the equations of motion of a large number of test-particle electrons in the time-dependent two-dimensional electric and magnetic fields generated from hybrid simulations (kinetic ions and fluid electron) using parameters typical of the solar flare plasma environment. The shock is produced by injecting plasma flow toward a rigid piston. Large-scale magnetic fluctuations-known to exist in plasmas and known to have important effects on the nonthermal electron acceleration at shocks-are also included in our simulations. For the parameters characteristic of the flaring region, our calculations suggest that the termination shock formed in the reconnection outflow region (above post-flare loops) could accelerate electrons to a kinetic energy of a few MeV within 100 ion cyclotron periods, which is of the order of a millisecond. Given a sufficient turbulence amplitude level ({delta}B{sup 2}/B 2{sub 0} {approx} 0.3), about 10% of thermal test-particle electrons are accelerated to more than 15 keV. We find that protons are also accelerated, but not to as high energy in the available time and the energy spectra are considerably steeper than that of the electrons for the parameters used in our simulations. Our results are qualitatively consistent with the observed hard X-ray emissions in solar flares.

79 ASTRONOMY AND ASTROPHYSICS↗

Neoclassical offset toroidal velocity and auxiliary ion heating in tokamaks

In conditions of ideal axisymmetry, for a magnetized plasma in a generic bounded domain, necessarily toroidal, the uniform absorption of external energy (e.g., RF or any isotropic auxiliary heating) cannot give rise to net forces or torques. Experimental evidence on contemporary tokamaks shows that the near central absorption of RF heating power (ICH and ECH) and current drive in presence of MHD activity drives a bulk plasma rotation in the co-I{sub p} direction, opposite to the initial one. Also the appearance of classical or neoclassical tearing modes provides a nonlinear magnetic braking that tends to clamp the rotation profile at the q-rational surfaces. The physical origin of the torque associated with P{sub RF} absorption could be due the effects of asymmetry in the equilibrium configuration or in power deposition, but here we point out also an effect of the response of the so-called neoclassical offset velocity to the power dependent heat flow increment. The neoclassical toroidal viscosity due to internal magnetic kink or tearing modes tends to relax the plasma rotation to this asymptotic speed, which in absence of auxiliary heating is of the order of the ion diamagnetic velocity. It can be shown by kinetic and fluid calculations, that the absorption of auxiliary power by ions modifies this offset proportionally to the injected power thereby forcing the plasma rotation in a direction opposite to the initial, to large values. The problem is discussed in the frame of the theoretical models of neoclassical toroidal viscosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a two-zone zero-dimensional mixing model for actively fueled prechamber engines

In this numerical work, a zero-dimensional, two-zone mixing model was developed and validated to capture the mixture preparation process of an actively fueled prechamber engine. The model was developed for ethanol-gasoline fuel blends but is applicable to various fuels of interest relevant to prechamber applications. The model was exercised to demonstrate the effects of boundary conditions, fuel composition, and residual content on the mixture stoichiometry predictions. The 0D model showed reasonable agreement in prechamber equivalence ratio predictions with higher fidelity computational fluid dynamics modeling, but the accuracy of model was constrained to the governing assumptions of the injection process. The modeling assumptions surrounding the non-dimensional injection of mass into the prechamber were scrutinized to identify key areas of model improvement. The results indicate heat transfer effects play a crucial role in the evaporation of the fuel within the confined auxiliary chamber, especially when considering fuels with high latent heat of vaporization. Additionally, the presence of fuel films that collect on the prechamber surfaces complicated the modeled evaporative cooling loses as the modeling suggests the films evaporate primarily by means of conduction from the boundary. Furthermore, the developed model serves as a simplified framework to generate calibration maps of prechamber fueling or run real time as an open loop control strategy such that desired prechamber mixture proportions are effectively targeted.

Zero-dimensional↗

EGS Collab Experiment 1: TOUGH2-CSM Simulation of Embedded Natural Fractures and Chemical Tracer Transport and Sorption

The EGS Collab SIGMA-V project is a multi-lab and university collaborative research project that is being undertaken at the Sanford Underground Research Facility (SURF) in South Dakota. The project consists of studying stimulation, fluid-flow, and heat transfer processes at a scale of 10-20 m, which is readily amenable to detailed characterization and monitoring. One objective of the project is to establish circulation from injector to producer by hydraulically fracturing the injector. Data generated during these experiments is to be compared with predictions from coupled thermal, hydrological, mechanical, and chemical simulators. One such a simulator, TOUGH2-CSM, has been enhanced in order to simulate EGS Collab SIGMA-V project experiments. These modifications include adding tracers, the capability to model tracer sorption, and an embedded fracture formulation. A set of example problems validate our conservative tracer transport and sorption formulations. We then simulated tracer transport and thermal breakthrough for the first EGS Collab SIGMA-V experiment. This dataset includes the TOUGH2-CSM input and output files associated with the thermal and tracer simulations. A conference paper is included for additional context.

15 GEOTHERMAL ENERGY↗

Reactive Transport Modeling of Anthropogenic Carbon Mineralization in Stacked Columbia River Basalt Reservoirs

Numerical simulation of CO2 storage in basalts and related reactive lithologies requires modeling complex, coupled hydrologic and chemical processes, including multi-phase flow and transport, partitioning of CO2 into the aqueous phase, and chemical interactions with aqueous fluids and rock minerals. We conducted reactive transport simulations of the Wallula pilot-scale CO2 injection into the flow tops of the Grande Ronde Basalt using our PNNL STOMP-CO2 simulator with the ECKEChem reactive module. Our mineralization simulation of the ~1,000 tons of injected CO2 into the interflow zones was based on the hydrologic transport model we previously developed. For this work, the simulations considered geochemical reactions involving the basalt components, precipitates, formation brine, and injected CO2. In our benchmark case, carbonate minerals precipitated, resulting in ~20% of the CO2 being mineralized in 10 years. Increasing the reaction rate of a single primary mineral phase (clinopyroxene) by an order of magnitude resulted in a carbon mineralization reaction extent of ~90% over the same time interval. Based on these initial sensitivity analysis results, it is clear that a thorough understanding of primary mineral dissolution rates is required for accurately predicting long-term fate and transport of injected CO2 into basalt formations. Our reactive transport numerical simulations will be key components of commercial-scale CO2 storage operation permitting, de-risking, and optimization in mafic and ultramafic reservoirs.

Cao, Ruoshi↗

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗

Monitoring and Numerical Modeling of Shallow CO 2 Injection, Greene County, Missouri

The project titled Monitoring and Numerical Modeling of Shallow CO 2 Injection, Greene County, Missouri provided training for three graduate students in areas related to carbon capture and storage. Numerical modeling of CO 2 injection into the St. Francois aquifer at the Southwest Power Plant Site in Greene County, Missouri indicates that up to 4.1 x 10 5 metric tons of CO 2 per year could be injected for 30 years without exceeding a 3 MPa differential injection pressure. The injected CO 2 would remain sequestered below the top of the overlying caprock (St. Francois confining unit) for more than 1000 years. Geochemical modeling indicates that portions of the injected CO 2 will react rapidly with trace minerals in the aquifer to form various solid carbonate mineral phases. These minerals would store significant portions of injected CO 2 over geologic time scales. Finally, a GIS data base on the pore-fluid chemistry of the overlying aquifer system in Missouri, the Ozark aquifer, was compiled from many sources. This data base could become useful in monitoring for leakage from future CO 2 sequestration sites.

54 ENVIRONMENTAL SCIENCES↗

Time-Domain Nuclear Magnetic Resonance Determination of Wettability Alteration: Analysis for Low-Salinity Water

Wettability has been shown to influence oil recovery. This property has become central to low-salinity (LSW) and smart (SWF) water flooding recovery mechanisms research. The challenge lies in the fact that oil recovery results from the combined effects of solid-liquid and liquid-liquid interactions. This demands methods that allow an independent interpretation of wettability alteration contributions. The primary objective of this work is to assess changes in wettability through the application of Time-Domain Nuclear Magnetic Resonance (TD-NMR) T 2 distribution and diffusion coefficient, starting with a well-controlled porous system, that is, glass beads, and then a model rock (Berea), in the presence of one phase, either oil or brine exclusively. Subsequently, two-phase fluid saturation was tested. For the glass beads, dimethyldichlorosilane was used to induce a hydrophobic response, as confirmed by contact angle experiments on slides of the same material. Sodium sulfate was used for its known positive influence on oil recovery during LSW and SWF. In cases where alteration of surface properties was expected, a leftward shift of the average T 2 distribution curve modes, accompanied by a reduction on the diffusion coefficient during the aging process was observed. The results of this work confirm that fluid-solid interactions during LSW and SWF, namely a shift in wettability, take place after the injection of low-salinity water.

42 ENGINEERING↗

Quantification of sauter mean diameter in diesel sprays using scattering-absorption extinction measurements

Quantitative measurements of the primary breakup process in diesel sprays are lacking due to a range of experimental and diagnostic challenges, including: high droplet number density environments, very small characteristic drop size scales (~1-10 μm), and high characteristic velocities in the primary breakup region (~600 m/s). Due to these challenges, existing measurement techniques have failed to resolve a sufficient range of the temporal and spatial scales involved and much remains unknown about the primary atomization process in practical diesel sprays. To gain a better insight into this process, we have developed a joint visible and x-ray extinction measurement technique to quantify axial and radial distributions of the path-integrated Sauter Mean Diameter (SMD) and Liquid Volume Fraction (LVF) for diesel-like sprays. This technique enables measurement of the SMD in regions of moderate droplet number density, enabling construction of the temporal history of drop size development within practical diesel sprays. The experimental campaign was conducted jointly at the Georgia Institute of Technology and Argonne National Laboratory using the Engine Combustion Network “Spray D” injector. X-ray radiography liquid absorption measurements, conducted at the Advanced Photon Source at Argonne, quantify the liquid-fuel mass and volume distribution in the spray. Diffused back-illumination liquid scattering measurements were conducted at Georgia Tech to quantify the optical thickness throughout the spray. By application of Mie-scatter equations, the ratio of the absorption and scattering extinction measurements is demonstrated to yield solutions for the SMD. This work introduces the newly developed scattering-absorption measurement technique and highlights the important considerations that must be taken into account when jointly processing these measurements to extract the SMD. These considerations include co-alignment of measurements taken at different institutions, identification of viable regions where the measurement ratio can be accurately interpreted, and uncertainty analysis in the measurement ratio and resulting SMD. Because the measurement technique provides the spatial history of the SMD development, it is expected to be especially informative to the diesel spray modeling community. Results from this work will aid in understanding the effect of ambient densities and injection pressures on primary breakup and help assess the appropriateness of spray submodels for engine computational fluid dynamics codes.

33 ADVANCED PROPULSION SYSTEMS↗

Transport studies in high-performance field reversed configuration plasmas

A significant improvement of field reversed configuration (FRC) lifetime and plasma confinement times in the C-2 plasma, called High Performance FRC regime, has been observed with neutral beam injection (NBI), improved edge stability, and better wall conditioning [Binderbauer et al., Phys. Plasmas 22, 056110 (2015)]. A Quasi-1D (Q1D) fluid transport code has been developed and employed to carry out transport analysis of such C-2 plasma conditions. The Q1D code is coupled to a Monte-Carlo code to incorporate the effect of fast ions, due to NBI, on the background FRC plasma. Numerically, the Q1D transport behavior with enhanced transport coefficients (but with otherwise classical parametric dependencies) such as 5 times classical resistive diffusion, classical thermal ion conductivity, 20 times classical electron thermal conductivity, and classical fast ion behavior fit with the experimentally measured time evolution of the excluded flux radius, line-integrated density, and electron/ion temperature. The numerical study shows near sustainment of poloidal flux for nearly 1 ms in the presence of NBI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Injection nozzle for a turbomachine

A turbomachine includes a compressor, a combustor operatively connected to the compressor, an end cover mounted to the combustor, and an injection nozzle assembly operatively connected to the combustor. The injection nozzle assembly includes a first end portion that extends to a second end portion, and a plurality of tube elements provided at the second end portion. Each of the plurality of tube elements defining a fluid passage includes a body having a first end section that extends to a second end section. The second end section projects beyond the second end portion of the injection nozzle assembly.

42 ENGINEERING↗

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗

Development and Validation of Physics-Based Sub-Models of High-Pressure Supercritical Fuel Injection at Diesel Conditions Insert

The project sought to develop and validate accurate, physics-based, numerical submodels and implement them in computational fluid dynamics (CFD) software codes to enable accurate predictions of diesel sprays at pressures higher than the current systems. The combined experimental-computational effort has four main objectives: (1) acquire spatially and temporally resolved scalar measurements of fuel-oxidizer mixing in the near field of the fuel jet, (2) develop software modules to calculate thermo-physical properties of diesel surrogates and their mixtures with oxidizer of interest in real-gas operating regime, (3) integrate the real-fluid model into a commercial and an open-source CFD code to simulate fuel-air mixing in high-pressure diesel sprays, and (4) assess the robustness, accuracy, and uncertainty of the integrated CFD solvers in computing diesel sprays. Project accomplishments are consistent with the original goals of the project. Although we acquired spatially and temporally resolved scalar measurements, they were limited to refractive index difference and not to thermodynamic properties such as temperature and equivalence ratio. However, we successfully developed the theoretical background needed to acquire such measurements at diesel conditions and demonstrated the capability of the RSD for scalar measurements in a canonical turbulent jet flow replicating features of the diesel jet. Additional work will be needed to implement these procedures at high-pressure, high-temperature engine-like conditions of the CPFR. Mainly, we did not anticipate the need to find the liquid boundary and its statistical features or the state relationships between refractive index and thermodynamic properties, assuming the widely known Dale-Gladstone equation would suffice. Still, efforts devoted to these tasks were essential to fulfil the goals of the project and serves as the foundation on which many of the conclusions and key contributions of this project are based. Further assessment and application of the developed CFD codes is recommended to fully realize the benefits of this research.

33 ADVANCED PROPULSION SYSTEMS↗

System for injecting fuel in a gas turbine combustor

A combustion system uses a fuel nozzle with an inner wall having a fuel inlet in fluid communication with a fuel outlet in a fuel cartridge. The inner wall defines a mounting location for inserting the fuel cartridge. A pair of annular lip seals around the cartridge outer wall on both sides of the fuel outlet seals the fuel passage between the fuel inlet and the fuel outlet.

33 ADVANCED PROPULSION SYSTEMS↗