Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random variables”

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 631 records · Page 35

Utilizing physics-based input features within a machine learning model to predict wind speed forecasting error

Machine learning is quickly becoming a commonly used technique for wind speed and power forecasting. Many machine learning methods utilize exogenous variables as input features, but there remains the question of which atmospheric variables are most beneficial for forecasting, especially in handling non-linearities that lead to forecasting error. This question is addressed via creation of a hybrid model that utilizes an autoregressive integrated moving-average (ARIMA) model to make an initial wind speed forecast followed by a random forest model that attempts to predict the ARIMA forecasting error using knowledge of exogenous atmospheric variables. Variables conveying information about atmospheric stability and turbulence as well as inertial forcing are found to be useful in dealing with non-linear error prediction. Streamwise wind speed, time of day, turbulence intensity, turbulent heat flux, vertical velocity, and wind direction are found to be particularly useful when used in unison for hourly and 3 h timescales. The prediction accuracy of the developed ARIMA–random forest hybrid model is compared to that of the persistence and bias-corrected ARIMA models. The ARIMA–random forest model is shown to improve upon the latter commonly employed modeling methods, reducing hourly forecasting error by up to 5 % below that of the bias-corrected ARIMA model and achieving an R 2 value of 0.84 with true wind speed.

17 WIND ENERGY↗

The Ability of GeoCarb to Constrain the Interannual Variability of Carbon Gases over the Amazon

We perform a number of idealized assimilation experiments with the GEOS constituent data assimilation system to test the ability of GeoCarb retrievals of CO, CO2, and CH4 to constrain the interannual variability of these gases over the Amazon. Retrievals for instruments on other satellites which observe in similar channels (e.g. MOPITT, GOSAT, and OCO-2) are limited due to persistent cloud coverage. Given its ability to sample the same location multiple times in one day, the expectation is that GeoCarb retrievals will return more soundings than those from previous missions. The goal of the assimilation experiments is to understand which scanning strategies lead to the best sounding densities and thus have the best chance of constraining interannual variability in the carbon species. The experiments each begin by picking a given year at random from a nature run (i.e., a model simulation meant to represent the truth). The model fields are sampled according to a given strategy and then screened to account for cloud coverage. Next, we pick another year at random and assimilate the synthetic GeoCarb samples into the GEOS model for that year. The output of the assimilation, 6-hourly, 3D fields of each constituent, is then directly comparable to the nature run. This comparison allows us to evaluate the ability of GeoCarb measurements to constrain the interannual variability of each gas.

Weir, B.↗

Dual origin of ferropericlase inclusions within super-deep diamonds

Ferropericlase [(Mg,Fe)O] is one of the major constituents of Earth's lower mantle and the most abundant mineral inclusion in sub-lithospheric diamonds. Although a lower mantle origin for ferropericlase inclusions has often been suggested, some studies have proposed that many of these inclusions may instead form at much shallower depths, in the deep upper mantle or transition zone. No straightforward method exists to discriminate ferropericlase of lower-mantle origin without characteristic mineral associations, such as co-existing former bridgmanite. To explore ferropericlase-diamond growth relationships, we have investigated the crystallographic orientation relationships (CORs), determined by single-crystal X-ray diffraction, between 57 ferropericlase inclusions and 37 diamonds from Juina (Brazil) and Kankan (Guinea). We show that ferropericlase inclusions can develop specific (16 inclusions in 12 diamonds), rotational statistical (9 inclusions in 7 diamonds) and random (32 inclusions in 25 diamond) CORs with respect to their diamond hosts. All measured inclusions showing a specific COR were found to be Fe-rich (X FeO > 0.30). Coexistence of non-randomly and randomly oriented ferropericlase inclusions within the same diamond indicates that their CORs may be variably affected by local growth conditions. However, the occurrence of specific CORs only for Fe-rich inclusions indicates that Fe-rich ferropericlases have a distinct genesis and are syngenetic with their host diamonds. In conclusion, this result provides strong support for a dual origin for ferropericlase in Earth's mantle, with Fe-rich compositions likely indicating redox growth in the upper mantle, while more Mg-rich compositions with random COR mostly representing ambient lower mantle trapped as protogenetic inclusions.

