## Abstract

With the increasing need of micro‐air‐vehicles (MAVs) and advances in MAV technology, Coandă MAVs offer new promises and challenges. In this context, Coandă MAVs capabilities are analyzed. As a baseline, a mathematical model for a spherical Coandă MAV in hover and translatory motion is developed and analyzed from first physical principles. A computational fluid dynamic (CFD) simulations for a Coandă MAV generic model are carried out to assess the theoretical prediction and obtaining further physical insight on the Coandă MAV flow physics. The mathematical model and performance measures are developed to assess the capability of the semi‐spherical Coandă MAV in performing effective flight as an aerial robotic platform, as indicated by the relationships between the relevant parameters of the mathematical model of the Coandă MAV to its system of flight forces.

### Keywords

- aerial robot
- Coandă micro air vehicle (CMAV)
- Coandă MAV aerodynamics
- Coandă MAV flight dynamics
- unmanned micro air vehicle (UMAV)

## 1. Introduction

The utilization of Coandă effect and Coandă jet [1, 2] for flight vehicles have been carried out for a range of aircrafts and vehicles for decades and for future aircrafts [3–14], due to its flow control effectiveness. Other applications are potential in wind‐turbines technology [15–20], and the Coanda effect have also been found and utilized in many other engineering applications and health science, such as in the upper respiratory system [20–24]. The science of flow control can be considered to originate from Prandtl [25] seminal paper in 1904, which has been hailed as a breakthrough in the science of fluid mechanics by the introduction of the boundary‐layer theory. Prandtl’s theory can be accounted for the physics of the separation phenomena and the boundary layers control.

The Coandă effect is based on the property of a jet flow to attach itself to a curved surface as it tangentially blows along and to remain attached until it leaves the surface or reach a critical stability breakdown. Coandă jet is one of the flow control techniques for modifying and enhancing aerodynamic performance, stability and control. Coandă jet has also been applied in the development of novel aircrafts for short take‐off and landing. With the progress of many enabling technologies and need for disaster mitigation, security, and environmental conservation, the relatively new application of Coandă effect in micro‐air‐vehicles (MAVs) as robotic platforms is very challenging. Capitalizing on the basic fundamental principles, the aero‐ and flight dynamics of a baseline Coandă MAV configuration are analyzed and developed to formulate the relationships between various essential parameters to the equation of motion in hover, and translatory motion, which should give further insight on and can be further elaborated into other Coandă MAV characteristics and performance related to maneuver and stability and control.

Coandă micro‐air‐vehicles have been developed in various configurations. To meet the desired Coandă MAV operation, mission, and design requirements, one should establish the basic working relationships among various significant design parameters and variables to the generated aerodynamics forces. For this purpose, several approaches can be followed. The most logical first approach is the analytical tool based on basic fundamental principles. With the impressive progress of computational fluid dynamics (CFD) which has brought with it flow visualization capabilities; CFD utilization will assist the identification to the problem and assessment of the analytical approach. Furthermore, CFD visualization facilitates better insight and identification of salient details. Another important approach is the use of experimental tools, which is the foundation of scientific observation through meticulous experimental design and has intricate relationships with analytical approach. For the conceptual and prototype design stages, analytical, and CFD computational and visualization tools will expedite the efforts to that end.

The theoretical analysis, mathematical modeling and synthesis or Coandă micro‐air‐vehicles as an aerial robot platform will start with the identification and definition of a Coandă jet as a relatively thin and slightly viscous jet flowing over and adhere on a smooth curved surface. The main physical parameters are the angle of separation, θ; slot width, *b*; radius of curvature, *R*; Reynolds number, *Re*; the Coandă jet velocity, pressure differential p_{s}−p_{∞}, Reynolds number and the Coandă jet momentum coefficient. Various platforms will be discussed to define the merit of Coandă MAV as a robotic platform. Hovering semi‐spherical Coandă MAV will be utilized in the baseline aerodynamic analysis. The mathematical model of Coandă MAV to generate the equation of motion starts with first principles and the basic flight dynamic equation. CFD simulation and demonstration experiment are elaborated to assess the functionality of Coandă MAV aerial robot platform.

