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,261 records · Page 70

Inferring Instantaneous, Multivariate and Nonlinear Sensitivities for the Analysis of Feedback Processes in a Dynamical System: Lorenz Model Case Study

A new approach is presented for the analysis of feedback processes in a nonlinear dynamical system by observing its variations. The new methodology consists of statistical estimates of the sensitivities between all pairs of variables in the system based on a neural network modeling of the dynamical system. The model can then be used to estimate the instantaneous, multivariate and nonlinear sensitivities, which are shown to be essential for the analysis of the feedbacks processes involved in the dynamical system. The method is described and tested on synthetic data from the low-order Lorenz circulation model where the correct sensitivities can be evaluated analytically.

Aires, Filipe↗

Sunspot Pattern Classification using PCA and Neural Networks (Poster)

The sunspot classification scheme presented in this paper is considered as a 2-D classification problem on archived datasets, and is not a real-time system. As a first step, it mirrors the Zuerich/McIntosh historical classification system and reproduces classification of sunspot patterns based on preprocessing and neural net training datasets. Ultimately, the project intends to move from more rudimentary schemes, to develop spatial-temporal-spectral classes derived by correlating spatial and temporal variations in various wavelengths to the brightness fluctuation spectrum of the sun in those wavelengths. Once the approach is generalized, then the focus will naturally move from a 2-D to an n-D classification, where "n" includes time and frequency. Here, the 2-D perspective refers both to the actual SOH0 Michelson Doppler Imager (MDI) images that are processed, but also refers to the fact that a 2-D matrix is created from each image during preprocessing. The 2-D matrix is the result of running Principal Component Analysis (PCA) over the selected dataset images, and the resulting matrices and their eigenvalues are the objects that are stored in a database, classified, and compared. These matrices are indexed according to the standard McIntosh classification scheme.

Rajkumar, T.↗

Thermophysical Properties Along Curiosity's Traverse in Gale Crater, Mars, Derived from the REMS Ground Temperature Sensor

