Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPUTER MODELING”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

An Overview of the NASA Lift+Cruise eVTOL Crash Test

Introduction – NASA RVLT Project Impact Dynamics / Crash Safety Task - Task Objective: “To improve the crashworthiness and impact safety of Urban Air Mobility (UAM) vehicle and provide data to simplify the certification process. Efforts will include development of validated computational models of these vehicles, as well as other impacting bodies such as birds and drones. Efforts will also focus on developing and evaluating energy absorbing and crush properties of emerging and non-traditional composite materials and processes. Finally, occupant protection will be addressed using computational models and physical assets as it pertains to all rotorcraft environments.” - Problem Statement: “There currently is a lack of data for requirements regarding the crashworthy performance of UAM vehicles and impact loads generated by a bird strike. To address this technology gap, NASA will develop test guidelines, adopt modeling methodologies demonstrating capability for ‘certification by analysis’, acquire vehicle and occupant data on full-scale representative vehicles, and provide data/guidance to consensus standards organizations and the UAM community.” - 4 Main focus points - The investigation of occupant injury using physical and computational assets - The development of energy absorbing technology - The generation of data from sub- and full-scale crash test data - The execution of advanced finite element modelling techniques

evtol↗
NASA NTRS1999

Damage Tolerance of Large Shell Structures

Progressive damage and fracture of large shell structures is investigated. A computer model is used for the assessment of structural response, progressive fracture resistance, and defect/damage tolerance characteristics. Critical locations of a stiffened conical shell segment are identified. Defective and defect-free computer models are simulated to evaluate structural damage/defect tolerance. Safe pressurization levels are assessed for the retention of structural integrity at the presence of damage/ defects. Damage initiation, growth, accumulation, and propagation to fracture are included in the simulations. Damage propagation and burst pressures for defective and defect-free shells are compared to evaluate damage tolerance. Design implications with regard to defect and damage tolerance of a large steel pressure vessel are examined.

Minnetyan, L.↗
DOE OSTI2016

Modeling and Qualification of a Modified Emission Unit for Radioactive Air Emissions Stack Sampling Compliance

A planned laboratory space and exhaust system modification to the Pacific Northwest National Laboratory Material Science and Technology Building indicated a new evaluation of the mixing at the air sampling system location would be required for compliance to ANSI/HPS N13.1-2011. The modified exhaust system would add a third fan thereby increasing the overall exhaust rate out the stack thus voiding the previous mixing study. Prior to modifying the radioactive air emissions exhaust system, a three-dimensional computational fluid dynamics computer model was used to evaluate the mixing at the sampling system location. Modeling of the new original three-fan system indicated that not all mixing criteria could be met. A second modeling effort was conducted with the addition of an air blender downstream of the confluence of the three fans which then showed satisfactory mixing results. The final installation included an air blender, and the exhaust system underwent full-scale tests to verify velocity, cyclonic flow, gas, and particulate uniformity. The modeling results and those of the full-scale tests show agreement between each of the evaluated criteria. The use of a computational fluid dynamics code was an effective aid in the design process and allowed the sampling system to remain in its original location while still meeting the requirements for sampling at a well-mixed location.

atmospheric emissions, air monitoring, environment↗
DOE OSTI2025

Computational Fluid Dynamics Modeling to Facilitate Qualification of Stack Sampling Probe Location

