DOI: 10.5281/zenodo.12085034
title: Virtual RHI lidar scans retrieved in high-fidelity wake vortex simulations of landing aircraft under turbulent crosswind conditions - LES Lidar Simulator (LLS)

description: This dataset contains two types of virtual measurements of multiple pulsed lidar integrated into high-fidelity hybrid RANS-LES wake vortex simulations of a landing Airbus A340 aircraft. Simulations have been performed for four different atmospheric conditions, varying in crosswind and therefore turbulence in the atmosphere.

The two lidar simulator types are:



LLS: LES Lidar Simulator (no noise, Signal to Noise Ratio >> 1). The analysis is based on the Range Gate Weighting function (RWF) [based on formulations in for example [1]]. Given the assumption of no instrument noise, no spectral analysis in the frequency domain is required.

LLSn: LES Lidar Simulator with noise (realisitic Signal to Noise Ratio). Based on formulations in [2], where non-linear low-pass spatial filters accurately model real field measurements. Both background noise and aerosol noise contained in real lidar measurements are modeled. 


Provided are the raw lidar scans in RHI format, with the position and strength of the wake vortices within each lidar scan given by a pressure-vorticity tracking algorithm from the wake vortex simulation (simulation truth, ST). Misidentifications have been removed from the provided dataset. In addition to the labels for each wake vortex lidar scan, also wind background scans are provided for the crosswind simulation cases. These background scans give a better understanding of the prevailing atmospheric condition within which the aircraft lands. 

Furthermore the evaluation of the LLS scans using the Method of Radial Velocities (RV method) [3], a state-of-the-art wake vortex characterization method for lidar scans is provided for a subset of the LLS dataset. For analysing the impact of the RWF for LLS scans, selected lidar scans are also given in a  'naive' fashion (lidar scans assuming point measurements are possible).

--------------------------------------------------------------------------------------------------------------------------------------------

Overview

In total each dataset (LLS and LLSn) contains 8 Aircraft landing simulation with associated virtual lidar RHI scans:



2x no wind

6x crosswind (specified at height b_0)

0.5w_0 from port direction

0.5w_0 from starboard direction

1.0w_0 from port direction

1.0w_0 from starboard direction

2.0w_0 from port direction

2.0w_0 from starboard direction


--------------------------------------------------------------------------------------------------------------------------------------------

The dataset has 13 folders:



POS1_POS2_POS8_scans  = Individual scans of various simulations as well as background wind scans (no wake vortices). 

POS1_POS2: 0_0, 0_5, 1_0, 2_0 corresponds to the wind, same as POS1_POS2 below

POS8: LLS or LLSn, corresponding to the type of lidar simulator (without or with instrument noise, respectively)


 

Scan naming convention: POS1_POS2_POS3_POS4_POS5_POS6_POS7.csv

 

Example: 0_5_D_248_8_161.2201878198302_168.4201878198237.csv

 





POS1_POS2: Together they form a factor which is multiplied with the initial descend speed of the wake vortex pair, w_0. The definition for w_0 can be found in [4]. It is common for the crosswind speed to be set according to the multiples of w_0. Due to the landing of the A340 aircraft, and the logarithmic nature of the wind simulation, we set the crosswind at the b_0 altitude of the simulation. For a definition of b_0, also see [4]. For the above example, 0.5w_0 is the crosswind speed. Note that every 26 lidar positions (POS4), the direction of the crosswind changes, if there is a crosswind.



          - Crosswind approaches from port side: POS4: 0-25, 52-77, 104-129, 156-181, 208-233, 260-285

          - Crosswind approaches from starboard side: POS4: 26-51, 78-103, 130-155, 182-207, 234-259, 286-311





POS3: Label indicating which part of the numerical simulation this scan belongs to (a full landing simulation consists of A-D). A: Hybrid RANS-LES, B,C,D: temporal LES (in order). For the above example, D indicates that the lidar scan was recorded during the last part of the simulation. W: Prior to the wake vortex simulations with crosswind, a background wind scan for each lidar position is recorded.







POS4: Specifies the number of the lidar, implications are the longitudinal position along the glide path of the aircraft. In the packs of lidar positions described in the description of POS1_POS2, the last lidar position is the closest to the touchdown point of the aircraft, smaller lidar positions within this pack are further away (in order) - also see the lid_plane_info directory. Note that the lidar position also adjusts the spectrum of the elevation angles used (and thus the size of the lidar scan). For that see the associated virtual lidar scan raw data. In the above example we have lidar position 248, thus crosswind approaching from starboard and rather mid-way of the longitudinal glide path direction.







POS5: Specifies the number of the scan for this lidar position (POS4) and simulation (POS1__POS2). In the above example this is scan number 8.







