SEARCH · Engineering Papers

Results for “Empirical characteristic function”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Machine Learning Prediction of H Adsorption Energies on Ag Alloys

Adsorption energies on surfaces are excellent descriptors of their chemical properties, including their catalytic performance. This has prompted a significant amount of work towards efficiently predicting adsorption energies without performing expensive quantum chemical calculations, based on semi-empirical or data scientific models. Here, we develop machine learning models to predict H adsorption energies on stepped Ag surface alloys using a database of energies obtained from first-principles electronic structure calculations based on density functional theory. Starting with a simple model and iteratively refining the features extracted from the calculations, we determine the key characteristics of successful models, and reduce the number of features required for accurate prediction of the adsorption energy. We find that nonlinear models are essential for accurate predictions. Moreover, detailed comparison of model predictions and calculations reveal a counterintuitive "selfpassivation" induced by some alloying elements, as well as unexpected discontinuous electronic structure effects. These findings are rationalized using tight-binding model hamiltonians and by analyzing the electronic density of states. Lastly, we show that we can estimate the reliability of a particular prediction, mitigating the effects of errors in the small fraction of cases where these errors are significant.

36 MATERIALS SCIENCE↗

NASA NTRS1978

Dual-loop model of the human controller

A dual-loop model of the human controller in single-axis compensatory tracking tasks is introduced. This model possesses an inner-loop closure that involves feeding back that portion of controlled element output rate that is due to control activity. A novel feature of the model is the explicit appearance of the human's internal representation of the manipulator-controlled element dynamics in the inner loop. The sensor inputs to the human controller are assumed to be system error and control force. The former can be sensed via visual, aural, or tactile displays, whereas the latter is assumed to be sensed in kinesthetic fashion. A set of general adaptive characteristics for the model is hypothesized, including a method for selecting simplified internal models of the manipulator-controlled element dynamics. It is demonstrated that the model can produce controller describing functions that closely approximate those measured in four laboratory tracking tasks in which the controlled element dynamics vary considerably in terms of ease of control. An empirically derived expression for the normalized injected error remnant spectrum is introduced.

Hess, R. A.↗

NASA NTRS1986

Nitrogen uptake and utilization by intact plants

The results of experiments support the proposed conceptual model that relates nitrogen uptake activity by plants as a balanced interdependence between the carbon-supplying function of the shoot and the nitrogen-supplying function of the roots. The data are being used to modify a dynamic simulation of plant growth, which presently describes carbon flows through the plant, to describe nitrogen uptake and assimilation within the plant system. Although several models have been proposed to predict nitrogen uptake and partitioning, they emphasize root characteristics affecting nutrient uptake and relay on empirical methods to describe the relationship between nitrogen and carbon flows within the plant. Researchers, on the other hand, propose to continue to attempt a mechanistic solution in which the effects of environment on nitrogen (as well as carbon) assimilation are incorporated through their direct effects on photosynthesis, respiration, and aging processes.

Raper, C. D., Jr.↗

NASA NTRS1993

Development and application of an empirical probability distribution for the prediction error of re-entry body maximum dynamic pressure

The relationship between actual and predicted re-entry maximum dynamic pressure is characterized using a probability density function and a cumulative distribution function derived from sounding rocket flight data. This paper explores the properties of this distribution and demonstrates applications of this data with observed sounding rocket re-entry body damage characteristics to assess probabilities of sustaining various levels of heating damage. The results from this paper effectively bridge the gap existing in sounding rocket reentry analysis between the known damage level/flight environment relationships and the predicted flight environment.

Lanzi, R. James↗

DOE OSTI2016

Particle Image Velocimetry Measurements in a Wire-Wrapped 61-Pin Hexagonal Fuel Bundle

