Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

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,387 records · Page 77

Machine-Learning Assisted Identification of Battery Life Models

Predictive battery life models are commonly utilized to extrapolate degradation trends observed during accelerated aging tests for simulation of degradation in real-world applications. Thus, fitting accelerated aging data as accurately as possible and with low uncertainty is crucial for making believable projections of battery lifetime, but it is challenging to identify algebraic expressions that accurately fit multivariate degradation trends. A review of models published in literature reveal some common expressions for fitting calendar aging data, which is only dependent on temperature and state-of-charge, but no consistency across many models for fitting cycle aging data, indicating the need for a statistically rigorous data driven approach for developing empirical models. This talk will describe a machine-learning assisted method for identification of predictive battery life models utilizing bilevel optimization and symbolic regression. Bilevel optimization with cross-validation is used to statistically determine cell- and stress-dependent model parameters, while symbolic regression identifies both linear and multiplicative candidate expressions to predict stress-dependent degradation rates by selecting low-order subsets of features from a generated feature library. Because model expressions are identified empirically, it is crucial to ensure resulting models behave according to physical expectations, so the stability of models for interpolation or extrapolation is interrogated qualitatively through simulation and quantitatively through cross-validation and uncertainty quantification via bootstrap resampling. This model identification approach substantially improves upon models identified purely using expert judgement in terms of both accuracy and uncertainty. Model simulation and validation is then conducted by deriving a state-equation form of the predictive model, enabling simulation of battery aging under dynamic stresses. This enables validation of the predictive battery model on lab-based tests with varying conditions or on drive-cycle or application-cycle testing protocols. Parameter uncertainty can be carried forward into model simulation, giving lifetime estimates and confidence windows for cell- or system-level lifetime. The financial impact of battery model uncertainty can be estimated by incorporating uncertainty into a technoeconomic model.

battery↗

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗

Unraveling TeV halos with the Cherenkov Telescope Array

Pulsars are observed to emit bright and spatially extended gamma-ray emission at multi-TeV energies. These so-called "TeV halos" are now understood to be a nearly universal feature of middle-aged pulsars. However, many of the key physical processes that govern these systems, particularly those affecting particle diffusion, remain poorly constrained. We aim to evaluate the ability of the Cherenkov Telescope Array (CTA) to probe the physical properties of TeV halos, with a focus on the nearby and well-studied case of the Geminga pulsar. We simulate gamma-ray emission from various TeV halo models, incorporating different assumptions for the injected electron spectrum, spin-down evolution, and energy-dependent diffusion. These models are then used to forecast CTA's sensitivity to spectral and spatial differences, based on realistic mock observations and instrument response simulations. We find that CTA will be able to distinguish between a wide range of TeV halo models that are currently consistent with existing data. In particular, CTA observations can constrain the normalization, energy dependence, and spatial extent of the diffusion coefficient surrounding Geminga, as well as the spectral shape of the injected electron population.

79 ASTRONOMY AND ASTROPHYSICS↗

Some effects of topography, soil moisture, and sea-surface temperature on continental precipitation as computed with the GISS coarse mesh climate model

The effects of terrain elevation, soil moisture, and zonal variations in sea/surface temperature on the mean daily precipitation rates over Australia, Africa, and South America in January were evaluated. It is suggested that evaporation of soil moisture may either increase or decrease the model generated precipitation, depending on the surface albedo. It was found that a flat, dry continent model best simulates the January rainfall over Australia and South America, while over Africa the simulation is improved by the inclusion of surface physics, specifically soil moisture and albedo variations.

Spar, J.↗

Tropospheric gas composition and cloud structure of the Jovian North Equatorial Belt

Voyager IRIS observations of the NEB reveal longitudinal variability of 5-micron brightness temperatures of order 100 C. An anisotropic multiple scattering radiative transfer model is used to calculate synthetic spectra for comparison with the IRIS observations. Mie theory is used to model the spectral dependence of cloud extinction from 180 to 2300/cm. Cloud base locations within the model vary with assumed gas abundances according to thermochemical equilibrium. It is found that spatial variations in the abundance profiles of the condensible species, parahydrogen profiles and cloud optical depths can be used as tracers of the local and large-scale dynamics. The variation of cloud opacity is strongly correlated with the variation of relative humidity, which suggests that dynamic depletion of water vapor above the cloud forming level is the most plausible model to explain the spatial variation in the water profile within the NEB.

