Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter space”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

An Advanced User Interface Approach for Complex Parameter Study Process Specification in the Information Power Grid

The creation of parameter study suites has recently become a more challenging problem as the parameter studies have now become multi-tiered and the computational environment has become a supercomputer grid. The parameter spaces are vast, the individual problem sizes are getting larger, and researchers are now seeking to combine several successive stages of parameterization and computation. Simultaneously, grid-based computing offers great resource opportunity but at the expense of great difficulty of use. We present an approach to this problem which stresses intuitive visual design tools for parameter study creation and complex process specification, and also offers programming-free access to grid-based supercomputer resources and process automation.

Yarrow, Maurice↗
NASA NTRS2019

Probing boson stars with extreme mass ratio inspirals

We propose to use gravitational waves from extreme mass ratio inspirals (EMRI), composed of a boson star and a supermassive black hole in the center of galaxies, as a new method to search for boson stars. Gravitational waves from EMRI have the advantage of being long-lasting within the frequency band of future space-based interferometer gravitational wave detectors and can accumulate large signal-to-noise ratio (SNR) for very sub-solar mass boson stars. Compared to gravitational waves from boson star binaries, which fall within the LIGO band, we find that much larger ranges of the mass and compactness of boson stars, as well as the underlying particle physics parameter space, can be probed by EMRI. We take tidal disruption of the boson stars into account and distinguish those which dissolve before the inner-most-stable-circular-orbit (ISCO) and those which dissolve after it. Due to the large number of cycles recorded, EMRIs can lead to a very precise mass determination of the boson star and distinguish it from standard astrophysical compact objects in the event of a signal. Tidal effects in inspiralling binary systems, as well as possible correlated electromagnetic signals, can also serve as potential discriminants.

Huai-Ke Guo↗
NASA NTRS2014

Bifurcation Analysis Using Rigorous Branch and Bound Methods

For the study of nonlinear dynamic systems, it is important to locate the equilibria and bifurcations occurring within a specified computational domain. This paper proposes a new approach for solving these problems and compares it to the numerical continuation method. The new approach is based upon branch and bound and utilizes rigorous enclosure techniques to yield outer bounding sets of both the equilibrium and local bifurcation manifolds. These sets, which comprise the union of hyper-rectangles, can be made to be as tight as desired. Sufficient conditions for the existence of equilibrium and bifurcation points taking the form of algebraic inequality constraints in the state-parameter space are used to calculate their enclosures directly. The enclosures for the bifurcation sets can be computed independently of the equilibrium manifold, and are guaranteed to contain all solutions within the computational domain. A further advantage of this method is the ability to compute a near-maximally sized hyper-rectangle of high dimension centered at a fixed parameter-state point whose elements are guaranteed to exclude all bifurcation points. This hyper-rectangle, which requires a global description of the bifurcation manifold within the computational domain, cannot be obtained otherwise. A test case, based on the dynamics of a UAV subject to uncertain center of gravity location, is used to illustrate the efficacy of the method by comparing it with numerical continuation and to evaluate its computational complexity.

Smith, Andrew P.↗
DOE OSTI2024

Search for a heavy neutral lepton that mixes predominantly with the tau neutrino

We report a search for a heavy neutral lepton (HNL) that mixes predominantly with ν τ . The search utilizes data collected with the Belle detector at the KEKB asymmetric energy e + e − collider. The data sample was collected at and just below the center-of-mass energies of the ϒ ( 4 S ) and ϒ ( 5 S ) resonances and has an integrated luminosity of 915 fb − 1 , corresponding to ( 836 ± 12 ) × 10 6 e + e − → τ + τ − events. We search for production of the HNL (denoted N ) in the decay τ − → π − N followed by its decay via N → μ + μ − ν τ . The search focuses on the parameter-space region in which the HNL is long-lived, so that the μ + μ − originate from a common vertex that is significantly displaced from the collision point of the KEKB beams. Consistent with the expected background yield, one event is observed in the data sample after application of all the event-selection criteria. We report limits on the mixing parameter of the HNL with the τ neutrino as a function of the HNL mass. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS1993

Chondrites, S asteroids, and space weathering: Thumping noises from the coffin?

