Search NASA⌕ Search

SEARCH · Search NASA

Results for “ORTHOGONAL 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.

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

At least 541 records · Page 30

Fitting a Two-Component Scattering Model to Polarimetric SAR Data from Forests

Two simple scattering mechanisms are fitted to polarimetric synthetic aperture radar (SAR) observations of forests. The mechanisms are canopy scatter from a reciprocal medium with azimuthal symmetry and a ground scatter term that can represent double-bounce scatter from a pair of orthogonal surfaces with different dielectric constants or Bragg scatter from a moderately rough surface, which is seen through a layer of vertically oriented scatterers. The model is shown to represent the behavior of polarimetric backscatter from a tropical forest and two temperate forest sites by applying it to data from the National Aeronautic and Space Agency/Jet Propulsion Laboratory's Airborne SAR (AIRSAR) system. Scattering contributions from the two basic scattering mechanisms are estimated for clusters of pixels in polarimetric SAR images. The solution involves the estimation of four parameters from four separate equations. This model fit approach is justified as a simplification of more complicated scattering models, which require many inputs to solve the forward scattering problem. The model is used to develop an understanding of the ground-trunk double-bounce scattering that is present in the data, which is seen to vary considerably as a function of incidence angle. Two parameters in the model fit appear to exhibit sensitivity to vegetation canopy structure, which is worth further exploration. Results from the model fit for the ground scattering term are compared with estimates from a forward model and shown to be in good agreement. The behavior of the scattering from the ground-trunk interaction is consistent with the presence of a pseudo-Brewster angle effect for the air-trunk scattering interface. If the Brewster angle is known, it is possible to directly estimate the real part of the dielectric constant of the trunks, a key variable in forward modeling of backscatter from forests. It is also shown how, with a priori knowledge of the forest height, an estimate for the attenuation coefficient of the canopy can be obtained directly from the multi-incidence-angle polarimetric observations. This attenuation coefficient is another key variable in forward models and is generally related to the canopy density.

temperate forests↗

JIGSAW-GEO (1.0): Locally Orthogonal Staggered Unstructured Grid Generation for General Circulation Modelling on the Sphere

An algorithm for the generation of non-uniform, locally orthogonal staggered unstructured spheroidal grids is described. This technique is designed to generate very high-quality staggered VoronoiDelaunay meshes appropriate for general circulation modelling on the sphere, including applications to atmospheric simulation, ocean-modelling and numerical weather prediction. Using a recently developed Frontal-Delaunay refinement technique, a method for the construction of high-quality unstructured spheroidal Delaunay triangulations is introduced. A locally orthogonal polygonal grid, derived from the associated Voronoi diagram, is computed as the staggered dual. It is shown that use of the Frontal-Delaunay refinement technique allows for the generation of very high-quality unstructured triangulations, satisfying a priori bounds on element size and shape. Grid quality is further improved through the application of hill-climbing-type optimisation techniques. Overall, the algorithm is shown to produce grids with very high element quality and smooth grading characteristics, while imposing relatively low computational expense. A selection of uniform and non-uniform spheroidal grids appropriate for high-resolution, multi-scale general circulation modelling are presented. These grids are shown to satisfy the geometric constraints associated with contemporary unstructured C-grid-type finite-volume models, including the Model for Prediction Across Scales (MPAS-O). The use of user-defined mesh-spacing functions to generate smoothly graded, non-uniform grids for multi-resolution-type studies is discussed in detail.

geophysical fluid dynamics↗

Analysis of radar images of the active volcanic zone at Krafla, Iceland: The effects of look azimuth biasing