In the model development and analysis of a semi‐spherical Coandă MAV, a propulsion system is utilized to introduce a Coandă jet blanket which is deflected downward due to the curved upper surface of the vehicle, to provide both hover as well as cruise capabilities. Such Coandă jet has been utilized for circulation enhancement for fixed‐wing aircraft in forward flight and turbine blades movement [4–9, 12–20]. Thus, lift will be generated [27–35], which can be utilized for hover and later for cruise propulsion [10, 35, 36]. The principle of Coandă MAV lift generation as well as the equation of motion for its translatory motion will be derived and elaborated.

Some results from the authors’ previous analysis [31–34] provide the establishment of Coandă MAV spherical configuration baseline. The equations derived for the baseline configuration are significant for establishing preliminary design configuration concepts, and can be readily modified for other configurations for assessment, comparative and improvement purposes. CFD visualization studies will be instrumental in gaining in‐depth understanding of the analytical and physical model development for further analysis.

Based on the physical and mathematical model conceived, some basic results were obtained to describe the physics of the flow field associated with the Coandă effect jet sheets. Having developed the basic fluid‐dynamic analysis of the Coandă MAV lift generation, the equation of motion for Coandă MAV in translatory motion can be derived.

Various CFD studies have also indicated that there are still significant discrepancies between CFD and experimental studies [5, 15], which necessitate the existence of specific baseline configurations for validation purposes. Coandă MAV is expected to be capable of maneuvering as illustrated in Figure 1 (a) and (b), and the analysis carried out here is associated with hovering condition, which should give further perception into its hovering performance.

The analysis carried out here is associated with hovering and cruising conditions, which should give further awareness on and can be further elaborated into Coandă MAV maneuvering performance. The baseline analysis on Coandă MAV spherical configuration carried out here refers to and confirms the authors’ previous analysis [31–35].

The generic and baseline physical and mathematical model developed have been utilized to obtain some basic results that are essential in describing the flow field and the physical phenomena related to the relevant surfaces influenced by the Coandă effect jet sheets. Comparison of the numerical computation results for some baseline cases with experimental data under similar conditions will be essential, while computational parametric study will be helpful as a preliminary effort for design and performance optimization. Some examples of Coandă effect applications and developments are illustrated in Figure 2.

It is with such objectives in mind that parametric studies may offer some clues on relevant parameters which may be utilized in a multivariable optimization (and to a larger scale, multidisciplinary optimization). The introduction of Coandă jet on both airfoils and aerodynamic surfaces results in enhanced L/D, which depends on the jet velocity or momentum coefficient. TE (Trailing-edge) rounding-off of in the introduction of the Coandă jet has been shown to be effective in increasing airfoil, L/D, as exemplified by recent applications of Coandă jet for STOL (Short Take-Off and Landing) or ESTOL (Extreme Short Take-Off and Landing), e.g., S809. The study on Coandă jet application to wind turbines that could produce maximum total energy output of Coandă configured airfoil exceeding that predicted by Betz limit could provide further insight on Coandă jet application for MAV.

The test and evaluation of unmanned systems presents special challenges, but these challenges are amplified when one moves into the realm of micro‐air‐vehicles. Further complications arise when these MAVs are fully autonomous [37]. Autonomous MAVs must be able to sense the environment and guide their movement through it. To find their targets, autonomous MAVs should be able to identify them using object recognition procedure, and move toward their goal based on a path planning algorithm. The associated sensors can be employed for an aerial robotic MAV. MAVs have been utilized for surveillance that uses a data link and line‐of‐sight control, as illustrated in Figure 3 (adapted from Davis et al. [37]); advanced MAVs can hover and navigate independently and carry multiple sensors. In addition, the production in large quantities of MAVs at low unit cost can be facilitated by exploiting microfabrication technology.

## 2. Coandă MAV Aerodynamics of hovering semi‐spherical Coandă MAV as a baseline configuration

The Coandă MAV mathematical model utilized here follows closely that of Djojodihardjo et al. [31, 32] and Ahmed et al [33, 34]. The model was conceived to be generic and basic, and allow the application of first principles. The detail is summarized and reproduced here as a baseline reference for developments that may follow. In the Coandă MAV configuration illustrated in Figure 4, an actuator rotor is located at the center of the upper part of the body. As a baseline, the dimension of the rotor can be designed to be small as required. As illustrated in the figure, the flow being drawn by the actuator is entirely utilized for lift, by producing radial flow on the surface of the body as a Coandă jet blanket. In addition, like in a helicopter, a part of the flow produced by the actuator can pass through a vertical stream tube to produce lift (or thrust).