Computational fluid dynamics (CFD) modeling was used to help evaluate modifications to a radiological effluent stack and assist with establishing a stack sampling location that met the mixing criteria for qualification. Requirements for stack sampling location are listed in the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-2021 standard. Modeling was performed to help develop a suitable design for increasing building ventilation for the radiological effluent stack. The ANSI/HPS N13.1-2021 criteria for the air monitoring probe location are that the coefficient of variation of velocity uniformity, gaseous tracer uniformity, and particulate tracer uniformity must be less than or equal to 20%. Furthermore, no point in the sampling location may have a gaseous tracer concentration that varies from the mean concentration by more than 30%. Additionally, the flow angle at the sampling location must not be more than 20 degrees. The ANSI/HPS N13.1-2021 standard allows for models (physical or computational) to be employed to perform the full suite of qualification tests, followed by a more limited set of verification tests on the actual stack to qualify the stack sampling location. Here, a series of computational model simulations were employed to evaluate the stack qualification criteria. Significant time and re-source savings are achieved using CFD modeling. CFD modeling demonstrated that the stack meets the criteria at the sample probe location. Verification tests were performed on the modified stack to measure the velocity uniformity and flow angle at the stack sampling location, and results demonstrated that the CFD model results may be used to support the qualification of the stack sampling location.

Air Monitoring↗
DOE OSTI2020

Promoting the cleavage of C-O bonds at the interface between a metal oxide cluster and a Co(0001) support

As a first step toward the rational design of Co-based catalysts with a higher activity and selectivity, we determine how one can activate a C-O bond at the interface between a metal oxide cluster and a Co(0001) support. The hypothesis here is that the metal ions in metal oxide clusters on a Co(0001) support enhance the adsorption of CO and weaken the C-O bond strength, which can then facilitate the dissociation of the CO reactant. To test this hypothesis, we developed three computational models of Ti4O8/Co(0001), Zr4O8/Co(0001), and Mn8O8/Co(0001). We quantify the CO adsorption behavior at the interface sites between an oxide cluster and the Co(0001) support as well as the corresponding IR spectra. We correlate the computed CO stretch frequencies with their CO adsorption energies, as well as the CO stretch frequency with the C-O bond length, and related these findings to the changes in the chemical bonding in the bound CO. The interface is the most favorable site for CO adsorption. Adsorption results in an increase of the C-O bond length and a decrease in its vibrational frequency. From a chemical bonding analysis, the bond order in CO at this site drops from 3 (in the gas phase) to 1. This decrease in bond order is a necessary precursor stage for CO dissociation. The experimental measurements of the corresponding FTIR spectra support this point. The favorability of CO adsorption at the interface sites is due to an electron transfer from the metal ion in the metal oxide cluster to the O atom in CO. We establish a linear relationship between C-O bond length and CO frequency and this relationship is found to be independent of the support, type of metal oxide cluster, or the adsorption site.

Fischer-Tropsch Synthesis, Heterogeneous Catalysis↗
DOE OSTI2016

Hot Rolling Scrap Reduction through Edge Cracking and Surface Defects Control

The design of future aircraft must address the combined demands for fuel efficiency, reduced emissions and lower operating costs. One contribution to these goals is weight savings through the development of new alloys and design techniques for airframe structures. This research contributes to the light-weighting through fabrication of monolithic components from advanced aluminum alloys by making a link between alloy processing history and in-service performance. Specifically, this research demonstrates the link between growing cracks with features of the alloy microstructure that follow from thermo-mechanical processing. This is achieved through a computer model of crack deviation. The model is validated against experimental data from production scale aluminum alloy plate, and demonstration of the effect of changes in processing history on crack growth is made. The model is cast in the open-source finite element code WARP3D, which is freely downloadable and well documented. This project provides benefit along several avenues. First, the technical contribution of the computer model offers the materials engineer a critical means of providing guidance both upstream, to process tuning to achieve optimal properties, and downstream, to enhance fault tolerance. Beyond the fuel savings and emissions reduction inherent in the light-weighting of aircraft structures, improved fault tolerance addresses demands for longer inspection intervals over baseline, and a lower life cycle cost.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRSDate not supplied

Testing of ACC-6 Carbon-Carbon at High Strain Rates and High Temperatures