The geomorphic expression of Mid-Ocean-Ridge (MOR) volcanism in a subaerial setting occurs uniquely on Earth in Iceland, and the most recent MOR eruptive activity has been concentrated in the Northeastern Volcanic Zone in an area known as Krafla. Within the Krafla region are many of the key morphologic elements of MOR-related basaltic volcanism, as well as volcanic explosion craters, subglacial lava shields, tectonic fissure swarms known as gjar, and basaltic-andesite flows with well developed ogives (pressure-ridges). The objective was to quantify the degree to which the basic volcanic and structural features can be mapped from directional SAR imagery as a function of the look azimuth. To accomplish this, the current expression of volcanic and tectonic constructs was independently mapped within the Krafla region on the E, W, and N-looking SAR images, as well as from SPOT Panchromatic imagery acquired in 1987. The initial observations of the E, W, and N images indicates that fresh a'a lava surfaces are extremely radar bright (rough at 3 cm to meter scales) independent of look direction; this suggests that these flows do not have strong flow direction related structures at meter and cm scales, which is consistent with typical Icelandic a'a lava surfaces in general. The basic impression from a preliminary analysis of the effects of look azimuth biasing on interpretation of the geology of an active MOR volcanic zone is that up to 30 percent of the diagnostic features can be missed at any given look direction, but that having two orthogonal look direction images is probably sufficient to prevent gross misinterpretation.

Garvin, J. B.↗

Mapping Aerosol Lidar Ratios Over Ocean for CALIPSO Using Modis AOD-Constrained Retrievals and A Global Aerosol Model

The CALIPSO aerosol algorithms currently assign one lidar ratio value for each of the seven tropospheric aerosol types. In a future data products release, the CALIPSO project aims to improve these algorithms by developing regional and seasonal lidar ratio climatologies. In this study, aerosol lidar ratios are inferred through Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) backscatter profiles constrained by collocated aerosol optical depth from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) datasets. The analysis is subsampled for only those profiles that are cloud-free and contain one aerosol type (based on CALIOP feature classification). In addition, the CALIOP profiles are collocated with aerosol volume fractions simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model. In this talk, the twelve-year (2006-2017) mean spatial distributions of inferred aerosol lidar ratios for CALIOP-classified marine aerosols and the corresponding modeled sea salt volume fractions (SSVF) will be shown. Model-assisted climatological lidar ratio maps on a seasonal scale will also be presented, developed from the empirical relationship found between the modeled SSVF and lidar ratios. A comparison with past studies will be provided as well as results of a sensitivity study regarding the variability of retrieved lidar ratios as a function of horizontal averaging resolution. While the majority of this talk will focus on lidar ratios for CALIOP-classified marine aerosols, preliminary results will be shown for other aerosol types over ocean, such as dust and elevated smoke.

Travis Toth↗

Spatial profiling of the interplay between cell type- and vision-dependent transcriptomic programs in the visual cortex

How early sensory experience during “critical periods” of postnatal life affects the organization of the mammalian neocortex at the resolution of neuronal cell types is poorly understood. We previously reported that the functional and molecular profiles of layer 2/3 (L2/3) cell types in the primary visual cortex (V1) are vision-dependent [S. Chenget al.,Cell185, 311–327.e24 (2022)]. Here, we characterize the spatial organization of L2/3 cell types with and without visual experience. Spatial transcriptomic profiling based on 500 genes recapitulates the zonation of L2/3 cell types along the pial–ventricular axis in V1. By applying multitasking theory, we suggest that the spatial zonation of L2/3 cell types is linked to the continuous nature of their gene expression profiles, which can be represented as a 2D manifold bounded by three archetypal cell types. By comparing normally reared and dark reared L2/3 cells, we show that visual deprivation-induced transcriptomic changes comprise two independent gene programs. The first, induced specifically in the visual cortex, includes immediate-early genes and genes associated with metabolic processes. It manifests as a change in cell state that is orthogonal to cell-type-specific gene expression programs. By contrast, the second program impacts L2/3 cell-type identity, regulating a subset of cell-type-specific genes and shifting the distribution of cells within the L2/3 cell-type manifold. Through an integrated analysis of spatial transcriptomics with single-nucleus RNA-seq data, we describe how vision patterns cortical L2/3 cell types during the critical period.

Science & Technology - Other Topics↗

The evaluation of the OSGLR algorithm for restructurable controls