In the analysis of Coandă effect, momentum conservation principle in the sense of Euler equation is applied in the analysis to find the relevant aerodynamic forces and to define the performance parameters of the Coandă MAV. This procedure, applied for a baseline and simplified semi‐spherical Coandă MAV configuration, is shown in Figure 5.

A rigorous analysis is then carried out to see how Coandă effect contributes to the Coandă MAV lift generation capability. For this purpose, the fundamental conservation analysis is applied on the control volume CV (dashed rectangle in Figure 5).

Applying the continuity equation (that is the mass conservation equation) given by

and

where *Vj‐R* is the jet flow velocity, *R* is the vehicle body radius, *h* is the jet slot thickness and *VJ‐out*, R_{o}, *ho*, *VJ‐in Ri*, and *hi* refer to velocity, radius and jet thickness at the outlet and inlet, respectively, of the Coandă jet blanket. As depicted in Figure 6, the, the Coandă jet blanket is exposed to the ambient pressure pa, which is then the prevailing static pressure on it. The momentum equation applied to the control volume in the z (vertical) direction can be differentiated into the contribution due to the Coandă Blanket momentum and the pressure difference on the body due to Coandă Blanket:

The contribution of lift from the momentum flux through the control volume CV in the y (vertical) direction is given by

with z and force positive in the upward direction.

It is noted that the momentum in the radial direction does not contribute to lift. Therefore, as exhibited in Figure 6, the momentum equation in the y direction for the control volume CV is

Before proceeding to analyze the contribution of the forces on the Coandă MAV due to the pressure difference across the control volume, the energy equation applied to the control volume will be considered. It reads

**Energy input to the control volume = Energy output from the control volume**

In the control volume analysis throughout this work, uniform properties across the sectional areas at the input and output of the Coandă jet blanket have been assumed. The entrainment energy exchange between the ambient air and the Coandă jet blanket as well as other energy losses have been ignored.

Hence, referring again to Figure 6, taking into consideration the mass conservation and assuming similar ambient pressure at the inlet and outlet, the application of the energy conservation equation within the control volume which only involves the fluid flow along the Coandă jet blanket results in

Therefore, upon substitution of Eq. (6) into Eq. (5), one obtains

Using the model as shown in Figure 7, the underbody represented by the outer part of the lower surface of the control volume as well as the wake of the MAV are at ambient pressure. The entire flow entering the rotor represented by the actuator disk can be utilized for the Coandă jet blanket after leaving the actuator disk with velocity V_{J‐in}, and then leaves the Coandă MAV after being ejected as *VJ‐out* and combined with the entrained fluid outside the Coandă jet blanket as the wake. If part of the flow is used to feed the Coandă jet blanket, then the rest of the flow leaving the actuator will form a stream tube leaving the Coandă MAV like a ducted fan and eventually will be combined with the exiting jet from the Coandă jet blanket.

Noting the specific characteristics of Coandă effect is attributed to the attachment of a jet to a curves surface, an analysis will be performed on the Coandă device control volume depicted in Figure 6 by focusing on a control volume enclosing radial slice of the Coandă blanket with thickness *h* at (*R*, *θ*). Considering a balance of force along the radial direction *R* in the Coandă jet slice control volume, centrifugal effects along *R* allows a radial pressure gradient to prevail perpendicular to Coandă jet; hence

where

so that

Here the ambient pressure *pa* acts on the outer surface of the Coandă jet blanket, while *pJ‐θ*_{ }acts on the lower surface of the Coandă jet blanket, thus on the upper surface of the Coandă MAV.

The contribution of the vertical forces on the Coandă MAV due to the pressure difference across the control volume can now be considered. This is then give by

Substituting Eq. (8c), the vertical force due to pressure difference is then given by

Since either the plenum chamber between the upper curved surface or the bottom part of the Coandă MAV is at the atmospheric pressure *pa*, then *pa* in Eq. (9b) will drop off. Hence

which is a downward force, or an upward force

Hence, the total lift due to Coandă jet blanket momentum and Coandă jet blanket pressure difference, or the total lift force due to the rate of momentum in the vertical direction and due to pressure difference across the control volume becomes

