Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Improved thermonuclear rate of 42 Ti( p , γ ) 43 V and its astrophysical implication in the rp process

Context. Accurate 42 Ti(p,γ) 43 V reaction rates are crucial for understanding the nucleosynthesis path of the rapid capture process (rp process) that occurs in X-ray bursts. Aims. We aim to improve the thermonuclear rates of 42 Ti(p,γ) 43 V based on more complete resonance information and a more accurate direct component, together with the recently released nuclear masses data. We also explore the impact of the newly obtained rates on the rp process. Methods. We reevaluated the reaction rate of 42 Ti(p,γ) 43 V by the sum of the isolated resonance contribution instead of the Hauser-Feshbach statistical model. We used a Monte Carlo method to derive the associated uncertainties of new rates. The nucleosynthesis simulations were performed via the NuGrid post-processing code ppn. Results. The new rates differ from previous estimations due to the use of a series of updated resonance parameters and a direct S factor. Compared with the previous results from the Hauser-Feshbach statistical model, which assumes compound nucleus 43 V with a sufficiently high-level density in the energy region of astrophysical interest, large differences exist over the entire temperature region of rp-process interest, up to two orders of magnitude. We consistently calculated the photodisintegration rate using our new nuclear masses via the detailed balance principle, and found the discrepancies among the different reverse rates are much larger than those for the forward rate, up to ten orders of magnitude at the temperature of 10 8 K. Using a trajectory with a peak temperature of 1.95×10 9 K, we performed the rp-process nucleosynthesis simulations to investigate the impact of the new rates. Our calculations show that the adoption of the new forward and reverse rates result in abundance variations for Sc and Ca of 128% and 49%, respectively, compared to the variations for the statistical model rates. On the other hand, the overall abundance pattern is not significantly affected. The results of using new rates also confirm that the rp-process path does not bypass the isotope 43 V. Conclusions. Our study found that the Hauser-Feshbach statistical model is inappropriate to the reaction rate evaluation for 42 Ti(p,γ) 43 V. The adoption of the new rates confirms that the reaction path of 42 Ti(p,γ) 43 V(p,γ) 44 Cr(β + ) 44 V is a key branch of the rp process in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Investigation into alternative statistical auroral oval models, stage 1

The statistical model sensitivity to different methods of binning the NOAA/TIROS data base was examined. Each partial auroral pass was analyzed to locate the latitude of peak electron energy flux. Then using this location as the final reference latitude, all other data was binned relative to it. The next line of attack was to bin according to where the latitudinal profile center of gravity lay. This procedure overcomes the problem of very narrow intense features controlling the binning. Again the same data sets were analyzed and a binning made. These profiles look better and indeed look similar to the original auroral model. The passes were then binned relative to the equatorward boundary. The problems associated with the equatorward edge were resolved.

Bowline, M.↗
NASA NTRS1997

Adaptation of a Fast Optimal Interpolation Algorithm to the Mapping of Oceangraphic Data

A fast, recently developed, multiscale optimal interpolation algorithm has been adapted to the mapping of hydrographic and other oceanographic data. This algorithm produces solution and error estimates which are consistent with those obtained from exact least squares methods, but at a small fraction of the computational cost. Problems whose solution would be completely impractical using exact least squares, that is, problems with tens or hundreds of thousands of measurements and estimation grid points, can easily be solved on a small workstation using the multiscale algorithm. In contrast to methods previously proposed for solving large least squares problems, our approach provides estimation error statistics while permitting long-range correlations, using all measurements, and permitting arbitrary measurement locations. The multiscale algorithm itself, published elsewhere, is not the focus of this paper. However, the algorithm requires statistical models having a very particular multiscale structure; it is the development of a class of multiscale statistical models, appropriate for oceanographic mapping problems, with which we concern ourselves in this paper. The approach is illustrated by mapping temperature in the northeastern Pacific. The number of hydrographic stations is kept deliberately small to show that multiscale and exact least squares results are comparable. A portion of the data were not used in the analysis; these data serve to test the multiscale estimates. A major advantage of the present approach is the ability to repeat the estimation procedure a large number of times for sensitivity studies, parameter estimation, and model testing. We have made available by anonymous Ftp a set of MATLAB-callable routines which implement the multiscale algorithm and the statistical models developed in this paper.

