Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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 1,477 records · Page 82

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

High-throughput combinatorial approach expedites the synthesis of a lead-free relaxor ferroelectric system

Developing novel lead-free ferroelectric materials is crucial for next-generation microelectronic technologies that are energy efficient and environment friendly. However, materials discovery and property optimization are typically time-consuming due to the limited throughput of traditional synthesis methods. In this work, we use a high-throughput combinatorial synthesis approach to fabricate lead-free ferroelectric superlattices and solid solutions of (Ba 0.7 Ca 0.3 )TiO 3 (BCT) and Ba(Zr 0.2 Ti 0.8 )O 3 (BZT) phases with continuous variation of composition and layer thickness. High-resolution x-ray diffraction (XRD) and analytical scanning transmission electron microscopy (STEM) demonstrate high film quality and well-controlled compositional gradients. Ferroelectric and dielectric property measurements identify the “optimal property point” achieved at the composition of 48BZT–52BCT. Displacement vector maps reveal that ferroelectric domain sizes are tunable by varying {BCT–BZT} N superlattice geometry. This high-throughput synthesis approach can be applied to many other material systems to expedite new materials discovery and properties optimization, allowing for the exploration of a large area of phase space within a single growth.

36 MATERIALS SCIENCE↗

New Global Electron Density Observations from GPS-RO in the D- and E-Region Ionosphere

A novel retrieval technique is developed for electron density (N(sub e)) in the D- and E-region (80-120 km) using the high-quality 50-Hz GPS radio occultation (GPS-RO) phase measurements. The new algorithm assumes a slow, linear variation in the F-region background when the GPS-RO passes through the D- and E-region, and extracts the N(sub e) profiles at 80-130 km from the phase advance signal caused by N(sub e). Unlike the conventional Abel function, the new approach produces a sharp N(sub e) weighting function in the lower ionosphere, and the N(sub e) retrievals are in good agreement with the IRI (International Reference Ionosphere) model in terms of monthly maps, zonal means and diurnal variations. The daytime GPS-RO N(sub e) profiles can be well characterized by the alpha-Chapman function of three parameters (N(sub mE), h(sub mE) and H), showing that the bottom of E-region is deepening and sharpening towards the summer pole. At high latitudes the monthly GPS-RO N(sub e) maps at 80-120 km reveal clear enhancement in the auroral zones, more prominent at night, as a result of energetic electron precipitation (EEP) from the outer radiation belt. The D-/E-region auroral N(sub e) is strongly correlated with K(sub p) on a daily basis. The new N(sub e) data allow further comprehensive analyses of the sporadic E (E(sub s)) phenomena in connection with the background N(sub e) in the E-region. The layered (2-10 km) and fluctuated (less than 2 km) E(sub s) components, namely N(sub e_Layer) than N(sub e_Pert), are extracted with respect to the background N( sub e_Region) on a profile-by-profile basis. The N(sub e_Layer) component has a strong but highly-refined peak at approximately 105 km, with an amplitude smaller than N(sub e_Region) approximately by an order of magnitude. The N(sub e_Pert) component, which was studied extensively in the past, is approximately 2 orders of magnitude weaker than N(sub e_Layer). Both N(sub e_Layer) and N(sub e_Pert) are subject to significant diurnal and semidiurnal variations, showing downward progression with local time in amplitude. The 11-year solar cycle dominates the N(sub e) interannual variations, showing larger N(sub e_Region) and N(sub e_Layer) but smaller N(sub e_Pert) amplitudes in the solar maximum years. Enhanced Ne profiles are often observed in the polar winter, showing good correlation with solar proton events (SPEs) and geomagnetic activity. The new methodology offers great potential for retrieving low N(sub e) in the D-region, where radio propagation and communication blackouts can occur due to enhanced ionization. For space weather applications it is recommended for GPSRO operations to raise the top of high-rate data acquisition to approximately 140 km in the future.

GSP-R↗

Computing Highly Accurate Spectroscopic Line Lists for Characterization of Planetary Atmospheres: CO2 and SO2 Line Lists Needed for Modeling Venus

