Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear mappings”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Internal Wave Study in the South China Sea Using SAR

Recently, the internal wave distribution maps in the China Seas have been compiled from hundreds of ERS-1/2, RADARSAT, and Space Shuttle SAR (Synthetic Aperture Radar) images from 1993 to 1999. Based on internal wave distribution map, most of internal waves in the northeast part of South China Sea were propagating westward. The wave crest can be as long as 200 km with amplitude of 100 m due to strong current from the Kuroshio branching out into the South China Sea. Based on the observations from drilling rigs near DongSha Island by Amoco Production Co., the solitons may be generated in a 4 km wide channel between Batan and Sabtang islands in Luzon Strait. The proposed generation mechanism is similar to the lee wave formation from a shallow topography. Both depression and elevation internal waves have been observed in the same RADARSAT ScanSAR image on May 4, 1998 near DongSha Island. Furthermore, depression and elevation internal waves have also been observed by SAR at the same location on the shelf in April and June, 1993 (in different seasons) respectively. Numerical models have been used to interpret their generation mechanism and evolution processes. Based on the SAR images, near DongSha Island, the westward propagating huge internal solitons are often encountered and diffracted/broken by the coral reefs on the shelf. After passing the island, the diffracted waves will re-merge or interact with each other. It has been observed that after the nonlinear wave-wave interaction, the phase of wave packet is shifted and wavelength is also changed. Examples of mesoscale features observed in SAR images, such as fronts, raincells, bathymetry, ship wakes, and oil spills will be presented. Recent mooring measurements in April 1999 near Dongsha Island, future field test ASIAEX (Asian Seas International Acoustics Experiment) planned for April 2001, and some pretest survey data will be discussed in this paper.

Liu, Antony K.↗
DOE OSTI2020

Stable Rate of Slip Along the Karakax Section of the Altyn Tagh Fault from Observation of Interglacial and Postglacial Offset Morphology and Surface Dating

Digital elevation maps obtained using TanDEM-X and Pleiades data combined with newly obtained surface age estimates using cosmogenic radionuclide (CRN) and optically simulated luminescence (OSL) methods are used to quantify the slip rate along the western section of the Altyn Tagh Fault in southern Xinjiang. The reconstruction of the conical shape of massive alluvial fans inferred to be from the Eemian (115±7 ka) from CRN dating shows consistent left-lateral offsets of 300±20 m, yielding a slip rate of 2.6±0.3 mm/year. Successive episodes of incision have left cut terraces inset in wide canyons, 10–25 m below the fans' surface. The incision was followed by the deposition of a broad terrace of early Holocene age, which is reincised by modern stream channels. Near the village of Shanxili, a 200 m wide valley is partially dammed by a shutter ridge displaced by the fault. A fill terrace deposited upstream from the ridge has an OSL age of 8.8±0.6 ka. The 23±2 m offset of the riser incising the terrace indicates a minimum postdepositional movement on the fault, yielding a Holocene rate of 2.6±0.5 mm/year, consistent with the 115-ka average slip rate. Furthermore, scarp degradation analysis using mass diffusion reveals a nonlinear relationship between fault displacement and degradation coefficient along the progressively exposed fault scarp, a pattern suggesting either seismic clustering or variable diffusion rate since the Eemian. Together with the Gozha Co-Longmu Co fault to the south, the Karakax section of the Altyn Tagh Fault contributes to the eastward movement of the western corner of Tibet.

58 GEOSCIENCES↗
DOE OSTI2026

Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems

In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.

Bayesian approximation error↗
NASA NTRS1987

Deconvolution of spectra for intimate mixtures

Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.

Mustard, John F.↗
DOE OSTI2025

Local charge-screening and polarization switching in a ferroelectric ionic-superconductor

For ferroelectric ionic-conductors, polarization switching is complicated by the interplay between ion mobility and charge screening effects. When the ionic charge carriers also play a key role in the domain reversal, such as in Rb-doped KTiOPO 4 (RKTP), a higher level of complexity is introduced. RKTP provides an ideal platform for investigating the relationship between ionic conductivity and polarization reversal because its highly anisotropic crystal properties allow selective modification of material characteristics through diffusive cation doping. Here, we use indiffused Ba/K doping to create a significant increase in the ionic conductivity. Time-of-flight secondary ion mass spectrometry is employed to map Ba/K doping within the RKTP crystal and correlate it to changes in ionic mobility and polarization switching characteristics under an external field applied to the nonpolar face. Using band-excitation piezoresponse force microscopy, we demonstrate a selective switching-inhibition mechanism driven by the enhanced charge screening.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Determining Key Factors for the Open-Loop Control of Molecular Fragmentation Using Shaped Strong Fields

