Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 181 records · Page 10

Numerical simulation of the subsolar magnetopause current layer in the sun-earth meridian plane

The formation and stability of the magnetopause current layer near the subsolar point in the sun-earth meridian plane are examined using a 2D electromagnetic particle simulation. For the case of zero IMF, the simulation results show that the current layer remains stable and is essentially the same as in the 1D simulation. The width of the current layer is given by the electron-ion hybrid gyroradius which is much smaller than the ion gyroradius. The current layer is found to remain stable for the northward IMF as well. As in the 1D simulation, the jump of the magnetic field at the current layer for the northward IMF remains small. For the southward IMF, collisionless magnetic reconnection is found to develop, leading to the formation of magnetic islands and density peaking within the current layer.

Okuda, H.↗

An Analysis of the Impact of Extreme Attitude Operation on a Turbofan Engine in a Regional Jet Aircraft

The objective of this work is to characterize the variability in engine performance at extreme attitudes, resulting from the use of different standard control variables, evaluated at different points in the engine’s lifecycle. This paper describes the procedure used to create a dynamic model of an aft-mounted jet engine on a T-tail regional jet aircraft. This model enables simulation of engine operation at extreme attitudes. The model was subsequently evaluated at various altitudes, Mach numbers, and power settings over a range of angles of attack and sideslip, and deterioration levels. Each case was simulated using both fan speed and engine pressure ratio as the engine control variable. The results show that the wing has a very large impact on the engine operation, overwhelming other sources of variation.

loss of control - inflight↗

A computer program for the simulation of heat and moisture flow in soils

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Camillo, P.↗

Shock Generation and Control Using DBD Plasma Actuators

This report is the final report of a NASA Phase I SBIR contract, with some revisions to remove company proprietary data. The Shock Boundary Layer Interaction (SBLI) phenomena in a supersonic inlet involve mutual interaction of oblique shocks with boundary layers, forcing the boundary layer to separate from the inlet wall. To improve the inlet efficiency, it is desired to prevent or delay shock-induced boundary layer separation. In this effort, Innovative Technology Applications Company (ITAC), LLC and the University of Notre Dame (UND) jointly investigated the use of dielectric-barrier-discharge (DBD) plasma actuators for control of SBLI in a supersonic inlet. The research investigated the potential for DBD plasma actuators to suppress flow separation caused by a shock in a turbulent boundary layer. The research involved both numerical and experimental investigations of plasma flow control for a few different SBLI configurations: (a) a 12 wedge flow test case at Mach 1.5 (numerical and experimental), (b) an impinging shock test case at Mach 1.5 using an airfoil as a shock generator (numerical and experimental), and (c) a Mach 2.0 nozzle flow case in a simulated 15 X 15 cm wind tunnel with a shock generator (numerical). Numerical studies were performed for all three test cases to examine the feasibility of plasma flow control concepts. These results were used to guide the wind tunnel experiments conducted on the Mach 1.5 12 degree wedge flow (case a) and the Mach 1.5 impinging shock test case (case b) which were at similar flow conditions as the corresponding numerical studies to obtain experimental evidence of plasma control effects for SBLI control. The experiments also generated data that were used in validating the numerical studies for the baseline cases (without plasma actuators). The experiments were conducted in a Mach 1.5 test section in the University of Notre Dame Hessert Laboratory. The simulation results from cases a and b indicated that multiple spanwise actuators in series and at a voltage of 75 kVp-p could fully suppress the flow separation downstream of the shock. The simulation results from case c showed that the streamwise plasma actuators are highly effective in creating pairs of counter-rotating vortices, much like the mechanical vortex generators, and could also potentially have beneficial effects for SBLI control. However, to achieve these effects, the positioning and the quantity of the DBD actuators used must be optimized. The wind tunnel experiments mapped the baseline flow with good agreement to the numerical simulations. The experimental results were conducted with spanwise actuators for cases a and b, but were limited by the inability to generate a sufficiently high voltage due to arcing in the wind-tunnel test-section. The static pressure in the tunnel was lower than the static pressure in an inlet at flight conditions, promoting arching and degrading the actuator performance.

Patel, Mehul P.↗

A simple lightning parameterization for calculating global lightning distributions

A simple parameterization has been developed to simulate global lightning distributions. Convective cloud top height is used as the variable in the parameterization, with different formulations for continental and marine thunderstorms. The parameterization has been validated using two lightning data sets: one global and one regional. In both cases the simulated lightning distributions and frequencies are in very good agreement with the observed lightning data. This parameterization could be used for global studies of lightning climatology; the earth's electric circuit; in general circulation models for modeling global lightning activity, atmospheric NO(x) concentrations, and perhaps forest fire distributions for both the present and future climate; and, possibly, even as a short-term forecasting aid.

