Search NASA⌕ Search

SEARCH · Search NASA

Results for “invariant imbeddings”

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.

52 records · Page 3

Number of terms required in the Fourier expansion of the reflection function for optically thick atmospheres

Computational results are presented for the separate terms in the Fourier expansion of the phase function and the reflection function of a semiinfinite, conservatively scattering atmosphere composed of cloud particles. The calculations involve successive applications of invariant imbedding, doubling, and asymptotic fitting methods to cover the range from very thin to very thick atmospheres. From the results, the ratio of the total reflection function to the first-order reflection function is determined as well as the number of terms required to describe the reflection function to an accuracy of 0.1 percent. The number of terms required depends strongly on the zenith angles of incidence and reflection as well as on details of the phase function. These results are compared with similar results obtained for a Henyey-Greenstein phase function having the same asymmetry factor as in the cloud model.

King, M. D.↗

Assessment of the Accuracy of the Conventional Ray-Tracing Technique: Implications in Remote Sensing and Radiative Transfer Involving Ice Clouds.

A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).

remote sensing↗

Ice Cloud Backscatter Study and Comparison with CALIPSO and MODIS Satellite Data

An invariant imbedding T-matrix (II-TM) method is used to calculate the single-scattering properties of 8-column aggregate ice crystals. The II-TM based backscatter values are compared with those calculated by the improved geometric-optics method (IGOM) to refine the backscattering properties of the ice cloud radiative model used in the MODIS Collection 6 cloud optical property product. The integrated attenuated backscatter-to-cloud optical depth (IAB-ICOD) relation is derived from simulations using a CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite) lidar simulator based on a Monte Carlo radiative transfer model. By comparing the simulation results and co-located CALIPSO and MODIS (Moderate Resolution Imaging Spectroradiometer) observations, the non-uniform zonal distribution of ice clouds over ocean is characterized in terms of a mixture of smooth and rough ice particles. The percentage of the smooth particles is approximately 6 percent and 9 percent for tropical and mid-latitude ice clouds, respectively.

MODIS↗

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗

On the Convergence of Numerical Computations for Both Exact and Approximate Solutions for Electromagnetic Scattering by Nonspherical Dielectric Particles

We summarize the size parameter range of the applicability of four lightscattering computational methods for nonspheric dielectric particles. These methods include two exact methods - the extended boundary condition method (EBCM) and the invariant imbedding T-matrix method (II-TM) and two approximate approaches - the physical-geometric optics method (PGOM) and the improved geometric optics method (IGOM). For spheroids, the single-scattering properties computed by EBCM and II-TM agree for size parameters up to 150, and the comparison gives us confidence in using IITM as a benchmark for size parameters up to 150 for other geometries (e.g., hexagonal columns) because the applicability of II-TM with respect to particle shape is generic, as demonstrated in our previous studies involving a complex aggregate. This study demonstrates the convergence of the exact II-TM and approximate PGOM solutions for the complete set of single-scattering properties of a nonspherical shape other than spheroids and circular cylinders with particle sizes of ~48λ, specifically a hexagonal column with a size parameter of length as kL=300, where k=2π/λ and L is the column length. IGOM is also quite accurate except near the exact 180º backscattering direction. This study demonstrates that a synergetic combination of the numerically-exact II-TM and the approximate PGOM can seamlessly cover the entire size parameter range of practical interest. To demonstrate the applicability of the approach, we compute the optical properties of dust particles with a downstream application to the retrieval of dust aerosol optical thickness and effective particle size from polarimetric observations.

Ping Yang↗

Massively Parallel Algorithms for Real-Time Wavefront Control of a Dense Adaptive Optics System

In this paper massively parallel algorithms and architectures for real-time wavefront control of a dense adaptive optic system (SELENE) are presented. We have already shown that the computation of a near optimal control algorithm for SELENE can be reduced to the solution of a discrete Poisson equation on a regular domain. Although this represents an optimal computation, due the large size of the system and the high sampling rate requirement, the implementation of this control algorithm poses a computationally challenging problem since it demands a sustained computational throughput of the order of 10 GFlops. We develop a novel algorithm, designated as Fast Invariant Imbedding algorithm, which offers a massive degree of parallelism with simple communication and synchronization requirements. Due to these features, our algorithm is significantly more efficient than other Fast Poisson Solvers for implementation on massively parallel architectures.

