Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel planes”

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 415 records · Page 23

Vector Radiative Transfer Code SORD: Performance Analysis and Quick Start Guide

We present a new open source polarized radiative transfer code SORD written in Fortran 9095. SORD numerically simulates propagation of monochromatic solar radiation in a plane-parallel atmosphere over a reflecting surface using the method of successive orders of scattering (hence the name). Thermal emission is ignored. We did not improve the method in any way, but report the accuracy and runtime in 52 benchmark scenarios. This paper also serves as a quick start users guide for the code available from ftp:maiac.gsfc.nasa.govpubskorkin, from the JQSRT website, or from the corresponding (first) author.

polarized radiative transfer↗

Controllable Planar Optical Focusing System

An optical device has a first metasurface. A high-contrast pattern of the first metasurface is operable for modifying, over a first phase profile, a phase front of an incident light beam. A second metasurface, is disposed over a plane parallel to the first metasurface with a second high-contrast pattern and operable for shaping, over a second phase profile, the modified phase front of the incident light beam into a converging spherical phase front. A spacer layer, in which the modified phase front of the incident light beam diffracts, is disposed in a controllably changeable separation between the first and second metasurfaces. Controllably changing the separation between the first and the second metasurfaces by a first distance correspondingly changes the position of the focus point of the converging spherical phase front by a second distance.

Arbabi, Amir↗

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jets used to suppress boundary layer separation and increase maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate stream wise direction. It also generate turbulent eddies at the oscillation frequency, which are typically large relative to boundary layer turbulence, and which augmenting mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from the oscillator, and the suppression of boundary layer separation by an array of sweep jets are performed. Simulation results are compared to data from a dedicated CFD validation experiment of a single oscillator and its sweep jet, and from a study of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets.2, 20 A critical step in the work is the development of realistic time-dependent sweep-jet in flow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the Over flow CFD solver, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗

Electromagnetic Scattering by Discrete Random Media. I: The Dispersion Equation and the Configuration-Averaged Exciting Field

We consider the scattering of a plane electromagnetic wave obliquely incident on a plane-parallel layer of discrete random medium with non-scattering boundaries. We solve the Lax integral equation for the conditional configuration-averaged exciting field coefficients by assuming a special-form solution, that is, by representing the conditional configuration-averaged exciting field coefficients as a linear combination of the coefficients corresponding to an up-going and a down-going wave. This solution representation is supposed to be valid within the whole domain occupied by the particles, even in the close proximity of the boundaries. By balancing the waves with different propagation directions and wavenumbers we derive two homogeneous systems of equations corresponding to the generalized Lorenz–Lorentz law and two inhomogeneous systems of equations corresponding to the generalized Ewald–Oseen extinction theorem. It is shown that (i) the two homogeneous systems of equations of the generalized Lorenz–Lorentz law reduce to a single homogeneous system of equations corresponding to a semi-infinite discrete random medium at normal incidence; (ii) the dispersion equation is direction and polarization independent; and (iii) the two inhomogeneous systems of equations of the generalized Ewald–Oseen extinction theorem can be reduced to two scalar equations by means of the addition theorem for vector spherical harmonics. It is also shown that the same dispersion equation can be obtained without assuming a special-form solution representation in the proximity of the boundaries.

Adrian Doicu↗

Matrix Exponential in C/C++ Version of Vector Radiative Transfer Code IPOL

We use only left eigenvectors to evaluate the matrix exponential in the method of discrete ordinates for the vector radiative transfer equation, which neglects circular polarization, in a plane-parallel atmosphere. This is contrary to a common practice of using the right eigenvectors to evaluate the matrix exponential combined with the left eigenvectors to avoid the inversion of the matrix of the right ones. Two numerical tests for Rayleigh and Aerosol scattering confirm our idea. For better explanation of our approach and for independent crosscheck of our results, we distribute an example in C/C++.

Sergey Korkin↗

Acoustic Testing of the Tiltrotor Test Rig in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel

The Tiltrotor Test Rig (TTR) was tested in the National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel from 2017 to 2018. The rotor system can be configured in airplane mode, with the rotor plane perpendicular to the wind flow, and in helicopter mode, with the rotor plane parallel to the wind flow. Four microphones were placed around the TTR: two on the wind tunnel floor and two on struts. The primary goal of the test was to understand the operational capabilities of the TTR, while also acquiring research data as available. Limited measurements of the blade vortex interaction (BVI) noise of the TTR rotor were taken to not only understand the acoustic testing capabilities of the TTR in the NFAC 40- by 80-Foot Wind Tunnel, but to also compare to previous tests and to be used for future validation studies. In particular, data will be compared to measurements of an XV-15 rotor previously acquired in the NFAC 80- by 120-Foot Wind Tunnel.