The REMS instrument onboard the Mars Science Laboratory rover, Curiosity, has measured ground temperature nearly continuously at hourly intervals for two Mars years. Coverage of the entire diurnal cycle at 1 Hz is available every few martian days. We compare these measurements with predictions of surface atmosphere thermal models to derive the apparent thermal inertia and thermally derived albedo along the rovers traverse after accounting for the radiative effects of atmospheric water ice during fall and winter, as is necessary to match the measured seasonal trend. The REMS measurements can distinguish between active sand, other loose materials, mudstone, and sandstone based on their thermophysical properties. However, the apparent thermal inertias of bedrock dominated surfaces [approx. 350-550 J m(exp. -2) K(exp. -1 s(exp. -1/2 )] are lower than expected. We use rover imagery and the detailed shape of the diurnal ground temperature curve to explore whether lateral or vertical heterogeneity in the surface materials within the sensor footprint might explain the low inertias. We find that the bedrock component of the surface can have a thermal inertia as high as 650-1700 J m(exp. -2) K(exp. -1) s(exp. -1/2) for mudstone sites and approx. 700 J m(exp. -2) K(exp. -1) s(exp. - 1/2) for sandstone sites in models runs that include lateral and vertical mixing. Although the results of our forward modeling approach may be non-unique, they demonstrate the potential to extract information about lateral and vertical variations in thermophysical properties from temporally resolved measurements of ground temperature.

Vasavada, Ashwin R.↗

The Past and Future Changes in Climate of the Rice-Wheat Cropping Zone in Punjab, Pakistan

Agriculture is ranked top among Pakistan economic sectors vulnerable to the potential impacts of climate change. The agricultural production system is directly affected by weather inputs (temperature, solar radiation, and rainfall) that are projected to change in future (following increases in carbon dioxide and other greenhouse gasses). Although Climatic extremes such as drought, floods and heat waves, are expected to increase with detrimental consequences for agriculture and livestock production, yet changes in mean climates also pose challenges to sustainable development. This study presents climate change results for five districts within the major rice-wheat productivity zone of Punjab province in Pakistan. The results are focused on RCP8.5 mid-century (2040-2069) scenarios derived from five global climate models (GCMs) output and the Kharif (June-October) and Rabi (November-April) seasons. Analysis of recent historical weather data of Sialkot and Sheikhupura districts shows an increase in minimum temperatures and maximum temperatures and a large variation in rainfall. This temperature change and variability in rainfall are expected to enhance further as we approach 2050s. The mean maximum temperature is projected to increase by 2-2.5 ºC during the rice growing (Kharif) season and 2.4-2.7ºC during the wheat growing (Rabi) season. Rainfall during the rice growing season is more uncertain, with projections indicating an increase of 25%-35 % in the study region, while a minimal change is expected during the Rabi season. Climax of the rice growing season corresponds with the projected increase in monsoon intensity, leaving no doubt about the crop water demand satisfaction. However, a simultaneous increase in the day and night temperature may affect the growth and development of some critical phenological stages and thus could impact final yield.

Climate Change↗

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of A Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of a Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗

Krylov spaces for truncated spectrum methodologies

We propose herein an extension of truncated spectrum methodologies, a nonperturbative numerical approach able to elucidate the low energy properties of quantum field theories. TSMs, in their various flavors, involve a division of a computational Hilbert space, H , into two parts, one part, H 1 that is “kept” for the numerical computations, and one part, H 2 , that is discarded or “truncated.” Even though H 2 is discarded, truncated spectrum methodologies will often try to incorporate the effects of H 2 in some effective way. In these terms, we propose to keep the dimension of H 1 small. We pair this choice of H 1 with a Krylov subspace iterative approach able to take into account the effects of H 2 . This iterative approach can be taken to arbitrarily high order and so offers the ability to compute quantities to arbitrary precision. In many cases it also offers the advantage of not needing an explicit UV cutoff. To compute the matrix elements that arise in the Krylov iterations, we employ a Feynman diagrammatic representation that is then evaluated with Monte Carlo techniques. Each order of the Krylov iteration is variational and is guaranteed to improve upon the previous iteration. The first Krylov iteration is akin to the next-to-leading order approach of Elias-Miró [NLO renormalization in the Hamiltonian truncation, ]. To demonstrate this approach, we focus on the ( 1 + 1 d )-dimensional ϕ 4 model and compute the bulk energy and mass gaps in both the Z 2 -broken and unbroken sectors. We estimate the critical ϕ 4 coupling in the broken phase to be g c = 0.2645 ± 0.002 . Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Autoregressive distributed lag-based dynamic uniformity modeling and monitoring approaches for superconductor manufacturing

High-temperature superconductors (HTS), known for their high efficiency and low energy loss, have found profound applications across various fields, driving the demand for long, uniformly performing tapes. However, ensuring uniform performance over extended lengths of HTS tapes, often characterized by the consistency of critical current, remains challenging due to fluctuations in growth conditions during manufacturing. Here, to elucidate the mechanisms underlying variations in tape uniformity and enable real-time monitoring of associated parameters, we propose an Autoregressive Distributed Lag (ADL)-based Dynamic Uniformity Modeling and Monitoring (ADUM2) approach. This method integrates uniformity measurement, the identification of critical process parameters and real-time monitoring within the manufacturing process. The ADUM2 approach is applied to the advanced metal organic chemical vapor deposition (A-MOCVD) process, a pilot-scale method for superconductor manufacturing. Our model demonstrates superior performance compared to benchmark methods, accounting for over 80% of the total variance in the data and identifying 13 key process parameters influencing the uniformity of HTS tapes. This study offers significant insights into the high-temperature superconductor manufacturing process and holds the potential to facilitate the production of cost-effective, uniformly performing long superconducting tapes in the future.

autoregressive distributed lag analysis↗

Non-Destructive Evaluation of Materials via Ultraviolet Spectroscopy

A document discusses the use of ultraviolet spectroscopy and imaging for the non-destructive evaluation of the degree of cure, aging, and other properties of resin-based composite materials. This method can be used in air, and is portable for field use. This method operates in reflectance, absorbance, and luminescence modes. The ultraviolet source is used to illuminate a composite surface of interest. In reflectance mode, the reflected response is acquired via the imaging system or via the spectrometer. The spectra are analyzed for organic compounds (conjugated organics) and inorganic compounds (semiconducting band-edge states; luminescing defect states such as silicates, used as adhesives for composite aerospace applications; and metal oxides commonly used as thermal coating paints on a wide range of spacecraft). The spectra are compared with a database for variation in conjugation, substitution, or length of molecule (in the case of organics) or band edge position (in the case of inorganics). This approach is useful in the understanding of material quality. It lacks the precision in defining the exact chemical structure that is found in other materials analysis techniques, but it is advantageous over methods such as nuclear magnetic resonance, infrared spectroscopy, and chromatography in that it can be used in the field to assess significant changes in chemical structure that may be linked to concerns associated with weaknesses or variations in structural integrity, without disassembly of or destruction to the structure of interest.

Pugel, Betsy↗

Lidar Measurements of Atmospheric CO2 From Regional to Global Scales

Atmospheric CO2 is a critical forcing for the Earth's climate and the knowledge on its distributions and variations influences predictions of the Earth's future climate. Large uncertainties in the predictions persist due to limited observations. This study uses the airborne Intensity-Modulated Continuous-Wave (IMCW) lidar developed at NASA Langley Research Center to measure regional atmospheric CO2 spatio-temporal variations. Further lidar development and demonstration will provide the capability of global atmospheric CO2 estimations from space, which will significantly advances our knowledge on atmospheric CO2 and reduce the uncertainties in the predictions of future climate. In this presentation, atmospheric CO2 column measurements from airborne flight campaigns and lidar system simulations for space missions will be discussed. A measurement precision of approx.0.3 ppmv for a 10-s average over desert and vegetated surfaces has been achieved. Data analysis also shows that airborne lidar CO2 column measurements over these surfaces agree well with in-situ measurements. Even when thin cirrus clouds present, consistent CO2 column measurements between clear and thin cirrus cloudy skies are obtained. Airborne flight campaigns have demonstrated that precise atmospheric column CO2 values can be measured from current IM-CW lidar systems, which will lead to use this airborne technique in monitoring CO2 sinks and sources in regional and continental scales as proposed by the NASA Atmospheric Carbon and Transport â€" America project. Furthermore, analyses of space CO2 measurements shows that applying the current IM-CW lidar technology and approach to space, the CO2 science goals of space missions will be achieved, and uncertainties in CO2 distributions and variations will be reduced.

Lin, Bing↗

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗

Evaluating the potential for bird‐habitat models to support biodiversity‐friendly urban planning

Abstract Urban expansion poses a major threat to wildlife populations. Biodiversity‐friendly urban landscapes could deliver benefits for both wildlife and people, by incorporating conservation and ecosystem services objectives. Well‐designed urban developments could also soften the ecological impacts of urbanization. However, delivering urban landscapes that integrate biodiversity requirements effectively remains challenging. Ecological models, designed to predict wildlife population responses to alternative urban designs, could prove effective in supporting the creation of biodiversity‐friendly urban landscapes. Here, we combine national‐scale bird abundance data with high resolution, spatially explicit habitat data to characterize relationships between bird densities and urban landscape form in Britain. From these analyses and cross validation, we evaluate the potential for well‐parameterized, species‐specific models to be used to predict bird densities in novel or modified urban areas. Our analyses indicate that responses of bird abundance to urban habitat are species‐specific and complex, with few variables consistently affecting a large proportion of species. However, contiguous areas of greenspace within urban sites are preferential for accommodating breeding birds, compared to a more fragmented arrangement of multiple, small greenspace patches. In combination, the bird‐habitat relationships identified could successfully predict observed variation in abundance for most bird species considered. Further evaluation of habitat descriptor variables, spatial scales of species' habitat use and analytical modelling approaches may be needed to improve the predictive ability of bird‐habitat models for certain species, particularly waterbirds and those observed less frequently in urban areas. Synthesis and applications . We modelled breeding bird abundance in built‐up areas with respect to the characteristics and contexts of urban environments. While most variables were important for multiple species, responses overall were species‐specific, so simple assemblage metrics, like diversity, will not describe the variation in bird communities well. However, the results illustrate the potential of an evidence‐based, spatially explicit evaluation of urban development impacts on biodiversity, by predicting the consequences for bird numbers. Subject to verification of predictive ability, practitioners can apply the models to compare, for example, land‐sparing and sharing within developments, or to quantify the biodiversity requirements for effective offsetting. This would be facilitated by incorporation into an online tool allowing user‐determined input scenarios.

Plummer, Kate E.↗

Flagging and Correction of Pattern Noise in the Kepler Focal Plane Array

In order for Kepler to achieve its required less than 20 PPM photometric precision for magnitude 12 and brighter stars, instrument-induced variations in the CCD readout bias pattern (our "2D black image"), which are either fixed or slowly varying in time, must be identified and the corresponding pixels either corrected or removed from further data processing. The two principle sources of these readout bias variations are crosstalk between the 84 science CCDs and the 4 fine guidance sensor (FGS) CCDs and a high frequency amplifier oscillation on less than 40% of the CCD readout channels. The crosstalk produces a synchronous pattern in the 2D black image with time-variation observed in less than 10% of individual pixel bias histories. We will describe a method of removing the crosstalk signal using continuously-collected data from masked and over-clocked image regions (our "collateral data"), and occasionally-collected full-frame images and reverse-clocked readout signals. We use this same set to detect regions affected by the oscillating amplifiers. The oscillations manifest as time-varying moir pattern and rolling bands in the affected channels. Because this effect reduces the performance in only a small fraction of the array at any given time, we have developed an approach for flagging suspect data. The flags will provide the necessary means to resolve any potential ambiguity between instrument-induced variations and real photometric variations in a target time series. We will also evaluate the effectiveness of these techniques using flight data from background and selected target pixels.

Kolodziejczak, Jeffery J.↗

Assessment of shear deformation theories for multilayered composite plates

A review is made of the different approaches used for modeling multilayered composite plates. Discussion focuses on different approaches for developing two-dimensional shear deformation theories; classification of two-dimensional theories based on introducing plausible displacement, strain and/or stress assumptions in the thickness direction; and first-order shear deformation theories based on linear displacement assumptions in the thickness coordinate. Extensive numerical results are presented showing the effects of variation in the lamination and geometric parameters of simply supported composite plates on the accuracy of the static and vibrational responses predicted by six different modeling approaches (based on two-dimensional shear deformation theories). The standard of comparison is taken to be the exact three-dimensional elasticity solutions. Some of the future directions for research on the modeling of multilayered composite plates are outlined.

Noor, Ahmed K.↗

A Genome-Based Model to Predict the Virulence of Pseudomonas aeruginosa Isolates

ABSTRACT: Variation in the genome of Pseudomonas aeruginosa , an important pathogen, can have dramatic impacts on the bacterium’s ability to cause disease. We therefore asked whether it was possible to predict the virulence of P. aeruginosa isolates based on their genomic content. We applied a machine learning approach to a genetically and phenotypically diverse collection of 115 clinical P. aeruginosa isolates using genomic information and corresponding virulence phenotypes in a mouse model of bacteremia. We defined the accessory genome of these isolates through the presence or absence of accessory genomic elements (AGEs), sequences present in some strains but not others. Machine learning models trained using AGEs were predictive of virulence, with a mean nested cross-validation accuracy of 75% using the random forest algorithm. However, individual AGEs did not have a large influence on the algorithm’s performance, suggesting instead that virulence predictions are derived from a diffuse genomic signature. These results were validated with an independent test set of 25 P. aeruginosa isolates whose virulence was predicted with 72% accuracy. Machine learning models trained using core genome single-nucleotide variants and whole-genome k-mers also predicted virulence. Our findings are a proof of concept for the use of bacterial genomes to predict pathogenicity in P. aeruginosa and highlight the potential of this approach for predicting patient outcomes. IMPORTANCE Pseudomonas aeruginosa is a clinically important Gram-negative opportunistic pathogen. P. aeruginosa shows a large degree of genomic heterogeneity both through variation in sequences found throughout the species (core genome) and through the presence or absence of sequences in different isolates (accessory genome). P. aeruginosa isolates also differ markedly in their ability to cause disease. In this study, we used machine learning to predict the virulence level of P. aeruginosa isolates in a mouse bacteremia model based on genomic content. We show that both the accessory and core genomes are predictive of virulence. This study provides a machine learning framework to investigate relationships between bacterial genomes and complex phenotypes such as virulence.

59 BASIC BIOLOGICAL SCIENCES↗

The Effect of Fiber Strength Stochastics and Local Fiber Volume Fraction on Multiscale Progressive Failure of Composites

Continuous fiber unidirectional polymer matrix composites (PMCs) can exhibit significant local variations in fiber volume fraction as a result of processing conditions that can lead to further local differences in material properties and failure behavior. In this work, the coupled effects of both local variations in fiber volume fraction and the empirically-based statistical distribution of fiber strengths on the predicted longitudinal modulus and local tensile strength of a unidirectional AS4 carbon fiber/ Hercules 3502 epoxy composite were investigated using the special purpose NASA Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC); local effective composite properties were obtained by homogenizing the material behavior over repeating units cells (RUCs). The predicted effective longitudinal modulus was relatively insensitive to small (~8%) variations in local fiber volume fraction. The composite tensile strength, however, was highly dependent on the local distribution in fiber strengths. The RUC-averaged constitutive response can be used to characterize lower length scale material behavior within a multiscale analysis framework that couples the NASA code FEAMAC and the ABAQUS finite element solver. Such an approach can be effectively used to analyze the progressive failure of PMC structures whose failure initiates at the RUC level. Consideration of the effect of local variations in constituent properties and morphologies on progressive failure of PMCs is a central aspect of the application of Integrated Computational Materials Engineering (ICME) principles for composite materials.