Carlson, Barbara E.↗

Recent observations and analyses of heliospheric 2-3 kHz radio emissions

Since the discovery of the 2-3 kHz heliospheric radio emissions in the Voyager data over twelve years ago, there have been two major events, the first in 1983-84 and the second in 1992-93, as well as several minor events. Strong evidence now exists that these radio emissions are generated by interactions that take place in the outer regions of the heliosphere, most likely in response to strong interplanetary shocks and associated disturbances propagating outward from the Sun. Just where these interactions take place is still a subject of controversy. The three main possibilities are at or near the termination shock, between the termination shock and the heliopause and at or near the heliopause. In this paper we will review the present state of knowledge concerning the 2-3 kHz radio emissions, including the most recent data from Voyagers 1 and 2. These observations will then be compared to models that have been proposed to explain the origin of the radio emissions. Depending on the model, various estimates can be made concerning the plasma densities in the region where the radio emission is produced, and the distance to the source. The constraints provided by these models on the size and structure of the heliosphere are discussed.

Gurnett, D. A.↗

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave transitionally in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave 'transitionally' in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗

A chemical model of the troposphere and stratosphere

A time-dependent, one-dimensional model of the coupled chemistry and vertical mixing of the atmosphere is used to compute the distribution for many atmospheric constituents of the troposphere and stratosphere. The model treats the photochemistry of the carbon-hydrogen-oxygen-nitrogen system and is (with the exception of H2O) self-consistent in the sense of requiring no assumptions regarding minor constituent distributions.

Stewart, R. W.↗

Model simulations of rainout and washout from a warm stratiform cloud

Calculations utilizing a one-dimensional, time-dependent stratiform cloud model with parameterized microphysics to study the washout and rainout of soluble trace gases are presented. It is shown that the in-cloud concentrations of relatively insoluble species remain essentially unaffected by the presence of cloudwater and rainwater, while the concentrations of highly soluble species are rapidly depleted within the cloud. Due to this depletion it is shown that for highly soluble species the below-cloud scavenging via washout can dominate over in-cloud scavenging via rainout; in-cloud, aqueous-phase cloudwater and rainwater concentrations can differ considerably and, near the cloudbase, cloudwater concentrations can exceed rainwater concentrations by factors of two to ten.

Xing, Lin↗

A two-dimensional kinematic dynamo model of the ionospheric magnetic field at Venus

The results of a high-resolution, two-dimensional, time dependent, kinematic dynamo model of the ionospheric magnetic field of Venus are presented. Various one-dimensional models are considered and the two-dimensional model is then detailed. In this model, the two-dimensional magnetic induction equation, the magnetic diffusion-convection equation, is numerically solved using specified plasma velocities. Origins of the vertical velocity profile and of the horizontal velocities are discussed. It is argued that the basic features of the vertical magnetic field profile remain unaltered by horizontal flow effects and also that horizontal plasma flow can strongly affect the magnetic field for altitudes above 300 km.

Cravens, T. E.↗

Designing for time-dependent material response in spacecraft structures