The detection and isolation of commercial aircraft control surface and actuator failures using the orthogonal series generalized likelihood ratio (OSGLR) test was evaluated. The OSGLR algorithm was chosen as the most promising algorithm based on a preliminary evaluation of three failure detection and isolation (FDI) algorithms (the detection filter, the generalized likelihood ratio test, and the OSGLR test) and a survey of the literature. One difficulty of analytic FDI techniques and the OSGLR algorithm in particular is their sensitivity to modeling errors. Therefore, methods of improving the robustness of the algorithm were examined with the incorporation of age-weighting into the algorithm being the most effective approach, significantly reducing the sensitivity of the algorithm to modeling errors. The steady-state implementation of the algorithm based on a single cruise linear model was evaluated using a nonlinear simulation of a C-130 aircraft. A number of off-nominal no-failure flight conditions including maneuvers, nonzero flap deflections, different turbulence levels and steady winds were tested. Based on the no-failure decision functions produced by off-nominal flight conditions, the failure detection performance at the nominal flight condition was determined. The extension of the algorithm to a wider flight envelope by scheduling the linear models used by the algorithm on dynamic pressure and flap deflection was also considered. Since simply scheduling the linear models over the entire flight envelope is unlikely to be adequate, scheduling of the steady-state implentation of the algorithm was briefly investigated.

Bonnice, W. F.↗

On-Orbit Calibration of Redundant Spacecraft Gyros by Optimal Reduction to Three Axes

The Aqua spacecraft will carry four single-axis gyros configured with three orthogonal axes and one skew axis. This redundancy presents a challenge for batch methods of on-orbit gyro calibration that use a spacecraft rotation model deterministically related to gyro data, in that sensor data can respond to at most three angular velocity components. When the number of gyros, N, is greater than 3, the 3xN matrix, G, that reduces the N gyro measurements to three body-frame angular-velocity components cannot be fully determined by such methods; there are many such matrices that produce essentially the same angular velocity history. In such a case, spacecraft operators require information about the Nx3 gyro linear response matrix, R, that relates gyro outputs to the body-frame angular velocities causing them. This matrix provides sufficient information to determine multiple reduced-dimension G-matrices for use in case of failure or degradation of one or more gyros, as well as to determine an optimal 3xN G for the fully-functional configuration. A viable proposal is to apply a 3xN pre-filter matrix, F, to the N gyro outputs before carrying out a conventional gyro calibration procedure. The angular-velocity history emerging from conventional calibration may then be used as input data, together with the same gyro data that generated it, to fit the alignment, scale-factor, and bias parameters of each gyro axis in turn. A difficulty of such a proposal is the arbitrariness in the choice of F. Due to gyro noise, different pre-filter matrices produce different calibrations. This paper presents a method of choosing F that is based on optimizing gyro consistency in the limit of infinite weight on gyro data, as compared to sensor data. The choice of F is independent of a priori alignment and is based on the gyro data alone. The method is applicable to any N of three or more, but reduces to conventional batch-estimation methodologies when N = 3. Results of computational comparison among calibration simulations using various choices of F will be presented for the Aqua gyro configuration with N = 4.

Radomski, M. S.↗

Supporting technology for the development of Controlled Ecological Life Support Systems (CELSS)

