Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete linearization method”

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 343 records · Page 19

Demonstration of Probabilistic Sensitivity Analyses Tools on the Structural Response of a Representative Inflatable Space Structure

This work provides an initial step toward demonstrating a probabilistic numerical simulation capability to support trade studies and the development of a certification plan for inflatable space habitats. This study concentrates on interpreting the results from probabilistic analysis and numerical simulation tools to identify parameter sensitivities for a novel inflatable airlock concept, specifically the Non‐Axisymmetric Inflatable Pressure Structure (NAIPS) that was designed and tested under NASA's Minimalistic Advanced Softgoods Hatch (MASH) Program. A brief overview of the finite element model is provided along with the probabilistic sensitivity analysis approach. The sensitivity studies required a model that was numerically stable and efficient enough that hundreds of simulations could be completed in the allotted time. Therefore, the existing full model was simplified by: extracting a quarter symmetry section of the dome; focusing on a single inflation pressure; and replacing the non‐linear material stress‐strain curves with linear, isotropic materials defined by elastic moduli. Responses of interest include the sensitivity of various structural component loads to material properties, cord lengths, inflation pressure and friction. Multiple sensitivity studies were completed and three are reported here. The first study focused on utilizing wide input parameter ranges to provide an opportunity to assess numerical robustness. The next two studies narrowed the parameter ranges to enable focus on understanding uncertainty at a fixed operating condition. The completion of the sensitivity studies improved understanding of the interdependence of multiple inputs on the responses. In addition, numerical stability of the simulations over wide parameter ranges, shows the feasibility of incorporating uncertainty‐based methods in the design and certification of inflatable space habitats. With the experience and trust gained, it is anticipated that these same methods will be applied to nonlinear, orthotropic models in the future.

Lyle, Karen H.↗

Calibration of a helium-cooled infrared spatial radiometer and grating spectrometer

Methods used by the Space Dynamics Laboratory of Utah State University (SDL/USU) to calibrate infrared sensors are described, using the Infrared Background Signature Survey (IBSS) spatial radiometer and grating spectrometer as examples. A calibration equation and a radiometric model are given for each sensor to describe their responsivity in terms of individual radiometric parameters. The calibration equation terms include dark offset, linearity, absolute responsivity, and measurement uncertainty, and the radiometric model domains include spatial, spectral, and temporal domains. A portable calibration facility, designed and fabricated by SDL/USU, provided collimated, extended, diffuse scatter, and Jones sources in a single cryogenic dewar. This multi-function calibrator allowed calibration personnel to complete a full calibration of the IBSS infrared radiometer and spectrometer in two 15-day periods. A calibration data system was developed to control and monitor the calibration facility, and to record and analyze sensor data.

Jacobsen, Larry↗

The application of 'strain range partitioning method' to torsional creep-fatigue interaction

The method of strain range partitioning was applied to a series of torsional fatigue tests conducted on tubular 304 stainless steel specimens at 1200 F (649 C). Creep strain was superimposed on cycling strain, and the resulting strain range was partitioned into four components; completely reversed plastic shear strain followed by creep strain, creep strain followed by plastic strain and completely reversed creep strain. Each strain component was related to the cyclic life of the material. The paper describes the experimental procedure used to achieve strain partitioning and the torsional test results are compared to those obtained from axial tests. The damaging effects of the individual strain components were expressed by a linear life fraction rule. The shear strain plastic component showed the least detrimental factor when compared to creep strain reversal by plastic strain. In the latter case, a reduction of torsional fatigue life in the order of magnitude of 1.5 was observed.

Zamrik, S. Y.↗

High Accuracy In-Flight Wavelength Calibration of Imaging Spectrometry Data

Accurate wavelength calibration of imaging spectrometer data is essential if proper atmospheric transmission corrections are to be made to obtain apparent surface reflectance. Accuracies of 0.1 nm are necessary for a 10 nm-sampling instrument in order to match the slopes of the deep atmospheric water vapor features that dominate the 0.7-2.3 micrometer regions. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is calibrated in the laboratory to determine the wavelength position and full-width-half-maximum (FWHM) response for each of the 224 channels. The accuracies are limited by the quality of the monochromator used as a source. The accuracies vary from plus or minus to plus or minus 1.5 nm depending on the wavelength region, in general decreasing with increasing wavelength. Green et al. make corrections to the wavelength calibrations by using the known positions of 14 atmospheric absorption features throughout the 0.4-2.5 micrometer wavelength region. These features, having varying width and intensity, were matched to the MODTRAN model with a non-linear least squares fitting algorithm. A complete calibration was developed for all 224 channels by interpolation. Instrument calibration cannot be assumed to be stable to 0.1 nm over a flight season given the rigors of thermal cycling and launch and landing loads. The upcoming sensor HYDICE will require a means for in-flight spectral calibration of each scene because the calibration is both temperature and pressure sensitive. In addition, any sensor using a two-dimensional array has the potential for systematic wavelength shifts as a function of cross-track position, commonly called 'smile'. Therefore, a rapid means for calibrating complete images will be required. The following describes a method for determining instrument wavelength calibration using atmospheric absorption features that is efficient enough to be used for entire images on workstations. This study shows the effect of the surface reflectance on the calibration accuracy and the calibration history for the AVIRIS B spectrometer over the 1992 flight season.