Over the last decade, it has become apparent that the most effective approach for determining highly accurate rotational and rovibrational line lists for molecules of interest in planetary atmospheres and other astrophysical environments is through a combination of highresolution laboratory experiments coupled with state-of-the art ab initio quantum chemistry methods. The approach involves computing the most accurate potential energy surface (PES) possible using state-of-the art electronic structure methods, followed by computing rotational and rovibrational energy levels using an exact variational method to solve the nuclear Schrödinger equation. Then, reliable experimental data from high-resolution experiments is used to refine the ab initio PES in order to improve the accuracy of the computed energy levels and transition energies. From the refinement step, we have been able to achieve an accuracy of approximately 0.015 cm-1 for rovibrational transition energies, and even better for purely rotational transitions. This combined "experiment / theory" approach allows for determination of essentially a complete line list, with hundreds of millions of transitions, and having the transition energies and intensities be highly accurate. Our group has successfully applied this approach to determine highly accurate line lists for NH3, CO2 and isotopologues, and SO2 and isotopologues. Here I will report our latest results for CO2 and SO2 including all isotopologues. Comparisons to the available data in HITRAN2012 and other available databases will be shown, though we note that our line lists for SO2 are significantly more complete than any other databases. Since it is important to span a large temperature range in order to model the spectral signature of Venus as well as exoplanets, we will demonstrate how the spectra change on going from low temperatures (100 K) to higher temperatures (500 K to 1500 K).

potential energy surface (PES)5405 PLANETARY SCIEN↗

Multi Model Monte Carlo with Python (MXMCPy)

Multi Model Monte Carlo with Python (\mxmc {}) is a software package developed as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Motivated by uncertainty propagation problems where classical Monte Carlo (MC) simulation is computationally intractable, various multi-model MC approaches have recently emerged that yield unbiased estimators with significantly reduced variance relative to MC for the same cost. These existing methods include multi-level Monte Carlo (MLMC), multi-fidelity Monte Carlo (MFMC), and approximate control variates (ACV). Given a fixed computational budget and a collection of models with varying cost/accuracy, each method seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. \mxmc {} is a versatile tool that enables convenient access to many existing multi-model MC approaches within one modular and extensible package. With \mxmc {}, users can easily compare existing methods to determine the best choice for their particular problem, while developers have a basis for implementing and sharing new variance reduction approaches. This report introduces the \mxmc {} software, providing a summary of the problem-solving workflow for users as well as a brief overview of the code layout for developers.

Geoffrey F Bomarito↗

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗

Topology control of human fibroblast cells monolayer by liquid crystal elastomer

Eukaryotic cells in living tissues form dynamic patterns with spatially varying orientational order that affects important physiological processes such as apoptosis and cell migration. The challenge is how to impart a predesigned map of orientational order onto a growing tissue. Here, we demonstrate an approach to produce cell monolayers of human dermal fibroblasts with predesigned orientational patterns and topological defects using a photoaligned liquid crystal elastomer (LCE) that swells anisotropically in an aqueous medium. The patterns inscribed into the LCE are replicated by the tissue monolayer and cause a strong spatial variation of cells phenotype, their surface density, and number density fluctuations. Unbinding dynamics of defect pairs intrinsic to active matter is suppressed by anisotropic surface anchoring allowing the estimation of the elastic characteristics of the tissues. The demonstrated patterned LCE approach has potential to control the collective behavior of cells in living tissues, cell differentiation, and tissue morphogenesis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Self‐Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces

The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as “Jacob's ladder.” This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBEα framework with mathematical equation being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A topological approach to computer-aided sensitivity analysis

Sensitivities of any arbitrary system are calculated using general purpose digital computer with available software packages for transfer function analysis. Sensitivity shows how element variation within system affects system performance. Signal flow graph illustrates topological system behavior and relationship among parameters in system.

Chan, S. P.↗

Assimilation of Geosat Altimetric Data in a Nonlinear Shallow-Water Model of the Indian Ocean by Adjoint Approach. Part II: Some Validation and Interpretation of the Assimilated Results: Some Validation and Interpretation of the Assimilated Results - Part 2

This paper examines the results of assimilating Geosat sea level variations relative to the November 1986-November 1988 mean reference, in a nonlinear reduced-gravity model of the Indian Ocean, Data have been assimilated during one year starting in November 1986 with the objective of optimizing the initial conditions and the yearly averaged reference surface. The thermocline slope simulated by the model with or without assimilation is validated by comparison with the signal, which can be derived from expandable bathythermograph measurements performed in the Indian Ocean at that time. The topography simulated with assimilation on November 1986 is in very good agreement with the hydrographic data. The slopes corresponding to the South Equatorial Current and to the South Equatorial Countercurrent are better reproduced with assimilation than without during the first nine months. The whole circulation of the cyclonic gyre south of the equator is then strongly intensified by assimilation. Another assimilation experiment is run over the following year starting in November 1987. The difference between the two yearly mean surfaces simulated with assimilation is in excellent agreement with Geosat. In the southeastern Indian Ocean, the correction to the yearly mean dynamic topography due to assimilation over the second year is negatively correlated to the one the year before. This correction is also in agreement with hydrographic data. It is likely that the signal corrected by assimilation is not only due to wind error, because simulations driven by various wind forcings present the same features over the two years. Model simulations run with a prescribed throughflow transport anomaly indicate that assimilation is rather correcting in the interior of the model domain for inadequate boundary conditions with the Pacific.