massively↗

Future ultra-speed tube-flight

Future long-link, ultra-speed, surface transport systems will require electromagnetically (EM) driven and restrained vehicles operating under reduced-atmosphere in very straight tubes. Such tube-flight trains will be safe, energy conservative, pollution-free, and in a protected environment. Hypersonic (and even hyperballistic) speeds are theoretically achievable. Ultimate system choices will represent tradeoffs between amoritized capital costs (ACC) and operating costs. For example, long coasting links might employ aerodynamic lift coupled with EM restraint and drag make-up. Optimized, combined EM lift, and thrust vectors could reduce energy costs but at increased ACC. (Repulsive levitation can produce lift-over-drag l/d ratios a decade greater than aerodynamic), Alternatively, vehicle-emanated, induced-mirror fields in a conducting (aluminum sheet) road bed could reduce ACC but at substantial energy costs. Ultra-speed tube flight will demand fast-acting, high-precision sensors and computerized magnetic shimming. This same control system can maintain a magnetic 'guide way' invariant in inertial space with inertial detectors imbedded in tube structures to sense and correct for earth tremors. Ultra-speed tube flight can complete with aircraft for transit time and can provide even greater passenger convenience by single-model connections with local subways and feeder lines. Although cargo transport generally will not need to be performed at ultra speeds, such speeds may well be desirable for high throughput to optimize channel costs. Thus, a large and expensive pipeline might be replaced with small EM-driven pallets at high speeds.

Salter, Robert M.↗

Manipulator interactive design with interconnected flexible elements

This paper describes the development of an analysis tool for the interactive design of control systems for manipulators and similar electro-mechanical systems amenable to representation as structures in a topological chain. The chain consists of a series of elastic bodies subject to small deformations and arbitrary displacements. The bodies are connected by hinges which permit kinematic constraints, control, or relative motion with six degrees of freedom. The equations of motion for the chain configuration are derived via Kane's method, extended for application to interconnected flexible bodies with time-varying boundary conditions. A corresponding set of modal coordinates has been selected. The motion equations are imbedded within a simulation that transforms the vector-dyadic equations into scalar form for numerical integration. The simulation also includes a linear, time-invariant controler specified in transfer function format and a set of sensors and actuators that interface between the structure and controller. The simulation is driven by an interactive set-up program resulting in an easy-to-use analysis tool.

Singh, R. P.↗

Forty-Year Calibrated Record of Earth-Surface Reflected Radiance from Landsat: A Review

Sensors on Landsat satellites have been collecting images of the Earth's surface for nearly 40 years. These images have been invaluable for characterizing and detecting changes in the land cover and land use of the world. Although initially conceived as primarily picture generating sensors, even the early sensors were radiometrically calibrated and spectrally characterized prior to launch and incorporated some capabilities to monitor their radiometric calibration once on orbit. Recently, as the focus of studies has shifted to monitoring Earth surface parameters over significant periods of time, serious attention has been focused toward bringing the data from all these sensors onto a common radiometric scale over this 40-year period. This effort started with the most recent systems and then was extended back in time. Landsat-7 ETM+, the best-characterized sensor of the series prior to launch and once on orbit, and the most stable system to date, was chosen to serve as the reference. The Landsat-7 project was the first of the series to build an image assessment system into its ground system, allowing systematic characterization of its sensors and data. Second, the Landsat-5 TM (still operating at the time of the Landsat-7 launch and continues to operate) calibration history was reconstructed based on its internal calibrator, vicarious calibrations, pseudo-invariant sites and a tie to Landsat-7 ETM+ at the time of the commissioning of Landsat-7. This process was performed in two iterations: the earlier one relied primarily on the TM internal calibrator. When this was found to have some deficiencies, a revised calibration was based more on pseudo-invariant sites, though the internal calibrator was still used to establish the short-term variations in response due to icing build up on the cold focal plane. As time progressed, a capability to monitor the Landsat-5 TM was added to the image assessment system. The Landsat-4 TM, which operated from 1982-1992, was the third system to which the radiometric scale was extended. The limited and broken use of the Landsat-4 TM made this analysis more difficult. Eight-day separated image pairs from Landsat-5 combined with analysis of pseudo invariant sites established this history. The fourth and most challenging effort was making the Landsat-1 to -5 MSS sensors' data internally radiometrically consistent. This effort was particularly complicated by the age of the MSS data, varying formats and processing levels in the archive, limited datasets, and limited documentation available. Ultimately, pseudo-invariant sites were identified in North America and used for this effort. Note that most of the Landsat-MSS archived data had already been calibrated using the MSS internal calibrators, so this processing was imbedded in the result. The final effort was developing an absolute scale for Landsat MSS similar to what was already established for the "TM" sensors. This was accomplished by using simultaneous data from Landsat-5 MSS and Landsat-5 TM, accounting for spectral differences between the sensors using EO-1 Hyperion data. The recalibrated history of the Landsat data and implications to users are discussed. The key result from this work is a consistently calibrated Landsat data archive that spans nearly 40 years with total uncertainties on the order of 10% or less for most sensors and bands.

