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 253 records · Page 14

Sensitivity Testing of the NSTAR Ion Thruster

During the Extended Life Test of the DS1 flight spare ion thruster, the engine was subjected to sensitvity testing in order to characterize the macroscopic dependence of discharge chamber sensitivity to a +\-3% vatiation in main flow, cathode flow and beam current, and to +\5% variation in beam and accelerator voltage, was determined for the minimum- (THO), half- (TH8) and full power (TH15) throttle levels. For each power level investigared, 16 high/low operating conditions were chosen to vary the flows, beam current, and grid voltages in in a matrix that mapped out the entire parameter space. The matrix of data generated was used to determine the partial derivative or senitivity of the dependent parameters--discharge voltage, discharge current, discharge loss, double-to-single-ion current ratio, and neutralizer-keeper voltage--to the variation in the independent parameters--main flow, cathode flow, beam current, and beam voltage. The sensititivities of each dependent parameter with respect to each independent parameter were determined using a least-square fit routine. Variation in these sensitivities with thruster runtime was recorded over the duration of the ELT, to detemine if discharge performance changed with thruster wear. Several key findings have been ascertained from the sensitivity testing. Discharge operation is most sensitve to changes in cathode flow and to a lesser degree main flow. The data also confirms that for the NSTAR configuration plasma production is limited by primary electron input due to the fixed neutral population. Key sensitivities along with their change with thruster wear (operating time) will be presented. In addition double ion content measurements with an ExB probe will also be presented to illustrate beam ion production and content sensitivity to the discharge chamber operating parameteres.

electric propulsion↗

Searching for Strongly-Interacting Dark Matter with the Heavy Photon Search Experiment