To study the influence on overall deformations of the time-dependent constitutive properties of fiber-reinforced polymeric matrix composite materials being considered for use in orbiting precision segmented reflectors, simple sandwich beam models are developed. The beam models include layers representing the face sheets, the core, and the adhesive bonding of the face sheets to the core. A three-layer model lumps the adhesive layers with the face sheets or core, while a five-layer model considers the adhesive layers explicitly. The deformation response of the three-layer and five-layer sandwich beam models to a midspan point load is studied. This elementary loading leads to a simple analysis, and it is easy to create this loading in the laboratory. Using the correspondence principle of viscoelasticity, the models representing the elastic behavior of the two beams are transformed into time-dependent models. Representative cases of time-dependent material behavior for the facesheet material, the core material, and the adhesive are used to evaluate the influence of these constituents being time-dependent on the deformations of the beam. As an example of the results presented, if it assumed that, as a worst case, the polymer-dominated shear properties of the core behave as a Maxwell fluid such that under constant shear stress the shear strain increases by a factor of 10 in 20 years, then it is shown that the beam deflection increases by a factor of 1.4 during that time. In addition to quantitative conclusions, several assumptions are discussed which simplify the analyses for use with more complicated material models. Finally, it is shown that the simpler three-layer model suffices in many situations.

Hyer, M. W.↗

An implicit adaptation algorithm for a linear model reference control system

This paper presents a stable implicit adaptation algorithm for model reference control. The constraints for stability are found using Lyapunov's second method and do not depend on perfect model following between the system and the reference model. Methods are proposed for satisfying these constraints without estimating the parameters on which the constraints depend.

Mabius, L.↗

Latitude-dependent sensitivity to stationary perturbations in simple climate models

The steady-state zonally averaged climate is perturbed by adding a latitude-dependent heat source to an energy balance equation of the simplified Budyko-Sellers type. The latitude of the ice edge, which is attached to an isotherm, becomes dependent on the strength of the perturbation. This dependence is given in terms of the well-known iceline-solar constant relation, and the latitude dependence of the perturbed temperature field is then uniquely determined. The exact analytical solution is linearized and expressed in terms of a superposition of line sources at various latitudes. The main features are: (1) The total temperature response is a sum of the direct effect of the perturbation and an indirect ice-albedo effect proportional to the solar ice-edge sensitivity; and (2) the indirect feedback effect produces an enhanced response in polar latitudes.

Salmun, H.↗

An efficient computational method for characterizing the effects of random surface errors on the average power pattern of reflectors

Based on the works of Ruze (1966) and Vu (1969), a novel mathematical model has been developed to determine efficiently the average power pattern degradations caused by random surface errors. In this model, both nonuniform root mean square (rms) surface errors and nonuniform illumination functions are employed. In addition, the model incorporates the dependence on F/D in the construction of the solution. The mathematical foundation of the model rests on the assumption that in each prescribed annular region of the antenna, the geometrical rms surface value is known. It is shown that closed-form expressions can then be derived, which result in a very efficient computational method for the average power pattern. Detailed parametric studies are performed with these expressions to determine the effects of different random errors and illumination tapers on parameters such as gain loss and sidelobe levels. The results clearly demonstrate that as sidelobe levels decrease, their dependence on the surface rms/wavelength becomes much stronger and, for a specified tolerance level, a considerably smaller rms/wavelength is required to maintain the low sidelobes within the required bounds.

Rahmat-Samii, Y.↗

Microwave Radiometric Signatures of Different Surface Types in Deserts

In arid environments, specific microwave signatures have been observed with the Special Sensor Microwave/Imager (SSM/I). For a given diurnal change in surface skin temperature, the corresponding change in the microwave brightness temperature is smaller than expected. With the help of a 1D, time-dependent heat conduction model, this behavior is explained by microwave radiation coming from different depths in the soil, depending on the soil type and on the microwave radiation frequency. Using the eight-times daily estimates of the surface skin temperature by the International Satellite Cloud Climatology Project (ISCCP) and a simple Fresnel model, collocated month-long time series of the SSM/I brightness temperatures and the surface skin temperatures give a consistent estimate of the effective microwave emissivity and penetration depth parameters. Results are presented and analyzed for the Sahara and the Arabian Peninsula, for July and November 1992. The case of the Australian desert is also briefly mentioned. Assuming a reasonable thermal diffusivity for the soil in desert areas, the microwave radiation is estimated to come from soil layers down to depths of at least five wavelengths in some locations. Regions where the microwave radiation comes from deeper soil layers also have large microwave emissivity polarization differences and large visible reflectances, suggesting that these areas correspond to sand dune fields.

Prigent, Catherine↗

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