To support the development of Controlled Ecological Life Support Systems (CELSS) in the space program, a metabolic simulator has been selected for use in a closed chamber to test functions of the CELSS. This metabolic simulator is a catalytic reactor which oxidizes the methyl acetate to produce carbon dioxide and water vapor. In this project, kinetic studies of catalytic oxidation of methyl acetate were conducted using monolithic and pellet catalysts with 0.5% (by weight) platinum (Pt) on aluminum oxide (Al2O3). The reaction was studied at a pressure of one atmosphere and at temperatures varying from 160 C to 420 C. By-products were identified at the exit of the preheater and reactor. For the kinetic study with the monolithic catalyst, a linear regression method was used to correlate the kinetic data with zero-order, first-order and Langmuir-Hinshelwood models. Results indicate that the first-order model represents the data adequately at low concentrations of methyl acetate. For higher concentrations of methyl acetate, the Langmuir-Hinshelwood model best represents the kinetic data. Both rate constant and adsorption equilibrium constants were estimated from the regression. A Taguchi orthogonal array (L(sub 9)) was used to investigate the effects of temperature, flow rate, and concentration on the catalytic oxidation of methyl acetate. For the monolithic catalyst, temperature exerts the most significant effect, followed by concentration of methyl acetate. For the pellet catalyst, reaction temperature is the most significant factor, followed by gas flow rate and methyl acetate concentration. Concentrations of either carbon dioxide or oxygen were seen to have insignificant effect on the methyl acetate conversion process. Experimental results indicate that the preheater with glass beads can accomplish thermal cracking and catalytic reaction of methyl acetate to produce acetic acid, methanol, methyl formate, and 1-propanol. The concentration of all by-products was measured in ppmv (parts per million by volume). At higher temperatures, greater amounts of these products are produced, as expected. In all cases, methanol was the predominant concentration detected, followed by methyl formate. At temperatures lower than 320 C for the P-type monolithic catalyst, methanol, acetic acid, and acetone were detected, whereas, for the E-type monolithic catalyst, only methanol was detected at 160 C. Both P and E types of the monolithic catalyst were specified with the same substrates (ceramic), washcoat (Al2O3), and promoter (Pt). However, the manufacturing and treatment procedures were quite different. It was therefore concluded that the performance of the E-type monolithic catalyst is superior to that of the P-type for oxidation of methyl acetate. At higher reaction temperatures, e.g., above 420 C, all reactants and byproducts were completely oxidized using these two types of monolithic catalyst to produce carbon dioxide and water vapor. A complex heterogenous catalytic reaction mechanism was proposed to explain the formation of the byproducts (methanol, acetic acid, and methyl formate) as the methyl acetate traveled through the preheater packed with glass beads. The by-product, 1-propanol, may be formed only through a homogeneous reaction, since it is difficult to develop a reasonable sequence of heterogeneous reaction steps to explain its formation. The homogeneous thermal decomposition of methyl acetate to form free radicals was proposed to explain the formation of 1-propanol, and also methanol, in the preheater. A dual-site catalytic reaction mechanism was proposed for the oxidation of methyl acetate over Pt/Al2O3 monolithic catalyst. The dual-site mechanism describes the chemisorption of oxygen molecules as well as a physical adsorption of methyl acetate on the active sites. On the active sites, methyl acetate is oxidized rapidly to form carbon dioxide and water vapor. A rate equation derived from this mechanism gives the Langmuir-Hinshelwood rate formula which has been observed from the experimental data obtained in this project for high methyl acetate concentration (greater than 1000 ppmv) over a monolithic catalyst. If the oxygen concentration is very high and methyl acetate concentration is very low, the reaction rate equation is then reduced to a first-order with respect to methyl acetate concentration. The first-order model has also been observed from the experimental data obtained in this project for low methyl acetate concentration (less than 1000 ppmv).

Li, Ku-Yen↗

Forward variable selection enables fast and accurate dynamic system identification with Karhunen-Loève decomposed Gaussian processes

A promising approach for scalable Gaussian processes (GPs) is the Karhunen-Loève (KL) decomposition, in which the GP kernel is represented by a set of basis functions which are the eigenfunctions of the kernel operator. Such decomposed kernels have the potential to be very fast, and do not depend on the selection of a reduced set of inducing points. However KL decompositions lead to high dimensionality, and variable selection thus becomes paramount. This paper reports a new method of forward variable selection, enabled by the ordered nature of the basis functions in the KL expansion of the Bayesian Smoothing Spline ANOVA kernel (BSS-ANOVA), coupled with fast Gibbs sampling in a fully Bayesian approach. It quickly and effectively limits the number of terms, yielding a method with competitive accuracies, training and inference times for tabular datasets of low feature set dimensionality. Theoretical computational complexities are O ( N P 2 ) in training and O ( P ) per point in inference, where N is the number of instances and P the number of expansion terms. The inference speed and accuracy makes the method especially useful for dynamic systems identification, by modeling the dynamics in the tangent space as a static problem, then integrating the learned dynamics using a high-order scheme. The methods are demonstrated on two dynamic datasets: a ‘Susceptible, Infected, Recovered’ (SIR) toy problem, along with the experimental ‘Cascaded Tanks’ benchmark dataset. Comparisons on the static prediction of time derivatives are made with a random forest (RF), a residual neural network (ResNet), and the Orthogonal Additive Kernel (OAK) inducing points scalable GP, while for the timeseries prediction comparisons are made with LSTM and GRU recurrent neural networks (RNNs) along with the SINDy package.