multiscale analysis↗

Radial Variations in the Io Plasma Torus during the Cassini Era

A radial scan through the midnight sector of the Io plasma torus was made by the Cassini Ultraviolet Imaging Spectrograph on 14 January 2001, shortly after closest approach to Jupiter. From these data, Steffl et al. (2004a) derived electron temperature, plasma composition (ion mixing ratios), and electron column density as a function of radius from L = 6 to 0 as well as the total luminosity. We have advanced our homogeneous model of torus physical chemistry (Delamere and Bagenal, 2003) to include latitudinal and radial variations in a manner similar to the two-dimensional model by Schreier et al. (1998). The model variables include: (1) neutral source rate, (2) radial transport coefficient, (3) the hot electron fraction, (4) hot electron temperature, and (5) the neutral O/S ratio. The radial variation of parameters 1-4 are described by simple power laws, making a total of nine parameters. We have explored the sensitivity of the model results to variations in these parameters and compared the best fit with previous Voyager era models (schreier et al., 1998), galileo data (Crary et al., 1998), and Cassini observations (steffl et al., 2004a). We find that radial variations during the Cassini era are consistent with a neutral source rate of 700-1200 kg/s, an integrated transport time from L = 6 to 9 of 100-200 days, and that the core electron temperature is largely determined by a spatially and temporally varying superthermal electron population.