58 GEOSCIENCES↗

The Chaotic Light Curves of Accreting Black Holes

We present model light curves for accreting Black Hole Candidates (BHC) based on a recently developed model of these sources. According to this model, the observed light curves and aperiodic variability of BHC are due to a series of soft photon injections at random (Poisson) intervals and the stochastic nature of the Comptonization process in converting these soft photons to the observed high energy radiation. The additional assumption of our model is that the Comptonization process takes place in an extended but non-uniform hot plasma corona surrounding the compact object. We compute the corresponding Power Spectral Densities (PSD), autocorrelation functions, time skewness of the light curves and time lags between the light curves of the sources at different photon energies and compare our results to observation. Our model reproduces the observed light curves well, in that it provides good fits to their overall morphology (as manifest by the autocorrelation and time skewness) and also to their PSDs and time lags, by producing most of the variability power at time scales 2 a few seconds, while at the same time allowing for shots of a few msec in duration, in accordance with observation. We suggest that refinement of this type of model along with spectral and phase lag information can be used to probe the structure of this class of high energy sources.

Kazanas, Demosthenes↗

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗

Metal Foam Morphology Affects the Run to Run Reproducibility of OER Using Nickel Catalysts

Nickel (Ni) foam-based electrodes are excellent catalysts for the oxygen evolution reaction; however, we found that the random pore-size distribution of Ni foams contributes to a significant variability in electrochemically active surface area, compromising experimental reproducibility. We provide insights into quantifying this critical material property, verified by four electrochemical laboratories.

25 ENERGY STORAGE↗

Probabilistic measures for biological adaptation and resilience

This paper introduces an approach to quantifying ecological resilience in biological systems, particularly focusing on noisy systems responding to episodic disturbances with sudden adaptations. Incorporating concepts from nonequilibrium statistical mechanics, we propose a measure termed “ecological resilience through adaptation,” specifically tailored to noisy, forced systems that undergo physiological adaptation in the face of stressful environmental changes. Randomness plays a key role, accounting for model uncertainty and the inherent variability in the dynamical response among components of biological systems. Our measure of resilience is rooted in the probabilistic description of states within these systems and is defined in terms of the dynamics of the ensemble average of a model-specific observable quantifying success or well-being. Our approach utilizes stochastic linear response theory to compute how the expected success of a system, originally in statistical equilibrium, dynamically changes in response to a environmental perturbation and a subsequent adaptation. Importantly, the resulting mathematical derivations allow for the estimation of resilience in terms of ensemble averages of simulated or experimental data. Finally, through a simple but clear conceptual example, we illustrate how our resilience measure can be interpreted and compared to other existing frameworks in the literature. The methodology is general but inspired by applications in plant systems, with the potential for broader application to complex biological processes.

60 APPLIED LIFE SCIENCES↗

Equatorial waves triggering extreme rainfall and floods in southwest Sulawesi, Indonesia

On the basis of detailed analysis of a case study and long-term climatology, it is shown that equatorial waves and their interactions serve as precursors for extreme rain and flood events in the central Maritime Continent region of southwest Sulawesi, Indonesia. Meteorological conditions on January 22, 2019, leading to heavy rainfall and devastating flooding in this area were studied. It is shown that a convectively coupled Kelvin wave (CCKW) and a convectively coupled Rossby wave (CCERW) embedded within the larger-scale envelope of the Madden-Julian Oscillation (MJO) enhanced convective phase, contributed to the onset of a mesoscale convective system which developed over the Java Sea. Low-Level convergence from the CCKW forced mesoscale convective organization and orographic ascent of moist air over the slopes of southwest Sulawesi. Climatological analysis shows that 92 % of December-January-February floods and 76% of extreme rain events in this region were immediately preceded by positive low-level westerly wind anomalies. It is estimated that both CCKWs and CCERWs propagating over Sulawesi double the chance of floods and extreme rain event development, which the probability of such hazardous events occurring during their combined activity is eight times greater than on a random day. While the MJO is a key component shaping tropical atmospheric variability, it is shown that its usefulness as a single factor for extreme weather-driven hazard prediction is limited.