Some materials, such as steel, titanium, and tungsten, have a yield stress that depends on the strain rate, and actually are stronger at high rates of load or deformation. The increase in strength is relevant; and if not taken into account, the computer models will not match the experimental behavior observed in typical shock loading, ordnance, or ballistic tests. For higher velocity (hypervelocity impacts), an equation of state that accounts for non-linearities and temperature in the pressure-volume relation, may also be needed. This rate-dependent behavior of the strength has also been observed in some composite materials. Southwest Research Institute® (SwRI®) has started an internal research effort to measure and better model the behavior of materials, in particular carbon/carbon composites like ACC-6, at ultra-high temperatures for both static and high rates of deformation. The ultimate objective is to obtain reliable computer models of carbon/carbon material under hypersonic flight conditions, hypervelocity impacts, shock-loading situations such as booster-vehicle separation, and munition operation events. SwRI® has successfully modeled in the past ballistic impact on other composites, namely S-2 glass/epoxy, Kevlar/resin, polyethylene fibers/resin, and carbon fiber composites. This presentation will show test results relevant to the strength properties of the ACC-6 C/C material which were measured at quasistatic, medium, and the moderately high strain rates (~103 s—1) achievable in a conventional Split Hopkinson bar apparatus. Some of the tests were performed at both high-rates and high temperatures typical of hypersonic flight, something that, to the authors’ knowledge, has never been published in the past for this particular material. The test method will be presented, as well as the comparison of the stress-strain curves for the different temperatures and strain-rates of interest. A discussion of how the results can be used in the typical hydrocodes of interest will be included.

carbon/carbon↗
NASA NTRS1979

Development of a model and computer code to describe solar grade silicon production processes

Two computer codes were developed for describing flow reactors in which high purity, solar grade silicon is produced via reduction of gaseous silicon halides. The first is the CHEMPART code, an axisymmetric, marching code which treats two phase flows with models describing detailed gas-phase chemical kinetics, particle formation, and particle growth. It can be used to described flow reactors in which reactants, mix, react, and form a particulate phase. Detailed radial gas-phase composition, temperature, velocity, and particle size distribution profiles are computed. Also, deposition of heat, momentum, and mass (either particulate or vapor) on reactor walls is described. The second code is a modified version of the GENMIX boundary layer code which is used to compute rates of heat, momentum, and mass transfer to the reactor walls. This code lacks the detailed chemical kinetics and particle handling features of the CHEMPART code but has the virtue of running much more rapidly than CHEMPART, while treating the phenomena occurring in the boundary layer in more detail.

Gould, R. K.↗
DOE OSTI2014

Confirmation of a realistic reactor model for BNCT dosimetry at the TRIGA Mainz

Purpose: In order to build up a reliable dose monitoring system for boron neutron capture therapy (BNCT) applications at the TRIGA reactor in Mainz, a computer model for the entire reactor was established, simulating the radiation field by means of the Monte Carlo method. The impact of different source definition techniques was compared and the model was validated by experimental fluence and dose determinations. Methods: The depletion calculation code ORIGEN2 was used to compute the burn-up and relevant material composition of each burned fuel element from the day of first reactor operation to its current core. The material composition of the current core was used in a MCNP5 model of the initial core developed earlier. To perform calculations for the region outside the reactor core, the model was expanded to include the thermal column and compared with the previously established ATTILA model. Subsequently, the computational model is simplified in order to reduce the calculation time. Both simulation models are validated by experiments with different setups using alanine dosimetry and gold activation measurements with two different types of phantoms. Results: The MCNP5 simulated neutron spectrum and source strength are found to be in good agreement with the previous ATTILA model whereas the photon production is much lower. Both MCNP5 simulation models predict all experimental dose values with an accuracy of about 5%. The simulations reveal that a Teflon environment favorably reduces the gamma dose component as compared to a polymethyl methacrylate phantom. Conclusions: A computer model for BNCT dosimetry was established, allowing the prediction of dosimetric quantities without further calibration and within a reasonable computation time for clinical applications. The good agreement between the MCNP5 simulations and experiments demonstrates that the ATTILA model overestimates the gamma dose contribution. The detailed model can be used for the planning of structural modifications in the thermal column irradiation channel or the use of different irradiation sites than the thermal column, e.g., the beam tubes.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

