Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “collocation”

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 415 records · Page 23

Accurate data-driven surrogates of dynamical systems for forward propagation of uncertainty

Stochastic collocation (SC) is a well-known non-intrusive method of constructing surrogate models for uncertainty quantification. In dynamical systems, SC is especially suited for full-field uncertainty propagation that characterizes the distributions of the high-dimensional solution fields of a model with stochastic input parameters. However, due to the highly nonlinear nature of the parameter-to-solution map in even the simplest dynamical systems, the constructed SC surrogates are often inaccurate. Here, this work presents an alternative approach, where we apply the SC approximation over the dynamics of the model, rather than the solution. By combining the data-driven sparse identification of nonlinear dynamics framework with SC, we construct dynamics surrogates and integrate them through time to construct the surrogate solutions. We demonstrate that the SC-over-dynamics framework leads to smaller errors, both in terms of the approximated system trajectories as well as the model state distributions, when compared against full-field SC applied to the solutions directly. We present numerical evidence of this improvement using three test problems: a chaotic ordinary differential equation, and two partial differential equations from solid mechanics.

42 ENGINEERING↗

A computationally efficient ductile damage model accounting for nucleation and micro-inertia at high triaxialities

The computational formulation of a micro-mechanical material model for the dynamic failure of ductile metals is presented in this paper. The statistical nature of porosity initiation is accounted for by introducing an arbitrary probability density function which describes the pores nucleation pressures. Each micropore within the representative volume element is modeled as a thick spherical shell made of plastically incompressible material. The treatment of porosity by a distribution of thick-walled spheres also allows for the inclusion of micro-inertia effects under conditions of shock and dynamic loading. The second order ordinary differential equation governing the microscopic porosity evolution is solved with a robust implicit procedure. A new Chebyshev collocation method is employed to approximate the porosity distribution and remapping is used to optimize memory usage. The adaptive approximation of the porosity distribution leads to a reduction of computational time and memory usage of up to two orders of magnitude. Moreover, the proposed model affords consistent performance: changing the nucleation pressure probability density function and/or the applied strain rate does not reduce accuracy or computational efficiency of the material model. The numerical performance of the model and algorithms presented is tested against three problems for high density tantalum: single void, one-dimensional uniaxial strain, and two-dimensional plate impact. Here, the results using the integration and algorithmic advances suggest a significant improvement in computational efficiency and accuracy over previous treatments for dynamic loading conditions.

42 ENGINEERING↗

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗

A bi-fidelity method for the multiscale Boltzmann equation with random parameters

In this paper, we study the multiscale Boltzmann equation with multi-dimensional random parameters by a bi-fidelity stochastic collocation (SC) method developed in [52], [70], [71]. By choosing the compressible Euler system as the low-fidelity model, we adapt the bi-fidelity SC method to combine computational efficiency of the low-fidelity model with high accuracy of the high-fidelity (Boltzmann) model. With only a small number of high-fidelity asymptotic-preserving solver runs for the Boltzmann equation, the bi-fidelity approximation can capture well the macroscopic quantities of the solution to the Boltzmann equation in the random space. Additionally, a priori estimate on the accuracy between the high- and bi-fidelity solutions together with a convergence analysis is established. Finally, we present extensive numerical experiments to verify the efficiency and accuracy of our proposed method.

97 MATHEMATICS AND COMPUTING↗

Quantifying UAS Observation Error Variance Used in Data Assimilation Systems and Its Impact on Predictive Skill

Observation error determines the weights of the observations and background state used in data assimilation to generate analyses. Quantifying observation error is critical for the optimal assimilation of observational data sets. Uncrewed Aircraft System (UAS) observations have shown potential benefits in filling observational gaps in the lower atmosphere; however, characterization of their error characteristics has been limited. To optimize the use of UAS observations in numerical weather prediction, UAS observation error is estimated based on the 3‐cornered hat diagnostic approach which uses three independent estimates of the atmospheric state. This approach is applied to data from the 2018 Lower Atmospheric Profiling Studies at Elevation‐a Remotely‐piloted Aircraft Team Experiment field campaign using collocated UAS and rawinsonde observations along with output from a set of convection‐permitting model simulations. The estimated observation error values for UAS temperature, wind, and relative humidity measurements were found to be only weakly dependent on height AGL with mean values equal to 0.5°C, 0.8 m s −1 , and 3%, respectively. Only the newly estimated observation error for temperature differed from that previously used to assimilate commercial aircraft observations into global models (1.0°C). However, using this reduced temperature observation error produced more accurate mesoscale analyses and forecasts of both terrain‐driven flows and convection initiation generated by colliding outflow boundaries within the San Luis Valley of Colorado.

