Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter continuation 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 145 records · Page 8

Adaptive model reduction for continuous systems via recursive rational interpolation

A method for adaptive identification of reduced-order models for continuous stable SISO and MIMO plants is presented. The method recursively finds a model whose transfer function (matrix) matches that of the plant on a set of frequencies chosen by the designer. The algorithm utilizes the Moving Discrete Fourier Transform (MDFT) to continuously monitor the frequency-domain profile of the system input and output signals. The MDFT is an efficient method of monitoring discrete points in the frequency domain of an evolving function of time. The model parameters are estimated from MDFT data using standard recursive parameter estimation techniques. The algorithm has been shown in simulations to be quite robust to additive noise in the inputs and outputs. A significant advantage of the method is that it enables a type of on-line model validation. This is accomplished by simultaneously identifying a number of models and comparing each with the plant in the frequency domain. Simulations of the method applied to an 8th-order SISO plant and a 10-state 2-input 2-output plant are presented. An example of on-line model validation applied to the SISO plant is also presented.

Lilly, John H.↗

The Fast Evolution of SN 2010bh Associated with XRF 100316D

The first observational evidence of a connection between supernovae (SNe) and gamma-ray bursts (GRBs) was found about a decade ago. Since then, only half a dozen spectroscopically confirmed associations have been discovered and XRF 1003160 associated with the type-Ic SN 20lObh is among the latest. Aims. We constrain the progenitor radius, the host-galaxy extinction, and the physical parameters of the explosion of XRF l00316D/SN 20lObh at z 0.059. We study the SN brightness and colours in the context of GRB-SNe. Methods. We began observations with the Gamma-Ray burst Optical and Near-infrared Detector (GROND) 12 hours after the GRB trigger and continued until 80 days after the burst. GROND provided excellent photometric data in six filter bands covering a wavelength range from approximately 350 to 1800 nm, significantly expanding the pre-existing data set for this event. Combining GROND and Swift data, the early broad-band spectral energy distribution (SED) is modelled with a blackbody and afterglow component attenuated by dust and gas absorption. The temperature and radius evolution of the thermal component are analysed and combined with earlier measurements available from the literature. Templates of SN 1998bw are fitted to the SN itself to directly compare the lightcurve properties. Finally, a two-component parametrized model is fitted to the quasi-bolometric light curve. which delivers physical parameters of the explosion. Results. The best-fit models to the broad-band SEDs imply moderate reddening along the line of sight through the host galaxy (A(sub v.host = 1.2 +/- 0.1 mag). Furthermore, the parameters of the blackbody component reveal a cooling envelope at an apparent initial radius of 7 x 10(exp 11) cm, which is compatible with a dense wind surrounding a Wolf-Rayet star. A multicolor comparison shows that SN 2010bh is 60 - 70% as bright as SN 1998bw. Reaching maximum brightness at 8 - 9 days after the burst in the blue bands, SN 20lObh proves to be the most rapidly evolving GRB-SNe to date. Modelling of the quasi-bolometric light curve yields M(sub Ni) = 0.21 +/- 0.03 solar M and M(sub ej) = 2.6 +/- 0.2 solar M, typical of values within the GRB-SN population. The kinetic energy is E(sub k) = (2.4 +/- 0.7) x 10(exp 52) erg, which is making this SN the second most energetic GRB-SN after SN 1998bw. Conclusions. This supernova has one of the earliest peaks ever recorded and thereafter fades more rapidly than other GRB-SNe, hypernovae, or typical type-Ic SNe. This could be explained by a thin envelope expanding at very high velocities, which is therefore unable to retain the gamma-rays that would prolong the duration of the SN event.

Olivares E., F.↗

Computational Methods for Stability and Control (COMSAC): The Time Has Come

Powerful computational fluid dynamics (CFD) tools have emerged that appear to offer significant benefits as an adjunct to the experimental methods used by the stability and control community to predict aerodynamic parameters. The decreasing costs for and increasing availability of computing hours are making these applications increasingly viable as time goes on and the cost of computing continues to drop. This paper summarizes the efforts of four organizations to utilize high-end computational fluid dynamics (CFD) tools to address the challenges of the stability and control arena. General motivation and the backdrop for these efforts will be summarized as well as examples of current applications.

Hall, Robert M.↗

Global Clear-Sky Surface Skin Temperature from Multiple Satellites Using a Single-Channel Algorithm with Angular Anisotropy Corrections

Surface skin temperature (T(sub s)) is an important parameter for characterizing the energy exchange at the ground/water-atmosphere interface. The Satellite ClOud and Radiation Property retrieval System (SatCORPS) employs a single-channel thermal-infrared (TIR) method to retrieve T(sub s) over clear-sky land and ocean surfaces from data taken by geostationary Earth orbit (GEO) and low Earth orbit (LEO) satellite imagers. GEO satellites can provide somewhat continuous estimates of T(sub s) over the diurnal cycle in non-polar regions, while polar T(sub s) retrievals from LEO imagers, such as the Advanced Very High Resolution Radiometer (AVHRR), can complement the GEO measurements. The combined global coverage of remotely sensed T(sub s), along with accompanying cloud and surface radiation parameters, produced in near-realtime and from historical satellite data, should be beneficial for both weather and climate applications. For example, near-realtime hourly T(sub s) observations can be assimilated in high-temporal-resolution numerical weather prediction models and historical observations can be used for validation or assimilation of climate models. Key drawbacks to the utility of TIR-derived T(sub s) data include the limitation to clear-sky conditions, the reliance on a particular set of analyses/reanalyses necessary for atmospheric corrections, and the dependence on viewing and illumination angles. Therefore, T(sub s) validation with established references is essential, as is proper evaluation of T(sub s) sensitivity to atmospheric correction source. This article presents improvements on the NASA Langley GEO satellite and AVHRR TIR-based T(sub s) product that is derived using a single-channel technique. The resulting clear-sky skin temperature values are validated with surface references and independent satellite products. Furthermore, an empirically adjusted theoretical model of satellite land surface temperature (LST) angular anisotropy is tested to improve satellite LST retrievals. Application of the anisotropic correction yields reduced mean bias and improved precision of GOES-13 LST relative to independent Moderate-resolution Imaging Spectroradiometer (MYD11_L2) LST and Atmospheric Radiation Measurement Program ground station measurements. It also significantly reduces inter-satellite differences between LSTs retrieved simultaneously from two different imagers. The implementation of these universal corrections into the SatCORPS product can yield significant improvement in near-global-scale, near-realtime, satellite-based LST measurements. The immediate availability and broad coverage of these skin temperature observations should prove valuable to modelers and climate researchers looking for improved forecasts and better understanding of the global climate model.

Scarino, Benjamin R.↗

Ensemble‐Based Spatially Distributed CLM5 Hydrological Parameter Estimation for the Continental United States

Abstract One of the major challenges in large‐domain hydrological modeling efforts lies in the estimation of spatially distributed hydrological parameters while simultaneously accounting for their associated uncertainties. Addressing this challenge is particularly difficult in ungauged locations. With growing societal demands for large‐scale streamflow projections to inform water resource management and long‐term planning, evaluating and constraining hydrological parameter uncertainty is increasingly vital. This study introduces a hybrid regionalization approach to enhance hydrological predictions of the Community Land Model version 5 (CLM5) across the Continental United States (CONUS), with a total of 50,629 1/8° grid cells. This hybrid method combines the strengths of two existing techniques: parameter regionalization and streamflow signature regionalization. It identifies ensemble behavioral parameters for each 1/8° grid cell across the CONUS domain, tailored to three distinct streamflow signatures focused on low flows, high flows, and annual water balance. Evaluating this hybrid method for 464 CAMELS (Catchment Attributes and Meteorology for Large‐sample Studies) basins demonstrates a significant improvement in CLM5 hydrological predictions, even in challenging arid regions. In CONUS applications, the derived spatially distributed parameter sets capture both spatial continuity and variation of parameters, highlighting their heterogeneous nature within specific regions. Overall, this hybrid regionalization approach offers a promising solution to the complex task of improving hydrological modeling over large domains for important hydrological applications.

CLM5↗

Discovery of Nine Gamma-Ray Pulsars in Fermi-Lat Data Using a New Blind Search Method

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs Jl803-2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6 x 10(exp 35) ergs per second and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J010622+3749, Jl620-4927, Jl746-3239, J2028+3332,J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b| greater than 10 degrees). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2x 10(exp 11)G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3 x l0(exp 33) erg per second) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars.