Menemenlis, Dimitris↗
DOE OSTI2023

Cross-section measurement of the 82 Kr ($p,γ$) 83 Rb reaction in inverse kinematics

The total cross section of the 82 Kr(p,γ) 83 Rb reaction was measured for the first time at effective center-of-mass energies between 2.4 and 3.0 MeV, within the relevant Gamow window for the astrophysical γ process. The experiment took place at the National Superconducting Cyclotron Laboratory at Michigan State University using the ReA facility. A 82 Kr beam was directed onto a hydrogen gas cell located at the center of the Summing NaI(Tl) (SuN) detector. Here, the obtained spectra were analyzed using the γ-summing technique and the extracted cross section was compared to standard statistical model calculations using the non-smoker and talys codes. The comparison indicates that standard statistical model calculations tend to overproduce the cross section of the 82 Kr(p,γ) 83 Rb reaction relative to the experimentally measured values. Furthermore, the experimental data were used to provide additional constraints on the nuclear level density and the γ-ray strength function used in the statistical model calculations.

59 ≤ A ≤ 89↗
DOE OSTI2023

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗
DOE OSTI2024

Statistical mechanical model for crack growth

Analytic relations that describe crack growth are vital for modeling experiments and building a theoretical understanding of fracture. Upon constructing an idealized model system for the crack and applying the principles of statistical thermodynamics, it is possible to formulate the rate of thermally activated crack growth as a function of load, but the result is analytically intractable. In this report an asymptotically correct theory is used to obtain analytic approximations of the crack growth rate from the fundamental theoretical formulation. These crack growth rate relations are compared to those that exist in the literature and are validated with respect to Monte Carlo calculations and experiments. The success of this approach is encouraging for future modeling endeavors that might consider more complicated fracture mechanisms, such as inhomogeneity or a reactive environment.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Comparison of transform coding methods with an optimal predictor for the data compression of digital elevation models

Statistical encoding techniques enable the reduction of the number of bits required to encode a set of symbols, and are derived from their probabilities. Huffman encoding is an example of statistical encoding that has been used for error-free data compression. The degree of compression given by Huffman encoding in this application can be improved by the use of prediction methods. These replace the set of elevations by a set of corrections that have a more advantageous probability distribution. In particular, the method of Lagrange Multipliers for minimization of the mean square error has been applied to local geometrical predictors. Using this technique, an 8-point predictor achieved about a 7 percent improvement over an existing simple triangular predictor.

Lewis, Michael↗
DOE OSTI2020

A Performance Comparison of Low-Cost Near-Infrared (NIR) Spectrometers to a Conventional Laboratory Spectrometer for Rapid Biomass Compositional Analysis

The performance of a conventional laboratory near-infrared (NIR) spectrometer and two NIR spectrometer prototypes (a Texas Instruments NIRSCAN Nano evaluation model (EVM) and an InnoSpectra NIR-M-R2 spectrometer) are compared by collecting reflectance spectra of 270 well-characterized herbaceous biomass samples, building calibration models using the partial least squares (PLS-2) algorithm to predict five constituents of the samples from the reflectance spectra, and comparing the resulting model statistics. The prediction models developed using spectra from the Foss XDS spectrometer were slightly better than the prediction models developed using spectra from either the TI NIRSCAN Nano EVM and the InnoSpectra NIR-M-R2 as measured by the root mean square error (RMSECV) and the correlation coefficient (R 2 _cv) for “leave-one-out” cross-validation (CV). The models built from the two prototype units were not statistically significantly different from each other (p = 0.05). The Foss spectrometer has a larger wavelength range (400–2500 nm) compared with the two prototypes (900–1700 nm). When the spectra from the Foss XDS spectrometer were truncated so their wavelength range matched the wavelength range of the two prototype units, the resulting model was not statistically significantly different from the models from either prototype.

47 OTHER INSTRUMENTATION↗
NASA NTRS1980

Statistical wind profile gust model