Price, Colin↗

Integrated Flight/Structural Mode Control for Very Flexible Aircraft Using L1 Adaptive Output Feedback Controller

This paper explores application of adaptive control architecture to a light, high-aspect ratio, flexible aircraft configuration that exhibits strong rigid body/flexible mode coupling. Specifically, an L(sub 1) adaptive output feedback controller is developed for a semi-span wind tunnel model capable of motion. The wind tunnel mount allows the semi-span model to translate vertically and pitch at the wing root, resulting in better simulation of an aircraft s rigid body motion. The control objective is to design a pitch control with altitude hold while suppressing body freedom flutter. The controller is an output feedback nominal controller (LQG) augmented by an L(sub 1) adaptive loop. A modification to the L(sub 1) output feedback is proposed to make it more suitable for flexible structures. The new control law relaxes the required bounds on the unmatched uncertainty and allows dependence on the state as well as time, i.e. a more general unmatched nonlinearity. The paper presents controller development and simulated performance responses. Simulation is conducted by using full state flexible wing models derived from test data at 10 different dynamic pressure conditions. An L(sub 1) adaptive output feedback controller is designed for a single test point and is then applied to all the test cases. The simulation results show that the L(sub 1) augmented controller can stabilize and meet the performance requirements for all 10 test conditions ranging from 30 psf to 130 psf dynamic pressure.

Che, Jiaxing↗

Effects of Mesoscale Eddies in the Active Mixed Layer: Test of the Parametrisation in Eddy Resolving Simulations

In eddy resolving simulations, we test a mixed layer mesoscale parametrisation, developed recently by Canuto and Dubovikov [Ocean Model., 2011, 39, 200-207]. With no adjustable parameters, the parametrisation yields the horizontal and vertical mesoscale fluxes in terms of coarse-resolution fields and eddy kinetic energy (EKE). We compare terms of the parametrisation diagnosed from coarse-grained fields with the eddy mesoscale fluxes diagnosed directly from the high resolution model. An expression for the EKE in terms of mean fields has also been found to get a closed parametrisation in terms of the mean fields only. In 40 numerical experiments we simulated two types of flows: idealised flows driven by baroclinic instabilities only, and more realistic flows, driven by wind and surface fluxes as well as by inflow-outflow. The diagnosed quasi-instantaneous horizontal and vertical mesoscale buoyancy fluxes (averaged over 1-2 degrees and 10 days) demonstrate a strong scatter typical for turbulent flows, however, the fluxes are positively correlated with the parametrisation with higher (0.5-0.74) correlations at the experiments with larger baroclinic radius Rossby. After being averaged over 3-4 months, diffusivities diagnosed from the eddy resolving simulations are consistent with the parametrisation for a broad range of parameters. Diagnosed vertical mesoscale fluxes restratify mixed layer and are in a good agreement with the parametrisation unless vertical turbulent mixing in the upper layer becomes strong enough in comparison with mesoscale advection. In the latter case, numerical simulations demonstrate that the deviation of the fluxes from the parametrisation is controlled by dimensionless parameter estimating the ratio of vertical turbulent mixing term to mesoscale advection. An analysis using a modified omega-equation reveals that the effects of the vertical mixing of vorticity is responsible for the two-three fold amplification of vertical mesoscale flux. Possible physical mechanisms, responsible for the amplification of vertical mesoscale flux are discussed.

simulation↗

A simulation investigation of cockpit display of aircraft traffic during curved, descending, decelerating approaches

The results of a simulation experiment involving the evaluation of cockpit display of aircraft traffic information are presented. The experiment was conducted using taped time dependent, noninteractive traffic in an approach to landing situation and two levels of pilot control models: 3-D automatic and computer augmented control. The tests involved two cases: the simulation aircraft flew approach paths which (1) followed another aircraft in between two other aircraft, and (2) merged between two other aircraft. Speed control via manual throttles was used in all tests (path stretching was not allowed for maintaining separation between aircraft). The approaches were conducted while the simulation aircraft was conducting a curved, descending, decelerating approach to landing. Performance data sets were examined, and subjective opinions regarding workload were gathered. Traffic positioning was varied to further evaluate the test subjects' monitoring performance.

Steinmetz, G. G.↗

Investigation of Transonic Wake Dynamics for Mechanically Deployable Entry Systems