Pulse shaping has long been employed for tailoring femtosecond laser pulses to study and control the fragmentation of polyatomic molecules. In many cases, a physical explanation connecting the properties of the field to the observed control is difficult to ascertain. We utilized 80 bit binary spectral phase functions to parametrize and map the search space, gaining insight into which pulse parameters most impact the ion yield and fragmentation pattern for the relatively large triethylamine [N(C 2 H 5 ) 3 ] molecule. Pulse structures used to control the m/z 86 branching ratio beyond a simple intensity dependence are identified and compared to pump–probe results. All of these findings are explained in terms of control via a dissociative Rydberg state in the neutral molecule. This methodology may be used to discover new control mechanisms and shed light onto which pulse parameters most influence the interaction between strong field lasers and matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Operator inference for non-intrusive model reduction with quadratic manifolds

This paper proposes a novel approach for learning a data-driven quadratic manifold from high-dimensional data, then employing this quadratic manifold to derive efficient physics-based reduced-order models. The key ingredient of the approach is a polynomial mapping between high-dimensional states and a low-dimensional embedding. This mapping consists of two parts: a representation in a linear subspace (computed in this work using the proper orthogonal decomposition) and a quadratic component. The approach can be viewed as a form of data-driven closure modeling, since the quadratic component introduces directions into the approximation that lie in the orthogonal complement of the linear subspace, but without introducing any additional degrees of freedom to the low-dimensional representation. Combining the quadratic manifold approximation with the operator inference method for projection-based model reduction leads to a scalable non-intrusive approach for learning reduced-order models of dynamical systems. Further, applying the new approach to transport-dominated systems of partial differential equations illustrates the gains in efficiency that can be achieved over approximation in a linear subspace.

42 ENGINEERING↗
DOE OSTI2021

Understanding Self-Assembly and the Stabilization of Liquid/Liquid Interfaces: The Importance of Ligand Tail Branching and Oil-Phase Solvation

Organophosphorus-based ligands represent a versatile set of solvent extraction reagents whose chemical makeup plays an important role in extraction mechanism. We hypothesize that the branching of the extractant hydrophobic tail and its oil-phase solvation affect the liquid/liquid interfacial structure. Understanding the structure mediated adsorption and interfacial ordering becomes key in designing ligands with enhanced selectivity and efficiency for targeted extractions. We employed vibrational sum frequency generation spectroscopy and interfacial tension measurements to extract thermodynamic adsorption energies, map interfacial ordering, and rationalize disparate behaviors of model di-(2-ethylhexyl) phosphoric acid and dioctyl phosphoric acid ligands at the hexadecane water interface. With increased surface loading, ligands with branched hydrophobic tails formed stable interfaces at much lower concentrations than those observed for ligands with linear alkyl tails. The lack of an oil phase and associated solvation results in markedly different interfacial properties, and thus measurements made at air/liquid surfaces cannot be assumed to correlate with the processes occurring at buried liquid/liquid interfaces. Here, we attribute these differences in the surface mediated self-assembly to key variations in hydrophobic interactions and tail solvation taking place in the oil phase demonstrating that interactions in both the polar and nonpolar phases are essential to understand self-assembly and function.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Domain wall motion across microstructural features in polycrystalline ferroelectric films

This paper describes the effect of microstructural features such as grain boundaries and triple points on the pinning of domain wall motion in perovskite Pb(Zr 0.52 Ti 0.48 )O 3 (PZT) films on polycrystalline SrRuO 3 /SrTiO 3 substrates. Here, spatial variability in the collective domain wall dynamics was assessed using non-linearity mapping via Band Excitation Piezoresponse Force Microscopy (BE-PFM). Collocating the non-linearity maps with triple point locations (as visualized by EBSD) allowed for exploration of the effects that local microstructure (e.g., grain boundary) have on domain wall motion. It was found that the extrinsic behavior varied with both the misorientation angle and the proximity to the grain boundary. The width of influence of individual grain boundaries on the motion of domain walls was a function of the character of the grain boundary; random grain boundaries exhibit deeper minima in α d /d 33,initial and larger widths of influence (up to 905 nm) compared to coincident site lattice (CSL) boundaries (up to 572 nm). Additionally, triple points containing larger numbers of random boundaries exhibited non-Rayleigh behavior to greater distances, suggesting that the triple point provides either a deep potential minimum or a region where domain wall motion is unfavorable.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Applications of Data Assimilation to Analysis of the Ocean on Large Scales