POS6_POS7: Specifies the simulation time within which the scan was measured during the aircraft landing simulation (POS1__POS2). In the above example this is 161.2201878198302 s to 168.4201878198237 s.



The scans purely with wind, an no wake vortices are stored in an additional directory within conv_scans. The case '0_0' (no wind) does not require wind scans. 

 

Within the scan files, we have a four (LLS) or five (LLSn) columns:





t: Simulation time

ELE: Elevation angle [deg] of the lidar beam

R: Range from lidar [m]

v_r: LOS velocity (radial velocity) along lidar beam

snr: Signal to Noise Ratio 


 

wind_POS8_scans = Simulated lidar scans of the background wind for LLS and LLSn.

POS8: LLS or LLSn, corresponding to the type of lidar simulator (without or with instrument noise, respectively)

Scans simulated here are in the same format as POS1_POS2_POS8_scans. Wind scans for point measurement scans are found in naive_scan_subset - the scans with POS4 = 0 should be used if more are available.

 

naive_scan_subset = Simulated lidar scans simulating velocity point measurements.




The format is the same as for POS1_POS2_POS8_scans, with the difference of being sorted first by wind subdirectories and then lidar number (LID).

 

labels = Labels of wake vortices for the above wake vortex scans. 


We have 4 main files, where each file represents the targets for one wind strength 0.0w_0, 0.5w_0, 1.0w_0, 2.0w_0.

Within the labels files, we have a multitude of columns with different data:





#Time: Simulation (not scan) time.

y_uw : lateral position in simulation domain of upwind (port) vortex.

z_uw : height position in simulation domain of upwind (port) vortex.

y_dw : lateral position in simulation domain of downwind (starboard) vortex.

z_dw : height position in simulation domain of downwind (starboard) vortex.

G_515_uw: Gamma 515 Circulation strength (see additional notes) of upwind (port) vortex [m^2/s].

G_515_dw: Gamma 515 Circulation strength (see additional notes) of downwind (starboard) vortex [m^2/s].

cr_uw: Core radius of upwind (port) vortex (in meters).

cr_dw: Core radius of downwind (starboard) vortex (in meters).

uw distance from lidar [m]: Horizontal distance of upwind (port) vortex from respective lidar.

dw distance from lidar [m]: Horizontal distance of downwind (starboard) vortex from respective lidar.

uw height [m]: Vertical distance of upwind (port) vortex from respective lidar (floor, as lidars are place at an altitude of 0 m).

dw height [m]: Vertical distance of downwind (starboard) vortex from respective lidar (floor, as lidars are place at an altitude of 0 m).

vortex age uw [s]: Age of upwind (port) vortex (after first generated by aircraft at respective measurement plane of lidar).

vortex age dw [s]: Age of downwind (starboard) vortex (after first generated by aircraft at respective measurement plane of lidar).

scan: Associated scan.

LID: Associated lidar number.


On top of the summarizing files for each simulation, in the subdirectory individual, the label for each scan is given in a separate file.

The above gives information on the simulation truth, furthermore the file labels_with_rv.csv can be found here, where the following extra label columns are given for a subset of scans: [in the following Conv and Naive refer to LLS scans and point measurement scans, respectively]



wind: Strength of crosswind, same as POS1_POS2 in  POS1_POS2_scans. The sign corresponds to the wind direction. Negative is a crosswind from the port side of the aircraft, positive is a crosswind from the starboard side of the aircraft.

Conv RV: 1 indicates the RV method has been evaluated for LLS scans, 0 if not.

Naive RV: 1 indicates the RV method has been evaluated for point measurement scans, 0 if not.

conv rv G_515_uw: Gamma 515 Circulation strength (see additional notes) of upwind (port) vortex computed using the RV method on LLS scans  [m^2/s].

conv rv G_515_dw: Gamma 515 Circulation strength (see additional notes) of downwind (starboard) vortex computed using the RV method on LLS scans [m^2/s].

conv rv uw distance from lidar [m]: Horizontal distance of upwind (port) vortex from respective lidar using the RV method on LLS scans.

conv rv dw distance from lidar [m]: Horizontal distance of downwind (starboard) vortex from respective lidar using the RV method on LLS scans.

conv rv uw height [m]: Vertical distance of upwind (port) vortex from respective lidar (floor, as lidars are place at an altitude of 0 m) using the RV method on LLS scans.

conv rv dw height [m]: Vertical distance of downwind (starboard) vortex from respective lidar (floor, as lidars are place at an altitude of 0 m) using the RV method on LLS scans.

naive rv G_515_uw: Gamma 515 Circulation strength (see additional notes) of upwind (port) vortex computed using the RV method on point measurement scans [m^2/s].

naive rv G_515_dw: Gamma 515 Circulation strength (see additional notes) of downwind (starboard) vortex computed using the RV method on point measurement scans [m^2/s].

