1. Introduction
UAV platforms have received increasing attention from both researchers and industries due to their low cost, efficiency, and reduced risk [
1]. Their wide-ranging applications have led to extensive research across multiple fields, as UAVs evolve from their original role in surveillance and reconnaissance operations to play essential roles in executing various operational tasks [
2]. For instance, in a near-future electronic engagement scenario, UAVs are estimated to perform quite similarly to how conventional military aircraft perform a variety of maneuvers while facing multiple threats. Likewise, UAVs and unmanned combat aerial vehicles (UCAVs) have become prevalent on the battlefield, valued for their affordability, low risk, and cost-effectiveness as carriers while also fulfilling critical operational functions [
3]. Central to these engagements is the interaction between land radar systems and these airborne platforms. Within this context, the radar system maintains a fixed position, continuously monitoring the aircraft’s aerial whereabouts. Nevertheless, the reflected power from the aircraft experiences temporal fluctuations, notably during specific maneuvers that modify the aspect angle discerned by the radar, thereby influencing its tracking and detection precision [
4]. This paper centers on modeling the potential radar cross-section (RCS) dynamics of a UAV/UCAV throughout a flight trajectory, with the tracking radar serving as the primary reference point. The literature presents numerous RCS simulation techniques, and the selection of the appropriate one hinges on factors such as the size of the object model and the wavelength of the radar. Subsequently, based on the frequency, the method for RCS calculation can be determined. For instance, in the case of an air superiority military aircraft, an L-band radar may be utilized [
5]. By assessing the ratio of the aircraft’s size (l) to the wavelength of the L-band radar (λ), since l is much bigger than λ, it becomes evident that the analysis should be conducted using high-frequency approximations [
6]. Consequently, methods suited for the optical region, such as geometrical optics (GO), physical optics (PO), and shooting and bouncing rays (SBRs) [
7,
8,
9] could be employed for accurate RCS analysis. The RCS of an aircraft is influenced by various factors including time, frequency, polarization, incident angle, physical geometry, and surface material [
10]. Traditionally, a single RCS value is considered as the time-averaged value of the aircraft at a particular polarization [
11]. However, for complex shapes such as jet fighters or UAV/UCAVs, the scattering geometry and scatterer units on the airplane vary significantly from simple shapes like spheres or cubes. Consequently, the physical characteristics of the aircraft play a crucial role in determining the scattered wave by influencing the incident radar wave. Therefore, they affect RCS characteristics of complex scatterers in both monostatic and bistatic scenarios. This paper focuses specifically on evaluating the monostatic case. The frequency of the emitted incident wave directly influences the phase variation of the scattered wave across various sections of the aircraft’s surfaces or scatterers, leading to fluctuations in RCS scintillation [
12]. While the importance of physical geometry and exterior features in the analysis is evident, changes in radar emission direction and frequency can significantly alter the aspect angle parameters and distinctive RCS characteristics of the targeted aircraft, respectively [
11].
One approach to mitigate radar threats during flight is through strategic flight path planning in relation to radar locations. To minimize the probability of detection by radars dispersed throughout a battlefield while optimizing fuel efficiency, aircraft employ various algorithms to determine the optimal path. A modified algorithm targeting, specifically, stealth drones navigating in an operational environment saturated with both fixed and dynamic radar systems, is presented in [
13]. The method relies on both flight path and attitude adjustments to enhance low observability (stealth) by optimizing the aircraft’s radar cross-section (RCS) profile, specifically, based on the RCS pattern of the F-117 platform. The authors focus on avoiding detection and demonstrate superior results with the proposed method. They do not consider scenarios where the aircraft has already been detected, so address different survivability challenges accordingly. Similarly, the authors in [
14] explore a comparable approach by analyzing the RCS of an ellipsoid-shaped glider, using the PO method combined with Gaussian filter post-processing. Aircraft attitude adjustments and flight trajectory optimization are key strategies to minimize detection by multiple radars while maintaining the shortest possible path.
The investigation conducted in [
12] examines aspect angle calculations within the context of “cone–cylinder rotated body aircraft”, which encompass commercial airliners and ballistic missiles. This analysis operates under the assumption that the RCS of such aircraft remains constant or isotropic in the roll dimension. Consequently, any roll maneuver is anticipated to have minimal to negligible impact on the RCS of the aircraft under consideration. In another work, detailed in [
15], the authors focused on aircraft tracking and classification utilizing passive radar technology and analyzing the radar cross-section of a hypothetical aircraft. They employed very high frequency (VHF) and ultra-high-frequency (UHF) bands, along with a method of moments (MoM) solver for the RCS analysis. Furthermore, the flight maneuvers were narrowed into three primary modes: constant velocity and altitude, climb or descent, and sustained turn at a constant speed. Subsequently, they formulated a Bayesian classification algorithm for the joint tracking and classification of aircraft by its RCS. In a further comprehensive investigation, detailed in [
16], the authors adopted an empirical approach to address the platform RCS and aspect angle analysis conundrum. They utilized three distinct datasets derived from different aircraft: the PA-28-181 Archer II, the Boeing 737, and the Inertial Navigation System (INS) units embedded within the aircraft. Upon obtaining the INS data, a simulation framework was employed to replicate the flight paths of the aircraft, incorporating a fixed-position radar to calculate aspect angles within the simulation environment. Although the degrees of freedom for the three aircraft did not exhibit significant variations, typically remaining within 10 degrees and often even less than 2 degrees due to predefined flight paths, a comprehensive environment was established for the analysis of target–radar interactions. In addition to the kinematic approaches of the dynamic RCS problem, research into cognitive radar applications is gaining momentum and demonstrating significant advancements in target tracking [
17]. The authors of [
18] focus on the detection of small drones, which pose a threat in civilian areas, using the YOLO algorithm. By utilizing data from an IRIS FMCW radar equipped with a rotating antenna on a moving platform, drone detection and classification were achieved with 99% accuracy in an environment containing various objects such as birds, wind turbines, and ground-based items. Traditional active decoys may be rendered obsolete by the superior processing and adaptive capabilities of cognitive tracking technologies on the battlefield. To effectively counter the feedback mechanisms that enable cognitive radars to sense, learn, and adapt to their environment, active decoy algorithms and capabilities must be enhanced [
19]. Hence, active decoys need to become more dynamic and adaptive. Integrating the RF properties of the aircraft deploying the decoy can provide critical inputs for refining the decoy’s jamming functions.
When employing radar for tracking an aircraft or platform, conventionally regarded as a mass entity, emphasis is typically placed solely on monitoring alterations in spatial location [
11]. However, it is imperative to acknowledge that the dynamic movement executed by the target encompasses not only translational motion but also variations in orientation [
14]. The radar cross-section depends on physical characteristics and external details of airplane and aspect direction of the tracking radar. Thus, in the evaluation of the RCS of an aircraft, due regard must be given to the fluctuation in RCS in relation to the diverse angles of attitude. Therefore, the correlation or fluctuation pattern between the RCS of the target and its orientation angles can be explained [
20]. In addition to the relative three-dimensional location of the aircraft with respect to the illuminating radar, the Euler angles, particularly, roll, pitch, and yaw, of the targeted aircraft establish a direct relation with the RCS [
21]. As mentioned, the Euler angles, which are comprised of roll, pitch, and yaw, serve to describe not only the orientation of a rigid body but also the attitude of a mobile frame of reference in 3-dimensional space. Throughout a complex flight trajectory, an aircraft undergoes continuous movement within both global coordinates by altering its own position, and local coordinates by adjusting its heading and attitude. The parameters of these coordinate systems operate independently, implying that the aircraft’s location and attitude are unrelated, each holding distinct significance for aspect angle calculations. These features are also known as six degrees of freedom (6DoF). An unclassified technical report [
22] presents a methodology for calculating the target aspect angle by leveraging aircraft attitude angles and the relative spherical position to the radar system. Consequently, the report delineates two angles that characterize the orientation of the radar beam’s incidence upon the aircraft. This approach delineates the optimal direction for aircraft pilots to locate radar systems, with calculation outcomes expressing the radar-to-platform vector direction in terms of platform spherical coordinates.
This study proposes an application domain for the aspect angle estimation method outlined in [
22], integrating aircraft aspect angles with corresponding RCS values at Euler angles. As illustrated in
Figure 1, an assessment of the verification of the aspect angle estimation approach is implemented by utilizing a controlled flight dataset, and then, the RCS simulation is conducted for an aircraft. Throughout this process, challenges encountered in implementing the method proposed in [
22] are addressed to acquire a usable dataset of aspect angles along a flight trajectory, aligned with the simulated RCS dataset. Furthermore, the commercial flight simulator software, DCS World [
23], is employed to generate an intricate and realistic flight time-series dataset. In brief, a highly maneuverable aircraft is flown in a realistic engagement scenario where air surveillance radar is present. Thanks to the freedom of controlling the aircraft in the simulator, numerous turns and roll maneuvers are conducted during the flight. Moreover, to be able to extract the time-series dataset with 6DoF from the flight conducted in the simulator, another commercially available simulation environment, TacView [
24], is utilized. The same software can be used for reviewing and monitoring the whole flight, from start to finish. Then, the acquired time series of 6DoF flight data is used as the input of the aspect angle estimation technique. The problems encountered in the verification of the proposed technique are considered, and some modifications are then applied. As a result, aircraft-to-radar aspect angle data are acquired for targeted flights created in the simulator. Ultimately, all acquired data are fused to establish a flight dynamic RCS profile specific to a designated radar location.
The contributions of the paper can be listed as follows:
We demonstrate the impact of aspect angle on the RCS fluctuations of an aircraft as observed by a ground-based radar threat. To this end, RCS estimation outputs from various flight scenarios were compared and analyzed.
We investigate the influence of flight dynamics on RCS fluctuations of a UAV-sized aircraft, utilizing readily available tools instead of traditional, costly flight tests or simulations. This demonstrates that practical data through affordable and accessible methods can be acquired.
We establish an operational framework that may help develop countermeasure strategies, including advanced active decoys.
We not only evaluate the effect of flight profiles on aspect angles but also thoroughly examine their impact on RCS fluctuations, offering a quantifiable evaluation.
Our dynamic RCS modeling can be generalized for all engagement scenarios with all operational actors.
The paper is structured as follows:
Section 2 describes the aspect angle procedure and verifies its results in the first two subheadings, then discusses the RCS simulation parameters and details of the 3D model used, and finally, details the flight simulation method.
Section 3 analyzes the results of the flight simulation characteristics, explains the aspect angle results of the subject flight, and evaluates the RCS fluctuation results obtained from integrating the spherical RCS data of the model and aspect angle data of the simulated flight.
Section 4 discusses the methods and results of the paper and outlines future work.
Section 5 concludes the paper.
3. Results
The RCS results were derived using a spherical scope of a 3D aircraft model, as opposed to the cylindrical model employed in [
12], where RCS was simplified for roll maneuvers. This method highlighted the roll angle’s significant impact on RCS fluctuations throughout the flight, enabling the simulation of aircraft movements to closely mimic real-world scenarios and not confining the analysis to a few specific flight modes as in [
15]. However, implementing the aspect angle calculations from [
22] for this study revealed several challenges in integrating the two datasets. To address the azimuth angle output issue—where the algorithm produced values only between 0° and 180°, and sometimes negative elevation values—assumptions of symmetry in the platform RCS were made. An orientation process was subsequently implemented to align the RCS values with the correct angles during data integration. Dynamic flight RCS data for the aircraft were gathered using accessible and cost-effective data acquisition methods, unlike in [
16], which used real INS data and uniform flight patterns across various aircraft, for which the attitude angles are minute. Essential characteristics of the flight data obtained from the flight simulator are detailed in
Figure 9, where angular degrees of freedom and altitude values throughout the flight are plotted. The sudden drop in the ζ value around the 30 s mark results from a 360° roll maneuver. Aside from that, the aircraft generally remained parallel to the ground most of the time. As seen in
Figure 10, the roll angle instances create a normal distribution centered around 0°. Sharp falls and rises in the heading (ξ) plot indicate transitions between 0° and 359° on the compass, rather than an unrealistic rotation on the compass plane. Between the 10 and 40 s marks, the nose of the aircraft generally pointed towards the west, resulting in 270° data values during that interval. Pitch (η) values are not sharply accumulated around the 0° mark like ζ, as the aircraft tends to point upwards, leading to a high number of instances around the 10° mark. Generally, samples are distributed between 0° and 20°. For a stable flight, all these parameters are supposed to peak at one point. However, in this case, the aircraft is quite maneuverable and utilizes this ability in an engagement scenario, thus the plots show variability in complex flight scenarios.
The azimuth aspect angle calculations resulted in a more dispersed data cluster compared to the data found in [
16], as shown in
Figure 11. Due to the radar’s position near the north-northeast of the aircraft, the azimuth aspect angle values accumulated in two intervals: 0–50° and 65–80°. The aircraft predominantly flew northward, leading to accumulation in the 0–50° range, while eastward flights caused the 65–80° interval to capture the remaining instances. The elevation values were predictable and comparable with [
16] because the aircraft’s bottom pointed down for most of the flight. Despite some extreme roll angle values, the roll angle peaked near 0°, as seen in
Figure 10. When roll and pitch angles were ignored and set to zero throughout the flight, resulting in a flight conducted solely with heading and altitude changes, the azimuth results did not change significantly, maintaining clear grouping. The elevation aspect angle values vary between −2° and −0.5°, contrasting with the −90° to 90° range seen in the fully disturbed data, demonstrating the significance of rotational degrees of freedom on aspect angles.
The effect of the flight profile on aspect angles is given in
Table 2. It is important to note that the cases in [
16] disabled all degrees of freedom, including heading, while in this paper, only roll and pitch are disabled. The rotation in heading significantly affects aspect angles in the case of this paper, so reaffirming this was deemed unnecessary. Additionally, the heading angle variation in [
16] was minimal, making the disabling of the heading variable closely align with real data. In contrast, the data used in this paper show significant changes in heading, and setting it to zero would result in unrealistic flight profiles, such as an aircraft traveling west while its nose points north, and overall, causing incomparability with the data provided in [
16]. To avoid this, heading was not altered primarily, allowing for direct comparison of elevation, which is independent of the heading parameter. Moreover, while in [
16] the azimuth scale includes negative values, this paper scales azimuth between 0° and 180°. This difference is manageable and does not significantly impact the analysis. The major difference in the maximum value of elevation is expected since it is largely dependent on roll and elevation. Regardless, if heading is set to zero, the gap between the maximum and median value of the azimuth aspect angle becomes narrow, indicating that the variation only depends on the location of the aircraft and relative distance to the radar, while the rotational effect plays an important role in data with full disturbance.
Following the acquisition of the two datasets, namely, the RCS simulation results and aspect angle calculations, they are finally integrated to elucidate the radar-specific RCS behavior of the aircraft, as seen in
Figure 12, visualizing the RCS scintillation level of the platform on the radar throughout the flight time, and some of these fluctuation behavior changes at the highlighted points are observed in the TacView [
24] flight monitor at the corresponding time. At point A, the aircraft makes a roll and pitch maneuver, resulting in only showing the bottom of the body with a skewed angle. At point B, another roll maneuver is made suddenly and the RCS is reduced drastically, which can also be seen at point C, where a 360° roll maneuver is conducted. After 40 s, the aircraft goes into a more stable flight profile, attempting less attitude changes. Lastly, at point D the plane turns its nose slightly towards the radar, resulting in a little peak in the flight dynamic RCS profile. Furthermore, the RCS fluctuation results derived from the flight data depicted in
Figure 11, where the roll and pitch angles were zeroed, are presented using the blue line. Beyond the annotated points, numerous instances exhibit significant deviations of varying magnitudes. Towards the end of the flight, a notable divergence is observed, with the blue line peaking near 15 dBsm, contrasting sharply with the red line, which shows a marked decrease to approximately −25 dBsm. Another distinct disparity appears around the 47 s mark, with an RCS difference of nearly 30 dBsm evident. Overall, the flight data reflecting realistic conditions with full degrees of freedom display a more dynamic RCS profile, while the constrained dataset results in a less volatile and less unstable RCS scintillation. Although a detailed analysis and evaluation of the major differences between the two RCS datasets may reveal more critical insights, the average difference between them is found to be 32.44% when using the full-disturbance RCS as the reference.
As illustrated in
Figure 13, fluctuations in the RCS are discernible along the flight path. Notably, these fluctuations are intricately tied to the specific location of the radar. It is imperative to note that even if the flight path remained unaltered, variations in the observer’s position would yield entirely different RCS fluctuations. In essence, the red plot delineates the exposure level of the radar throughout the duration of the flight, serving as a comprehensive representation of the radar’s interaction with the aircraft.
The results of the study show that the methodologic approach utilized near-realistic, accessible and repeatable simulations to depict RCS fluctuation behavior under dynamic flight conditions. The effects of the degrees of freedom are analyzed by comparing the statistics of different flight profiles and RCS scintillation results of the aspect angles.