Goetz, Alexander F. H.↗

Application of Synoptic Magnetograms for Prediction of Solar Activity Using Ensemble Kalman Filter

Solar activity predictions using the data assimilation approach have demonstrated great potential to build reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to a non-linear dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. These developments assume an empirical relationship between the mean toroidal magnetic field flux and the sunspot number. Estimated from the sunspot number series, variations of the toroidal field have been used to assimilate the data into the Parker-Kleeorin-Ruzmakin (PKR) dynamo model by applying the EnKF method. The dynamo model describes the evolution of the toroidal and poloidal components of the magnetic field and the magnetic helicity. Full-disk magnetograms provide more accurate and complete input data by constraining both the toroidal and poloidal global field components, but these data are available only for the last four solar cycles. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of limited time series observations on the accuracy of solar activity predictions. We present EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach.

Kitiashvili, Irina N.↗

Acquisition and analysis of accelerometer data

Acceleration data reduction must be undertaken with a complete understanding of the physical process, the means by which the data are acquired, and finally, the calculations necessary to put the data into a meaningful format. Discussed here are the acceleration sensor requirements dictated by the measurements desired. Sensor noise, dynamic range, and linearity will be determined from the physical parameters of the experiment. The digitizer requirements are discussed. Here the system from sensor to digital storage medium will be integrated, and rules of thumb for experiment duration, filter response, and number of bits are explained. Data reduction techniques after storage are also discussed. Time domain operations including decimating, digital filtering, and averaging are covered, as well as frequency domain methods, including windowing and the difference between power and amplitude spectra, and simple noise determination via coherence analysis. Finally, an example experiment using the Teledyne Geotech Model 44000 Seismometer to measure from 1 Hz to 10(exp -6) Hz is discussed. The sensor, data acquisition system, and example spectra are presented.

Verges, Keith R.↗

Aeroelastic stability and response of horizontal axis wind turbine blades

Coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine (HAWT) blade have been formulated. The analysis neglects blade-tower coupling. The final nonlinear equations have periodic coefficients. A new and convenient method of generating an appropriate time-dependent equilibrium position, required for the stability analysis, has been implemented and found to be computationally efficient. Steady-state response and stability boundaries for an existing (typical) HAWT blade are presented. Such stability boundaries have never been published in the literature. The results show that the isolated blade under study is basically stable. The tower shadow (wake) has a considerable effect on the out-of-plane response but leaves blade stability unchanged. Nonlinear terms can significantly affect linearized stability boundaries; however, they have a negligible effect on response, thus implying that a time-dependent equilibrium position (or steady-state response), based completely on the linear system, is appropriate for the type of HAWT blades under study.

Kottapalli, S. B. R.↗

Modeling and synthesis of multicomputer interconnection networks

The type of interconnection network employed has a profound effect on the performance of a multicomputer and multiprocessor design. Adequate models are needed to aid in the design and development of interconnection networks. A novel modeling approach using statistical and optimization techniques is described. This method represents an attempt to compare diverse interconnection network designs in a way that allows not only the best of existing designs to be identified but to suggest other, perhaps hybrid, networks that may offer better performance. Stepwise linear regression is used to develop a polynomial surface representation of performance in a (k+1) space with a total of k quantitative and qualitative independent variables describing graph-theoretic characteristics such as size, average degree, diameter, radius, girth, node-connectivity, edge-connectivity, minimum dominating set size, and maximum number of prime node and edge cutsets. Dependent variables used to measure performance are average message delay and the ratio of message completion rate to network connection cost. Response Surface Methodology (RSM) optimizes a response variable from a polynomial function of several independent variables. Steepest ascent path may also be used to approach optimum points.

Standley, Hilda M.↗

Application of Gaussian Mixture Regression for the Correction of Low Cost PM2.5 Monitoring Data in Accra, Ghana