Liquid metal fast reactor fuel assemblies using sodium as a coolant typically utilize a tightly packed triangular lattice of fuel pins enclosed in a hexagonal duct. The fuel pins are typically wrapped with a helical wire spacer. The primary function of the wire spacers is to maintain a gap between neighboring fuel pins. They also mitigate vortex-induced vibration and enhance sub-channel mixing to increase convective heat transfer. Many experiments have been performed to study thermal-hydraulic characteristics, primarily pressure drops, of the wire-wrapped fuel bundles. Novendstern developed a semi-empirical model to predict pressure losses in a hexagonal array of wire-wrapped pins in the turbulent flow regime. Rehme proposed pressure drop correlations for a wide range of geometrical parameters, such as pitch-to-diameter ratio P/D from 1.125 to 1.417, lead-to-wire diameter ratio H / dm from 6 to 45, and number of pins n from 7 to 61. Cheng and Todreas introduced their hydrodynamic models for subchannel friction factors and mixing parameters in a 37-pin fuel bundle, which covered the laminar, transition and turbulent flow regimes. Computational fluid dynamics (CFD) offers an opportunity to validate the above correlations and investigate the flow characteristics in the subchannels. Several numerical studies of the wire-wrapped fuel bundles have been performed and most of them are Reynolds-Averaged Navier-Stokes (RANS) based simulations with various turbulence models. Ahmad and Kim performed RANS simulations of 7- and 19-pin bundles using the k-omega SST turbulent model, while Gajapathy et al. performed RANS simulations on 7-pin bundles using the k-epsilons model. Recently, Fischer et al. carried out large-eddy simulation (LES) calculations on wire-wrapped fuel bundles to obtain higher-resolution numerical solutions of subchannel behavior. The geometrical complexity of the hexagonal wire-wrapped fuel bundle has limited experimental activities to pressure drop measurements. A database of pressure drops and flow-field measurements in a 61-pin wire-wrapped hexagonal fuel bundle is not available to benchmark the existing correlations and validate the CFD calculations. Texas A and M University is conducting isothermal flow experiments in a wire-wrapped 61-pin hexagonal fuel bundle to support the research on advanced nuclear fuel development sponsored by the Department of Energy (DOE). The purpose of these tests is to perform high spatial and temporal resolution measurements of the flow and pressure measurements at different locations in the fuel bundle. The experimental facility employs matched-index-of-refraction (MIR) techniques and laser diagnostic velocity measurement techniques such as particle image velocimetry (PIV) and laser Doppler velocimetry (LDV). The result of this experimental activity will be an experimental database of pressure and flow-field measurements that will be suitable for validating CFD computer codes. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2010

Global Patterns of Lightning Properties Derived by LIS

The Lightning Imaging Sensor LIS aboard the TRMM satellite provides unmatched empirical data of the global lightning distribution (up to approx.35deg S/N) since end of 1997. Climatological flash rate densities derived from LIS are standard references, e.g. for flash rate parameterizations used in GCMs. It is known that flash characteristics are quite variable, and that various quantities (like the flash energy or the NOx production per flash) vary considerably, statistically as well as systematically on regional and seasonal scales. LIS provides information beyond flash counts, in particular radiance and flash footprint. Here we present an analysis of global patterns of various lightning properties derived from LIS, in relation to the number of flashes. These normalized flash characteristics show consistent spatial patterns of regions with "strong" versus regions with "weak" lightning. Most striking is a clear land-ocean contrast, with oceanic flashes being "stronger" than continental flashes. But also over continents, flash strength shows systematic variations. Highest continental values are found over the US, while values over South America and India are quite low. These regional variations cannot be simply parameterized as function of latitude. Information on spatial patterns of mean flash "strength", though rather qualitative up to now, is potentially a valuable input for improving empirical parameterizations based on flash counts (like precipitation or lightning NOx). Further investigation is in progress to come to a more physical and quantitative understanding of the spatial patterns of the different LIS properties. In particular, it has to be checked how far they could be related to established lightning properties (like energy or the fraction of intra-cloud to cloud-to-ground flashes) or to meteorological quantities (like CAPE).

Beirle, S.↗

NASA NTRS2013

The Galex Time Domain Survey. I. Selection And Classification of Over a Thousand Ultraviolet Variable Sources

