Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesoscale”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-Scale Simulation of Wind Farm Performance during a Frontal Passage

Predicting the response of wind farms to changing flow conditions is necessary for optimal design and operation. Here in this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine–flow interaction, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with 10 m grid spacing. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed background velocity, turbine wake structure, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD model framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗
DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗
DOE OSTI2024

Focused Ion Beam Tomography of Alloy 617 Corroded in Molten Chloride Salt

Materials qualification of reactor structural materials is a critical step in rapid implementation of advanced nuclear reactor technologies, particularly to assess the corrosion performance in these designs. Accelerated qualification of reactor structural materials requires incorporating powerful computational toolsets, such as phase field modelling in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to predict the evolution of structural materials due to corrosion. Accordingly, computational toolsets will require experimental data generated at appropriate length scales to validate accuracy. Focused ion beam (FIB) provides a high degree of control over manipulation of materials for analytical purposes, including capturing data on the evolution in the microstructure and elemental composition of materials at the mesoscale, an appropriate length scale for phase field modelling of intergranular diffusion phenomena using the MOOSE framework. For instance, the FEI Helios G4 UX dual beam plasma FIB microscope at the Irradiated Materials Characterization Laboratory (IMCL) is capable of backscatter diffraction (EBSD) and energy-dispersive x-ray spectroscopy (EDS) documenting the evolution in the microstructure and elemental composition, respectively. The Helios can perform EDS and EBSD three-dimensionally (3D) using tomography, which is then combined using different software packages to visualize 3D volumes correlating elemental composition to microstructural data. The purpose of this investigation was to develop a streamlined characterization and data processing workflow for 3D tomography studies on the FEI Helios G4 plasma FIB. The investigation is segmented into three parts: 1) Optimizing the data collection workflow, 2) identifying appropriate data processing and visualization software (i.e. DREAM.3D, MIPAR, and VGStudioMax), and 3) establishing an infrastructure for public release. The optimization of the data collection workflow is in collaboration with members of the U220 department to setup formal training on the tomography operation of the G4, through ThermoFisher Scientific, and exploring DREAM.3D, MIPAR, and VGStudioMax data processing/visualization software packages. VGStudioMax currently demonstrates the most promise for future use. Optimization of the data collection and processing workflow is still ongoing. A collaboration with INL High Performance Computing (HPC) established an open-source license for expediting the public release of FIB tomography datasets through HPC. FIB tomography data generated by the G4 will provide comprehensive data for validating 3D phase field mesoscale modelling tools within the MOOSE framework for accelerated qualification of reactor structural materials.

Copeland-Johnson, Trishelle↗
DOE OSTI2021

Contrasting characteristics of open- and closed-cellular stratocumulus cloud in the eastern North Atlantic

Abstract. Extensive regions of marine boundary layer cloud impact the radiative balance through their significant shortwave albedo while having little impact on outgoing longwave radiation. Despite this importance, these cloud systems remain poorly represented in large-scale models due to difficulty in representing the processes that drive their life cycle and coverage. In particular, the mesoscale organization and cellular structure of marine boundary clouds have important implications for the subsequent cloud feedbacks. In this study, we use long-term (2013–2018) observations from the Atmospheric Radiation Measurement (ARM) Facility's Eastern North Atlantic (ENA) site on Graciosa Island, Azores, Portugal, to identify cloud cases with open- or closed-cellular organization. More than 500 h of each organization type are identified. The ARM observations are combined with reanalysis and satellite products to quantify the cloud, precipitation, aerosol, thermodynamic, and large-scale synoptic characteristics associated with these cloud types. Our analysis shows that both cloud organization populations occur during similar sea surface temperature conditions, but the open-cell cases are distinguished by stronger cold-air advection and large-scale subsidence compared to the closed-cell cases, consistent with their formation during cold-air outbreaks. We also find that the open-cell cases were associated with deeper boundary layers, stronger low-level winds, and higher rain rates compared to their closed-cell counterparts. Finally, raindrops with diameters larger than 1 mm were routinely recorded at the surface during both populations, with a higher number of large drops during the open-cellular cases. The similarities and differences noted herein provide important insights into the environmental and cloud characteristics during varying marine boundary layer cloud mesoscale organization and will be useful for the evaluation of model simulations for ENA marine clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Boundary layer moisture variability at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic observatory during marine conditions