Low-cost sensors (LCSs) for air quality monitoring have enormous potential to improve air quality data coverage in resource-limited parts of the world such as sub-Saharan Africa. LCSs, however, are affected by environment and source conditions. To establish high-quality data, LCSs must be collocated and calibrated with reference grade PM2.5 monitors. From March 2020, a low-cost PurpleAir PM2.5 monitor was collocated with a Met One Beta Attenuation Monitor 1020 in Accra, Ghana. While previous studies have shown that multiple linear regression (MLR) and random forest regression (RF) can improve accuracy and correlation between PurpleAir and reference data, MLR and RF yielded suboptimal improvement in the Accra collocation (R2 = 0.81 and R2 = 0.81, respectively). We present the first application of Gaussian mixture regression (GMR) to air quality data calibration and demonstrate improvement over traditional methods by increasing the collocated PM2.5 correlation and accuracy to R2 = 0.88 and MAE = 2.2 μg/cu. m. Gaussian mixture models (GMMs) are a probability density estimator and clustering method from which nonlinear regressions that tolerate missing inputs can be derived. We find that even when given missing inputs, GMR provides better correlation than MLR and RF performed with complete data. GMR also allows us to estimate calibration certainty. When evaluated, 95% confidence intervals agreed with reference PM2.5 data 96% of the time, suggesting that the model accurately assesses its own confidence. Additionally, clustering within the GMM is consistent with climate characteristics, providing confidence that the calibration approach can learn underlying relationships in data.

Sensors↗

Separated flow over bodies of revolution using an unsteady discrete-vorticity cross wake. Part 1: Theory and application

A method is developed to determine the flow field of a body of revolution in separated flow. The technique employed is the use of the computer to integrate various solutions and solution properties of the sub-flow fields which made up the entire flow field without resorting to a finite difference solution to the complete Navier-Stokes equations. The technique entails the use of the unsteady cross flow analogy and a new solution to the required two-dimensional unsteady separated flow problem based upon an unsteady, discrete-vorticity wake. Data for the forces and moments on aerodynamic bodies at low speeds and high angle of attack (outside the range of linear inviscid theories) such that the flow is substantially separated are produced which compare well with experimental data. In addition, three dimensional steady separation regions and wake vortex patterns are determined.

Marshall, F. J.↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗

Progress Towards Improved Analysis of TES X-ray Data Using Principal Component Analysis

The traditional method of applying a digital optimal filter to measure X-ray pulses from transition-edge sensor (TES) devices does not achieve the best energy resolution when the signals have a highly non-linear response to energy, or the noise is non-stationary during the pulse. We present an implementation of a method to analyze X-ray data from TESs, which is based upon principal component analysis (PCA). Our method separates the X-ray signal pulse into orthogonal components that have the largest variance. We typically recover pulse height, arrival time, differences in pulse shape, and the variation of pulse height with detector temperature. These components can then be combined to form a representation of pulse energy. An added value of this method is that by reporting information on more descriptive parameters (as opposed to a single number representing energy), we generate a much more complete picture of the pulse received. Here we report on progress in developing this technique for future implementation on X-ray telescopes. We used an 55Fe source to characterize Mo/Au TESs. On the same dataset, the PCA method recovers a spectral resolution that is better by a factor of two than achievable with digital optimal filters.

Low-temperature detectors A[center dot] Transition↗

Controls-structures-interaction dynamics during RCS control of the Orbiter/SRMS/SSF configuration

During the assembly flights of the Space Station Freedom (SSF), the Orbiter will either dock with the SSF and retract to the final berthed position, or will grapple the SSF using the Shuttle Remote Manipulator System (SRMS) and maneuver the SRMS coupled vehicles to their final berthed position. The SRMS method is expected to take approximately one to one and a half hours to complete and require periodic attitude corrections by either the Orbiter or the SSF reaction control system (RCS) or continuous control by a control moment gyro (CMG) system with RCS desaturation as required. Free drift of the attached vehicles is not currently thought to be acceptable because the desired system attitude will quickly deteriorate due to unbalanced gravity gradient and aerodynamic torques resulting in power generation problems, thermodynamic control problems, and communications problems. This paper deals with the simulation and control of the SRMS during trunnion/latch interaction dynamics and during RCS maneuvers. The SRMS servo drive joints have highly non-linear elastic characteristics which tend to degrade sensitive control strategies. In addition the system natural frequencies are extremely low and depend on the drive joint deflections and SRMS geometric position. The lowest mean period of oscillation for the Orbiter/SRMS/SSF(MB6) system in brakes hold mode positioned near the final berthed position is approximately 120 seconds. A detailed finite element model of the SRMS has been developed and used in a newly developed SRMS systems dynamics simulation to investigate the non-linear transient response dynamics of the Orbiter/SRMS/SSF systems. The present SRMS control strategy of brakes only recommended by the Charles Draper Labs is contrasted with a robust controller developed by the authors. The robust controller uses an optimal inear quadratic regulator (LQR) to optimally place the closed-loop poles of a multivariable continuous-time system within the common region of an open sector with the sector angle plus or minus 45 degrees from the negative real axis, and the left-hand side of a parallel to the imaginary axis in the complex s-plane. This guarantees that the critical damping ratio for the desired control modes is equal to or in excess of 0.707. The matrix sign function is used for solving the Riccati equations which appear in the controller design procedure. Fast and stable algorithms have recently been developed for the computation of the matrix sign function. Simulation results are given which demonstrate the potential CSI involvement for the current SRMS control system and the proposed control system.