We present the selection and classification of over a thousand ultraviolet (UV) variable sources discovered in approximately 40 deg(exp 2) of GALEX Time Domain Survey (TDS) NUV images observed with a cadence of 2 days and a baseline of observations of approximately 3 years. The GALEX TDS fields were designed to be in spatial and temporal coordination with the Pan-STARRS1 Medium Deep Survey, which provides deep optical imaging and simultaneous optical transient detections via image differencing. We characterize the GALEX photometric errors empirically as a function of mean magnitude, and select sources that vary at the 5 sigma level in at least one epoch. We measure the statistical properties of the UV variability, including the structure function on timescales of days and years. We report classifications for the GALEX TDS sample using a combination of optical host colors and morphology, UV light curve characteristics, and matches to archival X-ray, and spectroscopy catalogs. We classify 62% of the sources as active galaxies (358 quasars and 305 active galactic nuclei), and 10% as variable stars (including 37 RR Lyrae, 53 M dwarf flare stars, and 2 cataclysmic variables). We detect a large-amplitude tail in the UV variability distribution for M-dwarf flare stars and RR Lyrae, reaching up to absolute value(Δm) = 4.6 mag and 2.9 mag, respectively. The mean amplitude of the structure function for quasars on year timescales is five times larger than observed at optical wavelengths. The remaining unclassified sources include UV-bright extragalactic transients, two of which have been spectroscopically confirmed to be a young core-collapse supernova and a flare from the tidal disruption of a star by dormant supermassive black hole. We calculate a surface density for variable sources in the UV with NUV less than 23 mag and absolute value(Δm) greater than 0.2 mag of approximately 8.0, 7.7, and 1.8 deg(exp −2) for quasars, active galactic nuclei, and RR Lyrae stars, respectively. We also calculate a surface density rate in the UV for transient sources, using the effective survey time at the cadence appropriate to each class, of approximately 15 and 52 deg(exp −2 yr−1 for M dwarfs and extragalactic transients, respectively.

Gezari, S.↗

DOE OSTI2023

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗

NASA NTRS1981

Asymptotic form of the longitudinal correlation function for isotropic homogeneous turbulence

An asymptotic form is derived for the longitudinal correlation function for isotropic homogeneous turbulence in an incompressible fluid governed by the Navier-Stokes equation. The result is obtained from an analysis of the algebraic-differential structure of the two-point correlation tensor contained in the complex-valued Fourier transform of the probability measure over the turbulence ensemble, which reveals that the Hopf characteristic functional (the complex-valued Fourier transform of the probability measure) satisfies the Fourier interference inequality. Consequences of the expression obtained, in which the correlation function is a positive definite function of the inverse cube of the spatial coordinate as it approaches infinity, are shown to include the nonexistence of the Loitsianskii invariant, and the solution is shown to be consistent with empirical formulas.

Rosen, G.↗

NASA NTRS2011

Canopy Spectral Invariants: Application to Classification of Forest Types from Hyperspectral Data - Part 2

Many studies have been conducted to demonstrate the ability of hyperspectral data to discriminate plant dominant species. Most of them have employed the use of empirically based techniques, which are site specific, requires some initial training based on characteristics of known leaf and/or canopy spectra and therefore may not be extendable to operational use or adapted to changing or unknown land cover. In this paper we propose a physically based approach for separation of dominant forest type using hyperspectral data. The radiative transfer theory of canopy spectral invariants underlies the approach, which facilitates parameterization of the canopy reflectance in terms of the leaf spectral scattering and two spectrally invariant and structurally varying variables - recollision and directional escape probabilities. The methodology is based on the idea of retrieving spectrally invariant parameters from hyperspectral data first, and then relating their values to structural characteristics of three-dimensional canopy structure. Theoretical and empirical analyses of ground and airborne data acquired by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) over two sites in New England, USA, suggest that the canopy spectral invariants convey information about canopy structure at both the macro- and micro-scales. The total escape probability (one minus recollision probability) varies as a power function with the exponent related to the number of nested hierarchical levels present in the pixel. Its base is a geometrical mean of the local total escape probabilities and accounts for the cumulative effect of canopy structure over a wide range of scales. The ratio of the directional to the total escape probability becomes independent of the number of hierarchical levels and is a function of the canopy structure at the macro-scale such as tree spatial distribution, crown shape and size, within-crown foliage density and ground cover. These properties allow for the natural separation of dominant forest classes based on the location of points on the total escape probability vs the ratio log-log plane.

Schull, M. A.↗

NASA NTRS1975

Radar return from a continuous vegetation canopy

The radar backscatter coefficient, sigma deg, of alfalfa was investigated as a function of both radar parameters and the physical characteristics of the alfalfa canopy. Measurements were acquired with an 8-18 GHz FM-CW mobile radar over an angular range of 0 - 70 deg as measured from nadir. The experimental data indicates that the excursions of sigma deg at nadir cover a range of nearly 18 dB during one complete growing cycle. An empirical model for sigma deg was developed which accounts for its variability in terms of soil moisture, plant moisture and plant height.

Bush, T. F.↗

NASA NTRS1991

Macromolecular recognition: Structural aspects of the origin of the genetic system