naive rv uw distance from lidar [m]: Horizontal distance of upwind (port) vortex from respective lidar using the RV method on point measurement scans.

naive rv dw distance from lidar [m]: Horizontal distance of downwind (starboard) vortex from respective lidar using the RV method on point measurement scans.

naive rv uw height [m]: Vertical distance of upwind (port) vortex from respective lidar (floor, as lidars are place at an altitude of 0 m) using the RV method on point measurement scans.

naive rv dw height [m]: Vertical distance of downwind (starboard) vortex from respective lidar (floor, as lidars are place at an altitude of 0 m) using the RV method on point measurement scans.

uw phi [deg]: Elevation angle to the center of the upwind (port) vortex from the simulation truth.

uw range [m]: Range from the lidar to the center of the upwind (port) vortex from the simulation truth.

dw phi [deg]: Elevation angle to the center of the downwind (starboard) vortex from the simulation truth.

dw range [m]: Range from the lidar to the center of the downwind (starboard) vortex from the simulation truth.

conv rv uw phi [deg]: Elevation angle to the center of the upwind (port) vortex using the RV method on LLS scans.

conv rv uw range [m]: Range from the lidar to the center of the upwind (port) vortex using the RV method on LLS scans.

conv rv dw phi [deg]: Elevation angle to the center of the downwind (starboard) vortex using the RV method on LLS scans.

conv rv dw range [m]: Range from the lidar to the center of the downwind (starboard) vortex using the RV method on LLS scans.

naive rv uw phi [deg]: Elevation angle to the center of the upwind (port) vortex using the RV method on point measurement scans.

naive rv uw range [m]: Range from the lidar to the center of the upwind (port) vortex using the RV method on point measurement scans.

naive rv dw phi [deg]: Elevation angle to the center of the downwind (starboard) vortex using the RV method on point measurement scans.

naive rv dw range [m]: Range from the lidar to the center of the downwind (starboard) vortex using the RV method on point measurement scans.

y err conv uw [m]: Cartesian localization error (RV minus ST) in y-direction (lateral) between simulation truth and RV method of upwind (port) vortex center on LLS scans.

y err conv dw [m]: Cartesian localization error (RV minus ST) in y-direction (lateral) between simulation truth and RV method of downwind (starboard) vortex center on LLS scans.

z err conv uw [m]: Cartesian localization error (RV minus ST) in z-direction (vertical) between simulation truth and RV method of upwind (port) vortex center on LLS scans.

z err conv dw [m]: Cartesian localization error (RV minus ST) in z-direction (vertical) between simulation truth and RV method of downwind (starboard) vortex center on LLS scans.

y err naive uw [m]: Cartesian localization error (RV minus ST) in y-direction (lateral) between simulation truth and RV method of upwind (port) vortex center on point measurement scans.

y err naive dw [m]: Cartesian localization error (RV minus ST) in y-direction (lateral) between simulation truth and RV method of downwind (starboard) vortex center on point measurement scans.

z err naive uw [m]: Cartesian localization error (RV minus ST) in z-direction (vertical) between simulation truth and RV method of upwind (port) vortex center on point measurement scans.

z err naive dw [m]: Cartesian localization error (RV minus ST) in z-direction (vertical) between simulation truth and RV method of downwind (starboard) vortex center on point measurement scans.

D err conv uw [m]: Euclidean distance error for the Cartesian coordinates to the vortex center between simulation truth and RV method of upwind (port) vortex center on LLS scans.

D err conv dw [m]: Euclidean distance for the Cartesian coordinates to the vortex center between simulation truth and RV method of downwind (starboard) vortex center on LLS scans.

D err naive uw [m]: Euclidean distance error for the Cartesian coordinates to the vortex center between simulation truth and RV method of upwind (port) vortex center on point measurement scans.

D err naive dw [m]: Euclidean distance for the Cartesian coordinates to the vortex center between simulation truth and RV method of downwind (starboard) vortex center on point measurement scans.

plus_minus_phi: Indicates whether a RHI lidar scan features a positive lidar scanning rate (plus), or a negative one (negative).

conv err G_515_uw: Circulation error (RV minus ST) between simulation truth and RV method of upwind (port) on LLS scans [m^2/s]. 

conv err G_515_dw: Circulation error (RV minus ST) between simulation truth and RV method of downwind (starboard) on LLS scans [m^2/s]. 

naive err G_515_uw: Circulation error (RV minus ST) between simulation truth and RV method of upwind (port) on point measurement scans [m^2/s]. 

naive err G_515_dw: Circulation error (RV minus ST) between simulation truth and RV method of downwind (starboard) on point measurement scans [m^2/s]. 