In Silico Screen and Structural Analysis Identifies Bacterial Kinase Inhibitors which Act with β-Lactams To Inhibit Mycobacterial Growth

New tools and concepts are needed to combat antimicrobial resistance. Actinomycetes and firmicutes share several eukaryotic-like Ser/Thr kinases (eSTK) that offer antibiotic development opportunities, including PknB, an essential mycobacterial eSTK. Despite successful development of potent biochemical PknB inhibitors by many groups, clinically useful microbiologic activity has been elusive. Additionally, PknB kinetics are not fully described, nor are structures with specific inhibitors available to inform inhibitor design. We used computational modeling with available structural information to identify human kinase inhibitors predicted to bind PknB, and we selected hits based on drug-like characteristics intended to increase the likelihood of cell entry. The computational model suggested a family of inhibitors, the imidazopyridine aminofurazans (IPAs), bind PknB with high affinity. We performed an in-depth characterization of PknB and found that these inhibitors biochemically inhibit PknB, with potency roughly following the predicted models. A novel X-ray structure confirmed that the inhibitors bound as predicted and made favorable protein contacts with the target. These inhibitors also have antimicrobial activity toward mycobacteria and nocardia. Furthermore, we demonstrated that the inhibitors are uniquely potentiated by β-lactams but not antibiotics traditionally used to treat mycobacteria, consistent with PknB’s role in sensing cell wall stress. This is the first demonstration in the phylum actinobacteria that some β-lactam antibiotics could be more effective if paired with a PknB inhibitor. Collectively, our data show that in silico modeling can be used as a tool to discover promising drug leads, and the inhibitors we discovered can act with clinically relevant antibiotics to restore their efficacy against bacteria with limited treatment options.

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

Maximized Information Gain of Next Generation Pulsed Power Using Optimized Design of Z-Machine Experiments

This project develops a Bayesian optimization approach to extracting insights from Z Machine experimental data to determine if and how these insights can be used to extrapolate to a larger facility. The primary goal is to address the scientific challenge of informing how confidently experimental conditions can be predicted on a next generation facility, the design of which requires the reliable extrapolation of current high energy density technologies to regimes yet unobserved, except by costly high-fidelity computational models. Maximizing the use of presently available data and understanding how it informs future endeavors is critically important to enable transformative pulsed power and the science of extreme conditions. We explore a Bayesian optimization approach to experimental design which combines information theory, experimental data, and computational modeling to explore how information gain can be maximized.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

The spiral-compact galaxy pair AM 2208-251: Computer simulations versus observations

The system AM2208-251 is a roughly edge-on spiral extending east-west with a smaller round compact E system about 60 arcsec east of the spiral nucleus along the major axis of the spiral. Bertola, Huchtmeier, and Zeilinger (1990) have presented optical spectroscopic as well as single dish 21 cm observations of this system. Their spectroscopic data show, via emission lines lambda lambda 3727-29A, a rising rotation curve near the nucleus. These spectroscopic observations may indicate a tidal interaction in the system. In order to learn more about such pairs, the authors simulated the interaction using the computer model developed by Miller (1976 a,b, 1978) and modified by the authors (Byrd 1986, 1987, 1988). To do the simulation they need an idea of the mutual orbits of the two galaxies. Their computer model is a two-dimensional polar N-body program. It consists of a self-gravitating disk of particles, within an inert axially symmetric stabilizing halo potential. The particles are distributed in a 24(radial) by 36(azimuthal) polar grid. Self consistent calculations can be done only within the grid area. The disk is modeled with a finite Mestel disk, where all the particles initially move in circular orbits with constant tangential velocities (Mestel 1963), resulting in a flat rotation curve. The gas particles in the spiral's disk, which make up 30 percent of its mass, collide in the following manner. The number of particles in each bin of the polar grid is counted every time step. If it is greater than a given critical density, all the particles in the bin collide, obtaining in the result the same velocities, equal to the average for the bin. This process produces clumps of gas particles-the star formation sites. The authors suppress the collision in the inner part of the disk (within the circle r = 6) to represent the hole seen in the gas in the nuclear bulge of spirals. They thus avoid spurious effects due to collisions in that region. They also varied the size of the collisional bins, which did not affect their conclusions.