Greiner, Eric↗

Use of Smoothed Measured Winds to Predict and Assess Launch Environments

This paper describes an approach which was used to predict Ares I-X launch vehicle structural airloads and controllability prior to launch. This prediction along with the proper application of protection for variations in key environment and trajectory parameters was used to predict the potential for the launch trajectory to stay within the integrated vehicle certification envelopes. This data was applied by the launch team as a key element of the launch day go/no-go recommendation. NASA s Space Shuttle and the heavy lift versions of the Evolved Expendable Launch Vehicles assessments require the use of a high fidelity wind measurement which limits the number of wind measurement sources which can be used for a pre-launch go/no-go assessment. The approach described in this paper has the potential to allow the pre-launch assessment team to use wind measurements from a number of additional sources which will enhance launch availability and also has the potential to improve mission assurance and enhance safety

Cordova, Henry S.↗

Usefulness of heart measures in flight simulation

The results of three studies performed at the NASA Langley Research Center are presented to indicate the areas in which heart measures are useful for detecting differences in the workload state of subjects. Tasks that involve the arousal of the sympathetic nervous system, such as landing approaches, were excellent candidates for the use of average heart-rate and/or the increase in heart-rate during a task. The latter of these two measures was the better parameter because it removed the effects of diurnal variations in heart-rate and some of the intersubject variability. Tasks which differ in the amount of mental resources required are excellent candidates for heart-rate variability measures. Heart-rate variability measures based upon power spectral density techniques were responsive to the changing task demands of landing approach tasks, approach guidance options, and 2 versus 20 second interstimulus-intervals of a monitoring task. Heart-rate variability measures were especially sensitive to time-on-task when the task was characterized by minimal novelty, complexity, and uncertainty (i.e., heart-rate variability increases as a function of the subjects boredom).

Harris, Randall L., Sr.↗

Extending a 10‐node composite tetrahedral finite element for solid mechanics

Abstract We propose to extend the composite tetrahedral finite element first introduced by Thoutireddy et al. (2002) and recently reformulated by Ostien et al. (2016). We generalize the gradient operator and mass matrix to curved domains through analytical expressions weighted by subtetrahedra Jacobians. Optimal integration weights are constructed to increase the accuracy of Gaussian quadrature. We preserve the variational structure of the formulation through a new five‐field functional with additional, independent fields for the Jacobian and the pressure. This approach not only obviates volumetric locking but also yields symmetry. A deleterious soft mode, common to both the quadratic and composite tetrahedral element with constant pressure formulations is effectively stabilized through a novel convex energy penalty function. Numerous numerical examples spanning a patch test to the impact of a Taylor bar demonstrate the accuracy, robustness, and convergence of the extended composite tetrahedral element for application to structural metals.

Foulk III, James W.↗

Quantifying sensitivity in numerical weather prediction-modeled offshore wind speeds through an ensemble modeling approach

A decade of research has shown that numerical weather prediction (NWP)-modeled wind speeds can be highly sensitive to the inputs and setups within the NWP model. For wind resource characterization applications, this sensitivity is often addressed by constructing a range of setups and selecting the one that best validates against observations. However, this approach is not possible in areas that lack high-quality hub height observations, especially offshore wind areas. In such cases, techniques to quantify and disseminate confidence in NWP-modeled wind speeds in the absence of observations are needed. We address this need in the present study and propose best practices for quantifying the spread in NWP-modeled wind speeds. We implement an ensemble approach in which we consider 24 different setups to the Weather Research and Forecasting (WRF) model. We construct the ensemble by considering variations in WRF version, WRF namelist, atmospheric forcing, and sea surface temperature (SST) forcing. Our analysis finds that the standard deviation produces more consistent estimates compared to the interquartile range and tends to be the more conservative estimator for ensemble variability. We further find that model spread increases closer to the surface and on shorter time scales. In conclusion, we explore methods to attribute total ensemble variability to the different ensemble components (e.g., atmospheric forcing and SST product) and find that contributions by components also vary depending on time scale. We anticipate that the methods and results presented in this paper will provide a reasonable foundation for further research into ensemble-based wind resource data sets.