Latos, Beata↗

Internally and externally caused climate change

A numerical climate model is used to simulate climate change forced only by random fluctuations of the atmospheric heat transport. This short-term natural variability of the atmosphere is shown to be a possible 'cause' not only of the variability of the annual world average temperature about its mean, but also long-term excursions from the mean. Various external causes of climate change are also tested with the model and the results compared with observations for the past 100 years. Volcanic dust is shown to have been an important cause of climate change, while the effects of sunspot-related solar constant variation and anthropogenic forcing are not evident.

Robock, A.↗

Radio emission evolution of nonstationary sources in the Hedgehog model

Correlations are obtained for numerical calculation of flux F sub v and polarized radiation intensity of a cloud of arbitrary geometry, consisting of ultrarelativistic electrons that dissipate in a radial magnetic field of the nucleus at a random angle to the observer. It is possible that some of the variable extragalactic objects that were previously described by the Shklovskiy model are young formations in the examined model. Radio astronomical observations would permit a determination of their distance, age, and lifetime.

Kovalev, Y. A.↗

Modal methods in optimal control synthesis

Efficient algorithms for solving linear smoother-follower problems with quadratic criteria are presented. For time-invariant systems, the algorithm consists of one backward integration of a linear vector equation and one forward integration of another linear vector equation. Furthermore, the backward and forward Riccati matrices can be expressed in terms of the eigenvalues and eigenvectors of the Euler-Lagrange equations. Hence, the gains of the forward and backward Kalman-Bucy filters and of the optimal state-feedback regulator can be determined without integration of matrix Riccati equations. A computer program has been developed, based on this method of determining the gains, to synthesize the optimal time-invariant compensator in the presence of random disturbance inputs and random measurement errors. The program also computes the rms state and control variables of the optimal closed-loop system.

Bryson, A. E., Jr.↗

Reliability-Based Design Optimization of a Composite Airframe Component

A stochastic optimization methodology (SDO) has been developed to design airframe structural components made of metallic and composite materials. The design method accommodates uncertainties in load, strength, and material properties that are defined by distribution functions with mean values and standard deviations. A response parameter, like a failure mode, has become a function of reliability. The primitive variables like thermomechanical loads, material properties, and failure theories, as well as variables like depth of beam or thickness of a membrane, are considered random parameters with specified distribution functions defined by mean values and standard deviations.

Pai, Shantaram S.↗

Western Montana Ecological Forecasting: Modeling Habitat Suitability of Mustelid Species to Guide Detection Dog Surveys for Contaminants Monitoring, via Collected Scats in River Systems of Western Montana

Environmental contaminants are becoming increasingly prevalent in riverine ecosystems. The status of contaminants in western Montana’s relatively pristine river systems is largely unknown. Monitoring for heavy metals, brominated flame-retardants (BFRs), and pharmaceuticals is important due to their negative effects on ecosystems. Exposure to these contaminants can have significant endocrine, neurological, and reproductive effects. Contaminants easily travel up the food chain and bioaccumulate in apex predators. As predators with a largely aquatic diet, American mink (Mustela vison) and North American river otter (Lontra canadensis) serve as reliable indicator species of environmental health and the status of contaminants. Analysis of scat from these species is a noninvasive method to measure contaminant levels, and detection dogs from Working Dogs for Conservation (WD4C) have been used to locate these scat samples. To aid in the search of these samples, habitat suitability models were created for mink and otter for the years 2013-2020 and projected to 2040 using the random forest algorithm in the Software for Assisted Habitat Modeling (SAHM). Predictor variable data were acquired from Landsat 8 Operational Land Imager (OLI), Terra Moderate Resolution Imaging Spectroradiometer (MODIS), Global Precipitation Measurement Integrated Multi-satellite Retrievals for GPM (GPM IMERG), Shuttle Radar Topography Mission (SRTM), and Soil Moisture Active Passive (SMAP). Within these models, the most important variable for mink and otter habitat was distance to river. Suitable habitat also corresponded with emergent herbaceous land cover and deeper river locations. These habitat suitability models will inform sampling site section for further contaminant analysis.