It is commonplace to begin talks on this topic by noting that oceanographic data are too scarce and sparse to provide complete initial and boundary conditions for large-scale ocean models. Even considering the availability of remotely-sensed data such as radar altimetry from the TOPEX and ERS-1 satellites, a glance at a map of available subsurface data should convince most observers that this is still the case. Data are still too sparse for comprehensive treatment of interannual to interdecadal climate change through the use of models, since the new data sets have not been around for very long. In view of the dearth of data, we must note that the overall picture is changing rapidly. Recently, there have been a number of large scale ocean analysis and prediction efforts, some of which now run on an operational or at least quasi-operational basis, most notably the model based analyses of the tropical oceans. These programs are modeled on numerical weather prediction. Aside from the success of the global tide models, assimilation of data in the tropics, in support of prediction and analysis of seasonal to interannual climate change, is probably the area of large scale ocean modeling and data assimilation in which the most progress has been made. Climate change is a problem which is particularly suited to advanced data assimilation methods. Linear models are useful, and the linear theory can be exploited. For the most part, the data are sufficiently sparse that implementation of advanced methods is worthwhile. As an example of a large scale data assimilation experiment with a recent extensive data set, we present results of a tropical ocean experiment in which the Kalman filter was used to assimilate three years of altimetric data from Geosat into a coarsely resolved linearized long wave shallow water model. Since nonlinear processes dominate the local dynamic signal outside the tropics, subsurface dynamical quantities cannot be reliably inferred from surface height anomalies. Because of its potential for large scale synoptic coverage of the deep ocean, acoustic travel time data should be a natural complement to satellite altimetry. Satellite data give us vertical integrals associated with thermodynamic and dynamic processes.

Miller, Robert N.↗
DOE OSTI2026

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗
DOE OSTI2023

PCAfold 2.0—Novel tools and algorithms for low-dimensional manifold assessment and optimization

We describe an update to our open-source Python package, PCAfold, designed to help researchers generate, analyze and improve low-dimensional data manifolds. In the current version, PCAfold 2.0, we introduce novel tools and algorithms for assessing and optimizing low-dimensional manifolds. This includes a method that generates a “map” of local feature sizes that can help pinpoint researchers to problematic regions on a manifold. We introduce a novel cost function that characterizes the quality of a manifold topology with a single number. We develop two algorithms for feature selection based on principal component analysis (PCA) that use the cost function as an objective function to minimize. We introduce a quantity of interest (QoI)-aware dimensionality reduction strategy where data projections are computed using an artificial neural network and are directly optimized towards representing various projection-independent and projection-dependent QoIs. We also introduce an implementation of partition of unity networks (POUnets) for efficient reconstruction of QoIs from low-dimensional manifolds based on combining neural network classification with localized polynomial regression. Our software can be broadly applicable in all domains of science and engineering that aim to reduce data dimensionality, as well as in the fundamental research on representation learning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Functional Adaptation of LPS-affected Dentoalveolar Fibrous Joints in Rats

