Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 271 records · Page 15

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

Modeling Variability of Plasma Conditions in the Io Torus

Telescopic observations an in situ measurements of the Io plasma torus show the density, temperature and composition to vary over time, sometimes up to a factor of 2. While previous models of the physical and chemical processes in the Io plasma torus have reasonably modeled the conditions of the Voyager 1 era, their authors have not addressed the observed variability nor explored the sensitivity of torus conditions to input parameters. In this paper we present a homogeneous torus model parameterized by five variables (transport timescale, neutral source strength, ratio of oxygen sulfur to atoms in the source, fraction of superthermal electrons, temperature of these hot electrons). The model incorporates the most recent data for ionization, recombination, charge exchange and radiative energy losses for the major torus species (S, S(sup +), S(sup ++), S(sup +++), O, O(sup +), O(sup ++). We solve equations of conservation of mass and energy to find equilibrium conditions for a set of input parameters. We compare model plasma conditions with those observed by Voyager 1 Voyager 2, and Cassini. Furthermore, we explore the sensitivity of torus conditions to each parameter. We find that (1) torus conditions are distinctly different for the Voyager 1, Voyager 2 and Cassini eras, (2) unique torus input parameters for any given era are poorly constrained given the wide range of solution space that is consistent with the range of observed torus conditions, (3) ion composition is highly sensitive to the specification of a non-thermal electron distribution, (4) neutral O/S source ratio is highly variable with model values ranging between 1.7 for Cassini to 4.0 for Voyager 1 conditions, (5) transport times range between 23 days for Voyager 2 to 50 days for Voyager 1 and Cassini, (6) neutral source strengths range between 7 to 30 x 10(sup -4) cm (sup -3) s(sup -1) which corresponds to a net production of 0.4 to 1.3 tons/s for a torus volume of 1.4 x 10(sup 31) cm(sup 3), or 38 R(sub j)(sup 3).

Delamere, P. A.↗

LES ARM Symbiotic Simulation and Observation (LASSO) Implementation Strategy

This document illustrates the design of the Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) workflow to provide a routine, high-resolution modeling capability to augment the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility’s high-density observations. LASSO will create a powerful new capability for furthering ARM’s mission to advance understanding of cloud, radiation, aerosol, and land-surface processes. The combined observational and modeling elements will enable a new level of scientific inquiry by connecting processes and context to observations and providing needed statistics for details that cannot be measured. The result will be improved process understanding that facilitates concomitant improvements in climate model parameterizations. The initial LASSO implementation will be for ARM’s Southern Great Plains site in Oklahoma and will focus on shallow convection, which is poorly simulated by climate models due in part to clouds’ typically small spatial scale compared to model grid spacing, and because the convection involves complicated interactions of microphysical and boundary layer processes.

Gustafson Jr., WI↗

Recent Progresses in Incorporating Human Land-Water Management into Global Land Surface Models Toward Their Integration into Earth System Models

The global water cycle has been profoundly affected by human land-water management. As the changes in the water cycle on land can affect the functioning of a wide range of biophysical and biogeochemical processes of the Earth system, it is essential to represent human land-water management in Earth system models (ESMs). During the recent past, noteworthy progress has been made in large-scale modeling of human impacts on the water cycle but sufficient advancements have not yet been made in integrating the newly developed schemes into ESMs. This study reviews the progresses made in incorporating human factors in large-scale hydrological models and their integration into ESMs. The study focuses primarily on the recent advancements and existing challenges in incorporating human impacts in global land surface models (LSMs) as a way forward to the development of ESMs with humans as integral components, but a brief review of global hydrological models (GHMs) is also provided. The study begins with the general overview of human impacts on the water cycle. Then, the algorithms currently employed to represent irrigation, reservoir operation, and groundwater pumping are discussed. Next, methodological deficiencies in current modeling approaches and existing challenges are identified. Furthermore, light is shed on the sources of uncertainties associated with model parameterizations, grid resolution, and datasets used for forcing and validation. Finally, representing human land-water management in LSMs is highlighted as an important research direction toward developing integrated models using ESM frameworks for the holistic study of human-water interactions within the Earths system.

biogeochemistry↗