A numerical investigation of transonic flow around a mechanically deployable entry system being considered for a robotic mission to Venus has been performed, and preliminary results are reported. The flow around a conceptual representation of the vehicle geometry was simulated at discrete points along a ballistic trajectory using Detached Eddy Simulation (DES). The trajectory points selected span the low supersonic to transonic regimes with freestream Mach numbers from 1:5 to 0:8, and freestream Reynolds numbers (based on diameter) between 2:09 x 10(exp 6) and 2:93 x 10(exp 6). Additionally, the Mach 0:8 case was simulated at angles of attack between 0 and 5 . Static aerodynamic coefficients obtained from the data show qualitative agreement with data from 70deg sphere-cone wind tunnel tests performed for the Viking program. Finally, the effect of choices of models and numerical algorithms is addressed by comparing the DES results to those using a Reynolds Averaged Navier-Stokes (RANS) model, as well as to results using a more dissipative numerical scheme.

Stern, Eric↗

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗

Numerical simulation of cosmogenic nuclide production in lunar rocks

The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.

Reedy, R. C.↗

Introducing Multisensor Satellite Radiance-Based Evaluation for Regional Earth System Modeling

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.

Planetary Boundary Layer↗

Operating Small Sat Swarms as a Single Entity: Introducing SODA

Swarm concepts are a growing topic of interest in the small satellite community. Compared to a small satellite constellation, a swarm has the distinction of being multiple spacecraft in close proximity, in approximately the same orbit. Furthermore, we envision swarms to have capabilities for cross-link communication and station-keeping. Of particular interest is a means to maintain operator-specified geometry, alignment, and/or separation.From NASA's decadal survey, it is clear that simultaneous measurements from a 3D volume of space are desired for a variety of Earth scientific studies. As this mission concept is ultimately extended to deep space, some degree of local control for the swarm to self-correct its configuration is required. We claim that the practicality of ground commanding each individual satellite in the swarm is simply not a feasible concept of operations. In other words, the current state-of-practice does not scale to very large swarms (e.g. 100 spacecraft or more) without becoming cost prohibitive. To contain the operations costs and complexity, a new approach is required: the swarm must be operated as a unit, responding to high-level specifications for relative position and velocity.The Mission Design Division at NASA Ames Research Center is looking to the near future for opportunities to develop satellite swarm technology. As part of this effort, we are developing SODA (Swarm Orbital Dynamics Advisor), a tool that provides the orbital maneuvers required to achieve a desired type of relative swarm motion. The purpose of SODA is two-fold. First, it encompasses the algorithms and orbital dynamics model to enable the desired relative motion of the swarm satellites. The process starts with the user specifying the properties of a swarm configuration. This could be as simple as varying in-track spacing of the swarm in one orbit, or as complex as maintaining a specified 3D geometrical orientation. We presume that science objectives will drive this choice. Given these inputs, the tool provides the most efficient maneuver(s) to achieve the objective.Second, SODA provides a variety of visualization tools. We acknowledge that the relationship between a desired relative motion amongst the swarm, and the corresponding orbital parameters for each individual satellite may not be immediately apparent for ground controllers and mission planners. The purpose of SODA's visualization tools is to illustrate this concept clearly with a variety of graphics and animations. After computing the optimal orbital maneuvers to modify the swarm, these results are simulated to demonstrate successful swarm control.Our emphasis in this paper is on the importance of relating the desired motion of the swarm satellites relative to one another with the required orbital element changes. One cannot joystick a drifting swarm satellite back into position; the underlying orbital mechanics dictate the most efficient recovery maneuvers. To illustrate this point, results from several case study simulations are presented. We conclude with our forward work for ongoing SODA development and potential science applications.

small satellites↗

Demonstration of Initial GlennICE Relative Frame Capability: Axial-Flow Propeller

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

Aircraft Icing↗

Demonstration of Initial GlennICE Relative Frame Capability: Axial-Flow Propeller

Modifications have been implemented in the GlennICE software to accommodate a non-inertial reference frame. GlennICE accepts a flow solution from an external flow solver. It then introduces particles and tracks them through the flow field in a Lagrangian manner. Centrifugal and Coriolis terms were added to the GlennICE software to account for relative frame simulations. The objective of the present paper is twofold. First, to check that the new terms are implemented correctly and that the code still behaves as expected with respect to convergence. And second, to provide some initial insight into an upcoming propeller experiment in the NASA Icing Research Tunnel. The paper presents a description of the code modifications. In addition, results are presented for two operating conditions, and three particle sizes. Each case was simulated with four different grid densities to assess grid dependence.

Aircraft Icing↗

Development of a Novel Helmet Support Assembly for NASA Orion Crew Survival Suit