We report the functional interplay between cementum of the root and alveolar bone of the socket is tuned by a uniquely positioned 70–80 µm wide fibrous and lubricious ligament in a dentoalveolar joint (DAJ). In this study, structural and biomechanical properties of the DAJ, periodontal ligament space (PDL-space also known as the joint space), alveolar bone of the socket, and cementum of the tooth root that govern the biomechanics of a lipopolysaccharide (LPS)-affected DAJ were mapped both in space and time. The hemi-maxillae from 20 rats (4 control at 6 weeks of age, 4 control and 4 LPS-affected at 12 weeks of age, 4 control and 4 LPS-affected at 16 weeks of age) were investigated using a hybrid technique; micro-X-ray computed tomography (5 µm resolution) in combination with biomechanical testing in situ. Temporal variations in bone and cementum volume fractions were evaluated. Trends in mineral apposition rates (MAR) in additional six Sprague Dawley rats (3 controls, 3 LPS-affected) were revealed by transforming spatial fluorochrome signals to functional growth rates (linearity factor - RW) of bone, dentin, and cementum using a fast Fourier transform on fluorochrome signals from 100-µm hemi-maxillae sections. An overall change in LPS-affected DAJ biomechanics (a 2.5–4.5X increase in tooth displacement and 2X tooth rotation at 6 weeks, no increase in displacement and a 7X increase in rotation at 12 weeks; 27% increase in bone effective strain at 6 weeks and 11% at 12 weeks relative to control) was associated with structural changes in the coronal regions of the DAJ (15% increase in PDL-space from 0 to 6 weeks but only 5% from 6 to 12 weeks compared to control). A significant increase (p < 0.05) in PDL-space between ligated and age-matched control was observed. The bone fraction of ligated at 12 weeks was significantly lower than its age-matched control, and no significant differences (p > 0.05) between groups were observed at 6 weeks. Cementum in the apical regions grew faster but nonlinearly (11% and 20% increase in cementum fraction (CF) at 6 and 12 weeks) compared to control. Alveolar bone revealed site-specific nonlinear growth with an overall increase in MAR (108.5 µm/week to 126.7 µm/week after LPS treatment) compared to dentin (28.3 µm/week in control vs. 26.1 µm/week in LPS-affected) and cementum (126.5 µm/week in control vs. 119.9 µm/week in LPS-affected). A significant increase in CF (p < 0.05) in ligated specimens was observed at 6 weeks of age. Anatomy-specific responses of cementum and bone to the mechano-chemo stimuli, and their collective temporal contribution to observed changes in PDL-space were perpetuated by altered tooth movement. Data highlight the “resilience” of DAJ function through the predominance of nonlinear growth response of cementum, changes in PDL-space, and bone architecture. Despite the significant differences in bone and cementum architectures, data provided insights into the reactionary effects of cementum as a built-in compensatory mechanism to reestablish functional competence of the DAJ. The spatial shifts in architectures of alveolar bone and cementum, and consequently ligament space, highlight adaptations farther away from the site of insult, which also is another novel insight from this study. These adaptations when correlated within the context of joint function (biomechanics) illustrate that they are indeed necessary to sustain DAJ function albeit being pathological.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

GenMod: A generative modeling approach for spectral representation of PDEs with random inputs

Here, we propose a method for quantifying uncertainty in high-dimensional PDE systems with random parameters, where the number of solution evaluations is small. Parametric PDE solutions are often approximated using a spectral decomposition based on polynomial chaos expansions. For the class of systems we consider (i.e., high dimensional with limited solution evaluations) the coefficients are given by an underdetermined linear system in a regression formulation. This implies additional assumptions, such as sparsity of the coefficient vector, are needed to approximate the solution. Here, we present an approach where we assume the coefficients are close to the range of a generative model that maps from a low to a high dimensional space of coefficients. Our approach is inspired be recent work examining how generative models can be used for compressed sensing in systems with random Gaussian measurement matrices. Using results from PDE theory on coefficient decay rates, we construct an explicit generative model that predicts the polynomial chaos coefficient magnitudes. The algorithm we developed to find the coefficients, which we call GenMod, is composed of two main steps. First, we predict the coefficient signs using Orthogonal Matching Pursuit. Then, we assume the coefficients are within a sparse deviation from the range of a sign-adjusted generative model. This allows us to find the coefficients by solving a nonconvex optimization problem, over the input space of the generative model and the space of sparse vectors. We obtain theoretical recovery results for a Lipschitz continuous generative model and for a more specific generative model, based on coefficient decay rate bounds. We examine three high-dimensional problems and show that, for all three examples, the generative model approach outperforms sparsity promoting methods at small sample sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Land surface dynamics and meteorological forcings modulate land surface temperature characteristics