or

Defining Coandă jet input momentum rate as

The Coandă jet momentum coefficient (as initially defined by Poisson‐Quinton [38]) is conventionally defined as

where U and S should be referred to the ambient reference condition, which here can be chosen to be either the free stream velocity U_{∞} and the lifting area of the Coandă device (MAV), πR_{o}^{2},respectively, or the downward inlet velocity to the Coandă MAV jet inducing rotor, i.e., the velocity of the incoming downward flow to the rotor U_{R} (see Figure 5) and the area of the rotor A_{R}. Hence,

or

For rotor in‐flow completely dedicated for Coandă jet, from conservation of mass principle

and

Hence,

For convenience, in line with the rationale elaborated in [31–34], an indicator of the effectiveness of the Coandă jet device (MAV) to produce lift from rotor‐like momentum rate input from the ambient air, a Performance Measure PM can be defined. Following the previous development and different from [31–34], a new Performance Measure P_{M} can be defined as

Utilizing Eq. (16) and assuming *RJ‐in* ≈ *RR* (without loss of generalities), *C*_{μ‐Rotor} can be expressed as

The performance measure (PM) can alternatively be written in terms of momentum coefficient as

or

The above derivation is carried out on a Coandă MAV configuration with actuator (rotor) dedicated to produce Coandă Jet Blanket. It has tacitly been assumed that the Coandă MAV jet blanket covers the entire upper surface of the Coandă MAV (“perfect blanket”). Therefore, for this Coandă MAV semi‐spherical configuration design, which may represent a baseline of other similar Coandă MAV with moderately different upper surface curvature with MAV with (“perfect blanket”), the wake of the Coandă MAV is shed as a wake following the body of the Coandă MAV. Such situation can be verified at a later stage using CFD simulation for closely modelled flow situation. Then it will be reasonable to assume that the static pressure of the wake of the Coandă MAV the ambient pressure p_{a} prevails. For Coandă MAV with ducted‐fan‐line configuration, an actuator disk model should be added in the force system in addition to those elaborated above.

## 3. Lift generated by the rotor (actuator disk)

Figure7 depicts a MAV configuration consisting of a propeller at the center of the upper part with radially bled flow to generate Coandă effect along the upper surface of the spherical MAV. Here, UR is the axial inlet flow through the actuator along the vertical actuator shaft direction.

Therefore, the lift generated by the propeller can be calculated using actuator disk theory [39–42]

which should be added to Eq. (12), and modified to account for the actuator disk cylindrical part with assumed uniform radius *RI*. Therefore, the lift of the Coandă MAV due to the actuator and Coandă blanket becomes

Here, the mass flow rate into the actuator disk is given by

and one may assume that

and from actuator disk theory [40–43]

In this elaboration, the Betz limit [43], which is the theoretical maximum efficiency of the actuator disk (*η* = 16/27 ≈ 59.26%), is not considered.

A remark is in order here associated with the Coandă MAV as illustrated in Figures. 4, and 5 as represented and the top of Figure 7 above. To avoid the MAV body counter‐rotates due to the rotor, counter‐rotating rotors can be utilized as the actuator. Eqs (13)–(15) remain valid as the baseline equation, but the design configuration in the conceptual and detailed design of the Coandă MAV should be elaborated accordingly.

Several alternatives can be considered for the Coandă MAV with actuator configuration. First, if all the flow into the actuator is utilized for the Coandă MAV, one may obtain a comparison between the lift produced by Coandă MAV without actuator and ducted fan which utilizes actuator flow only. Hence

Here, the ratio between the mass flow rates injected for Coandă MAV jet _{Coandă‐mass}.

Eq. (26) establishes the benefit of Coandă MAV jet blanket lift compared to actuator (ducted fan, or helicopter) lift,

Eq. (22) can be written as

Then the performance measure (PM) for Lift can alternatively be written in terms of Coandă MAV mass ratio

If all the actuator flow is utilized for the Coandă MAV jet blanket, Eq. (28) reduces to

and the performance measure (PM) for Lift of Coandă MAV Eq. (29) reduces to

## 4. Flight dynamic analysis for semi‐spherical Coandă MAV in translatory motion

### 4.1. System of coordinates

