Search NASA⌕ Search

SEARCH · Search NASA

Results for “Randomized methods”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Delay-induced stochastic bifurcations in a bistable system under white noise

In this paper, the effects of noise and time delay on stochastic bifurcations are investigated theoretically and numerically in a time-delayed Duffing-Van der Pol oscillator subjected to white noise. Due to the time delay, the random response is not Markovian. Thereby, approximate methods have been adopted to obtain the Fokker-Planck-Kolmogorov equation and the stationary probability density function for amplitude of the response. Based on the knowledge that stochastic bifurcation is characterized by the qualitative properties of the steady-state probability distribution, it is found that time delay and feedback intensity as well as noise intensity will induce the appearance of stochastic P-bifurcation. Besides, results demonstrated that the effects of the strength of the delayed displacement feedback on stochastic bifurcation are accompanied by the sensitive dependence on time delay. Furthermore, the results from numerical simulations best confirm the effectiveness of the theoretical analyses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Long history of Cl-36 Assessment of Graphite Waste by EDF Engineering and the Latest Suggested Developments - 17043

Prior to 2005, the EDF initial version for the radiological inventory of graphite was based on the maximum values of the measures, namely a very pessimistic way. In 2008, a scientific method was developed by EDF to evaluate the inventory by reverse activation calculation, in order to limit the overestimation. The radiological inventory computation principles can be compared to those of all engineering studies. First, the calculation is simplified by taking margins, and if the results are not satisfactory enough, some simplifications are suppressed to reduce these margins even if computation becomes more complicated. On EDF piles, the gain obtained in 2008 represented a factor 50 relative to the initial version of Cl-36, even with a very penalizing multiplicative factor for uncertainty. Today, the accurate calculation of the uncertainty no longer justifies such a factor. A gain of about 150 now occurs compared to the initial version of Cl-36. It is more than 2 orders of magnitude, which is already enough to assess again the necessary type of disposal because Cl-36 is one of the key radionuclides. After global presentation of graphite waste management in France, the paper focuses on the inventory scientific process of EDF nuclear graphite, the purity of graphite and its consequences, its particular sampling, the lessons learned, the initial method description, the link with the demonstration of the completely random point process of chlorine in nuclear graphite and the last suggested precise uncertainty computation of inventory of Cl-36. The method confirmation for concentration assessment, the history of the EDF method development and its benefits are also emphasized. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Comparison of Experimental Methods for Estimating Matrix Diffusion Coefficients for Contaminant Transport Modeling

Diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2018

Comparison of experimental methods for estimating matrix diffusion coefficients for contaminant transport modeling

Here, diffusion cell and diffusion wafer experiments were conducted to compare methods for estimating effective matrix diffusion coefficients in rock core samples from Pahute Mesa at the Nevada Nuclear Security Site (NNSS). A diffusion wafer method, in which a solute diffuses out of a rock matrix that is pre-saturated with water containing the solute, is presented as a simpler alternative to the traditional through-diffusion (diffusion cell) method. Both methods yielded estimates of effective matrix diffusion coefficients that were within the range of values previously reported for NNSS volcanic rocks. The difference between the estimates of the two methods ranged from 14 to 30%, and there was no systematic high or low bias of one method relative to the other. From a transport modeling perspective, these differences are relatively minor when one considers that other variables (e.g., fracture apertures, fracture spacings) influence matrix diffusion to a greater degree and tend to have greater uncertainty than effective matrix diffusion coefficients. For the same relative random errors in concentration measurements, the diffusion cell method yields effective matrix diffusion coefficient estimates that have less uncertainty than the wafer method. However, the wafer method is easier and less costly to implement and yields estimates more quickly, thus allowing a greater number of samples to be analyzed for the same cost and time. Given the relatively good agreement between the methods, and the lack of any apparent bias between the methods, the diffusion wafer method appears to offer advantages over the diffusion cell method if better statistical representation of a given set of rock samples is desired.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Numerical models for the prediction of failure for multilayer fusion Al-alloy sheets