An evaluation of the E3SMv1 Arctic ocean and sea-ice regionally refined model

The Energy Exascale Earth System Model (E3SM) is a state-of-the-science Earth system model (ESM) with the ability to focus horizontal resolution of its multiple components in specific areas. Regionally refined global ESMs are motivated by the need to explicitly resolve, rather than parameterize, relevant physics within the regions of refined resolution, while offering significant computational cost savings relative to the respective cost of configurations with high-resolution (HR) everywhere on the globe. In this paper, we document results from the first Arctic regionally refined E3SM configuration for the ocean and sea-ice components (E3SM-Arctic-OSI), while employing data-based atmosphere, land, and hydrology components. Our aim is an improved representation of the Arctic coupled ocean and sea-ice state, its variability and trends, and the exchanges of mass and property fluxes between the Arctic and the sub-Arctic. We find that E3SM-Arctic-OSI increases the realism of simulated Arctic ocean and sea-ice conditions compared to a similar low-resolution E3SM simulation without the Arctic regional refinement in ocean and sea-ice components (E3SM-LR-OSI). In particular, exchanges through the main Arctic gateways are greatly improved with respect to E3SM-LR-OSI. Other aspects, such as the Arctic freshwater content variability and sea-ice trends, are also satisfactorily simulated. Yet, other features, such as the upper-ocean stratification and the sea-ice thickness distribution, need further improvements, involving either more advanced parameterizations, model tuning, or additional grid refinements. Overall, E3SM-Arctic-OSI offers an improved representation of the Arctic system relative to E3SM-LR-OSI, at a fraction (15 %) of the computational cost of comparable global high-resolution configurations, while permitting exchanges with the lower-latitude oceans that cannot be directly accounted for in Arctic regional models.

54 ENVIRONMENTAL SCIENCES↗

Venus lower thermosphere studies

Studies undertaken in this project have sought to understand lower thermospheric structure and dynamics (less than or equal to 145 km), particularly the processes responsible. This is a region just below the reach of in-situ instruments onboard PVO (Pioneer Venus Orbiter) during the first few diurnal cycles. PVO remote airglow observations (nitric oxide, O2, visible, O 1304A) have been coupled with ground-based observations (CO densities, winds, temperatures, O2 IR nightglow) to address the behavior of lower thermospheric winds and chemistry over 95 to 150 km. This interpretation of PVO and related data is accomplished by using the NCAR Venus thermospheric general circulation model (VTGCM) (Bougher et al., 1988; 1990). This model has been modified over the last two years to improve its ability to calculate O, CO, and O2 densities, temperatures, nightglow, and subsolar-to-antisolar and zonal winds over 95 to 150 km (Bougher and Borucki, 1992). Our VTGCM studies show that: (1) O2 visible and IR nightglow distributions can be used to trace lower thermosphere / upper mesosphere winds over 100-130 km. Typically, weak zonal winds (less than or equal to 25 m/sec) and nightglow maximum patches near 0100 LT prevail. Occasionally, strong zonal winds (30-60 m/sec) and airglow patches peaking near 0300 LT characterize the Venus lower thermosphere. (2) It is clear that the dynamics of the Venus 90-130 km region is highly variable on time scales as short as an hour. This is most likely due to the time variable nature of upward propagating gravity waves, which grow in amplitude and eventually break. The resulting turbulence gives rise to local time variable eddy diffusion and momentum drag, both of which strongly impact global density and nightglow distributions. (3) The oxygen chemistry (O, O2, etc.) over 90-120 km is strongly dependent on HO(x) and CLO(x) tracer species that must be properly included in any coupled chemical dynamical model. (4) The density profiles of light species (O, CO, N, He) are strongly affected by large-scale transport by the winds. Strong eddy diffusion is not a suitable model parameterization for approximating these light species, especially for extrapolation into the region below 140 km where PVO in-situ data is lacking. Instead, the fully coupled chemical dynamical VTGCM model should be used to improve estimates of densities within the VTS3 empirical model (Hedin et al., 1983) below 140 km.

Bougher, Stephen W.↗

Cloud-Based Model Calibration Using OpenStudio: Preprint