Tiltrotor Test Rig↗

A FDTD Solution of Scattering of Laser Beam with Orbital Angular Momentum by Dielectric Particles: Far-Field Characteristics

Electromagnetic (EM) beams with orbital angular momentum (OAM) may have great potential applications in communication technology and in remote sensing of the Earth-atmosphere system and outer planets. Study of their interaction with optical lenses and dielectric or metallic objects, or scattering of them by particles in the Earth-atmosphere system, is a necessary step to explore the advantage of the OAM EM beams. In this study, the 3-dimensional (3D) scattered-field (SF) finite-difference time domain (FDTD) technique with the convolutional perfectly matched layer (CPML) absorbing boundary conditions (ABC) is applied to calculate the scattering of the purely azimuthal (the radial mode number is assumed to be zero) Laguerre–Gaussian (LG) beams with the OAM by dielectric particles. We found that for OAM beam׳s interaction with dielectric particles, the forward-scattering peak in the conventional phase function (P11) disappears, and light scattering peak occurs at a scattering angle of ~15° to 45°. The disappearance of forward-scattering peak means that, in laser communications most of the particle-scattered noise cannot enter the receiver, thus the received light is optimally the original OAM-encoded signal. This feature of the OAM beam also implies that in lidar remote sensing of the atmospheric particulates, most of the multiple-scattering energy will be off lidar sensors, and this may result in an accurate profiling of particle layers in the atmosphere or in the oceans by lidar, or even in the ground when a ground penetration radar (GPR) with the OAM is applied. This far-field characteristics of the scattered OAM light also imply that the optical theorem, which is derived from plane-parallel wave scattering case and relates the forward scattering amplitude to the total cross section of the scatterer, is invalid for the scattering of OAM beams by dielectric particles.

Electromagnetic beams↗

TPSAS-NF1676L-17894-DND

Are domain-aggregated, long-term climate trends in retrieved reflected solar cloud properties detectable despite instrument related and algorithm-related aliasing errors? Science Question: By how much do changes in instrument calibration distort reflected solar cloud property decadal trends? Objective: Quantify how changes in calibration error feed into cloud property changes. Science Question: By how much do reflected solar cloud property retrieval assumptions distort trends in retrieved cloud properties? Objective: Quantify how the plane-parallel assumption in reflected solar satellite cloud retrievals distorts cloud optical thickness decadal trends by determining how and if the large-scale 3D optical thickness bias changes over time. - The three panel procedure below will be performed for each time (month), (zonal) region, set of solar and viewing geometries.

Yolanda Roberts↗

Global X-Ray Properties of the Vela and Puppis A Supernova Remnants

The Vela and Puppis A supernova remnants(SNRs)comprise a large emission region of∼8°diameter in the soft X-ray sky. The HaloSat CubeSat mission provides the first soft X-ray (0.4–7 keV) observation of the entire Vela SNR and Puppis A SNR region with a single pointing and moderate spectral resolution. HaloSat observations of the Vela SNR are best fit with a two-temperature thermal plasma model consisting of a cooler component with=-+kT0.1910.010.01keV in collisional ionization equilibrium and a hotter component with=-+kT1.0620.270.45 keV in nonequilibrium ionization. Observations of the Puppis A SNR are best fit with a single-component plane-parallel shocked plasma model with=-+kT0.860.050.06 keV in nonequilibrium ionization. For the first time, we find the total X-ray luminosities of both components of the Vela SNR spectrum in the 0.5–7 keV energy band to be= ́-+L4.410X1.41.434 erg s−1 for the cooler component and= ́-+L4.110X1.51.834 erg s−1for the hotter component. We find the total X-ray luminosities of the Vela and Puppis A SNRs to be= ́L8.4 10X34 erg s−1 and= ́-+L6.710X0.91.136 erg s−1

E. M. Silich↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Applying Orbital Multi-Angle Photopolarimetric Observations to Study Properties of Aerosols in the Earth's Atmosphere: Implications of Measurements in the 1.378 µm Spectral Channel to Retrieve Microphysical Characteristics and Composition of Stratospheric Aerosols

We analyze the possibilities of orbital photopolarimetric measurements to study properties of aerosols in the Earth's atmosphere. As an example, we consider the case when such measurements are performed within a narrow spectral channel centered at 1.378 µm that allows to retrieve microphysical characteristics of stratospheric aerosols separately from those of tropospheric aerosols. We consider the case of stratospheric aerosols caused by volcanic eruption, and adopt the model of the stratosphere in the form of a homogeneous plane-parallel layer composed of polydisperse spherical particles. We use numerically exact solutions of the vector radiative transfer equation to theoretically simulate measurements carried out at various numbers of scattering angles, including: (i) radiance measurements alone; (ii) polarization measurements alone; and (iii) radiance and polarization measurements together. The results of computations show that the simultaneous use of radiance and polarization measurements at a sufficiently large number of scattering angles enables one to retrieve the optical thickness, effective radius, and refractive index of aerosols with adequate accuracy. We demonstrate how the accuracy of the derived values of the optical parameters of aerosols depends on the accuracy of measurements of the intensity and polarization of the reflected light, optical thickness of aerosol layer itself, effective radius of aerosols, width of the particle size distribution, and number of viewing angles.