Initiation and propagation of cracks in monolithic and multi-layer aluminum alloys, called “Fusion”, is investigated. 2D plane strain finite element simulations are performed to model deformation due to bending and to predict failure. For this purpose, fracture strains are measured based on microscopic pictures of Nakajima specimens. In addition to, micro-structure of materials is taken into account by introducing a random grain distribution over the sheet thickness as well as a random distribution of the measured yield curve. It is shown that the performed experiments and the introduced FE-Model are appropriate methods to highlight the advantages of the Fusion material, especially for bending processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystallographic variant mapping using precession electron diffraction data

In this work, we developed three methods to map crystallographic variants of samples at the nanoscale by analyzing precession electron diffraction data using a high-temperature shape memory alloy and a VO2 thin film on sapphire as the model systems. The three methods are (I) a user-selecting-reference pattern approach, (II) an algorithm-selecting-reference-pattern approach, and (III) a k-means approach. In the first two approaches, Euclidean distance, Cosine, and Structural Similarity (SSIM) algorithms were assessed for the diffraction pattern similarity quantification. We demonstrated that the Euclidean distance and SSIM methods outperform the Cosine algorithm. We further revealed that the random noise in the diffraction data can dramatically affect similarity quantification. Denoising processes could improve the crystallographic mapping quality. With the three methods mentioned above, we were able to map the crystallographic variants in different materials systems, thus enabling fast variant number quantification and clear variant distribution visualization. The advantages and disadvantages of each approach are also discussed. We expect these methods to benefit researchers who work on martensitic materials, in which the variant information is critical to understand their properties and functionalities.

Crystallographic variant mapping↗
DOE OSTI2021

Development of a Random-Forest Cloud-Regime Classification Model Based on Surface Radiation and Cloud Products

Various methods have been developed to characterize cloud type, otherwise referred to as cloud regime. These include manual sky observations, combining radiative and cloud vertical properties observed from satellite, surface-based remote sensing, and digital processing of sky imagers. While each method has inherent advantages and disadvantages, none of these cloud-typing methods actually includes measurements of surface shortwave or longwave radiative fluxes. Here, a method that relies upon detailed, surface-based radiation and cloud measurements and derived data products to train a random-forest machine-learning cloud classification model is introduced. Measurements from five years of data from the ARM Southern Great Plains site were compiled to train and independently evaluate the model classification performance. A cloud-type accuracy of approximately 80% using the random-forest classifier reveals that the model is well suited to predict climatological cloud properties. Furthermore, an analysis of the cloud-type misclassifications is performed. While physical cloud types may be misreported, the shortwave radiative signatures are similar between misclassified cloud types. From this, we assert that the cloud-regime model has the capacity to successfully differentiate clouds with comparable cloud–radiative interactions. Therefore, we conclude that the model can provide useful cloud-property information for fundamental cloud studies, inform renewable energy studies, and be a tool for numerical model evaluation and parameterization improvement, among many other applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Stochastic many-body perturbation theory for anharmonic molecular vibrations

A new quantum Monte Carlo (QMC) method for anharmonic vibrational zero-point energies and transition frequencies is developed, which combines the diagrammatic vibrational many-body perturbation theory based on the Dyson equation with Monte Carlo integration. The infinite sums of the diagrammatic and thus size-consistent first- and second-order anharmonic corrections to the energy and self-energy are expressed as sums of a few m- or 2m-dimensional integrals of wave functions and a potential energy surface (PES) (m is the vibrational degrees of freedom). Each of these integrals is computed as the integrand (including the value of the PES) divided by the value of a judiciously chosen weight function evaluated on demand at geometries distributed randomly but according to the weight function via the Metropolis algorithm. In this way, the method completely avoids cumbersome evaluation and storage of high-order force constants necessary in the original formulation of the vibrational perturbation theory; it furthermore allows even higher-order force constants essentially up to an infinite order to be taken into account in a scalable, memory-efficient algorithm. The diagrammatic contributions to the frequency-dependent self-energies that are stochastically evaluated at discrete frequencies can be reliably interpolated, allowing the self-consistent solutions to the Dyson equation to be obtained. This method, therefore, can compute directly and stochastically the transition frequencies of fundamentals and overtones as well as their relative intensities as pole strengths, without fixed-node errors that plague some QMC. It is shown that, for an identical PES, the new method reproduces the correct deterministic values of the energies and frequencies within a few cm{sup −1} and pole strengths within a few thousandths. With the values of a PES evaluated on the fly at random geometries, the new method captures a noticeably greater proportion of anharmonic effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Uniform Subspace Correction Preconditioners for Discontinuous Galerkin Methods with hp-Refinement