17 WIND ENERGY↗

Three-body resonances in the φ 4 theory

We study the properties of three-body resonances using a lattice complex scalar φ4 theory with two scalars, with parameters chosen such that one heavy particle can decay into three light ones. We determine the two- and three-body spectra for several lattice volumes using variational techniques, and then analyze them with two versions of the three-particle finite-volume formalism: the Relativistic Field Theory approach and the Finite-Volume Unitarity approach. We find that both methods provide an equivalent description of the energy levels, and we are able to fit the spectra using simple parametrizations of the scattering quantities. By solving the integral equations of the corresponding three-particle formalisms, we determine the pole position of the resonance in the complex energy plane and thereby its mass and width. We find very good agreement between the two methods at different values of the coupling of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Interactions of Wind Energy Project Siting, Wind Resource Potential, and the Evolution of the U.S. Power System

The interactions of wind energy project siting, resource potential, and power system evolution are critical to understand given interest in high renewable energy systems, particularly as wind costs have fallen and grid operators implement solutions to manage variable generation. This study applies an integrative approach that combines spatially-explicit resource potential analysis with high spatial resolution U.S. electricity system modeling. Multiple wind supply curves, representing variations in siting regimes that account for interactions with built infrastructure, regulatory, physical, and social land use factors, are evaluated to determine how siting could influence prospective wind development. Different siting regimes lead to variations in future wind capacity, with the greatest impacts observed under scenarios with high demand for clean energy. With the tightest emissions limit modeled, 2050 onshore wind capacity varied by +7% (49 GW) in the least-restrictive siting regime to a decrease of 37% (270 GW) under the most-constrained case. More-stringent siting restrictions lead to higher electricity prices and emissions that should be weighed with local impacts. Under restrictive siting regimes, wind deployment is also sensitive to transmission availability and wind plant design. Overall, the findings highlight the importance of local land use considerations in regional and national power system planning.

17 WIND ENERGY↗

Strong Relativistic Effects in Lanthanide-Based Single-Molecule Magnets

Lanthanide-based single-molecule magnets (SMMs) are promising building blocks for quantum memory and spintronic devices. Designing lanthanide-based SMMs with long spin relaxation time requires a detailed understanding of their electronic structure, including the crucial role of the spin-orbit coupling (SOC). While traditional calculations of SOC using the perturbation theory applied to a solution of the non-relativistic Schrödinger equation are valid for light atoms, this approach is questionable for systems containing heavy elements such as lanthanides. We investigate the accuracy of the perturbation estimates of SOC by variationally solving the Dirac equation for the [DyO] + molecule, a prototype of a lanthanide-based SMM. We show that the energy splittings between the M J states involved in spin relaxation depend on the interplay between strong SOC and dynamic electron correlations. Here, we demonstrate that this interplay affects the resonances between the spin and vibrational transitions, and, therefore, the spin relaxation time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying the Value of Technology and Policy Innovation in Water Resource Portfolios

Abstract Recent water infrastructure planning models demonstrate that explicitly accounting for hydrological changes in long‐term water planning reduces costs and increases system robustness. However, current models do not consider the effects of other changes occurring in the system over long time horizons, namely technology and policy innovation. By leveraging state‐of‐the‐art dynamic control techniques, we propose a suite of tools that explicitly relate hydrologic change, technology innovation, and policy interventions to system‐level water costs. We design robust and cost‐optimal strategies for water supply expansion under hydrological variability and hundreds of plausible innovation states, for example, treatment cost and diversification improvements, deployment times, and demand reduction campaigns. Our approach circumvents the need to pre‐assign probabilities to innovation scenarios and, instead, directly calculates the effects of technology, policy, and hydrology variations on system cost and planning strategies to identify avenues of impactful innovation. These tools can support technology development hubs, urban, state, and federal water planners in anticipating cost implications of technology and policy innovation at the local level, and prioritizing high‐impact innovation goals based on quantitative targets. Results for the case study of Santa Barbara, California identify high‐value innovation attributes for the city as well as combinations of innovation attributes across technologies and policies, and quantify their potential utility cost reductions.

13 HYDRO ENERGY↗