Orion Spacecraft water landing environments have necessitated the design and testing of a novel Helmet Support Assembly (HSA) as part of the NASA Orion Crew Survival Systems (OCSS) launch and re-entry spacesuit. A series of development sled tests using 5 th female, 50th male, and 95th male Hybrid III Anthropometric Test Devices (ATDs) simulated worst case water landings and identified the need for an integrated suit countermeasure to control the helmet during dynamic loading. Initial prototype countermeasures rigidly restrained the helmet keeping the helmet in place as desired, but led to large force and moment transmission to the upper neck. Further development led to a flexible HSA which used high-strength steel wire contoured to the size-matched ATD thorax and attached to both the front and rear of the neck ring. The selected wire diameter and contoured shape provided some ability to flex under the expected worst-case dynamic loads. The existing shoulder harness was used to restrain the helmet motion during “eyeballs out/down” loading while also attempting to prevent head-to-helmet contact in “eyeballs down” spinal-axis loading. A new small prototype helmet was also developed and tested for the 5th and 50th ATDs with both HSA devices. Dynamic impact tests were conducted on the HIA (Horizontal Impulse Accelerator) at Wright Patterson Air Force Base (WPAFB) under simulated off-nominal landing conditions using a representative Orion seat and 5-point harness. For each ATD size, peak Neck Injury Criteria (Nij) were compared to equivalent sled tests with an unsuited ATD configuration. Adding a helmet without attempting to control motion increased peak Nij values for all ATD sizes. The use of the rigid countermeasure showed decreased peak Nij results for the 5th female, but increased peak Nij values for both 50th configurations and the 95th male ATDs. By comparison, the flexible HSA showed reduced peak Nij values for all ATD sizes. Overall, this series of dynamic impact tests identified a risk of increased crew neck loading when under- or over-restraining the helmet during simulated Orion water landings and also demonstrated an effective strategy to mitigate those loads with a compliant HSA design.

Safety↗

Experimental Demonstration of Microwave Signal/Electric Thruster Plasma Interaction Effects

An experiment was designed and conducted in the Electric Propulsion Laboratory of NASA Lewis Research Center to assess the impact of ion thruster exhaust plasma plume on electromagnetic signal propagation. A microwave transmission experiment was set up inside the propulsion test bed using a pair of broadband horn antennas and a 30 cm 2.3 kW ion thruster. Frequency of signal propagation covered from 6.5 to 18 GHz range. The stainless steel test bed when enclosed can be depressurized to simulate a near vacuum environment. A pulsed CW system with gating hardware was utilized to eliminate multiple chamber reflections from the test signal. Microwave signal was transmitted and received between the two hours when the thruster was operating at a given power level in such a way that the signal propagation path crossed directly through the plume volume. Signal attenuation and phase shift due to the plume was measured for the entire frequency band. Results for this worst case configuration simulation indicate that the effects of the ion thruster plume on microwave signals is a negligible attenuation (within 0.15 dB) and a small phase shift (within 8 deg.). This paper describes the detailed experiment and presents some of the results.

Zaman, Afroz J.↗

Scale Dependence of Cirrus Horizontal Heterogeneity Effects on TOA Measurements: MODIS Brightness Temperatures in the Thermal Infrared - Part I

This paper presents a study on the impact of cirrus cloud heterogeneities on MODIS simulated thermal infrared (TIR) brightness temperatures (BTs) at the top of the atmosphere (TOA) as a function of spatial resolution from 50 meters to 10 kilometers. A realistic 3-D (three-dimensional) cirrus field is generated by the 3DCLOUD model (average optical thickness of 1.4, cloudtop and base altitudes at 10 and 12 kilometers, respectively, consisting of aggregate column crystals of D (sub eff) equals 20 microns), and 3-D thermal infrared radiative transfer (RT) is simulated with the 3DMCPOL (3-D Monte Carlo Polarized) code. According to previous studies, differences between 3-D BT computed from a heterogenous pixel and 1-D (one-dimensional) RT computed from a homogeneous pixel are considered dependent at nadir on two effects: (i) the optical thickness horizontal heterogeneity leading to the plane-parallel homogeneous bias (PPHB); and the (ii) horizontal radiative transport (HRT) leading to the independent pixel approximation error (IPAE). A single but realistic cirrus case is simulated and, as expected, the PPHB mainly impacts the low-spatial resolution results (above approximately 250 meters), with averaged values of up to 5-7 K (thousand), while the IPAE mainly impacts the high-spatial resolution results (below approximately 250 meters) with average values of up to 1-2 K (thousand). A sensitivity study has been performed in order to extend these results to various cirrus optical thicknesses and heterogeneities by sampling the cirrus in several ranges of parameters. For four optical thickness classes and four optical heterogeneity classes, we have found that, for nadir observations, the spatial resolution at which the combination of PPHB and HRT effects is the smallest, falls between 100 and 250 meters. These spatial resolutions thus appear to be the best choice to retrieve cirrus optical properties with the smallest cloud heterogeneity-related total bias in the thermal infrared. For off-nadir observations, the average total effect is increased and the minimum is shifted to coarser spatial resolutions.

TOA↗