In this paper, we develop subspace correction preconditioners for discontinuous Galerkin (DG) discretizations of elliptic problems with hp-refinement. These preconditioners are based on the decomposition of the DG finite element space into a conforming subspace, and a set of small nonconforming edge spaces. The conforming subspace is preconditioned using a matrix-free low-order refined technique, which in this work, we extend to the hp-refinement context using a variational restriction approach. The condition number of the resulting linear system is independent of the granularity of the mesh h, and the degree of the polynomial approximation p. The method is amenable to use with meshes of any degree of irregularity and arbitrary distribution of polynomial degrees. Furthermore, numerical examples are shown on several test cases involving adaptively and randomly refined meshes, using both the symmetric interior penalty method and the second method of Bassi and Rebay (BR2).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Data from: Comparing traditional and Bayesian approaches to ecological meta-analysis

Despite the wide application of meta-analysis in ecology, some of the traditional methods used for meta-analysis may not perform well given the type of data characteristic of ecological meta-analyses. We reviewed published meta-analyses on the ecological impacts of global climate change, evaluating the number of replicates used in the primary studies (ni) and the number of studies or records (k) that were aggregated to calculate a mean effect size. We used the results of the review in a simulation experiment to assess the performance of conventional frequentist and Bayesian meta-analysis methods for estimating a mean effect size and its uncertainty interval. Our literature review showed that ni and k were highly variable, distributions were right-skewed, and were generally small (median ni =5, median k=44). Our simulations show that the choice of method for calculating uncertainty intervals was critical for obtaining appropriate coverage (close to the nominal value of 0.95). When k was low (<40), 95% coverage was achieved by a confidence interval based on the t-distribution that uses an adjusted standard error (the Hartung-Knapp-Sidik-Jonkman, HKSJ), or by a Bayesian credible interval, whereas bootstrap or z-distribution confidence intervals had lower coverage. Despite the importance of the method to calculate the uncertainty interval, 39% of the meta-analyses reviewed did not report the method used, and of the 61% that did, 94% used a potentially problematic method, which may be a consequence of software defaults. In general, for a simple random-effects meta-analysis, the performance of the best frequentist and Bayesian methods were similar for the same combinations of factors (k and mean replication), though the Bayesian approaches had higher than nominal (>95%) coverage for the mean effect when k was very low (k<15). Our literature review suggests that many meta-analyses that used z-distribution or bootstrapping confidence intervals may have over-estimated the statistical significance of their results when the number of studies was low; more appropriate methods need to be adopted in ecological meta-analyses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Complex Langevin simulation of a random matrix model at nonzero chemical potential

In this study we test the complex Langevin algorithm for numerical simulations of a random matrix model of QCD with a first order phase transition to a phase of finite baryon density. We observe that a naive implementation of the algorithm leads to phase quenched results, which were also derived analytically in this article. We test several fixes for the convergence issues of the algorithm, in particular the method of gauge cooling, the shifted representation, the deformation technique and reweighted complex Langevin, but only the latter method reproduces the correct analytical results in the region where the quark mass is inside the domain of the eigenvalues. In order to shed more light on the issues of the methods we also apply them to a similar random matrix model with a milder sign problem and no phase transition, and in that case gauge cooling solves the convergence problems as was shown before in the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A hybrid calibration approach to Hertz-type contact parameters for discrete element models