Klaric, Mario↗
NASA NTRS2000

Metabolic Adaptation to Muscle Ischemia

Although all tissues in the body can adapt to varying physiological/pathological conditions, muscle is the most adaptable. To understand the significance of cellular events and their role in controlling metabolic adaptations in complex physiological systems, it is necessary to link cellular and system levels by means of mechanistic computational models. The main objective of this work is to improve understanding of the regulation of energy metabolism during skeletal/cardiac muscle ischemia by combining in vivo experiments and quantitative models of metabolism. Our main focus is to investigate factors affecting lactate metabolism (e.g., NADH/NAD) and the inter-regulation between carbohydrate and fatty acid metabolism during a reduction in regional blood flow. A mechanistic mathematical model of energy metabolism has been developed to link cellular metabolic processes and their control mechanisms to tissue (skeletal muscle) and organ (heart) physiological responses. We applied this model to simulate the relationship between tissue oxygenation, redox state, and lactate metabolism in skeletal muscle. The model was validated using human data from published occlusion studies. Currently, we are investigating the difference in the responses to sudden vs. gradual onset ischemia in swine by combining in vivo experimental studies with computational models of myocardial energy metabolism during normal and ischemic conditions.

Cabrera, Marco E.↗
DOE OSTI2019

Calculating Dose Release Factors to Obtain the Worst-Case Maximally Exposed Location by Sector