54 ENVIRONMENTAL SCIENCES↗

An Efficient, Globally Convergent Method for Optimization Under Uncertainty Using Adaptive Model Reduction and Sparse Grids

We present a new approach to efficiently solve optimization problems constrained by partial differential equations (PDEs) with uncertain coefficients. The method leverages two sources of inexactness that trade accuracy for speed: (1) stochastic collocation based on dimension-adaptive sparse grids (SGs), which approximates the stochastic objective function with a limited number of quadrature nodes, and (2) projection-based reduced-order models (ROMs), which generate efficient approximations to PDE solutions. These two sources of inexactness lead to inexact objective function and gradient evaluations, which are managed by a trust-region method that guarantees global convergence by adaptively refining the SG and ROM until a proposed error indicator drops below a tolerance specified by trust-region convergence theory. A key feature of the proposed method is that the error indicator---which accounts for errors incurred by both the SG and ROM---must be only an asymptotic error bound, i.e., a bound that holds up to an arbitrary constant that need not be computed. This allows the approach to be applicable to a wide range of problems, including those where sharp, computable error bounds are not available; this distinguishes the proposed method from previous works. Numerical experiments performed on a model problem from optimal flow control under uncertainty verify global convergence of the method and demonstrate the method's ability to outperform previously proposed alternatives.

97 MATHEMATICS AND COMPUTING↗

Deterministic and Monte Carlo Nuclear Data Adjustment Methods [Slides]

For the Bayesian Monte Carlo methodology, a need to understand convergence of the posterior moments as a function of the number of parameter realizations is required. In high-dimensional systems, it can be very costly to sample entire parameter space and perform functional evaluation for every realization. Bayesian Monte Carlo allows one to relax the GLLS approximations of model linearity and prior/posterior PDF shape. The Bayesian Stochastic Collocation Method is a deterministic approach to “sample” the parameter space. It allows one to relax the GLLS approximations of model linearity and posterior PDF shape. Higher-order posterior moments (i.e., skewness, kurtosis, etc.) can be studied through polynomial expansion. Tensor product quadrature scales poorly and can use sparse grid quadrature methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Model Predictive Control to Improve Grid Resilience

The following article details a model predictive control (MPC) to improve grid resilience when faced with variable generation resources. This topic is of significant interest to utility power systems where distributed intermittent energy sources will increase significantly and be relied on for electric grid ancillary services. Previous work on MPCs has focused on narrowly targeted control applications such as improving electric vehicle (EV) charging infrastructure or reducing the cost of integrating Energy Storage Systems (ESSs) into the grid. In contrast, this article develops a comprehensive treatment of the construction of an MPC tailored to electric grids and then applies it integration of intermittent energy resources. To accomplish this, the following article includes a description of a reduced order model (ROM) of an electric power grid based on a circuit model, an optimization formulation that describes the MPC, a collocation method for solving linear time-dependent differential algebraic equations (DAEs) that result from the ROM, and an overall strategy for iteratively refining the behavior of the MPC. Next, the algorithm is validated using two separate numerical experiments. First, the algorithm is compared to an existing MPC code and the results are verified by a numerically precise simulation. It is shown that this algorithm produces a control comparable to existing algorithms and the behavior of the control carefully respects the bounds specified. Second, the MPC is applied to a small nine bus system that contains a mix of turbine-spinning-machine-based and intermittent generation in order to demonstrate the algorithm’s utility for resource planning and control of intermittent resources. This study demonstrates how the MPC can be tuned to change the behavior of the control, which can then assist with the integration of intermittent resources into the grid. The emphasis throughout the paper is to provide systematic treatment of the topic and produce a novel nonlinear control compatible design framework applicable to electric grids and the control of variable resources. This differs from the more targeted application-based focus in most presentations.

microgrid↗

Vegetation Monitoring Optimization with Normalized Difference Vegetation Index and Evapotranspiration Using Remote Sensing Measurements and Land Surface Models Over East Africa