This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. Here, the proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made.

Constitutive law↗
DOE OSTI2014

How accurate is the strongly orthogonal geminal theory in predicting excitation energies? Comparison of the extended random phase approximation and the linear response theory approaches

Performance of the antisymmetrized product of strongly orthogonal geminal (APSG) ansatz in describing ground states of molecules has been extensively explored in the recent years. Not much is known, however, about possibilities of obtaining excitation energies from methods that would rely on the APSG ansatz. In the paper we investigate the recently proposed extended random phase approximations, ERPA and ERPA2, that employ APSG reduced density matrices. We also propose a time-dependent linear response APSG method (TD-APSG). Its relation to the recently proposed phase including natural orbital theory is elucidated. The methods are applied to Li{sub 2}, BH, H{sub 2}O, and CH{sub 2}O molecules at equilibrium geometries and in the dissociating limits. It is shown that ERPA2 and TD-APSG perform better in describing double excitations than ERPA due to inclusion of the so-called diagonal double elements. Analysis of the potential energy curves of Li{sub 2}, BH, and H{sub 2}O reveals that ERPA2 and TD-APSG describe correctly excitation energies of dissociating molecules if orbitals involved in breaking bonds are involved. For single excitations of molecules at equilibrium geometries the accuracy of the APSG-based methods approaches that of the time-dependent Hartree-Fock method with the increase of the system size. A possibility of improving the accuracy of the TD-APSG method for single excitations by splitting the electron-electron interaction operator into the long- and short-range terms and employing density functionals to treat the latter is presented.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Xanthohumol Microbiome and Signature in Healthy Adults (the XMaS Trial): Safety and Tolerability Results of a Phase I Triple-Masked, Placebo-Controlled Clinical Trial

Scope Xanthohumol, a prenylflavonoid from hops, has been extensively studied preclinically but has undergone limited research in human subjects. A triple-masked, placebo-controlled phase I clinical trial was conducted to examine the safety and tolerability of xanthohumol. Methods and results Thirty healthy volunteers were randomized to 24?mg/day xanthohumol (99.8% pure) or placebo for eight weeks. Comprehensive metabolic panels, complete blood counts, body weight, vital signs, and health-related quality of life questionnaires were assessed every two weeks. Participants were interviewed for adverse events (AEs) throughout the trial. Thirteen of 14 (93%) and 14 of 16 (88%) participants completed the trial in Group A and B respectively. There were no withdrawals due to AEs. There were no clinically relevant, between-group differences in laboratory biomarkers, body weight, vital signs, or health-related quality of life. There were no severe or FDA-defined serious AEs, but non-serious AEs were documented in both Group A (n=42) and Group B (n=58). Conclusion Over an eight-week period, 24 mg daily xanthohumol was safe and well-tolerated by healthy adults.

Langley, Blake↗
DOE OSTI2013

Comparing Intra-Arterial Chemotherapy Combined With Intravesical Chemotherapy Versus Intravesical Chemotherapy Alone: A Randomised Prospective Pilot Study for T1G3 Bladder Transitional Cell Carcinoma After Bladder-Preserving Surgery