Anna Winter↗

The design of a turboshaft speed governor using modern control techniques

The objectives of this program were: to verify the model of off schedule compressor variable geometry in the T700 turboshaft engine nonlinear model; to evaluate the use of the pseudo-random binary noise (PRBN) technique for obtaining engine frequency response data; and to design a high performance power turbine speed governor using modern control methods. Reduction of T700 engine test data generated at NASA-Lewis indicated that the off schedule variable geometry effects were accurate as modeled. Analysis also showed that the PRBN technique combined with the maximum likelihood model identification method produced a Bode frequency response that was as accurate as the response obtained from standard sinewave testing methods. The frequency response verified the accuracy of linear models consisting of engine partial derivatives and used for design. A power turbine governor was designed using the Linear Quadratic Regulator (LQR) method of full state feedback control. A Kalman filter observer was used to estimate helicopter main rotor blade velocity. Compared to the baseline T700 power turbine speed governor, the LQR governor reduced droop up to 25 percent for a 490 shaft horsepower transient in 0.1 sec simulating a wind gust, and up to 85 percent for a 700 shaft horsepower transient in 0.5 sec simulating a large collective pitch angle transient.

Delosreyes, G.↗

Southern Rockies Western Slope Agriculture: Identifying Drivers of Rangeland Production for Drought Planning on the Western Slope of the Southern Rockies

Over the last decade, the southern Rocky Mountains of the United States experienced severe and variable drought. Local ranchers and landowners have reported strain on their operations, citing decreasing forage for their cattle and a need to adjust their business models. This study identified Major Land Resource Area-48 (MLRA-48) and northwestern Colorado as the key region for analysis. NASA DEVELOP partnered with the BLM Colorado River Field Office, Colorado State University Extension, USDA Forest Service, and the National Drought Mitigation Center to address concerns regarding the efficacy of remotely sensed rangeland production platforms and identify early warning climatic indicators of drought. The study identified two key platforms, The Rangeland Productivity Monitoring Service (RPMS) and Rangeland Analysis Platform (RAP), which use NASA Landsat 5 TM, Landsat 7 ETM+, Landsat 8 OLI, and Landsat 9 OLI-2 to estimate rangeland biomass. We regressed these with in situ biomass data to validate their efficacy and found that RAP was more effective than RPMS in estimating rangeland biomass, though it presents a tendency to overestimate. Our study performed a random forest analysis, comparing monthly RAP biomass estimates to a variety of climate variables, including mean precipitation, temperature, Palmer Drought Severity Index, snow water equivalent, snow persistence from Terra MODIS, wind speed and direction, and vapor pressure deficit. We determined that vapor pressure deficit and precipitation are key indicators in predicting forage production in MLRA-48. Our climate analysis provided our partners with greater understanding of the influence of various climate variables in determining rangeland production and allows them to assist land managers in drought mitigation.

remote sensing↗

Comparison of linear regression, k-nearest neighbour and random forest methods in airborne laser-scanning-based prediction of growing stock

Abstract In this study, for five sites around the world, we look at the effects of different model types and variable selection approaches on forest yield modelling performances in an area-based approach (ABA). We compared ordinary least squares regression (OLS), k-nearest neighbours (kNN) and random forest (RF). Our objective was to test if there are systematic differences in accuracy between OLS, kNN and RF in ABA predictions of growing stock volume. The analyses are based on a 5-fold cross-validation at five study sites: an eucalyptus plantation, a temperate forest and three different boreal forests. Two completely independent validation datasets were also available for two of the boreal sites. For the kNN, we evaluated multiple measures of distance including Euclidean, Mahalanobis, most similar neighbour (MSN) and an RF-based distance metric. The variable selection approaches we examined included a heuristic approach (for OLS, kNN and RF), exhaustive search among all combinations (OLS only) and all variables together (RF only). Performances varied by model type and variable selection approaches among sites. OLS and RF had similar accuracies and were more efficient than any of the kNN variants. Variable selection did not affect RF performance. Heuristic and exhaustive variable selection performed similarly for OLS. kNN fared the poorest amongst model types, and kNN with RF distance was prone to overfitting when compared with a validation dataset. Additional caution is therefore required when building kNN models for volume prediction though ABA, being preferable instead to opt for models based on OLS with some variable selection, or RF with all variables together.