Most of the spectral characteristics of ordinary chondrites and S-asteroids in the visible and infrared can be reduced to three numerical values. These values represent the depth of the absorption band resulting from octahedrally coordinated Fe(sup 2+), the reflectance at 0.56 microns and the slope of the continuum (as measured according to convention). By plotting these three characteristics, it is possible to immediately compare the spectral characteristics of large numbers of ordinary chondrites and S-asteroids. Commonality of spectral characteristics between these populations can thus be evaluated on the basis of overlap in position on three two-coordinate systems: albedo vs. band depth, band depth vs. slope, and slope vs. albedo. In order to establish identity, members of the two populations must overlap on all three of these independent parameter spaces. In this coordinate system, spectra of 23 ordinary chondrites (representing all metamorphic grades), and 39 S-asteroids were compared. It was found that there was no overlap between the two populations in terms of the slope vs. band depth parameters, nor were most chondrites identical to the S-asteroids with respect to the other criteria. However, the controversial question remains: Where are the parent bodies of the chondrites? Perhaps an even more critical question is: Where are our samples of the S-asteroids? Considering the geography of the asteroid belt and the theory that early solar-system electromagnetic induction heating differentiated protoasteroids in the inner portion of the main belt, it was suggested that although S-asteroids and ordinary chondrites have very similar mineralogy, the S-asteroids are mixtures of metallic nickel iron and silicates which resulted from magmatism induced by electromagnetic heating whereas chondrites were only slightly metamorphosed nebular condensates. In this scenario chondrites would have been derived from a population of bodies with thermal lag times so short that they were not subjected to melting during the phase of the electromagnetic induction heating event but only to various degrees of pervasive metamorphism. Furthermore, these objects would then have been too small to be observed and systematically included in the library of asteroidal spectra. It was also suggested that the parametric distribution of S-asteroid spectra could be reproduced by mixing various proportions of NiFe meteorite and achondritic materials. This has also been demonstrated in the laboratory.

Fanale, F. P.↗
NASA NTRS2000

Code IN Exhibits - Supercomputing 2000

The creation of parameter study suites has recently become a more challenging problem as the parameter studies have become multi-tiered and the computational environment has become a supercomputer grid. The parameter spaces are vast, the individual problem sizes are getting larger, and researchers are seeking to combine several successive stages of parameterization and computation. Simultaneously, grid-based computing offers immense resource opportunities but at the expense of great difficulty of use. We present ILab, an advanced graphical user interface approach to this problem. Our novel strategy stresses intuitive visual design tools for parameter study creation and complex process specification, and also offers programming-free access to grid-based supercomputer resources and process automation.

Yarrow, Maurice↗
NASA NTRS1992

SIMBAD quality-control

The astronomical database SIMBAD developed at the Centre de donnees astronomiques de Strasbourg presently contains 760,000 objects (stellar and non-stellar). It has the unique characteristic of being structured specifically for astronomical objects. All types of heterogeneous data (bibliographic references, measurements, and sets of identification) are connected with each object. The attributes that define quality of the database include the following. Reliability: cross-identification should not rely upon just exact values object coordinates. It also means that information attached to one simple object should be consistent. The existing data must be controlled in order to start with a reliable base and to cross-identify new data assuring the quality as data grows. Exhaustivity: delays between publication of new informations and their inclusion in the database should be as short as possible. The integrity of the database has to be maintained as data accumulates. Taking the amount of data into consideration and the rate of new data production, it is necessary to use automatic methods. One of the possibilities is to use multivariate data analysis. The factor-space is a n-dimensional relevancy space which is described by the n-axes representing a set of n subject matter headings; the words and phrases can be used to scale the axes and the documents are then a vector average of the terms within them. The application reported herein is based on the NASA-STI bibliographical database. The selected data concern astronomy, astrophysics, and space radiation (102,963 references from 1975 to 1991 included 8070 keywords). The F-space is built from this bibliographical data. By comparing the F-space position obtained from the NASA-STI keywords with the F-space position obtained from the SIMBAD references, the authors will be able to show whether it is possible to retrieve information with a restricted set of words only. If the comparison is valid, this will be a way to enter bibliographic information in the SIMBAD quality control process. Furthermore, it is possible to connect the physical measurements of stars from SIMBAD to literature concerning these stars from the NASA-STI abstracts. The physical properties of stars (e.g. UBV colors) are not randomly distributed. Stars are distributed among different clusters in a physical parameter space. The authors will show that there are some relations between this classification and the literature concerning these objects clusters in a factor space. They will investigate the nature of the relationship between the SIMBAD measurements and the bibliography. These would be new relationships that are not pre-established by an astronomer. In addition, the bibliography could be neutral information that can be used in combination with the measured parameters.

Lesteven, Soizick↗
NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗
NASA NTRS2015

Calibration of Predictor Models Using Multiple Validation Experiments