Theoretical simulation of prebiotic chemical processes is an invaluable tool for probing the phenomenon of evolution of life. Using computational and modeling techniques and guided by analogies from present day systems we, seek to understand the emergence of genetic apparatus, enzymatic catalysis and protein synthesis under prebiotic conditions. In one possible scenario, the RNA enzymatic reaction plays a key role in the emergence of the self-replicating and offers a clue to the onset of enzymatic catalysis prior to the existence of the protein biosynthetic machinery. Our ultimate goal is to propose a simple RNA segment which contains the specificity and catalytic activity of the contemporary RNA enzyme and which could emerge in a primordial chemical environment. To understand the mechanism of ribozyme catalyzed reactions, ab initio and semi-empirical (ZINDO) programs were used to investigate the reaction path of transphosphorylation. A special emphasis was placed on the possible catalytic and structural roles played by the coordinated magnesium cation. Both the inline and adjacent mechanisms of transphosphorylation have been studied. Another important aspect of this reaction is the identity of the functional groups which are essential for the acid base catalysis. The structural characteristics of the target helices, particularly a possible role of G center dot T pair, is under examination by molecular dynamics (MD) simulation technique. Modeling of the ancestral aminoacyl-tRNA synthetases (aRS) may provide important clues to the emergence of the genetic code and the protein synthetic machinery. Assuming that the catalytic function evolved before the elements of specific recognition of a particular amino acid, we are exploring the minimal structural requirements for the catalysis of tRNA aminoacylation. The molecular modeling system SYBYL was used for this study based on the high resolution crystallographic structures of the present day tyrosyl-adenylate:tyrRS and tRNA(Gln): ATP:glnRS complexes. The trinucleotide CCA of the 3'-end tRNA is placed into the active site pocket of tyrRS, based on the scheme of interaction between tRNA(Gln) and glnRS, and upon the stereochemistry of the tyrRS:tRNA:Tyr-AMP transition state. This provides a model of the non-specific recognition of a tRNA's 3'-end by an aRS, which might be similar to that of the ancestral aRS's. In the next step, modeling of the rest of the acceptor stem of tRNA (Tyr) with tyrRS is carried out.

Rein, Robert↗

NASA NTRS2001

Shock Formation of Slow Magnetosonic Waves in Coronal Plumes

We investigate the height of shock formation in coroner plumes for slow magnetosonic waves. The models take into account plume geometric spreading, heat conduction and radiative damping. The wave parameters as well as the spreading functions of the plumes and the base magnetic field strengths are given by empirical constraints mostly from Solar and Heliospheric Observatory/Ultraviolet Coronagraph Spectrometer (SOHO/UVCS). Our models show that shock formation occurs at low coronal heights, i.e., within 1.3 solar radius, depending on the model parameters. The shock formation is calculated using the well-established wave breaking condition given by the intersection of C+ characteristics in the space-time plane. Our models show that shock heating by slow magnetosonic waves is expected to be relevant at most heights in solar coronal plumes, although slow magnetosonic waves are most likely not a solely operating energy supply mechanism.

Cuntz, Manfred↗

DOE OSTI2023

A data-driven method for modelling dissipation rates in stratified turbulence

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$ . Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$ , where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$ , and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$ ), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

42 ENGINEERING↗

DOE OSTI2014

Spitzer ultra faint survey program (surfs up). I. An overview