Hayes, Kyle↗

A Novel Method for Characterizing Spacesuit Mobility through Metabolic Cost

Spacesuit mobility has historically been defined and characterized by a combination of range of motion and joint torque of the individual anatomical joints when performing isolated motions meant to drive that joint only in a given orthogonal plane. While this has been the standard approach for several decades, there are numerous shortcomings that suit designers and engineers would like to see rectified. First, the lack of a standardized method for collecting both range of motion and joint torque translates to many different test setups, procedures and methods of data analysis. Second, all of these previously used methods for data collection lack some degree of repeatability, even within the same test setup and the same conductor; in addition, attempts at higher fidelity data collection techniques require high overhead and cost with minimal improvement. Lastly, isolated motions in standard anatomical planes are not representative of real‐world tasks that a crewmember would be performing during an EVA, be it microgravity or surface exploration based. To address these shortcomings, options are being explored within the Space Suit and Crew Survival Systems Branch to ascertain the feasibility of an alternative approach to defining mobility - one that is more repeatable, lower overhead, and more tied to functional EVA tasks. This paper serves to document the first attempt at such an alternative option - one that looks at the metabolic energy‐cost of a spacesuit. In other words, can we objectively compare the mobility of a spacesuit by evaluating the metabolic cost of that suit to the wearer while performing a battery of functional EVA tasks?

McFarland, Shane↗

A One‐Pot Biocatalytic Cascade to Access Diverse l ‐Phenylalanine Derivatives from Aldehydes or Carboxylic Acids

Abstract Nonstandard amino acids (nsAAs) that are l ‐phenylalanine derivatives with aryl ring functionalization have long been harnessed in natural product synthesis, therapeutic peptide synthesis, and diverse applications of genetic code expansion. Yet, to date, these chiral molecules have often been the products of poorly enantioselective and environmentally harsh organic synthesis routes. Here, we reveal the broad specificity of multiple natural pyridoxal 5′‐phosphate (PLP)‐dependent enzymes, specifically an l ‐threonine transaldolase, a phenylserine dehydratase, and an aminotransferase, toward substrates that contain aryl side chains with diverse substitutions. We exploit this tolerance to construct a one‐pot biocatalytic cascade that achieves high‐yield synthesis of 18 diverse l ‐phenylalanine derivatives from aldehydes under mild aqueous reaction conditions. We demonstrate the addition of a carboxylic acid reductase module to this cascade to enable the biosynthesis of l ‐phenylalanine derivatives from carboxylic acids that may be less expensive or less reactive than the corresponding aldehydes. Finally, we investigate the scalability of the cascade by developing a lysate‐based route for preparative‐scale synthesis of 4‐formyl‐ l ‐phenylalanine, a nsAA with a bio‐orthogonal handle that is not readily market‐accessible. Overall, this work offers an efficient, versatile, and scalable route with the potential to lower manufacturing costs and democratize synthesis for many valuable nsAAs.