Aerosols↗

Analysis of simultaneous aerosol and ocean glint retrieval using multi-angle observations

Since early 2000, NASA's Multi-angle Imaging SpectroRadiometer (MISR) instrument has been performing remote sensing retrievals of aerosol optical properties from the polar-orbiting Terra spacecraft. A noteworthy aspect of MISR observations over the ocean is that, for much of the Earth, some of the multi-angle views have contributions from solar reflection by the ocean surface (glint, or glitter), while others do not. Aerosol retrieval algorithms often discard these glint-influenced observations because they can overwhelm the signal and are difficult to predict without knowledge of the (wind-speed-driven) ocean surface roughness. However, theoretical studies have shown that multi-angle observations of a location at geometries with and without reflected sun glint can be a rich source of information, sufficient to support simultaneous retrieval of both the aerosol state and the wind speed at the ocean surface. We are in the early stages of creating such an algorithm. In this paper, we describe our assessment of the appropriate level of parameterization for simultaneous aerosol and ocean surface property retrievals using sun glint. For this purpose, we use generalized nonlinear retrieval analysis (GENRA), an information content assessment (ICA) technique employing Bayesian inference, and simulations from the Ahmad–Fraser iterative radiative transfer code. We find that four parameters are suitable: aerosol optical depth (τ), particle size distribution (expressed as the fine mode fraction f of small particles in a bimodal size distribution), surface wind speed (w), and relative humidity (r, to define the aerosol water content and complex refractive index). None of these parameters define ocean optical properties, as we found that the aerosol state could be retrieved with the nine MISR near-infrared views alone, where the ocean body is strongly absorbing in the open ocean. We also found that retrieval capability varies with observation geometry and that as τ increases so does the ability to determine aerosol intensive optical properties (r and f, while it decreases for w). Increases in w decrease the ability to determine the true value of that parameter but have minimal impact on retrieval of aerosol properties. We explored the benefit of excluding the two most extreme MISR view angles for which radiative transfer with the plane-parallel approximation is less certain, but we found no advantage in doing so. Finally, the impact of treating wind speed as a scalar parameter, rather than as a two-parameter directional wind, was tested. While the simpler scalar model does contribute to overall aerosol uncertainty, it is not sufficiently large to justify the addition of another dimension to parameter space. An algorithm designed upon these principles is in development. It will be used to perform an atmospheric correction with MISR for coincident ocean color (OC) observations by the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument, also on the NASA Terra spacecraft. Unlike MISR, MODIS is a single-view-angle instrument, but it has a more complete set of spectral channels ideal for determination of optical ocean properties. The atmospheric correction of MODIS OC data can therefore benefit from MISR aerosol retrievals. Furthermore, higher-spatial-resolution data from coincident MISR observations may also improve glint screening.

Kirk Knobelspiesse↗

A Practical Guide to Writing a Radiative Transfer Code

Using our decades-long experience in radiative transfer (RT) code development for Earth science, we endeavor to reduce the knowledge gap of bringing RT from theory to code quickly. Despite numerous classic and recent literature, it is still hard to develop anRT code from scratch within a few weeks. It is equally hard to understand, not to mention modify, an existing “monster” RT code, for which the developer is either located remotely or has retired. Following the format of “Numerical Recipes” by Press et al., we collocate in this paper small pieces of necessary theory with corresponding small pieces of RT code. These are arranged in an order that is natural for code development, which is often opposite of the natural order for laying out the theoretical basis. We focus on the transfer of unpolarized monochromatic solar radiation in a plane-parallel atmosphere over a reflecting surface. Both the surface and the atmosphere are homogeneous (uniform) at all directions. The multiple scattering is numerically solved using the deterministic method of Gauss-Seidel iterations. Except for the presented Python-Numba open-source RT code gsit, the paper does not report any new scientific results, but rather serves as an academic demonstration. If development time is an issue or the reader is familiar with basic concepts of RT theory, we recommend proceeding directly to Sec.3 “RT code development.

multiple light scattering↗

Numerical Results for Polarized Light Scattering in a Spherical Atmosphere