The majority of people in East Africa rely on the agro-pastoral system for their livelihood, which is highly vulnerable to droughts and flooding. Agro-pastoral droughts are endemic to the region and are considered the main natural hazard that contributes to food insecurity. Drought begins with rainfall deficit, gradually leading to soil moisture deficit, higher land surface temperature, and finally impacts to vegetation growth. Therefore, monitoring vegetation conditions is essential in understanding the progression of drought, potential effects on food security, and providing early warning information needed for drought mitigation decisions. Because vegetation processes couple the land and atmosphere, monitoring of vegetation conditions requires consideration of both water provision and demand. While there is consensus in using either the Normalized Difference Vegetation Index (NDVI) or evapotranspiration(ET) for vegetation monitoring, a comprehensive assessment optimizing the use of both has not yet been done. Moreover, the evaluation methods for understanding the relationships between NDVI and ET for vegetation monitoring are also limited. Taking these gaps into account we have developed a framework to optimize vegetation monitoring using both NDVI and ET by identifying where they perform the best by using triple collocation and cross-correlation methods. We estimated the random error structure in Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI; ET from the Operational Simplified Surface Energy Balance (SSEBop) model; and ET from land surface models (LSMs). LSM ET and SSEBop ET have been found to be better indicators for vegetation monitoring during extreme drought events, while NDVI could provide better information on vegetation condition during wetter than normal conditions. The random error structures of these variables suggest that LSM ET is most likely to provide important information for vegetation monitoring over low and high ends of the vegetation fraction areas. Over moderate vegetative areas, any of these variables could provide important vegetation information for drought characterization and food security assessments. While this study provides a framework for optimizing vegetation monitoring for drought and food security assessments over East Africa, the framework can be adopted to optimize vegetation monitoring over any other drought and food insecure region of the world.

triple collocation↗

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗

Deducing Land-Atmosphere Coupling Regimes from SMAP Soil Moisture

In recent years, there has been a growing recognition of the significance of Land-Atmosphere (L-A) interactions and feedback mechanisms in understanding and predicting Earth’s water and energy cycles. Soil moisture plays a critical role in mediating the strength of L-A interactions and is important for understanding the complex and governing processes across this interface. This study aims to identify the significance of soil moisture in identifying L-A coupling strength within the Convective Triggering Potential (CTP) and Humidity Index (HI) framework. To address this, a consistent and reliable dataset of atmospheric profiles is created by merging CTP and HI using Triple Collocation (TC) with three reanalysis datasets. The merged CTP and HI product demonstrates enhanced performance globally as compared to the individual datasets when validated with radiosonde and satellite observations. This merged product of CTP and HI is then used to compare the L-A coupling strength based on Soil Moisture Active Passive Level 3 (SMAPL3) and SMAP Level 4 (SMAPL4) over two decades (2003-2022) where L-A coupling strength is defined as the persistence probability within the dry and wet coupling regimes. Results indicate that the persistency-based coupling strength is related to the ability of soil moisture to predict future atmospheric humidity and dry vs. wet coupling state. The coupling strength in SMAPL4 is consistently stronger than in SMAPL3 and is likely due to its reliance on a land surface model and reduced susceptibility to random noise. The difference in coupling strength based on the same CTP-HI underscores the importance of soil moisture data in estimating coupling strength within the CTP-HI framework. These findings lay the groundwork for understanding the role of L-A interactions and drought evolution due to soil moisture variations, by providing insight into the quantification of coupling strength and its role in drought monitoring and forecast efforts.

Land-atmosphere coupling↗

Optical properties and CCN activity of aerosols in a high-altitude Himalayan environment: Results from RAWEX-GVAX: CCN activity of aerosols over Himalayas

The seasonality and mutual dependence of aerosol optical properties and cloud condensation nuclei (CCN) activity under varying meteorological conditions at the high-altitude Nainital site (~2 km) in the Indo-Gangetic Plains were examined using nearly year-round measurements (June 2011 to March 2012) at the Atmospheric Radiation Measurement (ARM) mobile facility as part of the RAWEX-GVAX experiment of the Indian Space Research Organization and the U.S. Department of Energy. The results from collocated measurements provided enhanced aerosol scattering and absorption coefficients, CCN concentrations and total condensation nuclei (CN) concentrations during the dry autumn and winter months. The CCN concentration (at a supersaturation of 0.46) was higher during periods of high aerosol absorption (single-scattering albedo (SSA) < 0.80) than during periods of high aerosol scattering (SSA > 0.85), indicating that the aerosol composition seasonally changes and influences the CCN activity. The monthly mean CCN activation ratio (at a supersaturation of 0.46) was highest (> 0.7) in late autumn (November); this finding is attributed to the contribution of biomass-burning aerosols to CCN formation at high supersaturation conditions.