This study examines the effect of land cover, vegetation health, climatic forcings, elevation heat loads, and terrain characteristics (LVCET) on land surface temperature (LST) distribution over West Africa (WA). We employ fourteen machine-learning models, which preserve nonlinear relationships, to downscale LST and other predictands while preserving the geographical variability of WA. Our results showed that the random forest model performs best in downscaling predictands. This is important for the sub-region since it has limited access to mainframes to power multiplex machine-learning algorithms. In contrast to the northern regions, the southern regions consistently exhibit healthy vegetation. Also, areas with unhealthy vegetation coincide with hot LST clusters. The positive Normalized Difference Vegetation Index (NDVI) trends in the Sahel underscore rainfall recovery and subsequent Sahelian greening. The southwesterly winds cause the upwelling of cold waters, lowering LST in southern WA and highlighting the cooling influence of water bodies on LST. Identifying regions with elevated LST is paramount for prioritizing greening initiatives, and our study underscores the importance of considering LVCET factors in urban planning. Topographic slope-facing angles, heat loads, and diurnal anisotropic heat all contribute to variations in LST, emphasizing the need for a holistic approach when designing resilient and sustainable landscapes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2021

LOFTID Heat Flux Gauge Calibration: What is Truth?

The Low-Earth Orbit Flight Test of an Inflatable Decelerator (LOFTID) is a demonstration of Hypersonic Inflatable Aerodynamic Decelerator (HIAD) technology, which may enable the delivery of heavy payloads to Mars, Venus, and Titan, as well as return to Earth. Unlike rigid aeroshells that are constrained by the size of the rocket’s shroud, inflatable aeroshells can be deployed to a much larger scale, thus allowing a spacecraft to begin its deceleration earlier and experience less heating. On LOFTID, there will be 4 total heat flux gauges (HFG) with a range of 70 W/cm2 and 1 radiometer with a range of 3 W/cm2, arranged as shown in Fig. 1. Both the radiometer and total HFGs are Schmidt-Boelter gauges purchased from an external vendor. Radiative calibrations were performed in-house at NASA Ames’ Sensors and TPS Advanced Research Laboratories (STAR Labs) before and after environmental testing to investigate how the testing affected the sensors' response. Additional rounds of radiative calibration at STAR Labs were also performed in order to investigate the large uncertainties associated with these tests. For example, a survey of multiple calibration facilities concluded that the uncertainty within a given facility was +/-3% [1]. An additional NIST study that calibrated heat flux gauges at 7 different facilities also found the variation in calibration coefficients to be up to ~3% within a given facility, but up to 15% between facilities, suggesting systematic differences between test setups [2]. Finally, the response of heat flux gauges to radiative versus convective heat flux has shown to differ by up to 20% [3],[4]. Because the heat flux gauges on LOFTID will predominantly experience convective heat flux during flight, a convective calibration study was performed at Boeing's Large-Core Arc Tunnel (LCAT) facility. Radiative Calibration Procedure The calibrations performed at STAR Labs utilize a quartz lamp bank (QLB) that provides a maximum heat flux of 50 W/cm2, which bounds the expected LOFTID flight environment. The calibration involves exposing a water-cooled Gardon gauge (reference) and then the unit-under-test (UUT) to 5 different heat fluxes multiple times for 10 seconds each, and then calculating a linear fit. The test setup is shown in Fig. 2. The total HFGs were calibrated at STAR Labs 3 times, denoted as STAR 1 (before environmental testing), STAR 2 (after protoflight vibration and thermal-vacuum testing), and STAR 3 (no change from previous test). All 8 flight-lot total HFGs showed a decrease in full-scale output from STAR 1 to STAR 2 by between 0.5% and 10. The first portion of this investigation was to determine whether the change could be due to differences in temperature between the two calibration runs. A typical linear fit to the calibration data was performed using Eq. 1 where q’ is the heat flux in W/cm2, c is the calibration coefficient, and mV is the sensor output. To account for temperature, the data were fit to a nonlinear function that included both the sensor output (mV) and the temperature from the thermocouple embedded inside the HFG near the surface (T): q'=mV/(c1* T + c0}. The residuals between the fits and the actual data points were calculated for every point, and proven to be much smaller for the temperature-compensated fits than for the linear fits for all sensors. An example is shown in Fig. 3. When the temperature-compensated fits from STAR 1 were applied to the STAR 2 data, the residuals did not improve, suggesting that the change in sensitivity between these two calibration runs was not due to temperature. A third round of calibration (STAR 3) was conducted to further address the temperature dependence of the total HFGs, and the resulting sensitivities matched closely to STAR 2 (within 2%). Temperature-compensated calibration curves were once again fit to the data. In this case, when the temperature-compensated fits from STAR 3 were applied to STAR 2 data, the residuals between the fits and STAR 2 data were much lower than the residuals due to the linear fits. This suggests that the changes seen between STAR 1 and STAR 2 were likely due to actual changes in the sensors caused by the environmental testing between the two calibrations. A modification of the original calibration process, in which the UUT was exposed to each heat flux for just 3 seconds (instead of 10) to reduce the temperature increase during the test, was additionally performed on several of the HFGs. In general, the sensitivities were 1-1.5% lower than from the 10-second tests, but the temperatures were also significantly lower. When the temperature-compensated fits from the 10-second tests were applied to the 3-second test data, the residuals were greatly improved than when just using the linear fits, further suggesting that the temperature-compensated fits may lead to better accuracy than the linear fits in flight. Convective Calibration The second portion of this study was to create a mapping between the radiative and convective calibration coefficients. The majority of the heating during flight will be convective, so it is important to understand how the HFG response differs under these conditions. However, there are no standardized methods for convective calibration [5]. Because the TPS aerothermal response models were validated at LCAT, the same facility was chosen for convective calibration of two of the total HFGs (Fig. 4). Preliminary results showed that the full-scale output was 3% and 8% higher in convective heat flux as compared to radiative heat flux. However, tunnel variation may have contributed to noise and uncertainty in the measurements, and more testing and analysis remains to be done. Scope of Presentation The presentation will include an overview of the changes seen in HFG calibration before and after environmental testing, differences between radiative and convective calibrations, the modeling work done to aid in understanding the sensor response to varying environments, and recommended future work.