Anderson, Shelby R. [Department of Chemical and Bi↗

SnSe 1-x S x Alloys: Anisotropic Van der Waals Semiconductors with Tunable Bandgaps

Alloying is one of the main tools of bandgap engineering, allowing the tuning of crystal structure, lattice parameters, and electronic structure of 3D and 2D/layered semiconductors. Among the latter, it can play a key role in tailoring the properties of tin monochalcogenides, a class of van der Waals semiconductors of interest for optoelectronics, thermoelectrics, ferroelectrics, and valleytronics. Here, the study investigates the synthesis and properties of large flakes of the anion substitution alloys SnSe 1-x S x . Alloy flakes across a wide range of compositions are obtained systematically by repeated growth from the same mixed (SnS, SnSe) powder precursor. Combined experiment and theory show full miscibility for all compositions, along with tunable lattice constants, bandgaps, and vibrational modes. Atomic resolution imaging demonstrates the accumulation of S and Se in alternating layers in the SnSe 1-x S x unit cell, attributed to growth kinetics. Polarized Raman spectroscopy confirms anisotropic vibrational modes; the calculated and measured band structure shows systematic changes in the band edge energies and anisotropic electronic structure due to the anisotropic in-plane lattice of the monochalcogenides. Cathodoluminescence, finally, indicates that a unique configuration of two non-degenerate, direct valleys along orthogonal k-space directions persists all the way from SnS to SnSe, making SnSe 1-x S x alloys interesting for valleytronics.

36 MATERIALS SCIENCE↗

Transverse momentum-dependent heavy-quark fragmentation at next-to-leading order

The transverse momentum-dependent fragmentation functions (TMD FFs) of heavy (bottom and charm) quarks, which we recently introduced, are universal building blocks that enter predictions for a large number of observables involving final-state heavy quarks or hadrons. They enable the extension of fixed-order subtraction schemes to quasi-collinear limits, and are of particular interest in their own right as probes of the nonperturbative dynamics of hadronization. In this paper we calculate all TMD FFs involving heavy quarks and the associated TMD matrix element in heavy-quark effective theory (HQET) to next-to-leading order in the strong interaction. Our results confirm the renormalization properties, large-mass, and small-mass consistency relations predicted in our earlier work. We also derive and confirm a prediction for the large-z behavior of the heavy-quark TMD FF by extending, for the first time, the formalism of joint resummation to capture quark mass effects in heavy-quark fragmentation. Our final results in position space agree with those of a recent calculation by another group that used a highly orthogonal organization of singularities in the intermediate momentum-space steps, providing a strong independent cross check. As an immediate application, we present the complete quark mass dependence of the energy-energy correlator (EEC) in the back-to-back limit at $\mathcal{O}\left({\alpha}_s\right)$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Transportable Gravity Gradiometer Based on Atom Interferometry

A transportable atom interferometer-based gravity gradiometer has been developed at JPL to carry out measurements of Earth's gravity field at ever finer spatial resolutions, and to facilitate high-resolution monitoring of temporal variations in the gravity field from ground- and flight-based platforms. Existing satellite-based gravity missions such as CHAMP and GRACE measure the gravity field via precise monitoring of the motion of the satellites; i.e. the satellites themselves function as test masses. JPL's quantum gravity gradiometer employs a quantum phase measurement technique, similar to that employed in atomic clocks, made possible by recent advances in laser cooling and manipulation of atoms. This measurement technique is based on atomwave interferometry, and individual laser-cooled atoms are used as drag-free test masses. The quantum gravity gradiometer employs two identical atom interferometers as precision accelerometers to measure the difference in gravitational acceleration between two points (Figure 1). By using the same lasers for the manipulation of atoms in both interferometers, the accelerometers have a common reference frame and non-inertial accelerations are effectively rejected as common mode noise in the differential measurement of the gravity gradient. As a result, the dual atom interferometer-based gravity gradiometer allows gravity measurements on a moving platform, while achieving the same long-term stability of the best atomic clocks. In the laboratory-based prototype (Figure 2), the cesium atoms used in each atom interferometer are initially collected and cooled in two separate magneto-optic traps (MOTs). Each MOT, consisting of three orthogonal pairs of counter-propagating laser beams centered on a quadrupole magnetic field, collects up to 10(exp 9) atoms. These atoms are then launched vertically as in an atom fountain by switching off the magnetic field and introducing a slight frequency shift between pairs of lasers to create a moving rest frame for the trapped atoms. While still in this moving-frame molasses, the laser frequencies are further detuned from the atomic resonance (while maintaining this relative frequency shift) to cool the atom cloud's temperature to 2 K or below, corresponding to an rms velocity of less than 2 cm/s. After launch, the cold atoms undergo further state and velocity selection to prepare for atom interferometry. The atom interferometers are then realized using laser-induced stimulated Raman transitions to perform the necessary manipulations of each atom, and the resulting interferometer phase is measured using laser-induced fluorescence for state-normalized detection. More than 20 laser beams with independent controls of frequency, phase, and intensity are required for this measurement sequence. This instrument can facilitate the study of Earth's gravitational field from surface and air vehicles, as well as from space by allowing gravity mapping from a low-cost, single spacecraft mission. In addition, the operation of atom interferometer-based instruments in space offers greater sensitivity than is possible in terrestrial instruments due to the much longer interrogation times available in the microgravity environment. A space-based quantum gravity gradiometer has the potential to achieve sensitivities similar to the GRACE mission at long spatial wavelengths, and will also have resolution similar to GOCE for measurement at shorter length scales.

Yu, Nan↗

Target-Tracking Camera for a Metrology System

An analog electronic camera that is part of a metrology system measures the varying direction to a light-emitting diode that serves as a bright point target. In the original application for which the camera was developed, the metrological system is used to determine the varying relative positions of radiating elements of an airborne synthetic aperture-radar (SAR) antenna as the airplane flexes during flight; precise knowledge of the relative positions as a function of time is needed for processing SAR readings. It has been common metrology system practice to measure the varying direction to a bright target by use of an electronic camera of the charge-coupled-device or active-pixel-sensor type. A major disadvantage of this practice arises from the necessity of reading out and digitizing the outputs from a large number of pixels and processing the resulting digital values in a computer to determine the centroid of a target: Because of the time taken by the readout, digitization, and computation, the update rate is limited to tens of hertz. In contrast, the analog nature of the present camera makes it possible to achieve an update rate of hundreds of hertz, and no computer is needed to determine the centroid. The camera is based on a position-sensitive detector (PSD), which is a rectangular photodiode with output contacts at opposite ends. PSDs are usually used in triangulation for measuring small distances. PSDs are manufactured in both one- and two-dimensional versions. Because it is very difficult to calibrate two-dimensional PSDs accurately, the focal-plane sensors used in this camera are two orthogonally mounted one-dimensional PSDs.

Liebe, Carl↗

Microgravity-Driven Optic Nerve/Sheath Biomechanics Simulations

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a concern for long-duration space flight. Current thinking suggests that the ocular changes observed in VIIP syndrome are related to cephalad fluid shifts resulting in altered fluid pressures [1]. In particular, we hypothesize that increased intracranial pressure (ICP) drives connective tissue remodeling of the posterior eye and optic nerve sheath (ONS). We describe here finite element (FE) modeling designed to understand how altered pressures, particularly altered ICP, affect the tissues of the posterior eye and optic nerve sheath (ONS) in VIIP. METHODS: Additional description of the modeling methodology is provided in the companion IWS abstract by Feola et al. In brief, a geometric model of the posterior eye and optic nerve, including the ONS, was created and the effects of fluid pressures on tissue deformations were simulated. We considered three ICP scenarios: an elevated ICP assumed to occur in chronic microgravity, and ICP in the upright and supine positions on earth. Within each scenario we used Latin hypercube sampling (LHS) to consider a range of ICPs, ONH tissue mechanical properties, intraocular pressures (IOPs) and mean arterial pressures (MAPs). The outcome measures were biomechanical strains in the lamina cribrosa, optic nerve and retina; here we focus on peak values of these strains, since elevated strain alters cell phenotype and induce tissue remodeling. In 3D, the strain field can be decomposed into three orthogonal components, denoted as first, second and third principal strains. RESULTS AND CONCLUSIONS: For baseline material properties, increasing ICP from 0 to 20 mmHg significantly changed strains within the posterior eye and ONS (Fig. 1), indicating that elevated ICP affects ocular tissue biomechanics. Notably, strains in the lamina cribrosa and retina became less extreme as ICP increased; however, within the optic nerve, the occurrence of such extreme strains greatly increased as ICP was elevated (Fig. 2). In particular, c. 48 of simulations in the elevated ICP condition showed peak strains in the optic nerve that exceeded the strains expected on earth. Such extreme strains are likely important, since they represent a larger signal for mechano-responsive resident cells [2]. The models predicted little to no anterior motion of the prelaminar neural tissue (optic nerve swelling, or papilledema, secondary to axoplasmic stasis), typically seen with elevated ICP. Specialized FE models to capture axoplasmic stasis would be required to study papilledema. These results suggest that the most notable effect of elevated ICP may occur via direct optic nerve loading, rather than through connective tissue deformation. These FE models can inform the design of future studies designed to bridge the gap between biomechanics and pathophysiological function in VIIP.

biomechanics↗

Parameterizing Coefficients of a POD-Based Dynamical System

A method of parameterizing the coefficients of a dynamical system based of a proper orthogonal decomposition (POD) representing the flow dynamics of a viscous fluid has been introduced. (A brief description of POD is presented in the immediately preceding article.) The present parameterization method is intended to enable construction of the dynamical system to accurately represent the temporal evolution of the flow dynamics over a range of Reynolds numbers. The need for this or a similar method arises as follows: A procedure that includes direct numerical simulation followed by POD, followed by Galerkin projection to a dynamical system has been proven to enable representation of flow dynamics by a low-dimensional model at the Reynolds number of the simulation. However, a more difficult task is to obtain models that are valid over a range of Reynolds numbers. Extrapolation of low-dimensional models by use of straightforward Reynolds-number-based parameter continuation has proven to be inadequate for successful prediction of flows. A key part of the problem of constructing a dynamical system to accurately represent the temporal evolution of the flow dynamics over a range of Reynolds numbers is the problem of understanding and providing for the variation of the coefficients of the dynamical system with the Reynolds number. Prior methods do not enable capture of temporal dynamics over ranges of Reynolds numbers in low-dimensional models, and are not even satisfactory when large numbers of modes are used. The basic idea of the present method is to solve the problem through a suitable parameterization of the coefficients of the dynamical system. The parameterization computations involve utilization of the transfer of kinetic energy between modes as a function of Reynolds number. The thus-parameterized dynamical system accurately predicts the flow dynamics and is applicable to a range of flow problems in the dynamical regime around the Hopf bifurcation. Parameter-continuation software can be used on the parameterized dynamical system to derive a bifurcation diagram that accurately predicts the temporal flow behavior.

Kalb, Virginia L.↗

Evaluating Nighttime CALIOP 0.532 micron Aerosol Optical Depth and Extinction Coefficient Retrievals

NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) Version 3.01 5-km nighttime 0.532 micron aerosol optical depth (AOD) datasets from 2007 are screened, averaged and evaluated at 1 deg X 1 deg resolution versus corresponding/co-incident 0.550 micron AOD derived using the US Navy Aerosol Analysis and Prediction System (NAAPS), featuring two-dimensional variational assimilation of quality-assured NASA Moderate Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) AOD. In the absence of sunlight, since passive radiometric AOD retrievals rely overwhelmingly on scattered radiances, the model represents one of the few practical global estimates available from which to attempt such a validation. Daytime comparisons, though, provide useful context. Regional-mean CALIOP vertical profiles of night/day 0.532 micron extinction coefficient are compared with 0.523/0.532 micron ground-based lidar measurements to investigate representativeness and diurnal variability. In this analysis, mean nighttime CALIOP AOD are mostly lower than daytime (0.121 vs. 0.126 for all aggregated data points, and 0.099 vs. 0.102 when averaged globally per normalised 1 deg. X 1 deg. bin), though the relationship is reversed over land and coastal regions when the data are averaged per normalised bin (0.134/0.108 vs. 0140/0.112, respectively). Offsets assessed within single bins alone approach +/- 20 %. CALIOP AOD, both day and night, are higher than NAAPS over land (0.137 vs. 0.124) and equal over water (0.082 vs. 0.083) when averaged globally per normalised bin. However, for all data points inclusive, NAAPS exceeds CALIOP over land, coast and ocean, both day and night. Again, differences assessed within single bins approach 50% in extreme cases. Correlation between CALIOP and NAAPS AOD is comparable during both day and night. Higher correlation is found nearest the equator, both as a function of sample size and relative signal magnitudes inherent at these latitudes. Root mean square deviation between CALIOP and NAAPS varies between 0.1 and 0.3 globally during both day/night. Averaging of CALIOP along-track AOD data points within a single NAAPS grid bin improves correlation and RMSD, though day/night and land/ocean biases persist and are believed systematic. Vertical profiles of extinction coefficient derived in the Caribbean compare well with ground-based lidar observations, though potentially anomalous selection of a priori lidar ratios for CALIOP retrievals is likely inducing some discrepancies. Mean effective aerosol layer top heights are stable between day and night, indicating consistent layer-identification diurnally, which is noteworthy considering the potential limiting effects of ambient solar noise during day.

aerosol optical depth↗