The Heavy Photon Search Experiment (HPS) is a fixed-target experiment at Jefferson Lab’s Hall B, designed to explore a hidden sector (HS) of particles containing dark matter and a new force mediator known as the “heavy photon” (A'). The A' is a massive spin-1 gauge boson associated with a new U (1)D symmetry in the HS that kinetically mixes with the Standard Model () photon with a weak coupling strength parameterized by ¿, with ¿2 ~ 10-2 -10-10. HPS utilizes a high-intensity electron beam on a thin tungsten target to produce heavy photons in the MeV-GeV mass range via “dark bremsstrahlung,” a process analogous to SM bremsstrahlung but suppressed by ¿2. The A' can decay resonantly to SM leptons, allowing HPS to conduct both mass resonance searches for prompt decays (large ¿) and displaced vertex searches for long-lived particles (small ¿). In addition to the minimal A' model, HPS probes more complex extensions such as the QCD-like strongly-interacting massive particles (SIMPs) HS containing “dark” pions (pD) and vector mesons (VD), with pD as dark matter candidates. These particles introduce new thermal dark matter freeze- out scenarios and visible signals through long-lived VD decays to SM leptons, which are accessible to HPS. This analysis conducted a displaced vertex search for VD ¿ e-e+ in the mass range 30 MeV to 124 MeV and ¿ between 10-6 < ¿ < 10-2 using data from the 2016 Engineering Run (10.753 nb-1) at 2.3 GeV. Unlike the minimal A' search, SIMP signal kinematics required new approaches to signal normalization and SM background rejection. The strongest signal evidence was a local p-value of 0.01317 for mVD = 119 MeV, corresponding to a global significance of 0.9s. Although no signal was found, this search excluded a region of the SIMP parameter space at 90 % confidence. This work demonstrates HPS’s competitive capability to probe SIMP sectors within cosmologically significant parameters and introduces a new method for HPS displaced vertex searches using track vertical impact parameter cuts

Spellman, Alic [Univ. of California, Santa Cruz, C↗

The dependence of the shape and stability of captive rotating drops on multiple parameters

Asymptotic and numerical techniques in bifurcation theory are applied to the Young-Laplace equation governing meniscus shape in order to analyze the dependence of the shape and stability of rigidly rotating drops held captive between corotating solid faces on multiple parameters. Asymptotic analysis of the evolution of drop shape from the cylindrical as a function of distance between the solid faces, drop volume, rotational Bond number and gravitational Bond number shows that some shape bifurcations from cylinders to wavy, axisymmetric menisci are ruptured by small changes in drop volume or gravity. Computer calculations of axisymmetric drop shapes based on a finite element representation of the interface and numerical algorithms for tracking shape families and singular points are then used to map drop stability for the four-dimensional parameter space. The results of the asymptotic and numerical analyses are shown to agree well within the limited range of parameters where the asymptotic analysis is valid.

Ungar, L. H.↗

Switching LPV Control for High Performance Tactical Aircraft

This paper examines a switching Linear Parameter-Varying (LPV) control approach to determine if it is practical to use for flight control designs within a wide angle of attack region. The approach is based on multiple parameter-dependent Lyapunov functions. The full parameter space is partitioned into overlapping subspaces and a family of LPV controllers are designed, each suitable for a specific parameter subspace. The hysteresis switching logic is used to accomplish the transition among different parameter subspaces. The proposed switching LPV control scheme is applied to an F-16 aircraft model with different actuator dynamics in low and high angle of attack regions. The nonlinear simulation results show that the aircraft performs well when switching among different angle of attack regions.

Lu, Bei↗

Adjoint Methods for Adjusting Three-Dimensional Atmosphere and Surface Properties to Fit Multi-Angle/Multi-Pixel Polarimetric Measurements

This paper derives an efficient procedure for using the three-dimensional (3D) vector radiative transfer equation (VRTE) to adjust atmosphere and surface properties and improve their fit with multi-angle/multi-pixel radiometric and polarimetric measurements of scattered sunlight. The proposed adjoint method uses the 3D VRTE to compute the measurement misfit function and the adjoint 3D VRTE to compute its gradient with respect to all unknown parameters. In the remote sensing problems of interest, the scalar-valued misfit function quantifies agreement with data as a function of atmosphere and surface properties, and its gradient guides the search through this parameter space. Remote sensing of the atmosphere and surface in a three-dimensional region may require thousands of unknown parameters and millions of data points. Many approaches would require calls to the 3D VRTE solver in proportion to the number of unknown parameters or measurements. To avoid this issue of scale, we focus on computing the gradient of the misfit function as an alternative to the Jacobian of the measurement operator. The resulting adjoint method provides a way to adjust 3D atmosphere and surface properties with only two calls to the 3D VRTE solver for each spectral channel, regardless of the number of retrieval parameters, measurement view angles or pixels. This gives a procedure for adjusting atmosphere and surface parameters that will scale to the large problems of 3D remote sensing. For certain types of multi-angle/multi-pixel polarimetric measurements, this encourages the development of a new class of three-dimensional retrieval algorithms with more flexible parametrizations of spatial heterogeneity, less reliance on data screening procedures, and improved coverage in terms of the resolved physical processes in the Earth's atmosphere.

linearization↗

Adjoint Methods for Adjusting Three-Dimensional Atmosphere and Surface Properties to Fit Multi-Angle Multi-Pixel Polarimetric Measurements

This paper derives an efficient procedure for using the three-dimensional (3D) vector radiative transfer equation (VRTE) to adjust atmosphere and surface properties and improve their fit with multi-angle/multi-pixel radiometric and polarimetric measurements of scattered sunlight. The proposed adjoint method uses the 3D VRTE to compute the measurement misfit function and the adjoint 3D VRTE to compute its gradient with respect to all unknown parameters. In the remote sensing problems of interest, the scalar-valued misfit function quantifies agreement with data as a function of atmosphere and surface properties, and its gradient guides the search through this parameter space. Remote sensing of the atmosphere and surface in a three-dimensional region may require thousands of unknown parameters and millions of data points. Many approaches would require calls to the 3D VRTE solver in proportion to the number of unknown parameters or measurements. To avoid this issue of scale, we focus on computing the gradient of the misfit function as an alternative to the Jacobian of the measurement operator. The resulting adjoint method provides a way to adjust 3D atmosphere and surface properties with only two calls to the 3D VRTE solver for each spectral channel, regardless of the number of retrieval parameters, measurement view angles or pixels. This gives a procedure for adjusting atmosphere and surface parameters that will scale to the large problems of 3D remote sensing. For certain types of multi-angle/multi-pixel polarimetric measurements, this encourages the development of a new class of three-dimensional retrieval algorithms with more flexible parametrizations of spatial heterogeneity, less reliance on data screening procedures, and improved coverage in terms of the resolved physical processes in the Earth׳s atmosphere.

adjoint methods↗

CONSTELL: NASA's Satellite Constellation Model

The CONSTELL program represents an initial effort by the orbital debris modeling group at NASA/JSC to address the particular issues and problems raised by the presence of LEO satellite constellations. It was designed to help NASA better understand the potential orbital debris consequences of having satellite constellations operating in the future in LEO. However, it could also be used by constellation planners to evaluate architecture or design alternatives that might lessen debris consequences for their constellation or lessen the debris effects on other users of space. CONSTELL is designed to perform debris environment projections rapidly so it can support parametric assessments involving either the constellations themselves or the background environment which represents non-constellation users of the space. The projections need to be calculated quickly because a number of projections are often required to adequately span the parameter space of interest. To this end CONSTELL uses the outputs of other NASA debris environment models as inputs, thus doing away with the need for time consuming upfront calculations. Specifically, CONSTELL uses EVOLVE or ORDEM96 debris spatial density results as its background environment, debris cloud snapshot templates to simulate debris cloud propagation, and time dependent orbit profiles of the intact non- functional constellation spacecraft and upper stages. In this paper the environmental consequences of the deployment of particular LEO satellite constellations using the CONSTELL model will be evaluated. Constellations that will undergo a parametric assessment will reflect realistic parameter values. Among other results the increase in loss rate of non-constellation spacecraft, the number of collisions involving constellation elements, and the replacement rate of constellation satellites as a result of debris impact will be presented.

Theall, Jeffrey R.↗

Reduced-dimension Bayesian optimization for model calibration of transient vapor compression cycles

Development and calibration of first-principles dynamic models of vapor compression cycles (VCCs) is of critical importance for applications that include control design and fault detection and diagnostics. Nevertheless, the inherent complexity of models that are represented by large systems of differential–algebraic equations leads to significant challenges for model calibration processes that utilize classical gradient-based methods. Bayesian optimization (BO) is a sample-efficient and gradient-free approach using a probabilistic surrogate model and optimal search over a feasible parameter space. Despite the benefits of BO in reducing computational costs, challenges remain in dealing with a high-dimensional calibration task resulting from a large set of parameters that have significant impacts on system behavior and need to be calibrated simultaneously. This paper presents a reduced-dimension BO framework for calibrating transient VCCs models where the calibration space is projected to a low-dimensional subspace for accelerating convergence of the solution algorithm and consequently reducing the number of transient simulations. The proposed approach was demonstrated via two case studies associated with different VCC applications where 10 parameters were calibrated in each case using laboratory measurements. The reduced-dimension BO framework only required 1 / 8 th of the iterations associated with a standard BO method that deals with high-dimensional calibration parameters for converged solutions and yielded comparable accuracy. Furthermore, both calibrated models revealed significant accuracy improvements compared to uncalibrated models.

Ma, Jiacheng↗

Conditional diffusion machine-learning framework for mapping valence electron distribution from convergent beam electron diffraction

Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Discovering Planetary Nebula Geometries: Explorations with a Hierarchy of Models

Astronomical objects known as planetary nebulae (PNe) consist of a shell of gas expelled by an aging medium-sized star as it makes its transition from a red giant to a white dwarf. In many cases this gas shell can be approximately described as a prolate ellipsoid. Knowledge of the physics of ionization processes in this gaseous shell enables us to construct a model in three dimensions (3D) called the Ionization-Bounded Prolate Ellipsoidal Shell model (IBPES model). Using this model we can generate synthetic nebular images, which can be used in conjunction with Hubble Space Telescope (HST) images of actual PNe to perform Bayesian model estimation. Since the IBPES model is characterized by thirteen parameters, model estimation requires the search of a 13-dimensional parameter space. The 'curse of dimensionality,' compounded by a computationally intense forward problem, makes forward searches extremely time-consuming and frequently causes them to become trapped in local solutions. We find that both the speed and of the search can be improved by judiciously reducing the dimensionality of the search space. Our basic approach employs a hierarchy of models of increasing complexity that converges to the IBPES model. Earlier studies establish that a hierarchical sequence converges more quickly, and to a better solution, than a search relying only on the most complex model. Here we report results for a hierarchy of five models. The first three models treat the nebula as a 2D image, while the last two models explore its characteristics as a 3D object and enable us to characterize the physics of the nebula. This five-model hierarchy is applied to HST images of ellipsoidal PNe to estimate their geometric properties and gas density profiles.

Huyser, Karen A.↗

Astrophysical Model Selection in Gravitational Wave Astronomy

Theoretical studies in gravitational wave astronomy have mostly focused on the information that can be extracted from individual detections, such as the mass of a binary system and its location in space. Here we consider how the information from multiple detections can be used to constrain astrophysical population models. This seemingly simple problem is made challenging by the high dimensionality and high degree of correlation in the parameter spaces that describe the signals, and by the complexity of the astrophysical models, which can also depend on a large number of parameters, some of which might not be directly constrained by the observations. We present a method for constraining population models using a hierarchical Bayesian modeling approach which simultaneously infers the source parameters and population model and provides the joint probability distributions for both. We illustrate this approach by considering the constraints that can be placed on population models for galactic white dwarf binaries using a future space-based gravitational wave detector. We find that a mission that is able to resolve approximately 5000 of the shortest period binaries will be able to constrain the population model parameters, including the chirp mass distribution and a characteristic galaxy disk radius to within a few percent. This compares favorably to existing bounds, where electromagnetic observations of stars in the galaxy constrain disk radii to within 20%.

hierarchical↗

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗

Computational methods for estimation of parameters in hyperbolic systems

Approximation techniques for estimating spatially varying coefficients and unknown boundary parameters in second order hyperbolic systems are discussed. Methods for state approximation (cubic splines, tau-Legendre) and approximation of function space parameters (interpolatory splines) are outlined and numerical findings for use of the resulting schemes in model "one dimensional seismic inversion' problems are summarized.

Banks, H. T.↗

Aircraft ride quality controller design using new robust root clustering theory for linear uncertain systems

The aspect of controller design for improving the ride quality of aircraft in terms of damping ratio and natural frequency specifications on the short period dynamics is addressed. The controller is designed to be robust with respect to uncertainties in the real parameters of the control design model such as uncertainties in the dimensional stability derivatives, imperfections in actuator/sensor locations and possibly variations in flight conditions, etc. The design is based on a new robust root clustering theory developed by the author by extending the nominal root clustering theory of Gutman and Jury to perturbed matrices. The proposed methodology allows to get an explicit relationship between the parameters of the root clustering region and the uncertainty radius of the parameter space. The current literature available for robust stability becomes a special case of this unified theory. The bounds derived on the parameter perturbation for robust root clustering are then used in selecting the robust controller.

Yedavalli, R. K.↗

Gravitational-Wave Data Analysis with Spinning Merger-Ringdown Waveforms

The recent availability of high-quality, gravitational merger-ringdown waveforms from spinning black-hole systems has made possible the development of multi-mode GW templates for use in data-analysis studies of current and proposed interferometric GW detectors. We report on recent work at NASA Goddard, analyzing the most significant modes from aligned-spin black-hole-binary mergers. From these, we have developed time-domain merger-ringdown GW templates covering the aligned-spin portion of parameter space. We also discuss how using the full information content of aligned-spin mergers can significantly reduce uncertainties in some parameters, emphasizing the significant gains possible in the last stages of merger, inaccessible to inspiral-only post-Newtonian templates.

Kelly Bernard J.↗

Development of a Data Fusion Methodology for Lineload Aerodynamic Databases for a Launch Vehicle during Liftoff and Transition

The need for databases for the distributed loading on launch vehicles during the early portion of flight necessitates the use of expensive computational flows in regimes where wake effects dominate. While also being expensive, this is a regime that computational tools tend to historically have problems simulating accurately. To help tackle this problem, a method of data fusion to combine computational results to wind tunnel derived force and moment data is developed. Using this method, significant reduction in computational costs and increases in confidence of the final product is possible and has been used to generate several databases for the Space Launch System (SLS) at NASA. While the full details of database generation are not part of this work, the crucial method at its core is developed here. Two SLS geometries are used throughout the work to demonstrate the techniques. These are two of the larger geometries and represent both planned crewed missions to the Moon as well as potential cargo missions to deep space. The method uses principal component analysis (PCA) to generate a reduced ordered model (ROM) to help fill in the full parameter space. Other similar techniques are explored, but were not found to have a significant result on the predictions of the ROM. Because the full number of components are kept to generate the model, this lack of difference is expected. This method is then extended to ensure that predicted surfaces match trusted force and moment data derived from wind tunnel testing. This extension is done by setting up a constrained optimization problem in order to minimize the deviation from the surface resolved computational data while still integrating to the desired values. When generating the constrained optimization problem, a weighting factor to balance these competing needs is introduced. The work compares previously introduced weighting terms from similar work to the proposed terms and shows that the previously used terms do not have as desirable behavior in this flow regime. This method is then expanded by developing a technique to incorporate uncertainty quantification into the developed data fusion methodology. This expansion takes a two pronged approach. One examines transferring the uncertainties in the force and moment database and characterizes how those adjustments change the predicted lineloads. The second looks at model form error and looks how rebuilding the model using slightly different data changes the predictions. These two terms are then combined in order to create an uncertainty model that takes both effects into account. The limitations of the proposed methods is then discussed as well as possible techniques to address these shortcomings.

Launch Vehicles↗

CSD 2374310: Experimental Crystal Structure Determination

An entry from the Inorganic Crystal Structure Database, the world’s repository for inorganic crystal structures. The entry contains experimental data from a crystal diffraction study. The deposited dataset for this entry is freely available from the joint CCDC and FIZ Karlsruhe Access Structures service and typically includes 3D coordinates, cell parameters, space group, experimental conditions and quality measures.

Cell Parameters↗

Viscosity estimates for the crust and upper mantle from patterns of lacustrine shoreline deformation in the Eastern Great Basin

The deformed shorelines of Lake Bonneville constitute a classic source of information on lithospheric elastic thickness and upper mantle viscosity. We describe and apply a new model to a recently augmented data set. New data better constrain both the complex spatio-temporal pattern of the lake load and the crustal deformation response to that load. The history of lake level fluctuations has been significantly refined and somewhat modified. This is due to both more radiocarbon dates from within the Bonneville basin and to an improved calibration of the radiocarbon timescale itself. The data which constrain the crustal deformation pattern consist of ages and shoreline elevations from several hundred points which sample three major levels of Lake Bonneville and corresponding elevations from the high stands of three smaller lakes situated to the west of Lake Bonneville. The geometry of the Earth model incorporates an arbitrary number of layers overlying a half-space, and the rheology of each level can accommodate an arbitrary number of Maxwell viscoelastic elements in parallel. The inverse modeling comprises three complementary approaches: for the simplest configurations, we performed a direct search of the parameter space and delineated the irregular boundary of the subspace of acceptable models. For more complex configurations, we constrained the elastic parameters to their seismically determined values and then solved for viscosity versus depth profiles by either expressing the log(viscosity) versus log(depth) profile as a series of specially constructed orhtogonal polynomials, or by allowing each of 8-10 layers (plus the half-space) to have an independently determined viscosity. We found that the data do not strongly support (nor can they conclusively exclude) a more complex rheology than simple Maxwell viscoelasticity. The orthogonal polynomial solution exhibits an essentially monotonic decrease in viscosity with depth.

Bills, Bruce G.↗