Cosenza, Diogo N.↗

Modeling freight mode choice using machine learning classifiers: a comparative study using Commodity Flow Survey (CFS) data

This study explores the usefulness of machine learning classifiers for modeling freight mode choice. We investigate eight commonly used machine learning classifiers, namely Naïve Bayes, Support Vector Machine, Artificial Neural Network, K-Nearest Neighbors, Classification and Regression Tree, Random Forest, Boosting and Bagging, along with the classical Multinomial Logit model. US 2012 Commodity Flow Survey data are used as the primary data source; we augment it with spatial attributes from secondary data sources. The performance of the classifiers is compared based on prediction accuracy results. The current research also examines the role of sample size and training-testing data split ratios on the predictive ability of the various approaches. In addition, the importance of variables is estimated to determine how the variables influence freight mode choice. The results show that the tree-based ensemble classifiers perform the best. Specifically, Random Forest produces the most accurate predictions, closely followed by Boosting and Bagging. With regard to variable importance, shipment characteristics, such as shipment distance, industry classification of the shipper and shipment size, are the most significant factors for freight mode choice decisions.

42 ENGINEERING↗

Recursive Branching Simulated Annealing Algorithm

This innovation is a variation of a simulated-annealing optimization algorithm that uses a recursive-branching structure to parallelize the search of a parameter space for the globally optimal solution to an objective. The algorithm has been demonstrated to be more effective at searching a parameter space than traditional simulated-annealing methods for a particular problem of interest, and it can readily be applied to a wide variety of optimization problems, including those with a parameter space having both discrete-value parameters (combinatorial) and continuous-variable parameters. It can take the place of a conventional simulated- annealing, Monte-Carlo, or random- walk algorithm. In a conventional simulated-annealing (SA) algorithm, a starting configuration is randomly selected within the parameter space. The algorithm randomly selects another configuration from the parameter space and evaluates the objective function for that configuration. If the objective function value is better than the previous value, the new configuration is adopted as the new point of interest in the parameter space. If the objective function value is worse than the previous value, the new configuration may be adopted, with a probability determined by a temperature parameter, used in analogy to annealing in metals. As the optimization continues, the region of the parameter space from which new configurations can be selected shrinks, and in conjunction with lowering the annealing temperature (and thus lowering the probability for adopting configurations in parameter space with worse objective functions), the algorithm can converge on the globally optimal configuration. The Recursive Branching Simulated Annealing (RBSA) algorithm shares some features with the SA algorithm, notably including the basic principles that a starting configuration is randomly selected from within the parameter space, the algorithm tests other configurations with the goal of finding the globally optimal solution, and the region from which new configurations can be selected shrinks as the search continues. The key difference between these algorithms is that in the SA algorithm, a single path, or trajectory, is taken in parameter space, from the starting point to the globally optimal solution, while in the RBSA algorithm, many trajectories are taken; by exploring multiple regions of the parameter space simultaneously, the algorithm has been shown to converge on the globally optimal solution about an order of magnitude faster than when using conventional algorithms. Novel features of the RBSA algorithm include: 1. More efficient searching of the parameter space due to the branching structure, in which multiple random configurations are generated and multiple promising regions of the parameter space are explored; 2. The implementation of a trust region for each parameter in the parameter space, which provides a natural way of enforcing upper- and lower-bound constraints on the parameters; and 3. The optional use of a constrained gradient- search optimization, performed on the continuous variables around each branch s configuration in parameter space to improve search efficiency by allowing for fast fine-tuning of the continuous variables within the trust region at that configuration point.

Bolcar, Matthew↗