The aim of this project was to update the SRS Dose Release Factors (DRF's) for the Maximally Exposed Individual (MEI) for 20 area release points onsite. To update, potential MEIs now include onsite locations and modeled stack heights were updated to match onsite operations. The DRF's and MEI locations are used to ensure DOE and EPA compliance is met for dose standards. The Environmental Protection Agency (EPA) requires the use of the computer model CAP88 to estimate doses for demonstrating compliance with 40 CFR 61, the National Emission Standard for Hazardous Air Pollutants (NESHAP). The location of the maximally exposed individual (MEI) by sector and distance is required for execution of CAP88. The worst-case sector is the sector and distance where an offsite individual would receive the highest dose due to a potential release from a designated release location at a facility. The worst-case sector analysis addresses potential ground-level and stack-height releases from 20 areas of Savannah River Site (SRS). Weather patterns were analyzed at each of these areas relative to sixteen compass points. The analysis makes use of recently updated distance data and meteorological data. The computer model, CAP88 is EPA's software system used by the Department of Energy (DOE) to calculate the receptor dose due to routine atmospheric releases via a gaussian plume at the Savannah River Site (SRS). The dose-release factors (DRFs) are used to estimate the offsite and onsite dose from facility-specific releases and represent the dose to a receptor exposed to 1 Ci of a specified radionuclide being released into the atmosphere. The DRFs are simply a ratio of the TED to the release rate and, therefore, are directly scalable. The DRF values listed in the Appendix can be applied to expected release values from each area to estimate the potential dose to an individual adult at each location. NESHAP DRFs have been supplied to Environmental Compliance and Area Closure Projects (EC and ACP) for many years and are periodically updated to include changes in Cap88. estimate the potential dose to an individual adult at each location DRF in mrem Ci{sup -1} are estimated by modeling the total effective dose (TED) in mrem yr{sup -1} assuming an annual atmospheric release rate of 1 Ci yr-1 of the associated radionuclides. TED = Effective Dose (ED) + Committed Effective Dose (CED). TED = DRF/RR. Approach and Methods: Determine possible MEI locations Onsite/Offsite including buildings and farms; Determine the distances from the release area to the possible MEI location and the direction from the release to the MEI location; Using the most recent meteorological data (2007-2011) and parameters, run CAP88 for each release area for each compass sector and previously calculated distances to MEI locations using a 1 Ci release rate for tritium; Use the calculated the air concentration values to determine the MEI for each sector from CAP88. The relationship between air concentration and received dose is directly proportional and can be used to determine which locations to calculate DRF's and a dose for an MEI; The DRFs are found for each compass sector from 20 SRS release areas using CAP88-PC Version 4.0; Each location is calculated with ground level releases evaluated for each area for a list of radionuclides, some with and some without progenies; Compare the DRFs of the potential onsite MEIs to offsite MEIs. Based on the 2018 MEI Report, the Site worst case sector MEI is located 593 mW of B-Area. This MEI location represents a point along the main road that travels through a tract of land leased to the SC National Guard. While this location is estimated to have the highest exposure, minimal radioactivity is released from B-Area. In 2017, SRS requested and was approved by the EPA to change the location for modeling site releases from the Center of Site to H-Area, and to use six stack heights. H-Area emission account for approximately 90% of all radioactive operations onsite. The MEI for H-Area was found to be 11,458 m NNE which locates the Savannah River Ecology Lab Conference Center. Comparisons for 2018 data and 2014 shows the new methods of calculation are more conservative yet still remain well below the DOE 100 mrem limit.

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

Modeling Interfacial Glass-Water Reactions: Recent Advances and Current Limitations

Describing the reactions that occur at the glass-water interface and control the development of the altered layer constitutes one of the main scientific challenges impeding existing models from providing accurate radionuclide release estimates. Radionuclide release estimates are a critical component of the safety basis for geologic repositories. The altered layer (i.e., amorphous hydrated surface layer and crystalline reaction products) represents a complex region, both physically and chemically, sandwiched between two distinct boundaries pristine glass surface at the inner most interface and aqueous solution at the outer most interface. Computational models, spanning different length and time-scales, are currently being developed to improve our understanding of this complex and dynamic process with the goal of accurately describing the pore-scale changes that occur as the system evolves. These modeling approaches include geochemical simulations [i.e., classical reaction path simulations and glass reactivity in allowance for alteration layer (GRAAL) simulations], Monte Carlo simulations, and Molecular Dynamics methods. Finally, in this manuscript, we discuss the advances and limitations of each modeling approach placed in the context of the glass-water reaction and how collectively these approaches provide insights into the mechanisms that control the formation and evolution of altered layers.

36 MATERIALS SCIENCE↗
NASA NTRS1991

A Scheduling Algorithm for Replicated Real-Time Tasks

We present an algorithm for scheduling real-time periodic tasks on a multiprocessor system under fault-tolerant requirement. Our approach incorporates both the redundancy and masking technique and the imprecise computation model. Since the tasks in hard real-time systems have stringent timing constraints, the redundancy and masking technique are more appropriate than the rollback techniques which usually require extra time for error recovery. The imprecise computation model provides flexible functionality by trading off the quality of the result produced by a task with the amount of processing time required to produce it. It therefore permits the performance of a real-time system to degrade gracefully. We evaluate the algorithm by stochastic analysis and Monte Carlo simulations. The results show that the algorithm is resilient under hardware failures.

Yu, Albert C.↗
NASA NTRS2018

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗
DOE OSTI2016

UQ Toolkit v. 3.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterization and propagation of uncertainties in computational models. This library provides Bayesian inference tools to infer uncertain model parameters, Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, methods for constructing Karhunen-Loève representations of stochastic processes, and global sensitivity analysis tools used to compute Sobol indices in order to characterize the importance of uncertain inputs or the interactions between them. The basis for many of these methods relies on representing random variables with Polynomial Chaos expansions (PCEs). This library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For the propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Sargsyan, Khachik↗