For setting the equation of motion of the Coandă MAV, the system coordinates that are required will be defined. Four main reference frames can be identified for this purpose; these are the inertial coordinate system, the local horizon reference frame, the body reference frame, and the wind reference frame. These reference frames are shown in Figure 8.

For simplicity and instructiveness, only two‐dimensional coordinate systems in the plane of symmetry of the Coandă MAV will be defined. These four frame of reference are:

**Inertial coordinate system**, which is used for defining and the application of the Newton’s law of motion. In two dimensions, this coordinate system is depicted by the O_{E}xyz.**The local horizon coordinate system Oxhyhzh**, which is fixed at the center of mass of the vehicle, and is parallel with the earth Inertial Coordinate system.**The body coordinate system Oxbybzb**, which is fixed to the vehicle and follows a conveniently chosen axis of the vehicle.**The wind axes system Oxwywzw**, which moves with the vehicle and the*xw*axis coincides with the flight path of the vehicle, hence its velocity vector.

The wind axes are oriented to the flight path angle γ relative the local horizon axes, as shown in Figure 8a, and by the angle of attack α relative to the body axes. Figure 8b exhibits the two dimensional configuration in the vertical plane perpendicular to the earth surface.

### 4.2. Flight dynamic analysis

Without loss of generalities, the equation of motion of the flight dynamics of Coandă MAV in translatory motion is developed in the vertical plane perpendicular to the earth motion (two‐dimensional planar motion). The hovering state, illustrated in Figure 9, is employed as a reference. In addition, in the present analysis, only the Coandă MAV jet is incorporated in the derivation. The influence of the actuator can be linearly superposed as demonstrated in the hover case, Eqs. (13)–(15).

Balance of forces in the free‐body diagram exhibited in Figure 10 implies

The use of the mass conservation equation

and the use of the Coandă lift from Eq. (11b) then leads to:

Further, simplifying assumptions will be made for the development of the equation of motion during translatory motion. This equation could be further refined to incorporate more realistic cases.

These comprise

During hover, the thrust of the Coandă MAV will not be influenced by the inclination of the Coandă MAV. This assumption can be later modified for more elaborate analytical model, which could be assisted by CFD simulation.

The thrust generated by the Coandă MAV during translatory motion is assumed to stay the same; therefore:

Using the basic flight dynamic equation for flight vehicle, the equation of motion of the Coandă MAV in translational motion in the two‐dimensional vertical plane can be written by referring to Figure 10

Kinematic:

Coandă MAV Fluid Dynamics:

Coandă MAV Flight Dynamics:

Weight‐Fuel Consumption:

where C is the specific fuel consumption, weight per unit thrust [Newton/Newton].

The above relationships are derived for balance of forces and Newton’s equation for semi‐spherical Coandă MAV, which is treated as a point mass moving in a plane perpendicular to the planar Earth.

## 5. Performance measure

### 5.1. Performance measure during hover

The feasibility of using Coandă techniques to enhance aerodynamic performance of Coandă MAV can be outlined through some non‐dimensional quantities of performance measure. The most logical measure that has commonly been utilized is the Coandă jet momentum coefficient, *Cμ*(Poisson‐Quinton [38]), which are defined as (14a) or (14b), as appropriate:

or

or

One can define another Performance Measure to evaluate the aerodynamic performance of spherical Coandă MAV considered here during hover and lift‐off, by assessing the output lift as compared to the input momentum rate, following consideration discussed in [31, 33–34]. These are given by Eqs. (18) and (20) following a different derivation elaborated in previous section.

The Performance Measure may be defined in a different way by referring to the rate of the energy of the Coandă jet ejected at its ideal inlet velocity *Vj‐out* at its peripheral outlet, and comparing it to the rate of momentum influx of the Coandă jet times its inlet velocity *Vj‐in* (non‐dimensional):

based on Coandă MAV geometry and kinematics, or

based on Coandă MAV geometry, kinematics and momentum coefficient. Although Eq. (18) is based on rate of momentum flux while Eq. (45) on energy input, both equations are similar.