H S Alpert↗
NASA NTRS2011

Project Management Using Modern Guidance, Navigation and Control Theory

Implementing guidance, navigation, and control (GN&C) theory principles and applying them to the human element of project management and control is not a new concept. As both the literature on the subject and the real-world applications are neither readily available nor comprehensive with regard to how such principles might be applied, this paper has been written to educate the project manager on the "laws of physics" of his or her project (not to teach a GN&C engineer how to become a project manager) and to provide an intuitive, mathematical explanation as to the control and behavior of projects. This paper will also address how the fundamental principles of modern GN&C were applied to the National Aeronautics and Space Administration's (NASA) Constellation Program (CxP) space suit project, ensuring the project was managed within cost, schedule, and budget. A project that is akin to a physical system can be modeled and managed using the same over arching principles of GN&C that would be used if that project were a complex vehicle, a complex system(s), or complex software with time-varying processes (at times nonlinear) containing multiple data inputs of varying accuracy and a range of operating points. The classic GN&C theory approach could thus be applied to small, well-defined projects; yet when working with larger, multiyear projects necessitating multiple organizational structures, numerous external influences, and a multitude of diverse resources, modern GN&C principles are required to model and manage the project. The fundamental principles of a GN&C system incorporate these basic concepts: State, Behavior, Feedback Control, Navigation, Guidance and Planning Logic systems. The State of a system defines the aspects of the system that can change over time; e.g., position, velocity, acceleration, coordinate-based attitude, and temperature, etc. The Behavior of the system focuses more on what changes are possible within the system; this is denoted in the state of the system. The behavior of a system, as captured in the system modeling, when properly done will aid in accurately predicting future system performance. The Feedback Control system understands the state and behavior of the system and uses feedback to adjust control inputs into the system. The feedback, which is the right arm of the Control system, allows change to be affected in the overall system; it therefore is important to not only correctly identify the system feedback inputs, but also the system response to the feedback inputs. The Navigation system takes multiple data inputs and based on a priori knowledge of the inputs, develops a statistically based weighting of the inputs and measurements to determine the system's state. Guidance and Planning Logic of the system, complete with an understanding of where the system is (provided by the Navigation system), will in turn determine where the system needs to be and how to get it there. With any system/project, it is critical that the objective of the system/project be clearly defined -- not only to plan but to measure performance and to aid in guiding the system or the project. The system principles discussed above, which can be and have been applied to the current CxP space suit development project, can also be mapped to real-world constituents, thus allowing project managers to apply systems theories that are well defined in engineering and mathematics to a discipline (i.e., Project Management) that historically has been based in personal experience and intuition. This mapping of GN&C theory to Project Management will, in turn, permit a direct, methodical approach to Project Management, planning and control providing a tool to help predict (and guide) performance and an understanding of the project constraints, how the project can be controlled, and the impacts to external influences and inputs. This approach, to a project manager, flows down to the three bottom-line variables of cost, schedule, and scope ando the needed control of these three variables to successfully perform and complete a project.

Hill, Terry R.↗