OpenStudio is a free, open source Software Development Kit (SDK) and application suite for performing building energy modeling and analysis. The OpenStudio Parametric Analysis Tool has been extended to allow cloud-based simulation of multiple OpenStudio models parametrically related to a baseline model. This paper describes the new cloud-based simulation functionality and presents a model cali-bration case study. Calibration is initiated by entering actual monthly utility bill data into the baseline model. Multiple parameters are then varied over multiple iterations to reduce the difference between actual energy consumption and model simulation results, as calculated and visualized by billing period and by fuel type. Simulations are per-formed in parallel using the Amazon Elastic Cloud service. This paper highlights model parameterizations (measures) used for calibration, but the same multi-nodal computing architecture is available for other purposes, for example, recommending combinations of retrofit energy saving measures using the calibrated model as the new baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Anthropogenic enhancements to production of highly oxygenated molecules from autoxidation

As particulate matter abatement strategies reduce anthropogenic nitrogen oxide (NOx) emissions, the chemical regime of an urban atmosphere may change to one in which autoxidation becomes increasingly important. Developed countries such as the United States have already experienced large declines in nitrogen oxides because of controls on emissions, and the possibility for increased aerosol exists as autoxidation results in highly oxidized molecules (HOM) which efficiently form particulate matter. Using a-pinene reaction with hydroxyl radicals as an archetype system, a molecular representation of organic aerosol was created and evaluated against laboratory experiments. The applica-tion of parameters constrained by these experiments to the ambient atmosphere indicated autoxidation is becoming more competitive in urban areas of the United States while at the same time HOM production rates are declining due to decreases in oxidant abundance. For locations such as Guangzhou, China, where emissions of NOx generally increased from the 1990s to late 2000s, HOM production was also predicted to increase. These results demonstrate that anthro-pogenic emissions can enhance HOM and the resulting biogenic organic aerosol formation through oxidant abundance. Ambient particulate matter is responsible for a significant fraction of the global burden of disease (Lim et al., 2012) and affects climate (Myhre et al., 2013). The organic portion of ambient aerosol largely forms due to the condensation of low vapor pressure (Jimenez et al., 2009) or highly soluble (Saxena and Hildemann, 1996) gases in the atmosphere . Oxidation of monoterpenes from vegetation, results in significant formation of particle mass (Pandis et al., 1991), and is likely the dominant source of organic aerosol in the southeastern United States (Zhang et al., 2018; Xu et al., 2018). While monoterpenes are included as a secondary organic aerosol (SOA) source in most chemical transport models, many current regional to global scale model parameterizations (Pye et al., 2010, 2017; Koo et al., 2014) lack a mechanistic dependence of monoterpene SOA on NOx and oxidant identity (e.g. OH vs. ozone). This limits the ability of models to predict how changes in emissions will affect ambient concentrations and new particle formation events. Building a mechanistic understanding of formation pathways is critical for being able to determine historic drivers of ambient particulate matter as well as inform predictions of future levels

Pye, Havala O.↗

Sandia Unstructured Triangle Table Generator

The software generates data tables for thermodynamic and transport properties of materials as described by a set of input models. For each input model parameterization, an associated table is created on an unstructured triangular grid. These grids all conform to the same topology. A statistical accuracy guarantee is provided for the tabular representation of the model. Details of the model and table specification are given in a XML input deck.

Carpenter, JohnH.↗

Application of aerosol-ice nucleating particle closure to establish the leading parameters governing ice crystal number concentration under commonly observed mixed-phase cloud conditions

This is the final technical report for DOE DE-SC0021034 award entitled “Application of aerosol-ice nucleating particle closure to establish the leading parameters governing ice crystal number concentration under commonly observed mixed-phase cloud conditions”. The overarching objective of this project was to gain a predictive understanding of the processes that lead from aerosols to ice-nucleating particles (INPs) to ice crystal number concentrations (ICNCs) in stratiform mixed-phase clouds where supercooled water droplets and ice crystals coexist, and to inform climate model parameterizations. During this four-year project we addressed the three main research goals of this objective: 1) Study the aerosol-INP link; 2) Develop an aerosol-cloud column model (AC-1D); and 3) Study the INP-ICNC link. Those activities significantly advanced our predictive understanding of how to represent INPs in could models and its implications on ice crystal number concentrations. The results have been published in five peer-reviewed publications with three more in review and preparation. It allowed for collaboration resulting in additional fourteen peer-reviewed publications that are concerned with the objective of this research project. This research has been communicated to the scientific community at national and international meetings and conferences by numerous talks, seminars, and posters including fifteen invited presentations. This project allowed training of four graduate students and one postdoctoral research associate resulting in one M.S. thesis and advancement to Ph.D. candidacy.