A statistical wind profile gust model for the Space Transportation Operations and Trade Studies is developed by using 1800 Jimsphere wind profile data collected at Cape Kennedy during 1965 to 1972. Wind profiles from the surface to 20 km in component form, i.e., zonal and meridional are processed through the digital filters of different wave length ranges bases on the Martin-Graham cosine rolloff model. The residuals obtained from the filtering processes for the data base for the statistical analysis. For each wind component the gust and gust length at a specified reference altitude in a residual profile are defined. A two parameter gamma probability marginal distribution seems to fit the component gust amplitude and the gust length when redefined. The problem of finding an appropriate bivariate joint distribution of the gust amplitude and length remains to be solved. The probability distribution of the modulus of the gust amplitudes was derived under the assumption that they are independently distributed as gamma variates.

Doss, D. C.↗
NASA NTRS1977

Chain Pooling modeling selection as developed for the statistical analysis of a rotor burst protection experiment

As many as three iterated statistical model deletion procedures were considered for an experiment. Population model coefficients were chosen to simulate a saturated 2 to the 4th power experiment having an unfavorable distribution of parameter values. Using random number studies, three model selection strategies were developed, namely, (1) a strategy to be used in anticipation of large coefficients of variation, approximately 65 percent, (2) a strategy to be sued in anticipation of small coefficients of variation, 4 percent or less, and (3) a security regret strategy to be used in the absence of such prior knowledge.

Holms, A. G.↗
NASA NTRS1987

Development of a two-dimensional zonally averaged statistical-dynamical model. II - The role of eddy momentum fluxes in the general circulation and their parameterization

The role of eddy momentum fluxes in the general circulation was investigated using a two-dimensional zonally averaged statistical-dynamical model described by Yao and Stone (1987), which is almost two orders of magnitude faster than the three-dimensional climate model of Hansen et al. (1983). Results show that the vertical structure of the meridional eddy flux has relatively little impact on the general circulation, presumably because the vertical structure is strongly constrained by the thermal wind relation and surface friction. On the other hand, it was found that, in order to simulate accurately the general circulation and its response to climate changes, parameterization of the vertically integrated meridional eddy flux of angular momentum is necessary. A new parameterization of this eddy momentum transport was carried out, which is intended to represent the transport due to large-scale transient eddies arising from baroclinic instability.

Stone, Peter H.↗
DOE OSTI2021

Data Science Meets Physical Organic Chemistry

At the heart of synthetic chemistry is the holy grail of predictable catalyst design. In particular, researchers involved in reaction development in asymmetric catalysis have pursued a variety of strategies toward this goal. This is driven by both the pragmatic need to achieve high selectivities and the inability to readily identify why a certain catalyst is effective for a given reaction. While empiricism and intuition have dominated the field of asymmetric catalysis since its inception, enantioselectivity offers a mechanistically rich platform to interrogate catalyst-structure response patterns that explain the performance of a particular catalyst or substrate. In the early stages of an asymmetric reaction development campaign, the overarching mechanism of the reaction, catalyst speciation, the turnover limiting step, and many other details are unknown or posited based on related reactions. Considering the unclear details leading to a successful reaction, initial enantioselectivity data are often used to intuitively guide the ultimate direction of optimization. However, if the conditions of the Curtin-Hammett principle are satisfied, then measured enantioselectivity can be directly connected to the ensemble of diastereomeric transition states (TSs) that lead to the enantiomeric products, and the associated free energy difference between competing TSs (ΔΔ G ‡ = - RT ln[( S )/( R )], where ( S ) and ( R ) represent the concentrations of the enantiomeric products). We, and others, speculated that this important piece of information can be leveraged to guide reaction optimization in a quantitative way. Although traditional linear free energy relationships (LFERs), such as Hammett plots, have been used to illuminate important mechanistic features, we sought to develop data science derived tools to expand the power of LFERs in order to describe complex reactions frequently encountered in modern asymmetric catalysis. Specifically, we investigated whether enantioselectivity data from a reaction can be quantitatively connected to the attributes of reaction components, such as catalyst and substrate structural features, to harness data for asymmetric catalyst design. In this context, we developed a workflow to relate computationally derived features of reaction components to enantioselectivity using data science tools. The mathematical representation of molecules can incorporate many aspects of a transformation, such as molecular features from substrate, product, catalyst, and proposed transition states. Statistical models relating these features to reaction outputs can be used for various tasks, such as performance prediction of untested molecules. Perhaps most importantly, statistical models can guide the generation of mechanistic hypotheses that are embedded within complex patterns of reaction responses. Overall, merging traditional physical organic experiments with statistical modeling techniques creates a feedback loop that enables both evaluation of multiple mechanistic hypotheses and future catalyst design. In this Account, we highlight the evolution and application of this approach in the context of a collaborative program based on chiral phosphoric acid catalysts (CPAs) in asymmetric catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Modeling the dosimetry of organ-at-risk in head and neck IMRT planning: An intertechnique and interinstitutional study