conv err uw phi: Elevation angle error (RV minus ST) between simulation truth and RV method of upwind (port) vortex center on LLS scans [deg].

conv err dw phi: Elevation angle error (RV minus ST) between simulation truth and RV method of downwind (starboard) vortex center on LLS scans [deg].

naive err uw phi: Elevation angle error (RV minus ST) between simulation truth and RV method of upwind (port) vortex center on LLS scans [deg].

naive err dw phi: Elevation angle error (RV minus ST) between simulation truth and RV method of downwind (starboard) vortex center on point measurement scans [deg].

conv err uw range: Range from lidar error (RV minus ST) between simulation truth and RV method of upwind (port) vortex center on LLS scans [m].

conv err dw range: Range from lidar error (RV minus ST) between simulation truth and RV method of downwind (starboard) vortex center on LLS scans [m].

naive err uw range: Range from lidar error (RV minus ST) between simulation truth and RV method of upwind (port) vortex center on point measurement scans [m].

naive err dw range: Range from lidar error (RV minus ST) between simulation truth and RV method of downwind (starboard) vortex center on point measurement scans [m].


0_0_movement.csv gives insight on the vertical movement on selected vortices in selected lidar scans of the 0_0 wind case. 

The following data columns exist within the file (note csv file with ; separator):



Str_dw: Vertical movement of the downwind (starboard) vortex, 0 is neutral, -1 is downward, and 1 is upward.

Prt_uw: Vertical movement of the upwind (port) vortex, 0 is neutral, -1 is downward, and 1 is upward.

scan_code: Identfies associated lidar scan with POS3_POS4_POS5 as in POS1_POS2_scans.


 

lid_plane_info = Additional guidance on the simulated lidars.

Each simulation case has a separate file due to minimal ground speed changes of the aircraft. 

The following data columns exist within each file:



Index: Gives the lidar number corresponding to POS4 and LID.

Plane Pos [x]: Gives the longitudinal position of the lidar with respect to the glide path of the aircraft in meters.

time plane passed [s]: Gives the simulation time when the aircraft first passes the measurement plane (related to the previous column.

height plane passed [m]: Aircraft altitude at the respective lidar positon.

lateral position lidar from GP [m]: Lateral position of the lidar from the glide path of the aircraft.


--------------------------------------------------------------------------------------------------------------------------------------------


Additional Notes: 





Gamma 515 circulation [4]: The Gamma 515 circulation is the averaged circulation of a vortex evaluated at radii 5-15m from the vortex center. 









Lidar scans may contain 2 or 1 vortex. In the case of 1 vortex, the other has G_515_?? set to nan. Note that the position may still available in the labels data sets, but only vortices with available G_515_?? values should be used.

The upwind and downwind definitions should purely be taken as names, rather than physical meaning. The names are derived from the simulation, due to mirroring and other post-processing these can be misleading however.




--------------------------------------------------------------------------------------------------------------------------------------------



Funding: This dataset was generated within the mFUND KIWI project funded by the Federal Ministry for Digital and Transportation Germany and the DLR undertaking "Wetter und Disruptive Ereignisse".

--------------------------------------------------------------------------------------------------------------------------------------------

Acknowledgements:

The necessary RANS simulations were performed as part of the EU-funded AWIATOR by DLR's Institute of Aerodynamics and Flow Technology.

We acknowledge Airbus for the allowance to use the RANS data.

The LES was performed with the incompressible Navier-Stokes code MGLET [5] and kindly provided by the Technical University of Munich, Hydromechanics.

The wake vortex simulations were computed on the high performance computer SuperMUC-NG by Leibniz-Rechenzentrum (LRZ).


--------------------------------------------------------------------------------------------------------------------------------------------

 

References: 

 

[1] Robey, Rachel, and Julie K. Lundquist. "Behavior and mechanisms of Doppler wind lidar error in varying stability regimes." Atmospheric Measurement Techniques 15.15 (2022): 4585-4622.

 

[2] Stephan, Anton, Norman Wildmann, and Igor Smalikho. "Effectiveness of the MFAS Method for Determining the Wind Velocity Vector from Windcube 200s Lidar Measurements." Atmospheric and Oceanic Optics 32.5 (2019): 555-563.

 

[3] Smalikho, Igor., et al. "Method of radial velocities for the estimation of aircraft wake vortex parameters from data measured by coherent Doppler lidar." Optics Express      23.19 (2015): A1194-A1207.

 

[4] Gerz, Thomas, Frank Holzäpfel, and Denis Darracq. "Commercial aircraft wake vortices." Progress in Aerospace Sciences 38.3 (2002): 181-208. 

 

[5] Manhart, Michael. "A zonal grid algorithm for DNS of turbulent boundary layers." Computers & fluids 33.3 (2004): 435-461.