This paper presents a framework for calibrating computational models using data from several and possibly dissimilar validation experiments. The offset between model predictions and observations, which might be caused by measurement noise, model-form uncertainty, and numerical error, drives the process by which uncertainty in the models parameters is characterized. The resulting description of uncertainty along with the computational model constitute a predictor model. Two types of predictor models are studied: Interval Predictor Models (IPMs) and Random Predictor Models (RPMs). IPMs use sets to characterize uncertainty, whereas RPMs use random vectors. The propagation of a set through a model makes the response an interval valued function of the state, whereas the propagation of a random vector yields a random process. Optimization-based strategies for calculating both types of predictor models are proposed. Whereas the formulations used to calculate IPMs target solutions leading to the interval value function of minimal spread containing all observations, those for RPMs seek to maximize the models' ability to reproduce the distribution of observations. Regarding RPMs, we choose a structure for the random vector (i.e., the assignment of probability to points in the parameter space) solely dependent on the prediction error. As such, the probabilistic description of uncertainty is not a subjective assignment of belief, nor is it expected to asymptotically converge to a fixed value, but instead it casts the model's ability to reproduce the experimental data. This framework enables evaluating the spread and distribution of the predicted response of target applications depending on the same parameters beyond the validation domain.

Crespo, Luis G.↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Maximum likelihood algorithm using an efficient scheme for computing sensitivities and parameter confidence intervals

Improved techniques for estimating airplane stability and control derivatives and their standard errors are presented. A maximum likelihood estimation algorithm is developed which relies on an optimization scheme referred to as a modified Newton-Raphson scheme with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort compared to integrating the analytically-determined sensitivity equations or using a finite difference scheme. An aircraft estimation problem is solved using real flight data to compare MNRES with the commonly used modified Newton-Raphson technique; MNRES is found to be faster and more generally applicable. Parameter standard errors are determined using a random search technique. The confidence intervals obtained are compared with Cramer-Rao lower bounds at the same confidence level. It is observed that the nonlinearity of the cost function is an important factor in the relationship between Cramer-Rao bounds and the error bounds determined by the search technique.

Murphy, P. C.↗
NASA NTRS1986

Influence of transient baroclinic eddies on planetary-scale waves

Low-order models with a mean zonal flow forced by zonally symmetric heating and one planetary-scale wave forced by topography have multiple flow equilibria. Two stable equilibria are characterized by strong zonal flow with low wave amplitude (high index circulation) and weak zonal flow with a high wave amplitude (low index circulation) fixed in phase with the topography. These two states presumably represent normal zonal circulation and a blocking configuration, respectively. When a shorter, baroclinically unstable wave is dded to this low-order model, the planetary-scale wave no longer stays in stable equilibrium states. Instead, the long wave remains in weather regimes or preferred regions in phase space. These regimes are, in general, different from the equilibria of the model with the planetary-scale wave only. Thus, the short unstable wave adds some randomness to planetary-scale circulation and changes its position with respect to the planetary-scale topography. The Reinhold and Pierrehumbert model is further explored here by adding a long-wave in the thermal forcing and exploring wider parameter space. When the symmetric thermal forcing is weak so that the short wave is stable, the planetary scale wave is fixed in amplitude and phase by the asymmetric forcing. However, when the symmetric forcing is increased to an unstable level, the planetary-scale wave becomes less organized. With asymmetric thermal forcing only (i.e., no topography) the long-wave is randomly distributed through phase space. Thermal forcing seems to be less effective in organizing weather regimes than topographic forcing.

Branscome, L. E.↗
DOE OSTI2025

Leveraging differentiable programming in the inverse problem of neutron stars

Neutron stars (NSs) probe the high-density regime of the nuclear equation of state (EOS). However, inferring the EOS from observations of NSs is a computationally challenging task. Here, in this work, we efficiently solve this inverse problem by leveraging differential programming in two ways. First, we enable full Bayesian inference in under one hour of wall time on a GPU by using gradient-based samplers, without requiring pretrained machine learning emulators. Moreover, we demonstrate efficient scaling to high-dimensional parameter spaces. Second, we introduce a novel gradient-based optimization scheme that recovers the EOS of a given NS mass-radius curve. We demonstrate how our framework can reveal consistencies or tensions between nuclear physics and astrophysics. First, we show how the breakdown density of a metamodel description of the EOS can be determined from NS observations. Second, we demonstrate how degeneracies in EOS modeling using nuclear empirical parameters can influence the inverse problem during gradient-based optimization. Looking ahead, our approach opens up new theoretical studies of the relation between NS properties and the EOS, while effectively tackling the data analysis challenges brought by future detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

TCAD-Machine Learning Enabled TID Compact Model Development for Commercial SiC MOSFET