Purpose: To build a statistical model to quantitatively correlate the anatomic features of structures and the corresponding dose-volume histogram (DVH) of head and neck (HN) Tomotherapy (Tomo) plans. To study if the model built upon one intensity modulated radiation therapy (IMRT) technique (such as conventional Linac) can be used to predict anticipated organs-at-risk (OAR) DVH of patients treated with a different IMRT technique (such as Tomo). To study if the model built upon the clinical experience of one institution can be used to aid IMRT planning for another institution. Methods: Forty-four Tomotherapy intensity modulate radiotherapy plans of HN cases (Tomo-IMRT) from Institution A were included in the study. A different patient group of 53 HN fixed gantry IMRT (FG-IMRT) plans was selected from Institution B. The analyzed OARs included the parotid, larynx, spinal cord, brainstem, and submandibular gland. Two major groups of anatomical features were considered: the volumetric information and the spatial information. The volume information includes the volume of target, OAR, and overlapped volume between target and OAR. The spatial information of OARs relative to PTVs was represented by the distance-to-target histogram (DTH). Important anatomical and dosimetric features were extracted from DTH and DVH by principal component analysis. Two regression models, one for Tomotherapy plan and one for IMRT plan, were built independently. The accuracy of intratreatment-modality model prediction was validated by a leave one out cross-validation method. The intertechnique and interinstitution validations were performed by using the FG-IMRT model to predict the OAR dosimetry of Tomo-IMRT plans. The dosimetry of OARs, under the same and different institutional preferences, was analyzed to examine the correlation between the model prediction and planning protocol. Results: Significant patient anatomical factors contributing to OAR dose sparing in HN Tomotherapy plans have been analyzed and identified. For all the OARs, the discrepancies of dose indices between the model predicted values and the actual plan values were within 2.1%. Similar results were obtained from the modeling of FG-IMRT plans. The parotid gland was spared in a comparable fashion during the treatment planning of two institutions. The model based on FG-IMRT plans was found to predict the median dose of the parotid of Tomotherapy plans quite well, with a mean error of 2.6%. Predictions from the FG-IMRT model suggested the median dose of the larynx, median dose of the brainstem and D2 of the brainstem could be reduced by 10.5%, 12.8%, and 20.4%, respectively, in the Tomo-IMRT plans. This was found to be correlated to the institutional differences in OAR constraint settings. Re-planning of six Tomotherapy patients confirmed the potential of optimization improvement predicted by the FG-IMRT model was correct. Conclusions: The authors established a mathematical model to correlate the anatomical features and dosimetric indexes of OARs of HN patients in Tomotherapy plans. The model can be used for the setup of patient-specific OAR dose sparing goals and quality control of planning results. The institutional clinical experience was incorporated into the model which allows the model from one institution to generate a reference plan for another institution, or another IMRT technique.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1987

Development of a two-dimensional zonally averaged statistical-dynamical model. I - The parameterization of moist convection and its role in the general circulation

The moist convection parameterization used in the GISS 3-D GCM is adapted for use in a two-dimensional (2-D) zonally averaged statistical-dynamical model. Experiments with different versions of the parameterization show that its impact on the general circulation in the 2-D model does not parallel its impact in the 3-D model unless the effect of zonal variations is parameterized in the moist convection calculations. A parameterization of the variations in moist static energy is introduced in which the temperature variations are calculated from baroclinic stability theory, and the relative humidity is assumed to be constant. Inclusion of the zonal variations of moist static energy in the 2-D moist convection parameterization allows just a fraction of a latitude circle to be unstable and enhances the amount of deep convection. This leads to a 2-D simulation of the general circulation very similar to that in the 3-D model. The experiments show that the general circulation is sensitive to the parameterized amount of deep convection in the subsident branch of the Hadley cell. The more there is, the weaker are the Hadley cell circulations and the westerly jets. The experiments also confirm the effects of momentum mixing associated with moist convection found by earlier investigators and, in addition, show that the momentum mixing weakens the Ferrel cell. An experiment in which the moist convection was removed while the hydrological cycle was retained and the eddy forcing was held fixed shows that moist convection by itself stabilizes the tropics, reduces the Hadley circulation, and reduces the maximum speeds in the westerly jets.