For the baseline semi‐spherical Coandă MAV case considered, *VJ‐out*/*VJ‐in* will be the same, and hence Eqs. (45) and (46) reduces to Eqs. (18) and (20). This notion may not be met in more general cases. Since the theoretical analysis does not consider viscosity and other losses, more information is required to insure uniqueness. Further realistic assumption can be drawn by carrying out computational fluid dynamic (CFD) simulation and visualization, which incorporate viscous effects and turbulence, to add further understanding into the problem and introduce further refinement in the theoretical modeling and analysis.

### 5.2. Performance measure during translation

For further reference and as a baseline, Performance Measure can be defined for two cases during translation: for level flight and for climbing.

For level flight (by referring to Figure 11):

For purely level flight, *γ* = 0 (by referring to Figures 10 and 11). Taking into account that the prevailing cross‐sectional area of the Coandă MAV perpendicular to the flow is depicted in Figure 11(b), Eq. (47) reduces to:

Alternatively, another related Performance Measure can be defined as follows:

Here *CD→x* is the Drag coefficient of the Coandă MAV in the x direction, and *α* is its inclination angle with respect to the flight direction as depicted as depicted in Figure 11.

The general expression for the Performance Measure for climbing is:

Here, all variables are taken for the same flow and flight condition.

For purely climbing flight, as schematically depicted in Figure 12, *γ* = *π*/2, and Eq. (45) reduces to

Considering the fuel consumed, Eq. (46) becomes

If we define that, the thrust T of the Coandă MAV should be able to produce climbing motion for a period of θ seconds, straight forward arithmetic calculation using Eq. (47) will require that the propulsion of the Coandă MAV should have a power capacity exceeding

corresponding to providing thrust of

where

These expressions comprise further development of those given in [31, 33–34]. The Performance Measure P_{MCoandă‐MAV Level Flight} can be expressed as a function of C_{D,}R_{o, m˙}*hi, α, V, VJ‐in*. These basic relationships can further be altered to address other flight conditions.

## 6. Results and examples

### 6.1. Results and examples for hover

The theoretical prediction for the performance measures *PM* as given by Eqs. (41) and (42) are depicted in Figures. 13–16. Comparison of the analytical model and the CFD simulation results for another Performance Measure *PM2* defined as the Coandă MAV lift per Coandă jet inlet velocity as suggested in [11, 29] is shown in Figure 17 and elaborated in [31, 33–34]. The Performance Measure P_{M2} is calculated for various Jet slot thickness, and serve as an assessment and validation of the theoretical model. Certainly the CFD model for the same Coandă MAV configuration is based on Reynolds averaged Navier‐Stokes equations which take into account viscous effects, which is not considered in the theoretical model and analysis. Thus, some discrepancies will appear, and the comparison should be judged as such.

Comparison of analytical and CFD [32] models for the Performance Measure *PM2* defined by total lift *L*_{MAV}/ Inlet Velocity *Vj‐in* for spherical Coandă is depicted in Figure 17. The figure clearly exhibits the noticeable influence of the ratio of the jet slot thickness *h* to the reference radius *Ro*, to the Performance Measure *PM2*. The influence of the jet inlet velocity on lift enables designers to appropriately select the right propulsion system specifications for the Coandă MAV.

Figure 18 exemplifies the effect of viscosity for two different locations of the Coandă jet introduction that can be revealed by CFD simulation. The CFD images in Figure 18 (a) and (b), which simulate the velocity magnitude contours for similar configurations considered by Mirkov and Rasuo [11, 29] and Ghozali [30], show the velocity magnitude contours for each of these configurations and exhibit the influence of the inlet configuration. The resulting lift as a function of input jet velocity diagram for each of these configurations are also shown. The corresponding computational grid for each of the configurations utilized by Mirkov and Rasuo, and Ghozali, is also shown as inset in each figure.

For these CFD simulations, which were performed using steady RANS (Reynolds Averaged Navier-Stokes), two equations k‐ω shear stress transport (SST) turbulence model is employed. In the solution procedure followed, simple pressure‐velocity coupling scheme with the least squares cell base as discretization gradient has been applied; second order upwind for the momentum equation and the turbulent kinetic energy were also imposed.