We report numerical results for polarized light reflection from the top of a Rayleigh scattering spherical atmosphere with height-dependent single scattering albedo over a dark surface. Michael Mishchenko considered this scenario back in the 1990’s, for a plane-parallel atmosphere of unit optical thickness (OT), for which radiance errors arising from neglecting polarization reaches their highest values. To further extend Mishchenko’s results, we consider a value of OT = 0.25, for which the effect of atmospheric curvature is pronounced. New results are generated using three state-of-the art radiative transfer (RT) codes. These are: the MYSTIC and MCSSA models, which simulate light scattering in a true- spherical atmosphere using Monte Carlo methods; and the discrete ordinate code VLIDORT, operating with a new multiple-scatter spherical correction designed to deliver reasonable approximations to spherical-medium scattering. In this work, we report results for both single and multiple scattering; this will help to support the validation of existing and future polarized spherical RT codes, especially those using approximative methods to deal with sphericity.

S Korkin↗

3D Cloud Tomography and Droplet Size Retrieval from Multi-Angle Polarimetric Imaging of Scattered Sunlight from Above

Tomography aims to recover a three-dimensional (3D) density map of a medium or an object. In medical imaging,it is extensively used for diagnostics via X-ray computed tomography (CT). We define and derive a tomographyof cloud droplet distributions via passive remote sensing. We use multi-view polarimetric images to fit a 3Dpolarized radiative transfer (RT) forward model. Our motivation is 3D volumetric probing of vertically-developedconvectively-driven clouds that are ill-served by current methods in operational passive remote sensing. Currenttechniques are indeed based on strictly 1D RT modeling and applied to a single cloudy pixel, where cloud geometrydefaults to that of a plane-parallel slab. Incident unpolarized sunlight, once scattered by cloud droplets, changesits polarization state according to droplet size. Therefore, polarimetric measurements in the rainbow and gloryangular regions can be used to infer the droplet size distribution. This work defines and derives a framework for afull 3D tomography of cloud droplets for both their mass concentration in space and their distribution across arange of sizes. This gridded 3D retrieval of key microphysical properties is made tractable by our novel approachthat involves a restructuring and partial linearization of an open-source polarized 3D RT code to accommodate aspecial two-step iterative optimization technique. Physically-realistic synthetic clouds are used to demonstrate themethodology with rigorous uncertainty quantification, while a real-world cloud imaged by AirMSPI is processedto illustrate the new remote sensing capability

Schechner, Yoav Y.↗

Radiative interaction of atmosphere and surface: write up with elements of code

In passive satellite remote sensing of the Earth, separation of the path radiance (atmosphere-only contribution) from the surface reflection remains a “significant challenge”. Recent literature names it among the gaps in radiative transfer (RT) topics that “require continued research in the near future”. The challenge comes from multiple reflections (bouncing) between the atmosphere and surface – radiative interaction. In this paper we use a known RT technique, the matrix-operator method (MOM), and a new modification of the monochromatic vector RT (vRT) code IPOL (Intensity and POLarization) to simulate the interaction of a plane-parallel atmosphere and a few widely used surface reflection models. Following the idea of the Green’s function method, IPOL no longer takes the surface model parameters on input. Instead, it provides the path radiance, and the atmospheric reflection and transmission matrices as output. Despite many RT codes use the MOM formalism, this output does not seem common. The surface reflection matrix is computed externally. Therefore, this paper extends the Green’s function atmospheric correction technique to the case of polarized light. Aiming clarity rather than performance, we explain in Python the structure of the surface matrices for the isotropic (Lambertian), directional unpolarized, and polarized ocean reflection models. We then combine these surface matrices and the precomputed IPOL output to get numerically accurate signal at the top of atmosphere (TOA) and test it vs. published benchmarks. Then, for each benchmark scenario we show how to get the surface from the TOA signal, i.e. perform the RT-based atmospheric correction.

radiative transfer↗

Hyper Illumination of Exoplanets: Analytical and Numerical Approaches

This work describes the illumination of exoplanets whose orbits are close enough to their host star that the finite angular size of their host star causes hyper illumination, in which more than 50% of the planet receives light. Such exoplanets include the hot Jupiters KELT-9 b (64.5% illuminated) and Kepler-91 b (69.6% illuminated). We describe the geometry of three primary illumination zones: the fully illuminated zone, penumbral zone, and unilluminated zone. The integrals required to determine the incident radiation as a function of position from the substellar point on the exoplanet are explained and derived, and the analytical solution is presented within the fully illuminated zone. We find that the illumination predicted by our model is greater at the substellar point than the typical plane-parallel ray model used would suggest. In addition, it is greater within the region of the penumbral zone extending into the antistellar side of the exoplanet. Finally, we compare our model to that used in starry, an open-source software package used to create albedo maps. It appears that starry may be overestimating the illumination of closely orbiting exoplanets because the foreshortening of the area element of the host star is not included in its calculation.

Exoplanets↗