Spitzer UltRa Faint SUrvey Program is a joint Spitzer and Hubble Space Telescope Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z ≳ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest-frame optical light, which only Spitzer can probe for sources at z ≳ 7, for a large enough sample of typical galaxies. Our program consists of 550 hr of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ∼30 hr exposure time in both 3.6 μm and 4.5 μm in the central 4' × 4' field and ∼15 hr in the flanking fields. This results in 3σ sensitivity limits of ∼26.6 and ∼26.2 AB magnitudes for the central field in the IRAC 3.6 and 4.5 μm bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z ≳ 7 sources (using a z = 9.5 galaxy behind MACS J1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical point-spread function models. In the future we plan to release these data products for the entire survey.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2017

Modeling for Battery Prognostics

For any battery-powered vehicles (be it unmanned aerial vehicles, small passenger aircraft, or assets in exoplanetary operations) to operate at maximum efficiency and reliability, it is critical to monitor battery health as well performance and to predict end of discharge (EOD) and end of useful life (EOL). To fulfil these needs, it is important to capture the battery's inherent characteristics as well as operational knowledge in the form of models that can be used by monitoring, diagnostic, and prognostic algorithms. Several battery modeling methodologies have been developed in last few years as the understanding of underlying electrochemical mechanics has been advancing. The models can generally be classified as empirical models, electrochemical engineering models, multi-physics models, and molecular/atomist. Empirical models are based on fitting certain functions to past experimental data, without making use of any physicochemical principles. Electrical circuit equivalent models are an example of such empirical models. Electrochemical engineering models are typically continuum models that include electrochemical kinetics and transport phenomena. Each model has its advantages and disadvantages. The former type of model has the advantage of being computationally efficient, but has limited accuracy and robustness, due to the approximations used in developed model, and as a result of such approximations, cannot represent aging well. The latter type of model has the advantage of being very accurate, but is often computationally inefficient, having to solve complex sets of partial differential equations, and thus not suited well for online prognostic applications. In addition both multi-physics and atomist models are computationally expensive hence are even less suited to online application An electrochemistry-based model of Li-ion batteries has been developed, that captures crucial electrochemical processes, captures effects of aging, is computationally efficient, and is of suitable accuracy for reliable EOD prediction in a variety of operational profiles. The model can be considered an electrochemical engineering model, but unlike most such models found in the literature, certain approximations are done that allow to retain computational efficiency for online implementation of the model. Although the focus here is on Li-ion batteries, the model is quite general and can be applied to different chemistries through a change of model parameter values. Progress on model development, providing model validation results and EOD prediction results is being presented.

Prognostics↗

DOE OSTI2014

Three-region characteristic temperature in p-doped quantum dot lasers

We have investigated the temperature dependence of threshold in p-doped 1.3 μm InAs/GaAs quantum dot (QD) lasers with ten layers of QDs in the active region. It is found that the dependence of threshold current density on the temperature within the temperature range from 10 to 90 °C can be divided into three regions by its characteristic temperature (T{sub 0}): negative, infinite, and positive T{sub 0} regions. Furthermore, the T{sub 0} region width is dependent on the cavity length: the longer cavity length of the QD lasers correspondingly the wider T{sub 0} region. Additionally, for the broad area laser, the threshold modal gains of the lasers with different cavity lengths can be fitted by an empirical expression as a function of the threshold current density, when at the temperatures of 30, 50, and 70 °C. We find that the transparency current density (J{sub tr}) remains almost unchanged under different temperatures according to the extracted parameters from these fitted results, which indicates that J{sub tr} plays an important role in balancing the T{sub 0} between negative region and positive one.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS1999

The Prospect for Remote Sensing of Cirrus Clouds with a Submillimeter-Wave Spectrometer

Given the substantial radiative effects of cirrus clouds and the need to validate cirrus cloud mass in climate models, it is important to measure the global distribution of cirrus properties with satellite remote sensing. Existing cirrus remote sensing techniques, such as solar reflectance methods, measure cirrus ice water path (IWP) rather indirectly and with limited accuracy. Submillimeter/wave radiometry is an independent method of cirrus remote sensing based on ice particles scattering the upwelling radiance emitted by the lower atmosphere. A new aircraft instrument, the Far Infrared Sensor for Cirrus (FIRSC), is described. The FIRSC employs a Fourier Transform Spectrometer (FTS). which measures the upwelling radiance across the whole submillimeter region (0.1 1.0-mm wavelength). This wide spectral coverage gives high sensitivity to most cirrus particle sizes and allows accurate determination of the characteristic particle size. Radiative transfer modeling is performed to analyze the capabilities of the submillimeter FTS technique. A linear inversion analysis is done to show that cirrus IWP, particle size, and upper-tropospheric temperature and water vapor may be accurately measured, A nonlinear statistical algorithm is developed using a database of 20000 spectra simulated by randomly varying most relevant cirrus and atmospheric parameters. An empirical orthogonal function analysis reduces the 500-point spectrum (20 - 70/cm) to 15 "pseudo-channels" that are then input to a neural network to retrieve cirrus IWP and median particle diameter. A Monte Carlo accuracy study is performed with simulated spectra having realistic noise. The retrieval errors are low for IWP (rms less than a factor of 1.5) and for particle sizes (rins less than 30%) for IWP greater than 5 g/sq m and a wide range of median particle sizes. This detailed modeling indicates that there is good potential to accurately measure cirrus properties with a submillimeter FTS.

Evans, K. Franklin↗