Figure 18 also exhibits structured numerical mesh used in the flow domain. Through meticulous attempts and grid sensitivity studies the size of the mesh cells was generated small to enable visualizations of the flow around the whole body with best details. The CFD simulations are carried out on two configuration versions: in Figure 18 (a), the flow configuration and grid mesh are closely similar to those utilized by Mirkov and Rasuo [11, 29], and in Figure 18 (b), the flow configuration and grid mesh are closely similar to those utilized by Ghozali [30]. The inset of each figure shows the mesh details adopted by Mirkov and Rasuo [29], and Ghozali [30], respectively. Moderately small Reynolds number, *Re* = 68458, based on mean velocity and jet inlet height has been employed in the numerical study. Additionally, across the inlet which has a thickness of 0.05 m, the velocity is st to be uniform (*Vj‐in* = 20 m/s). The CFD simulations utilize computational grid with 52,830 mesh cells. Assessment of the results and referring to the procedure guidelines assure that the mesh quality was acceptable. The CFD computational results indicate that orthogonality quality was 4.96943e‐01. Again, from the ANSYS orthogonality quality requirements, these figures should be acceptable. Further details are elaborated by Djojodihardjo et al. [31, 32–35] and Ahmed et al. [33, 34–36].

Grid sensitivity [18–19] and iteration convergence CFD simulation [44] studies were carried out and exhibited in Figure 19. The mesh quality and comparable agreements between the present CFD simulation results exhibited in Figure 18, as well as the supporting grid mesh and convergence studies exemplified by Figure 19 demonstrates that the generation of mesh in the present CFD simulations meet the present quality requirements as well as those stipulated in the ANSYS™© guidelines; therefore these can be considered to be acceptable.

The CFD simulations are also carried out to investigate the influence of the jet inlet radius on the air vehicle performance (lift). The CFD simulations were carried out for two inlet jet radius, *Ri* = 5 mm and *Ri* = 50 mm for the baseline semi‐spherical Coandă MAV configurations studied. The results are presented as velocity magnitude contours, which exhibit the influence of the inlet radius on the flow field. Figures 20 and 21 exhibits the pressure contours along semi‐spherical Coandă MAV surface obtained by CFD simulation results and the static pressure associated with it in color‐coded contour. Figure 22 compares the CFD simulation results for the pressure coefficient of on the Coandă MAV surface jet blanket obtained from the present CFD simulation and those obtained by Mirkov and Rasuo [29]. Overall, these results serve to verify the plausibility of the analytical as well as the CFD simulation efforts.

### 6.2. Results and examples for translation

For translatory flight, some examples of Performance Measures simulation studies for Coandă MAV in purely level are exhibited by Figures 23 and 24, obtained using Eqs. (44) and (48), respectively.

Figure 23 is worked out for the Coandă Jet inlet Radius R_{J‐in} = 0.02 m, Coandă jet slot thickness *hJ‐In* = 0.01 m, and Coandă jet inlet velocity *VJ‐n*= 5 m/sec.. Figure 24 shows the Thrust ratio *Γ* as a function of versus flight velocity for pure climbing, following Eq. (49) with equal sign, assuming C_{D} ≈ 2.0 and specific fuel consumption C = 0.000139 and period = 1 h. The example is also worked out for the Coandă Jet inlet Radius *RJ‐in* = 0.02 m, Coandă jet slot thickness *hJ‐in* = 0.01 m, and Coandă jet inlet velocity *VJ‐in* = 5 m/sec.

Parametric study results exhibited by these figures can be utilized to give insight on the performance of the semi‐spherical Coandă MAV considered as well as for preliminary sizing purposes.

## 7. Concluding remarks

A comprehensive and meticulous analysis has been made to analyze the potential utilization of Coandă MAV as an aerial robotic platform (vehicle) and deriving the Coandă MAV governing equations. Utilizing first principles, i.e., conservation principles for control volume analysis as baseline fluid‐dynamic analysis for the Coandă MAV in hover and free‐body diagram and Newton’s law of motion as baseline flight‐dynamic analysis for in‐plane translatory motion, the analysis establishes the basic working relationships among various relevant variables and parameters governing the aerodynamics forces and performance measures of Coandă MAV in hover and translatory motion.

To gain insight in the development of the analytical model, and to assess and validate the analytical model subject to the theoretical assumptions, CFD simulations for the Coandă MAV generic models are carried out. Comparison of the theoretical analysis and CFD simulation results indicates the plausibility of the theoretical approach. Performance measures are derived to elaborate the significance of Coandă jet momentum coefficient and to assess the effectiveness of the Coandă MAV as an aerial robot platform.