54 ENVIRONMENTAL SCIENCES↗

Insights into riming and aggregation processes as revealed by aircraft, radar, and disdrometer observations for a 27 April 2011 widespread precipitation event: Insights into Riming and Aggregation

Our study presents aircraft spiral ascent and descent observations intercepting a transition to riming processes during widespread stratiform precipitation. The sequence is documented using collocated scanning and profiling radar, including longer-wavelength dual polarization measurements and shorter-wavelength Doppler spectra. Riming regions are supported using aircraft measurements recording elevated liquid water concentrations, spherical particle shapes, and saturation with respect to water. Profiling cloud radar observations indicate riming regions during the event as having increasing particle fall speeds, rapid time-height changes, and bimodalities in Doppler spectra. These particular riming signatures are coupled to scanning dual polarization radar observations of higher differential reflectivity (ZDR) aloft. Moreover, reduced melting layer enhancements and delayed radar bright-band signatures in the column are also observed during riming periods, most notably with the profiling radar observations. The bimodal cloud radar Doppler spectra captured near riming zones indicate two time-height spectral ice peaks, one rimed particle peak, and one peak associated with pristine ice needle generation and/or growth between -4°C and -7°C also sampled by aircraft probes. We observe this pristine needle population near the rimed particle region which gives a partial explanation for the enhanced ZDR. The riming signatures aloft and radar measurements within the melting level are weakly lag correlated (r~0.6) with smaller median drop sizes at the surface, as compared with later times when aggregation of larger particle sizes was believed dominant.

54 ENVIRONMENTAL SCIENCES↗

On the unified estimation of turbulence eddy dissipation rate using Doppler cloud radars and lidars: Radar and Lidar Turbulence Estimation

Coincident profiling observations from Doppler lidars and radars are used to estimate the turbulence energy dissipation rate (ε) using three different data sources: (i) Doppler radar velocity (DRV), (ii) Doppler lidar velocity (DLV), and (iii) Doppler radar spectrum width (DRW) measurements. Likewise, the agreement between the derived ε estimates is examined at the cloud base height of stratiform warm clouds. Collocated ε estimates based on power spectra analysis of DRV and DLV measurements show good agreement (correlation coefficient of 0.86 and 0.78 for both cases analyzed here) during both drizzling and nondrizzling conditions. This suggests that unified (below and above cloud base) time-height estimates of ε in cloud-topped boundary layer conditions can be produced. This also suggests that eddy dissipation rate can be estimated throughout the cloud layer without the constraint that clouds need to be nonprecipitating. Eddy dissipation rate estimates based on DRW measurements compare well with the estimates based on Doppler velocity but their performance deteriorates as precipitation size particles are introduced in the radar volume and broaden the DRW values. And, based on this finding, a methodology to estimate the Doppler spectra broadening due to the spread of the drop size distribution is presented. Furthermore, the uncertainties in ε introduced by signal-to-noise conditions, the estimation of the horizontal wind, the selection of the averaging time window, and the presence of precipitation are discussed in detail.

54 ENVIRONMENTAL SCIENCES↗

Retrievals of ice cloud microphysical properties of deep convective systems using radar measurements: Convective Cloud Microphysical Retrieval