Markham, Brian↗

Calculation of three-dimensional turbulent subsonic flows in transition ducts

A method for computing three-dimensional turbulent subsonic flow in curved ducts is being developed. A set of tube-like surface oriented coordinates is employed for a general class of geometries applicable to subsonic diffusers with offset bends. The geometric formulation is complex and no previous treatment of this class of viscous flow problems is known to the authors. The duct centerline is a space curve specified by piecewise polynomials. A Frenet frame is located on the centerline at each axial location. The cross sections are described by superellipses imbedded in the Frenet frame. Duct surfaces are also coordinate surfaces, which greatly simplifies the boundary conditions. The resulting coordinates are nonorthogonal. An approximate set of governing equations is employed for viscous flows having strong flow in a primary flow direction. The derivation is coordinate invariant and the resulting equations are expressed in tensor form. These equations are solved by an efficient alternating direction implicit (ADI) method. This numerical method is generally stable and permits solution in difficult geometries using the general tensor formulation.

Levy, R.↗

On the equation of state for an electron gas in an intense magnetic field

In this paper we derive the equation of state for a relativistic electron gas imbedded in a static homogeneous magnetic field of arbitrary strength. The derivation is based on the evaluation of the energy-momentum tensor and the use of Dirac's equation for such a problem. Contrary to a derivation presented several years ago, the present derivation is completely gauge-invariant. We also show how to recover, in an exact manner, the perfect gas law for the case of weak magnetic fields.

Canuto, V.↗

Modern control concepts in hydrology

Two approaches to an identification problem in hydrology are presented based upon concepts from modern control and estimation theory. The first approach treats the identification of unknown parameters in a hydrologic system subject to noisy inputs as an adaptive linear stochastic control problem; the second approach alters the model equation to account for the random part in the inputs, and then uses a nonlinear estimation scheme to estimate the unknown parameters. Both approaches use state-space concepts. The identification schemes are sequential and adaptive and can handle either time invariant or time dependent parameters. They are used to identify parameters in the Prasad model of rainfall-runoff. The results obtained are encouraging and conform with results from two previous studies; the first using numerical integration of the model equation along with a trial-and-error procedure, and the second, by using a quasi-linearization technique. The proposed approaches offer a systematic way of analyzing the rainfall-runoff process when the input data are imbedded in noise.

Duong, N.↗

Modern control concepts in hydrology

Two approaches to an identification problem in hydrology are presented, based upon concepts from modern control and estimation theory. The first approach treats the identification of unknown parameters in a hydrologic system subject to noisy inputs as an adaptive linear stochastic control problem; the second approach alters the model equation to account for the random part in the inputs, and then uses a nonlinear estimation scheme to estimate the unknown parameters. Both approaches use state-space concepts. The identification schemes are sequential and adaptive and can handle either time-invariant or time-dependent parameters. They are used to identify parameters in the Prasad model of rainfall-runoff. The results obtained are encouraging and confirm the results from two previous studies; the first using numerical integration of the model equation along with a trial-and-error procedure, and the second using a quasi-linearization technique. The proposed approaches offer a systematic way of analyzing the rainfall-runoff process when the input data are imbedded in noise.

Duong, N.↗