54 ENVIRONMENTAL SCIENCES↗

Climate Predicts Geographic and Temporal Variation in Mosquito-Borne Disease Dynamics on Two Continents

Climate drives population dynamics through multiple mechanisms, which can lead to seemingly context-dependent effects of climate on natural populations. For climate-sensitive diseases such as dengue, chikungunya, and Zika, climate appears to have opposing effects in different contexts. Here we show that a model, parameterized with laboratory measured climate-driven mosquito physiology, captures three key epidemic characteristics across ecologically and culturally distinct settings in Ecuador and Kenya: the number, timing, and duration of outbreaks. The model generates a range of disease dynamics consistent with observed Aedes aegypti abundances and laboratory-confirmed arboviral incidence with variable accuracy (28–85% for vectors, 44–88%for incidence). The model predicted vector dynamics better in sites with a smaller proportion of young children in the population, lower mean temperature, and homes with piped water and made of cement. Models with limited calibration that robustly capture climate-virus relationships can help guide intervention efforts and climate change disease projections.

Climate-change ecology↗

Incorporating Valve Switching Losses Into a Static Optimal Control Algorithm for the Hybrid Hydraulic-Electric Architecture (HHEA)

The Hybrid Hydraulic-Electric Architecture (HHEA) has, in recent years, been proposed as an energy efficient alternative to conventional load-sensing architectures in mobile machines such as excavators and wheel-loaders. HHEA leverages the advantages of hydraulic power and electric power to eliminate throttling valves while also improving the energy and control performance of the system. The architecture utilizes a set of common pressure rails to provide a majority of power and and a small electric motor driven pump to modulate this power to meet the exact demand. Previous work has developed a computationally efficient Lagrange Multiplier approach for determining the optimal pressure rail selections that minimizes the energy losses in the system. The static model used considers only the energy use for each pressure rail selection but not the losses associated with the valves during the transition. This paper presents an approach to include the switching losses in the model and in the optimization procedure. To capture the switching losses, switching events between different rails and at various input and output flow rates were simulated with consideration of valve spool dynamics. A parameterized model that summarizes the losses is then obtained, allowing switching losses to be added to the previous energy analysis. The performance of the switching loss model was compared with reference data obtained from a high-fidelity simulation model. To incorporate the switching losses into optimal control algorithm, an efficient dynamic programming approach that prevents frequent switching is adopted in place of the Lagrange multiplier approach. The overall effect of switching losses on the energy consumption and optimal control decisions is presented. In general, switching losses contribute to about 9–10% of input energy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Understanding the Effect of Aerosol Properties on Cloud Droplet Formation during TCAP Field Campaign Report

The formation of clouds is an essential element in understanding the Earth’s radiative budget. Liquid water clouds form when the relative humidity exceeds saturation and condensedphase water nucleates on atmospheric particulate matter. The effect of aerosol properties such as size, morphology, and composition on cloud droplet formation has been studied theoretically as well as in the laboratory and field. Almost without exception these studies have been limited to parallel measurements of aerosol properties and cloud formation or collection of material after the cloud has formed, at which point nucleation information has been lost. Studies of this sort are adequate when a large fraction of the aerosol activates, but correlations and resulting model parameterizations are much more uncertain at lower supersaturations and activated fractions.

54 ENVIRONMENTAL SCIENCES↗

Polynomial chaos expansions for dependent random variables