Abstract. Boundary layer moisture variability at the Eastern North Atlantic (ENA) site during marine conditions is examined at monthly and daily timescales using 5 years of ground-based observations and output from the European Center for Medium range Weather Forecast (ECMWF) reanalysis model. The annual cycle of the mixed-layer total water budgets is presented to estimate the relative contribution of large-scale advection, local moisture tendency, entrainment, and precipitation to balance the moistening due to surface latent heat flux on monthly timescales. When marine conditions prevail, advection of colder and dry air from the north acts as an important moisture sink (∼ 50 % of the overall budget) during fall and winter driving the seasonality of the budget. Entrainment and precipitation contribute to the drying of the boundary layer (∼ 25 % and ∼ 15 % respectively), and the local change in moisture contributes to a small residual part. On a daily temporal scale, moist and dry mesoscale columns of vapor (∼ 10 km) are analyzed during 10 selected days of precipitating stratocumulus clouds. Adjacent moist and dry columns present distinct mesoscale features that are strongly correlated with clouds and precipitation. Dry columns adjacent to moist columns have more frequent and stronger downdrafts immediately below the cloud base. Moist columns have more frequent updrafts, stronger cloud-top cooling, and higher liquid water path and precipitation compared to the dry columns. This study highlights the complex interaction between large-scale and local processes controlling the boundary layer moisture and the importance of spatial distribution of vapor to support convection and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Identifying synoptic controls on boundary layer thermodynamic and cloud properties in a regional forecast model

Although most of our understanding of boundary layer cloudiness is based on idealized, subtropical, barotropic marine environments, boundary layer clouds exist across a range of conditions. In this study, we use the Naval Research Laboratory's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) and an automated cold-front-relative analysis framework to explore the boundary layer structure associated with low clouds across a transect through the cold front of a midlatitude synoptic cyclone. The model credibly captures boundary layer structure in line with conceptual models of midlatitude cyclones and ground-based observations at Graciosa Island in the Azores. The warm sector is conditionally unstable, with clouds that are too shallow and with too little liquid water, compared to cloud property retrievals from satellite and surface-based instruments. The cold-frontal region exhibits convection associated with weak stability and ascent. Northwest of the cold front, the boundary layer is well mixed, deeper, and capped by a strong inversion maintained by large-scale subsidence. Simulated clouds in frontal and post-frontal regions are mostly too thick, with too much liquid water and too little cloud base drizzle. The post-frontal clouds are associated with grid-scale updrafts, which appear to be the model's attempt to represent mesoscale organization of cellular convection typically observed in the cold sector of midlatitude cyclones. The deep, well-mixed post-frontal boundary layers and cloudiness are maintained by strong surface fluxes, as in cold air outbreaks. Our analysis framework serves as a unique approach to model verification, and our results offer insights into differences in boundary layer cloud behavior between subtropical and synoptic cold-sector regimes.