Delamere, P. A.↗

VpROM: a novel variational autoencoder-boosted reduced order model for the treatment of parametric dependencies in nonlinear systems

Reduced Order Models (ROMs) are of considerable importance in many areas of engineering in which computational time presents difficulties. Established approaches employ projection-based reduction, such as Proper Orthogonal Decomposition. The limitation of the linear nature of such operators is typically tackled via a library of local reduction subspaces, which requires the assembly of numerous local ROMs to address parametric dependencies. Our work attempts to define a more generalisable mapping between parametric inputs and reduced bases for the purpose of generative modeling. We propose the use of Variational Autoencoders (VAEs) in place of the typically utilised clustering or interpolation operations, for inferring the fundamental vectors, termed as modes, which approximate the manifold of the model response for any and each parametric input state. The derived ROM still relies on projection bases, built on the basis of full-order model simulations, thus retaining the imprinted physical connotation. However, it additionally exploits a matrix of coefficients that relates each local sample response and dynamics to the global phenomena across the parametric input domain. The VAE scheme is utilised for approximating these coefficients for any input state. This coupling leads to a high-precision low-order representation, which is particularly suited for problems where model dependencies or excitation traits cause the dynamic behavior to span multiple response regimes. Moreover, the probabilistic treatment of the VAE representation allows for uncertainty quantification on the reduction bases, which may then be propagated to the ROM response. The performance of the proposed approach is validated on an open-source simulation benchmark featuring hysteresis and multi-parametric dependencies, and on a large-scale wind turbine tower characterised by nonlinear material behavior and model uncertainty.

Conditional VAEs↗