Schliesing, J. A.↗

Carotid Intima Media Thickness in the Astronaut Corps: Association to Spacecraft

Background: Carotid Intima Media Thickness (CIMT) has been demonstrated to be predictive of future cardiovascular events. Within various populations, radiation exposure, stress, and physical confinement have all been linked to an increased CIMT. Recent research discovered CIMT was significantly increased in ten long duration astronauts from pre-flight to four days post flight. The relationship between spaceflight and CIMT is not understood and trends in CIMT within the larger astronaut population are unknown. Methods: In 2010, CIMT was offered as part of the astronaut annual exam at the JSC Flight Medicine Clinic using a standardized CIMT screening protocol and professional sonographers. Between 2010 and 2016, CIMT measurements were collected on 213 NASA astronauts and payload specialists. The values used in this retrospective chart review are the mean of the CIMT from the right and left. Spaceflight exposure was categorized based on the total number of days spent in space at the time of the ground-based ultrasound (0, 1-29, 30-100, 101-200, ≥200). Linear regression with generalized estimating equations were used to estimate the association between spaceflight exposures and CIMT. Results: 530 studies were completed among 213 astronauts with a mean of 2.5 studies (range 1-6) per astronaut over the six year period. As in other populations, CIMT was significantly associated with age; however, gender was not. While there was no significant direct correlation between total spaceflight exposure and CIMT found, astronauts with 30-100 spaceflight days and astronauts with greater than 100 spaceflight days had significantly increased CIMT over astronauts who had never flown (p=0.002 and p=<0.0001 respectively) after adjustment for age. Conclusion: Further work is needed to fully understand CIMT and its association to spaceflight. Current occupational surveillance activities are under way to study CIMT values in conjunction with other cardiovascular risk factors among astronauts as compared to the general population.

Suffredini, John↗

Design of satellite flexibility experiments

A preliminary study has been completed to begin development of a flight experiment to measure spacecraft control/flexible structure interaction. The work reported consists of two phases: identification of appropriate structural parameters which can be associated with flexibility phenomena, and suggestions for the development of an experiment for a satellite configuration typical of near-future vehicles which are sensitive to such effects. Recommendations are made with respect to the type of data to be collected and instrumentation associated with these data. The approach consists of developing the equations of motion for a vehicle possessing a flexible solar array, then linearizing about some nominal motion of the craft. A set of solutions are assumed for array deflection using a continuous normal mode method and important parameters are exposed. Inflight and ground based measurements are distinguished. Interrelationships between these parameters, measurement techniques, and input requirements are discussed which assure minimization of special vehicle maneuvers and optimization of data to be obtained during the normal flight sequence.

Kaplan, M. H.↗

The Application of Strainrange Partitioning Method to Multiaxial Creep-Fatigue Interaction

The method of strain range partitioning was applied to a series of torsional fatigue tests conducted on tubular 304 stainless steel specimens at 1200 F (649 C). Creep strain was superimposed on cycling strain, and the resulting strain range was partitioned into four components: (1) completely reversed plastic shear strain, (2) plastic shear strain followed by creep strain, (3) creep strain followed by plastic strain, and (4) completely reversed creep strain. Each strain component was related to the cyclic life of the material. The experimental procedure used to achieve strain partitioning is described, and the torsional test results are compared to those obtained from axial tests. The damaging effects of the individual strain components were expressed by a linear life fraction rule.

Zamrik, S. Y.↗

Final Report - Subcontract B629114

During my visits to LLNL August 20-24, 2018, and September 10-14, 2018, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from $hp$ adaptive finite element calculations, and is implemented in the PLTMG software package, version 13. This preconditioner typically requires 3 – 20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive definite elliptic equations, convection dominated convection-diffusion equations, and highly indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree $p$ as $p → ∞$, but but remains strong for $p ≤$ 9, which is the highest polynomial degree allowed in PLTMG , due to limitations of the numerical quadrature rules implemented in the software package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.

97 MATHEMATICS AND COMPUTING↗

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