This study presents new algorithms for retrieving ice cloud microphysical properties (ice water content (IWC) and median mass diameter (Dm)) for the stratiform and thick anvil regions of Deep Convective Systems (DCSs) using Next-Generation Radar (NEXRAD) reflectivity and recently developed empirical relationships from aircraft in situ measurements during the Midlatitude Continental Convective Clouds Experiment (MC3E). A classic DCS case on 20 May 2011 is used to compare the retrieved IWC profiles with other retrieval and cloud-resolving model simulations. The mean values of each retrieved and simulated IWC fall within one standard derivation of the other two. The statistical results from six selected cases during MC3E show that the aircraft in situ derived IWC and Dm are 0.47 ± 0.29 g m-3 and 2.02 ± 1.3 mm, while the mean values of retrievals have a positive bias of 0.16 g m-3 (34%) and a negative bias of 0.39 mm (19%). To validate the newly developed retrieval algorithms from this study, IWC and Dm are performed with other DCS cases during Bow Echo and Mesoscale Convective Vortex Experiment (BAMEX) field campaign using composite gridded NEXRAD reflectivity and compared with in situ IWC and Dm from aircraft. A total of 64 1-min collocated aircraft and radar samples are available for comparisons, and the averages of radar retrieved and aircraft in situ measured IWCs are 1.22 g m-3 and 1.26 g m-3 with a correlation of 0.5, and their averaged Dm values are 2.15 and 1.80 mm. These comparisons have shown that the retrieval algorithms 45 developed during MC3E can retrieve similar ice cloud microphysical properties of DCS to aircraft in situ measurements during BAMEX with median errors of ~40% and ~25% for IWC and Dm retrievals, respectively. This is indicating our retrieval algorithms are suitable for other midlatitude continental DCS ice clouds, especially at stratiform rain and thick anvil regions. In addition, based on the averaged IWC and Dm values during MC3E and BAMEX, the DCS IWC values over midlatitude are significantly different, while their Dm values are close to each other. On the other hand, these DCS IWC and Dm values are 1-2 orders of magnitude larger than those of single-layered cirrus clouds over midlatitudes.

Tian, Jingjing↗

Low-CCN concentration air masses over the eastern North Atlantic: Seasonality, meteorology, and drivers

In this paper, a 20 month cloud condensation nucleus concentration (N CCN ) data set from Graciosa Island (39°N, 28°W) in the remote North Atlantic is used to characterize air masses with low cloud condensation nuclei (CCN) concentrations. Low-CCN events are defined as 6 h periods with mean N CCN <20 cm -3 (0.1% supersaturation). A total of 47 low-CCN events are identified. Surface, satellite, and reanalysis data are used to explore the meteorological and cloud context for low-CCN air masses. Low-CCN events occur in all seasons, but their frequency was 3 times higher in December–May than during June–November. Composites show that many of the low-CCN events had a common meteorological basis that involves southerly low-level flow and rather low wind speeds at Graciosa. Anomalously low pressure is situated to the west of Graciosa during these events, but back trajectories and lagged SLP composites indicate that low-CCN air masses often originate as cold air outbreaks to the north and west of Graciosa. Finally, low-CCN events were associated with low cloud droplet concentrations (N d ) at Graciosa, but liquid water path (LWP) during low-CCN events was not systematically different from that at other times. Satellite N d and LWP estimates from MODIS collocated with Lagrangian back trajectories show systematically lower N d and higher LWP several days prior to arrival at Graciosa, consistent with the hypothesis that observed low-CCN air masses are often formed by coalescence scavenging in thick warm clouds, often in cold air outbreaks.

54 ENVIRONMENTAL SCIENCES↗

Influence of Emissions and Aqueous Processing on Particles Containing Black Carbon in a Polluted Urban Environment: Insights From a Soot Particle-Aerosol Mass Spectrometer

Inorganic and organic coatings on black carbon (BC) particles can enhance light absorption and affect atmospheric lifetimes of BC-containing particles and thus have significant implications for climate. To study the physical and chemical characteristics of atmospheric BC and BC-associated coatings, a soot particle-aerosol mass spectrometer was deployed during the winter of 2014–2015 in Fresno, a city located in the San Joaquin Valley of California, to selectively analyze BC-containing particles. Comparing soot particle-aerosol mass spectrometer measurements to those from the collocated single-particle soot photometer (SP2) and high-resolution aerosol mass spectrometer, we found that 17% of total submicrometer aerosol mass was associated with BC-containing particles, suggesting that a majority of the fine particles in Fresno contained no BC. Most BC-containing particles appeared to be associated with residential wood burning and vehicular traffic. These particles typically had a bulk-average mass ratio of coating to BC (R coat/rBC ) less than 2. However, during periods of persistent fog larger R coat/rBC values were observed, with the coatings primarily composed of secondary inorganic and organic components that likely resulted from aqueous-phase processing. Specifically, compared to periods with less fog, the BC coating increased in concentration and contained a larger fraction of nitrate and oxidized organic matter. The size distributions of BC and associated organic coating were generally centered around 300 nm in vacuum aerodynamic diameter. However, during foggy periods BC had an additional peak at ~400 nm and organics and nitrate displayed a prominent mode in the accumulation size range.

54 ENVIRONMENTAL SCIENCES↗

A minimally invasive, efficient method for propagation of full-field uncertainty in solid dynamics

In this work, we present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

42 ENGINEERING↗