We propose a TCAD (Technology Computer Aided Design)-machine learning coupled approach that combines a TCAD tool (Charon), optimization/uncertainty quantification tool (Dakota), surrogate models, and Bayesian learning capabilities. The coupling approach is used for accurate modeling and calibration of total ionizing dose (TID) induced threshold voltage (V th ) shifts in Commercial-Off-The-Shelf (COTS) semiconductor devices and to develop physics-informed TID compact models. This versatile approach is applied to model the TID effect in an exemplar COTS 3.3 kV SiC power MOSFET (Metal-Oxide-Semiconductor Field-Effect Transistor). With the Charon-Dakota coupling, we can determine key device geometry and doping values based on device physics, which are difficult to obtain or not available for COTS devices but important for TCAD simulation; additionally, we can efficiently generate thousands of simulation results in a large parameter space, which makes it possible to develop data-driven surrogate models and perform Bayesian calibration. Utilizing the full tool-coupling approach, we achieve calibrated TCAD simulation models that accurately capture the average TID-induced V th shifts behavior with total doses and V th shifts saturation at high doses as observed in experimental data. More importantly, the calibrated TCAD simulations are obtained with determined TID model parameters (e.g., hole trap density and capture cross section) values that contain well quantified uncertainties. Furthermore, we can isolate and quantify the noises that are not captured by the TCAD models but exist in the measured data due to measurements and devices variabilities. Lastly, the calibrated surrogate models are used to develop physics-informed TID compact models. The method is generalizable to other devices and/or radiation conditions with few modifications and can provide well-determined uncertainties.

COTS↗
NASA NTRS2000

Trajectories in Close Proximity to Asteroids

Spacecraft motion in close proximity to irregularly shaped, rotating bodies such as asteroids presents a unique dynamical environment as compared to most space missions. There are several fundamental novelties in this environment that spacecraft must deal with. These include the possibility of orbital instabilities that can act over very short time spans (on the order of hours for some systems), possible non-uniform rotation of the central gravity field, divergence of traditional gravity field representations when close to the asteroid surface, dominance of perturbing forces, an extremely large asteroid model parameter space that must be prepared for in the absence of reliable information, and the possibility of employing new and novel trajectory control techniques such as hovering and repeated landings on the asteroid surface. An overview of how these novelties impact the space of feasible close proximity operations and how different asteroid model properties will affect their implementation is given. In so doing, four fundamental types of close proximity operations will be defined. Listed in order of increasing technical difficulty these are: (1) close, stable orbits; (2) low-altitude flyovers; (3) landing trajectories; and (4) hovering trajectories. The feasibility and difficulty of implementing these operations will vary as a function of the asteroid shape, size, density, and rotation properties, and as a function of the spacecraft navigation capability. Additional information is contained in the original extended abstract.

Scheeres, D. J.↗
DOE OSTI2025

Scalable Bayesian Physics-Informed Kolmogorov-Arnold Networks

Uncertainty quantification (UQ) plays a pivotal role in scientific machine learning, especially when surrogate models are used to approximate complex systems. Although multilayer perceptions (MLPs) are commonly employed as surrogates, they often suffer from overfitting due to their large number of parameters. Kolmogorov-Arnold networks (KANs) offer an alternative solution with fewer parameters. However, gradient-based inference methods, such as Hamiltonian Monte Carlo (HMC), may result in computational inefficiency when applied to KANs, especially for large-scale datasets, due to the high cost of back-propagation. To address these challenges, we propose a novel approach, combining the dropout Tikhonov ensemble Kalman inversion (DTEKI) with Chebyshev KANs. This gradient-free method effectively mitigates overfitting and enhances numerical stability. In addition, we incorporate the active subspace method to reduce the parameter-space dimensionality, allowing us to improve the accuracy of predictions and obtain more reliable uncertainty estimates. Extensive experiments demonstrate the efficacy of our approach in various test cases, including scenarios with large datasets and high noise levels. Our results show that the new method achieves comparable or better accuracy, much higher efficiency as well as stability compared to HMC, in addition to scalability. Moreover, by leveraging the low-dimensional parameter subspace, our method preserves prediction accuracy while substantially reducing further the computational cost.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

A methodology for airplane parameter estimation and confidence interval determination in nonlinear estimation problems

An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. With the fitted surface, sensitivity information can be updated at each iteration with less computational effort than that required by either a finite-difference method or integration of the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, and thus provides flexibility to use model equations in any convenient format. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. The degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels and to predict the degree of agreement between CR bounds and search estimates.

Murphy, P. C.↗
NASA NTRS1994

Numerical simulations of upstream blocking, columnar disturbances, and bores in stably stratified shear flows over an obstacle

A two-dimensional, nonhydrostatic, elastic numerical model has been used to study the generation of gravity waves for a stably stratified shear flow over an obstacle. When a low-level wind shear is included in the simulation, we find that the predictions for noticeable upstream effects based on Froude number for a uniform flow are no longer accurate. Upstream effects are encountered in the form of upstream propagating columnar disturbances and internal bores away from the obstacle. The limited parameter space studies conducted in this study suggest that the ratio of the shear depth to the obstacle heigh (d/H), the obstacle aspect ratio (H/L), and the Froude number (U/NH) are instrumental in determining the strength and the existence of these upstream disturbances. Thus, the present theoretical and empirical understanding of the importance of the Froude number for determining the nature of upstream effects should be modified substantially to include additional nondimensional parameters when shear is present.

Chen, Chaing↗