Purpose: To compare the efficacy of intra-arterial chemotherapy combined with intravesical chemotherapy versus intravesical chemotherapy alone for T1G3 bladder transitional cell carcinoma (BTCC) followed by bladder-preserving surgery. Materials and Methods: Sixty patients with T1G3 BTCC were randomly divided into two groups. After bladder-preserving surgery, 29 patients (age 30-80 years, 24 male and 5 female) received intra-arterial chemotherapy in combination with intravesical chemotherapy (group A), whereas 31 patients (age 29-83 years, 26 male and 5 female) were treated with intravesical chemotherapy alone (group B). Twenty-nine patients were treated with intra-arterial epirubicin (50 mg/m{sup 2}) + cisplatin (60 mg/m{sup 2}) chemotherapy 2-3 weeks after bladder-preserving surgery once every 4-6 weeks. All of the patients received the same intravesical chemotherapy: An immediate prophylactic was administered in the first 6 h. After that, therapy was administered one time per week for 8 weeks and then one time per month for 8 months. The instillation drug was epirubicin (50 mg/m{sup 2}) and lasted for 30-40 min each time. The end points were tumour recurrence (stage Ta, T1), tumour progression (to T2 or greater), and disease-specific survival. During median follow-up of 22 months, the overall survival rate, tumour-specific death rate, recurrence rate, progression rate, time to first recurrence, and adverse reactions were compared between groups. Results: The recurrence rates were 10.3 % (3 of 29) in group A and 45.2 % (14 of 31) in group B, and the progression rates were 0 % (0 of 29) in group A and 22.6 % (7 of 31) in group B. There was a significant difference between the two groups regarding recurrence (p = 0.004) and progression rates (p = 0.011). Median times to first recurrence in the two groups were 15 and 6.5 months, respectively. The overall survival rates were 96.6 and 87.1 %, and the tumour-specific death rates were 0 % (0 of 29) and 13.5 % (4 of 31) in groups A and B, respectively. During the intra-arterial chemotherapy cycle, although more than 50 % patients experienced some toxicities, most were minor and reversible [grade 1-2 (46.7 %) vs. grade 1-2 (6.9 %)]. Conclusion: These findings suggest that combining intra-arterial chemotherapy with intravesical chemotherapy could delay tumour recurrence and progression compared with intravesical chemotherapy alone and this type treatment is relatively safe.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Balloon-Expandable Biodegradable Stents Versus Self-Expandable Metallic Stents: A Comparison Study of Stent-Induced Tissue Hyperplasia in the Rat Urethra

Purpose: To compare the degrees of stent-induced tissue hyperplasia of balloon-expandable, biodegradable stents (BEBSs) with those of self-expandable metallic stents (SEMSs) in a rat urethral model. Materials and Methods: A total of 20 rats were randomized into two groups. The BEBS group (n = 10) received a poly-l-lactic acid (PLLA) biodegradable stent. The SEMS group (n = 10) received a nitinol bare stent. All rats were killed eight weeks after stent placement. The degree of stent-induced tissue hyperplasia was assessed by comparing the results of retrograde urethrography and histologic examination between the two groups. Results: Stent placement was technically successful in all rats. Two rats in the BEBS group were excluded due to procedure-related death. The mean luminal diameter of stented urethra on urethrograms was not significantly different at 4 and 8 weeks between the two groups. On histologic analysis, the percentage of granulation tissue area (p < 0.001) and the thickness of papillary projection (p < 0.001) were significantly higher in the BEBS group compared with the SEMS group. The inflammatory cell infiltration showed a clear tendency to significance (p = 0.050). There were no statistical differences in the number of epithelial layers and the thickness of submucosal fibrosis between the two groups. Conclusion: Formation of stent-induced tissue hyperplasia was significantly evident in the rat urethra with similar degrees between the BEBS and the SEMS. The BEBS was associated with a thicker papillary projection and larger granulation tissue area resulting from higher inflammation compared with the SEMS.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Finite deformation cohesive polygonal finite elements for modeling pervasive fracture

In this work, we introduce a framework for modeling dynamic fracture problems using cohesive polygonal finite elements. Random polygonal meshes provide a robust, efficient method for generating an unbiased network of fracture surfaces. Further, these meshes have more facets per element than standard triangle or quadrilateral meshes, providing more possible facets per element to insert cohesive surfaces. This property of polygonal meshes is advantageous for the modeling of pervasive fracture. We use both Wachspress and maximum entropy shape functions to form a finite element basis over the polygons. Fracture surfaces are captured through dynamically inserted cohesive zone elements at facets between the polygons in the mesh. Contact is enforced through a penalty method that is applied to both closed cohesive surfaces and general interpenetration of two polygonal elements. Finally, several numerical examples are presented that illustrate the capabilities of the method and demonstrate convergence of solutions.

97 MATHEMATICS AND COMPUTING↗