Eissner, Jordan M. [Univ. of Kansas, Lawrence, KS ↗
DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The remarkable inefficiency of stratocumulus

Marine stratocumulus clouds play a central role in Earth's climate system by reflecting incoming solar radiation and exerting a strong cooling effect. Their organization into open and closed mesoscale cellular morphologies can be thought of as an example of bistable dynamics driven by aerosol–cloud interactions and mesoscale processes. From the perspective of non-equilibrium thermodynamics, these structures are an example of a far-from-equilibrium open system that continuously produces and exports entropy. While entropy production has been studied in idealized deep convective systems, it has not yet been quantified for shallow clouds. Here, we compute and decompose the internal entropy production of open- and closed-cell stratocumulus using an ensemble of large-eddy simulations. We show that the overall entropy production of stratocumulus is low, reflecting the limited vertical extent and corresponding reduced ability to utilize the energy fluxes at the system's boundaries. Moist processes dominate the overall irreversibility, which, combined with their low entropy production, leads to a mechanical efficiency about an order of magnitude smaller than in deep convective systems. Although the dominant irreversible processes differ between open- and closed-cell regimes, the distributions of total entropy production largely overlap across the ensemble, limiting the ability to distinguish the dynamics of individual cases based solely on total entropy production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

PyFLEXTRKR: a flexible feature tracking Python software for convective cloud analysis

Abstract. This paper describes the new open-source framework PyFLEXTRKR (Python FLEXible object TRacKeR), a flexible atmospheric feature tracking software package with specific capabilities to track convective clouds from a variety of observations and model simulations. This software can track any atmospheric 2D objects and handle merging and splitting explicitly. The package has a collection of multi-object identification algorithms, scalable parallelization options, and has been optimized for large datasets including global high-resolution data. We demonstrate applications of PyFLEXTRKR on tracking individual deep convective cells and mesoscale convective systems from observations and model simulations ranging from large-eddy resolving (∼100s m) to mesoscale (∼10s km) resolutions. Visualization, post-processing, and statistical analysis tools are included in the package. New Lagrangian analyses of convective clouds produced by PyFLEXTRKR applicable to a wide range of datasets and scales facilitate advanced model evaluation and development efforts as well as scientific discovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗
DOE OSTI2021

Assessing boundary condition and parametric uncertainty in numerical-weather-prediction-modeled, long-term offshore wind speed through machine learning and analog ensemble

To accurately plan and manage wind power plants, not only does the time-varying wind resource at the site of interest need to be assessed but also the uncertainty connected to this estimate. Numerical weather prediction (NWP) models at the mesoscale represent a valuable way to characterize the wind resource offshore, given the challenges connected with measuring hub-height wind speed. The boundary condition and parametric uncertainty associated with modeled wind speed is often estimated by running a model ensemble. However, creating an NWP ensemble of long-term wind resource data over a large region represents a computational challenge. Here, we propose two approaches to temporally extrapolate wind speed boundary condition and parametric uncertainty using a more convenient setup in which a mesoscale ensemble is run over a short-term period (1 year), and only a single model covers the desired long-term period (20 year). We quantify hub-height wind speed boundary condition and parametric uncertainty from the short-term model ensemble as its normalized across-ensemble standard deviation. Then, we develop and apply a gradient-boosting model and an analog ensemble approach to temporally extrapolate such uncertainty to the full 20-year period, for which only a single model run is available. As a test case, we consider offshore wind resource characterization in the California Outer Continental Shelf. Both of the proposed approaches provide accurate estimates of the long-term wind speed boundary condition and parametric uncertainty across the region (R 2 >0.75), with the gradient-boosting model slightly outperforming the analog ensemble in terms of bias and centered root-mean-square error. At the three offshore wind energy lease areas in the region, we find a long-term median hourly uncertainty between 10 % and 14 % of the mean hub-height wind speed values. Finally, we assess the physical variability in the uncertainty estimates. In general, we find that the wind speed uncertainty increases closer to land. Also, neutral conditions have smaller uncertainty than the stable and unstable cases, and the modeled wind speed in winter has less boundary condition and parametric sensitivity than summer.

17 WIND ENERGY↗
DOE OSTI2020

Radiation Damage Studies on Titanium Alloys as High Intensity Proton Accelerator Beam Window Materials

A high-strength dual a+ß phase titanium alloy Ti–6Al–4V is utilized as a material for beam windows in several accelerator target facilities. However, relatively little is known about how material properties of this alloy are affected by high-intensity proton beam irradiation. With plans to upgrade neutrino facilities at J-PARC and Fermilab to over 1 MW beam power, the radiation damage in the window material will reach a few displacements per atom (dpa) per year, significantly above the ~0.3 dpa level of existing data. The international RaDIATE collaboration, Radiation Damage In Accelerator Target Environments, has conducted a high intensity proton beam irradiation of various target and window material specimens at Brookhaven Linac Isotope Producer (BLIP) facility, including a variety of titanium alloys. Post-Irradiation Examination (PIE) of the specimens in the 1st capsule, irradiated at up to 0.25 dpa, is in progress. Tensile tests in a hot cell at Pacific Northwest National Laboratory (PNNL) exhibited a clear signature of radiation hardening and loss of ductility for Ti–6Al–4V, while Ti–3Al–2.5V, with less ß phase, exhibited less severe hardening. Microstructural investigations will follow to study the cause of the difference in tensile behavior between these alloys. High-cycle fatigue (HCF) performance is critical to the lifetime estimation of beam windows exposed to a periodic thermal stress from a pulsed proton beam. The first HCF data on irradiated titanium alloys are to be obtained by a conventional bend fatigue test at Fermilab and by an ultrasonic mesoscale fatigue test at Culham Laboratory. Specimens in the 2nd capsule, irradiated at up to ~1 dpa, cover typical titanium alloy grades, including possible radiation-resistant candidates. These systematic studies on the effects of radiation damage of titanium alloy materials are intended to enable us not only to predict realistic lifetimes of current beam windows made of Ti–6Al–4V, but also to extend the lifetime by choosing a more radiation and thermal shock tolerant alloy with a preferable heat treatment, or even by developing new materials.

Titanium Alloy, Beam Window, Proton Beam, Radiatio↗
NASA NTRS1974

Results of the US contribution to the joint US/USSR Bering Sea experiment

The atmospheric circulation which occurred during the Bering Sea Experiment, 15 February to 10 March 1973, in and around the experiment area is analyzed and related to the macroscale morphology and dynamics of the sea ice cover. The ice cover was very complex in structure, being made up of five ice types, and underwent strong dynamic activity. Synoptic analyses show that an optimum variety of weather situations occurred during the experiment: an initial strong anticyclonic period (6 days), followed by a period of strong cyclonic activity (6 days), followed by weak anticyclonic activity (3 days), and finally a period of weak cyclonic activity (4 days). The data of the mesoscale test areas observed on the four sea ice option flights, and ship weather, and drift data give a detailed description of mesoscale ice dynamics which correlates well with the macroscale view: anticyclonic activity advects the ice southward with strong ice divergence and a regular lead and polynya pattern; cyclonic activity advects the ice northward with ice convergence, or slight divergence, and a random lead and polynya pattern.

Campbell, W. J.↗
NASA NTRS1976

Three-dimensional interplanetary stream magnetism and energetic particle motion

Cosmic rays interact with mesoscale configurations of the interplanetary magnetic field. A technique is presented for calculating such configurations in the inner solar system, which are due to streams and source conditions near the sun, and maps of magnetic field are constructed for some plausible stream and source conditions. One effect of these mesoscale configurations on galactic cosmic rays is shown to be an out-of-the-ecliptic gradient drift sufficient to explain Forbush decreases. The effects on solar energetic particles include small polar drifts due to the field gradients and a possibly large modification of the time-intensity profiles and anisotropy characteristics due to the formation of mirror configurations in space. If a diffusion model is applicable to solar particles, the true diffusion coefficient will be masked by the effects of streams. A conceptual model which incorporates these ideas and those of several other models is presented.

Barouch, E.↗
NASA NTRS1976

Advanced atmospheric sounder and imaging radiometer /AASIR/ for STORMSAT

The principal mission of the three-axis stabilized STORMSAT spacecraft is to provide the necessary meteorological data for tracking, studying the detailed structure, and modeling mesoscale weather phenomena. In the area of mesoscale events, the following meteorological objectives are indicated: high-quality imagery, visible and infrared; wind velocity from cloud tracers (1 m/sec), atmospheric temperature profiles (1 K), and atmospheric humidity sounding. These objectives are reflected in the functional characteristics of the AASIR, which is a second generation meteorological sensor based on the Visible Infrared Spin-Scan Radiometer (VISSR) and the Atmospheric Sounder (VAS). The AASIR design and interface constraints with the STORMSAT spacecraft is discussed.

Chase, S. C.↗
NASA NTRS1979

NASA's participation in the AVE-SESAME '79 program

NASA's Marshall Space Flight Center participated with its AVE (Atmospheric Variability Experiment) in a large interagency mesoscale and severe storms experiment identified herein as AVE-SESAME '79 (Atmospheric Variability Experiment-Severe Environmental Storms and Mesoscale Experiment 1979). A primary objective of NASA was to support an effort to acquire carefully edited sets of rawinsonde data during selected severe weather events for use in correlative and diagnostic studies with satellite and radar data obtained at approximately the same times. Data were acquired during six individual 24-h experiments on both the regional and storm scales over a network in the central United States that utilized approximately 20 supplemental rawinsonde sites meshed among 23 standard National Weather Service sites. Included among the six experiments are data obtained between 1200 GMT on April 10 and 1200 GMT on April 11, encompassing the formation and development period for the tornado-producing systems that devastated Wichita Falls, Texas, and other sections of Oklahoma and Texas. The other dates for which data sets are available are April 19-20 and 25-26, May 9-10 and 20-21, and June 7-8, 1979.

Hill, K.↗
NASA NTRS1981

Case studies of clear air turbulence using the diagnostic Richardson Number Tendency formulation

The results of four case studies of clear air turbulence (CAT) using the diagnostic Richardson number tendency (DRT) formulation are highlighted. The performance of this technique in resolving regions of documented CAT encounters is encouraging. Its operational adaptability appears particularly attractive in that input data can be supplied by the currently operational rawinsonde system. Two CAT indices are calculated deterministically, sensing synoptic scale changes in static stability and vertical wind shear conductive for supporting mesoscale CAT layers. These two indices reveal volumes of the troposphere which act as source regions for patches of CAT. The first, the time to reach the critical Richardson number necessary for the initiation of turbulent conditions, was first used by Oard (1974). A second index is devised which is an adaptation of Roach's (1970) work relating synoptic scale and mesoscale energetical coupling and gives more information on the relative intensity of these source regions. The output from the DRT computer module highlights specific regions of the atmosphere which can be interpreted operationally in terms of a CAT encounter probability.

Keller, J. L.↗