Polynomial chaos expansions (PCE) are well-suited to quantifying uncertainty in models parameterized by independent random variables. The assumption of independence leads to simple strategies for building multivariate orthonormal bases and for sampling strategies to evaluate PCE coefficients. In contrast, the application of PCE to models of dependent variables is much more challenging. Three approaches can be used to construct PCE of models of dependent variables. The first approach uses mapping methods where measure transformations, such as the Nataf and Rosenblatt transformation, can be used to map dependent random variables to independent ones; however we show that this can significantly degrade performance since the Jacobian of the map must be approximated. A second strategy is the class of dominating support methods. In these approaches a PCE is built using independent random variables whose distributional support dominates the support of the true dependent joint density; we provide evidence that this approach appears to produce approximations with suboptimal accuracy. A third approach, the novel method proposed here, uses Gram-Schmidt orthogonalization (GSO) to numerically compute orthonormal polynomials for the dependent random variables. This approach has been used successfully when solving differential equations using the intrusive stochastic Galerkin method, and in this paper we use GSO to build PCE using a non-intrusive stochastic collocation method. The stochastic collocation method treats the model as a black box and builds approximations of the input-output map from a set of samples. Building PCE from samples can introduce ill-conditioning which does not plague stochastic Galerkin methods. To mitigate this ill-conditioning we generate weighted Leja sequences, which are nested sample sets, to build accurate polynomial interpolants. We show that our proposed approach, GSO with weighted Leja sequences, produces PCE which are orders of magnitude more accurate than PCE constructed using mapping or dominating support methods.

97 MATHEMATICS AND COMPUTING↗

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗

Development and Validation of an Initial Electric Motor Rotor Vibration Model

An overview of the development of a low fidelity prediction tool for electric motor rotor resonant modes and frequencies. These predictions are useful for the prediction of UAM motor noise during the design phase. Finite Element Analysis of motor rotors has led to the development of a parameterized model to predict these modes shapes. This has in turn allowed to a curve-fit based model from parametric sweeps. These curves can predict resonant frequencies and shapes for 1-5kW class motors without running a simulation.

acoustics↗

A fully Bayesian method for jointly fitting instrumental calibration and X-ray spectral models

Owing to a lack of robust principled methods, systematic instrumental uncertainties have generally been ignored in astrophysical data analysis despite wide recognition of the importance of including them. Ignoring calibration uncertainty can cause bias in the estimation of source model parameters and can lead to underestimation of the variance of these estimates. We previously introduced a pragmatic Bayesian method to address this problem. The method is 'pragmatic' in that it introduced an ad hoc technique that simplified computation by neglecting the potential information in the data for narrowing the uncertainty for the calibration product. Following that work, we use a principal component analysis to efficiently represent the uncertainty of the effective area of an X-ray (or γ-ray) telescope. Here, however, we leverage this representation to enable a principled, fully Bayesian method that coherently accounts for the calibration uncertainty in high-energy spectral analysis. In this setting, the method is compared with standard analysis techniques and the pragmatic Bayesian method. The advantage of the fully Bayesian method is that it allows the data to provide information not only for estimation of the source parameters but also for the calibration product—here the effective area, conditional on the adopted spectral model. In this way, it can yield more accurate and efficient estimates of the source parameters along with valid estimates of their uncertainty. Provided that the source spectrum can be accurately described by a parameterized model, this method allows rigorous inference about the effective area by quantifying which possible curves are most consistent with the data.

79 ASTRONOMY AND ASTROPHYSICS↗

Global Nonlinear Parametric Modeling with Application to F-16 Aerodynamics

A global nonlinear parametric modeling technique is described and demonstrated. The technique uses multivariate orthogonal modeling functions generated from the data to determine nonlinear model structure, then expands each retained modeling function into an ordinary multivariate polynomial. The final model form is a finite multivariate power series expansion for the dependent variable in terms of the independent variables. Partial derivatives of the identified models can be used to assemble globally valid linear parameter varying models. The technique is demonstrated by identifying global nonlinear parametric models for nondimensional aerodynamic force and moment coefficients from a subsonic wind tunnel database for the F-16 fighter aircraft. Results show less than 10% difference between wind tunnel aerodynamic data and the nonlinear parameterized model for a simulated doublet maneuver at moderate angle of attack. Analysis indicated that the global nonlinear parametric models adequately captured the multivariate nonlinear aerodynamic functional dependence.

Morelli, Eugene A.↗