Yao, Mao-Sung↗
DOE OSTI2017

Statistical mechanical model of gas adsorption in porous crystals with dynamic moieties

Some nanoporous, crystalline materials possess dynamic constituents, for example, rotatable moieties. These moieties can undergo a conformation change in response to the adsorption of guest molecules, which qualitatively impacts adsorption behavior. Here, we pose and solve a statistical mechanical model of gas adsorption in a porous crystal whose cages share a common ligand that can adopt two distinct rotational conformations. Guest molecules incentivize the ligands to adopt a different rotational configuration than maintained in the empty host. Our model captures inflections, steps, and hysteresis that can arise in the adsorption isotherm as a signature of the rotating ligands. The insights disclosed by our simple model contribute a more intimate understanding of the response and consequence of rotating ligands integrated into porous materials to harness them for gas storage and separations, chemical sensing, drug delivery, catalysis, and nanoscale devices. Particularly, our model reveals design strategies to exploit these moving constituents and engineer improved adsorbents with intrinsic thermal management for pressure-swing adsorption processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Extinction Dynamics of the Foot-and-Mouth Disease Virus Carrier State Under Natural Conditions

Foot-and-mouth disease (FMD) is one of the most economically important livestock diseases worldwide. Following the clinical phase of FMD, a large proportion of ruminants remain persistently infected for extended periods. Although extinction of this carrier state occurs continuously at the animal and population levels, studies vary widely in their estimates of the duration of persistent infection. There is a need for robust statistical models to capture the dynamics of persistent infection for the sake of guiding FMD control and trade policies. The goal of the current study was to develop and assess statistical models to describe the extinction of FMD virus (FMDV) persistent infection using data from primary longitudinal studies of naturally infected cattle and Asian buffalo in Vietnam and India. Specifically, accelerated failure time (AFT) models and generalized linear mixed models (GLMM) were developed to predict the probability of persistent infection in seropositive animals and identified carriers at the individual animal level at sequential time points after outbreaks. The primary studies were analyzed by country and combined using an individual-participant data meta-analysis approach. The models estimated similar trends in the duration of persistent infection for the study/species groups included in the analyses, however the significance of the trends differed between the models. The overall probabilities of persistent infection were similar as predicted by the AFT and GLMM models: 6 months: 99% (AFT) /80% (GLMM), 12 months: 51% (AFT) /32% (GLMM), 18 months: 6% (AFT) /5% (GLMM), 24 months: 0.8% (AFT) /0.6% (GLMM). These models utilizing diverse and robust data sets predict higher probabilities of persistence than previously published, suggesting greater endurance of carriers subsequent to an outbreak. This study demonstrates the utility of statistical models to investigate the dynamics of persistent infection and the importance of large datasets, which can be achieved by combining data from several smaller studies in meta-analyses. Results of this study enhance current knowledge of the FMDV carrier state and may inform policy decisions regarding FMDV persistent infection.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1975

UV stars and the interstellar medium - A statistical time-dependent model

A statistical time-dependent model for the effects on the interstellar gas of the ionizing radiation from UV stars is presented. The radiation from the stars is assumed to be emitted as monochromatic bursts at random times and locations. The thermal and ionization history of the gas is followed in a Monte Carlo simulation. Probability densities as well as mean values for the ionization fraction, temperature, and other observable quantities are derived. The fluctuations in observed quantities due to the stochastic nature of the model are estimated. Models of the interstellar medium perpendicular to the galactic plane are developed and compared with observations. The action of the UV stars can explain the bulk properties of the interstellar gas toward the galactic poles and is consistent with the Copernicus satellite observations toward lambda Scorpii.

Lyon, J.↗