Celik-Tinmaz, Ozlem↗

Spin Diagnostics Using Mott Polarimeter for the EIC Pre-Injector

Polarization measurement and front-end commissioning are critical steps in the EIC preinjector program. We designed a spin-diagnostics beamline at around 4 MeV for front-end commissioning. This paper describes a high-energy Mott polarimeter, including Mott scattering, the Sherman function, and extrapolation toward single elastic scattering in finite-thickness foils. We summarize the key equations used to calculate essential design parameters, including detection time, effective Sherman function, differential cross section, and Mott polarimeter figure of merit. We also present the diagnostics beamline lattice design and kicker design for continuous spin-polarization monitoring. These notes provide concepts, formulas, and methods for diagnostic-beamline design and commissioning measurements.

43 PARTICLE ACCELERATORS↗

A Comparative Study of Probability Collectives Based Multi-agent Systems and Genetic Algorithms

We compare Genetic Algorithms (GA's) with Probability Collectives (PC), a new framework for distributed optimization and control. In contrast to GA's, PC-based methods do not update populations of solutions. Instead they update an explicitly parameterized probability distribution p over the space of solutions. That updating of p arises as the optimization of a functional of p. The functional is chosen so that any p that optimizes it should be p peaked about good solutions. The PC approach works in both continuous and discrete problems. It does not suffer from the resolution limitation of the finite bit length encoding of parameters into GA alleles. It also has deep connections with both game theory and statistical physics. We review the PC approach using its motivation as the information theoretic formulation of bounded rationality for multi-agent systems. It is then compared with GA's on a diverse set of problems. To handle high dimensional surfaces, in the PC method investigated here p is restricted to a product distribution. Each distribution in that product is controlled by a separate agent. The test functions were selected for their difficulty using either traditional gradient descent or genetic algorithms. On those functions the PC-based approach significantly outperforms traditional GA's in both rate of descent, trapping in false minima, and long term optimization.

Huang, Chien-Feng↗

Broad Search for Unstable Resonant Orbits in the Planar Circular Restricted Three-Body Problem

Unstable resonant orbits in the circular restricted three-body problem have increasingly been used for trajectory design using optimization and invariant manifold techniques.In this study, several methods for computing these unstable resonant orbits are explored including flyby maps, continuation from two-body models, and grid searches. Families of orbits are computed focusing on the Jupiter-Europa system, and their characteristics are explored. Different parameters such as period and stability are examined for each set of resonantor bits, and the continuation of several specific orbits is explored in more detail.

unstable resonant orbits↗

Polynomial chaos expansions on principal geodesic Grassmannian submanifolds for surrogate modeling and uncertainty quantification

In this work we introduce a manifold learning-based surrogate modeling framework for uncertainty quantification in high-dimensional stochastic systems. Our first goal is to perform data mining on the available simulation data to identify a set of low-dimensional (latent) descriptors that efficiently parameterize the response of the high-dimensional computational model. To this end, we employ Principal Geodesic Analysis on the Grassmann manifold of the response to identify a set of disjoint principal geodesic submanifolds, of possibly different dimension, that captures the variation in the data. Since operations on the Grassmann require the data to be concentrated, we propose an adaptive algorithm based on Riemannian K-means and the minimization of the sample Fréchet variance on the Grassmann manifold to identify “local” principal geodesic submanifolds that represent different system behavior across the parameter space. Polynomial chaos expansion is then used to construct a mapping between the random input parameters and the projection of the response on these local principal geodesic submanifolds. Here, the method is demonstrated on four test cases, a toy-example that involves points on a hypersphere, a Lotka-Volterra dynamical system, a continuous-flow stirred-tank chemical reactor system, and a two-dimensional Rayleigh-Bénard convection problem.

42 ENGINEERING↗

Analysis of the surface load and radiated sound of a vibrating airfoil with application to the experiment of Brooks

A method is developed for calculating the surface load and radiated sound from a vibrating surface in a compressible viscous fluid. The method is applied to a thin two-dimensional elliptic cross-section. For large values of the viscous diffusion parameter, the surface load tends to an elliptic distribution in agreement with the results of inviscid theory when edge pressure continuity is enforced. For thin surfaces, the surface load is insensitive to variations in the thickness ratio. A three-dimensional spectral technique is developed to calculate the inviscid surface load and radiated sound from a thin vibrating airfoil. The inviscid theory predicts the correct form of the far field sound pressure and its phase. The actual levels are somewhat sensitive to the choice of theoretical spanwise surface pressure mode but are in better agreement with the experiment than the surface pressure. The comparison of theoretical and experimental surface pressure indicates that the viscous theory, used to validate the inviscid theory, is either inadequate or there is a source of experimental error.

Yates, J. E.↗

Development of Airport Surface Required Navigation Performance (RNP)

The U.S. and international aviation communities have adopted the Required Navigation Performance (RNP) process for defining aircraft performance when operating the en-route, approach and landing phases of flight. RNP consists primarily of the following key parameters - accuracy, integrity, continuity, and availability. The processes and analytical techniques employed to define en-route, approach and landing RNP have been applied in the development of RNP for the airport surface. To validate the proposed RNP requirements several methods were used. Operational and flight demonstration data were analyzed for conformance with proposed requirements, as were several aircraft flight simulation studies. The pilot failure risk component was analyzed through several hypothetical scenarios. Additional simulator studies are recommended to better quantify crew reactions to failures as well as additional simulator and field testing to validate achieved accuracy performance, This research was performed in support of the NASA Low Visibility Landing and Surface Operations Programs.

Cassell, Rick↗

Discovery Planetary Mission Operations Concepts

The NASA Discovery Program of small planetary missions will provide opportunities to continue scientific exploration of the solar system in today's cost-constrained environment. Using a multidisciplinary team, JPL has developed plans to provide mission operations within the financial parameters established by the Discovery Program. This paper describes experiences and methods that show promise of allowing the Discovery Missions to operate within the program cost constraints while maintaining low mission risk, high data quality, and reponsive operations.

low cost planetary exploration Discovery mission o↗

An analysis of the threshold method for measuring area-average rainfall

Experimental evidence shows that the area-average rain rate and the fractional area covered by rain rate exceeding a fixed threshold are highly correlated; that is, are highly linearly related. A precise theoretical explanation of this fact is given. The explanation is based on the observation that rain rate has a mixed distribution, one that is a mixture of a discrete distribution and a continuous distribution. Under a homogeneity assumption, the slope of the linear relationship depends only on the continuous part of the distribution and as such is found to be markedly immune to parameter changes. This is illustrated by certain slope surfaces obtained from three specific distributions. The threshold level can be chosen in an optimal way by minimizing a certain distance function defined over the threshold range. In general, the threshold level should be not too far from the mean rain rate conditional on rain. The so-called threshold method advocates measuring rainfall from fractional area exploiting the observed linear relationship of the later with the area average rain rate. The method is potentially useful for the estimation of rainfall from space via satellites.

Kedem, Benjamin↗

Approaches to the utilization of space for CO2 research

A system study of the potential of space technology to monitor climate changes and improve the understanding of the coupling beteen CO2 and cloud cover is summarized. The basis for the study is scientific data requirements pertinent to the U.S. Department of Energy's CO2 Research Program. The capabilities of space-based sensor systems are matched to meet these requirements. New sensor system concepts are identified, including a Space Shuttle-launched recalibration package to provide for continuity of measurement and recalibration between satellites, a high-orbit radiation budget satellite, or parallax sensor to measure cloud altitude, a passive method for the direct measurements of CO2 and a high-altitude powered platform to monitor select regional parameters. Space-based sensor systems that could be the development focus for the time frames of 0-5, 5-10 and 10-20 years are recommended.

Glaser, P.↗

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Numerical solutions of jet-induced mixing in a partially full cryogenic tank are presented. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear-free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to solve the nondimensional form of steady state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of above stated parameters on the condensation Nusselt and Stanton numbers which characterize the steady-state interface condensation process are investigated. Detailed analysis to gain a better understanding of the fundamentals of fluid mixing and interface condensation is performed.

Lin, Chin-Shun↗

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Numerical solutions of jet-induced mixing in a partially full cryogenic tank are presented. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear-free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to solve the nondimensional form of steady state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of above stated parameters on the condensation Nusselt and Stanton numbers which characterize the steady-state interface condensation process are investigated. Detailed analysis to gain a better understanding of the fundamentals of fluid mixing and interface condensation is performed.

Lin, Chin-Shun↗

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