Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Predictive analytics”

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 325 records · Page 18

Physiological fidelity or model parsimony? The relative performance of reverse-toxicokinetic modeling approaches

Physiologically-based toxicokinetic (PBTK) models are often developed to facilitate in vitro to in vivo extrapolation (IVIVE) using a top-down, compartmental approach, favoring architectural simplicity over physiological fidelity despite the lack of general guidelines relating model design to dynamical predictions. Here we explore the impact of design choice (high vs. low fidelity) on chemical distribution throughout an animal’s organ system. Results We contrast transient dynamics and steady states of three previously proposed PBTK models of varying complexity in response to chemical exposure. The steady states for each model were determined analytically to predict exposure conditions from tissue measurements. Steady state whole-body concentrations differ between models, despite identical environmental conditions, which originates from varying levels of physiological fidelity captured by the models. These differences affect the relative predictive accuracy of the inverted models used in exposure reconstruction to link effects-based exposure data with whole-organism response thresholds obtained from in vitro assay measurements. Conclusions Our results demonstrate how disregarding physiological fideltiy in favor of simpler models affects the internal dynamics and steady state estimates for chemical accumulation within tissues, which, in turn, poses significant challenges for the exposure reconstruction efforts that underlie many IVIVE methods. Developing standardized systems-level models for ecological organisms would not only ensure predictive consistency among future modeling studies, but also ensure pragmatic extrapolation of in vivo effects from in vitro data or modeling exposure-response relationships.

59 BASIC BIOLOGICAL SCIENCES↗

FY2022 Progress on Imbibition Testing in Containment Science

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as for vapor-phase contaminant transport and/or containment of noble gases in the subsurface. In this second progress report of FY22, we present two ongoing activities related to imbibition testing on volcanic rock samples. We present the development of a new analytical solution predicting the temperature response observed during imbibition into dry samples, as discussed in our previous first progress report for FY22. We also illustrate the use of a multi-modal capillary pressure distribution to simulate both early- and late-time imbibition data collected on tuff core that can exhibit multiple pore types. These FY22 imbibition tests were conducted for an extended period (i.e., far beyond the time required for the wetting front to reach the top of the sample), which is necessary for parameter estimation and characterization of two different pore types within the samples.

42 ENGINEERING↗

Data analytics for intermodal freight transportation applications

With the growth of intermodal freight transportation, it is important that transportation planners and decision-makers are knowledgeable about freight flow data to make informed decisions. This is particularly true with Intelligent Transportation Systems (ITS) offering new capabilities for intermodal freight transportation. Specifically, ITS enables access to multiple different data sources, but they have different formats, resolutions, and time scales. Thus, knowledge of data science is essential to be successful in future ITS-enabled intermodal freight transportation systems. This chapter discusses the commonly used descriptive and predictive data analytic techniques in intermodal freight transportation applications. These techniques cover the entire spectrum of univariate, bivariate, and multivariate analyses. In addition to illustrating how to apply these techniques manually, this chapter will also show how to apply them using the statistical software R. Additional exercises are provided for those who wish to apply the described techniques to more complex problems.

Huynh, Nathan↗

Reconstructing a f ( R ) theory from the α-Attractors

We show an analogy at high curvature between a f ( R ) = R + aR {sup n} {sup −} {sup 1} + bR {sup 2} theory and the α-Attractors. We calculate the expressions of the parameters a , b and n as functions of α and the predictions of the model f ( R ) = R + aR {sup n} {sup −} {sup 1} + bR {sup 2} on the scalar spectral index n {sub s} and the tensor-to-scalar ratio r . We find that the power law correction R {sup n} {sup −} {sup 1} allows for a production of gravitational waves enhanced with respect to the one in the Starobinsky model, while maintaining a viable prediction on n {sub s}. We numerically reconstruct the full α-Attractors class of models testing the goodness of our high-energy approximation f ( R ) = R + aR {sup n} {sup −} {sup 1} + bR {sup 2}. Moreover, we also investigate the case of a single power law f ( R ) = γ R {sup 2} {sup −} {sup δ} theory, with γ and δ free parameters. We calculate analytically the predictions of this model on the scalar spectral index n {sub s} and the tensor-to-scalar ratio r and the values of δ which are allowed from the current observational results. We find that −0.015 < δ < 0.016, confirming once again the excellent agreement between the Starobinsky model and observation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Color Confinement, Hadron Dynamics, and Hadron Spectroscopy from Light-Front Holography and Superconformal Algebra

The QCD light-front Hamiltonian equation H L F Ψ = M 2 Ψ derived from quantization at fixed LF time τ = t + z / c provides a causal, frame-independent method for computing hadron spectroscopy as well as dynamical observables such as structure functions, transverse momentum distributions, and distribution amplitudes. The QCD Lagrangian with zero quark mass has no explicit mass scale. de Alfaro, Fubini, and Furlan (dAFF) have made an important observation that a mass scale can appear in the equations of motion without affecting the conformal invariance of the action if one adds a term to the Hamiltonian proportional to the dilatation operator or the special conformal operator. If one applies the dAFF procedure to the QCD light-front Hamiltonian, it leads to a color-confining potential κ 4 ζ 2 for mesons, where ζ 2 is the LF radial variable conjugate to the q q ¯ invariant mass squared. The same result, including spin terms, is obtained using light-front holography, the duality between light-front dynamics and A d S 5 , if one modifies the A d S 5 action by the dilaton e κ 2 z 2 in the fifth dimension z . When one generalizes this procedure using superconformal algebra, the resulting light-front eigensolutions provide a unified Regge spectroscopy of meson, baryon, and tetraquarks, including remarkable supersymmetric relations between the masses of mesons and baryons and a universal Regge slope. The pion q q ¯ eigenstate has zero mass at m q = 0 . The superconformal relations also can be extended to heavy-light quark mesons and baryons. This approach also leads to insights into the physics underlying hadronization at the amplitude level. I will also discuss the remarkable features of the Poincaré invariant, causal vacuum defined by light-front quantization and its impact on the interpretation of the cosmological constant. AdS/QCD also predicts the analytic form of the nonperturbative running coupling α s ( Q 2 ) ∝ e - Q 2 / 4 κ 2 . The mass scale κ underlying hadron masses can be connected to the parameter Λ M S ¯ in the QCD running coupling by matching the nonperturbative dynamics to the perturbative QCD regime. The result is an effective coupling α s ( Q 2 ) defined at all momenta. One obtains empirically viable predictions for spacelike and timelike hadronic form factors, structure functions, distribution amplitudes, and transverse momentum distributions. Finally, I address the interesting question of whether the momentum sum rule is valid for nuclear structure functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards developing a compact model for magnetization switching in straintronics magnetic random access memory devices

Strain-mediated magnetization switching in a magnetic tunneling junction (MTJ) by exploiting a combination of piezoelectricity and magnetostriction has been proposed as an energy efficient alternative to spin transfer torque (STT) and field induced magnetization switching methods in MTJ-based magnetic random access memories (MRAM). Theoretical studies have shown the inherent advantages of strain-assisted switching, and the dynamic response of the magnetization has been modeled using the Landau-Lifshitz-Gilbert (LLG) equation. However, an attempt to use LLG for simulating dynamics of individual elements in large-scale simulations of multi-megabyte straintronics MRAM leads to extremely time-consuming calculations. Hence, a compact analytical solution, predicting the flipping delay of the magnetization vector in the nanomagnet under stress, combined with a liberal approximation of the LLG dynamics in the straintronics MTJ, can lead to a simplified model of the device suited for fast large-scale simulations of multi-megabyte straintronics MRAMs. In this work, a tensor-based approach is developed to study the dynamic behavior of the stressed nanomagnet. First, using the developed method, the effect of stress on the switching behavior of the magnetization is investigated to realize the margins between the underdamped and overdamped regimes. The latter helps the designer realize the oscillatory behavior of the magnetization when settling along the minor axis, and the dependency of oscillations on the stress level and the damping factor. Next, a theoretical model to predict the flipping delay of the magnetization vector is developed and tested against LLG-based numerical simulations to confirm the accuracy of findings. Lastly, the obtained delay is incorporated into the approximate solutions of the LLG dynamics, in order to create a compact model to liberally and quickly simulate the magnetization dynamics of the MTJ under stress. Using the developed delay equation, the efficiency of the straintronics switching over the STT method is highlighted by analytically investigating the energy-delay trade-off of both methodologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting boron coordination in multicomponent borate and borosilicate glasses using analytical models and machine learning

Accurate prediction of boron coordination in multicomponent glasses is critical in glass science and technology as it strongly affects the properties of borate and borosilicate glasses. We have collected a dataset containing 657 glasses from literature with boron coordination values and developed models using analytical functions based on the well accepted Dell, Xiao and Bray model. Good prediction of boron coordination with a R 2 value higher than 0.8 was obtained. The large variation of boron coordination from experiments, originated from sample preparations and characterizations, led to difficulties in obtaining models with better prediction performance. Various machine learning (ML) algorithms were evaluated and slightly better prediction performance was observed; however, interpretation of the ML models is less straight forward. In conclusion, this study developed various models capable of providing quantitative boron coordination predictions, providing insights into its structural roles in multi-component glasses, and suggesting fruitful areas for future research.

36 MATERIALS SCIENCE↗

Spatial correlation of proton irradiation-induced activity and dose in polymer gel phantoms for PET/CT delivery verification studies

Purpose: This work demonstrates a novel application of BANG3-Pro2 polymer gel dosimeter as a dosimetric phantom able to accurately capture both dose and induced activity. Methods: BANG3-Pro2 dosimeters were irradiated with a clinical proton beam using an unmodulated beam and a spread-out Bragg peak (SOBP) modulation, the latter with a Lucite compensator to introduce a range offset in one quadrant of the circular field. The dosimeters were imaged in a nearby positron emission tomography/computed tomography (PET/CT) unit starting within 5 min of beam-off. Induced positron emission (PE) activity along the central axis of the beam was compared to analytical calculations. Dose distributions were read out using an optical CT scanner and were validated against ion chamber measurements and the treatment plan. The offset between the distal fall-off of dose and activity (50% level) was determined over the entire irradiated field. Lateral profiles of PE were correlated to measured dose for the unmodulated beam delivery. Results: Measured profiles of PE activity along the central beam axis were found to be within 10% of the predictions of analytical calculations. The depth-dose profiles agreed with the reference values (ion chamber or treatment plan) within 3%. The offset between the depth profiles of dose and activity for the unmodulated beam was 8.4 {+-} 1.4 mm. For the compensator-based SOBP delivery, the distribution of offsets throughout the field was found to be bimodal, with the mean of 8.9 {+-} 2.8 mm for the thinner region of the compensator and 4.3 {+-} 2.5 mm for the thicker region. For the pristine beam delivery, lateral profiles of dose and activity were found to exhibit fair spatial correlation throughout the beam range, with the mean 2D gamma index of 0.42 and 91% of the evaluated pixels passing the test. Conclusions: This work presents the first demonstration of simultaneous and accurate experimental measurement of three-dimensional distributions of dose and induced activity and lays the groundwork for further investigations using BANG3-Pro2 as a dosimetric phantom in PET/CT delivery verification studies.

61 RADIATION PROTECTION AND DOSIMETRY↗

Validation of SPH code Spheral to model interacting solid bodies in a supersonic flow

Contemporary discussions of planetary defense involve analyzing the risks posed by smaller sized, 20 to 200 m diameter, asteroids which are capable of breaking up in the atmosphere and generating a blast wave. Consequence assessments for this size class of asteroids are performed through fast-running analytic or semi-analytic models which are informed by high-fidelity hydrocode simulations of asteroid entry and breakup. However, insufficient historical data necessitates validating the independent physical processes which dominate airburst events. Here, the Fluid Solid Interface Smoothed Particle Hydrodynamics solver was previously used by Pearl et al. in 2023 to model the Chelyabinsk airburst and is used here to perform a series of validation simulations. The first effort involves modeling a cylinder in a hypersonic flow and comparing the bow shock geometry to that predicted by analytic theory. The second effort involves modeling the separation of two spherical bodies in supersonic flow and validating against experimental footage. Combined, these exercises demonstrate the ability of the code to model the flight-path of interacting solid bodies in a hypersonic flow.

Airburst↗

Structural short-range forces between solid-melt interfaces

We predict the structural interaction of crystalline solid-melt interfaces using amplitude equations which are derived from classical density functional theory or phase field crystal modeling. The solid ordering decays exponentially on the scale of the interface thickness at solid-melt interfaces; the overlap of two such profiles leads to a short-range interaction, which is mainly carried by the longest-range density waves, which can facilitate grain boundary premelting. Here, we calculate the tail of these interactions, depending on the relative translation of the two crystals, fully analytically and predict the interaction potential, and compare it to numerical simulations. For grain boundaries the interaction is predicted to decay twice as fast as for two crystals without misorientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

MODULATION OF THE GALACTIC LOW-ENERGY PROTON SPECTRUM IN THE INNER HELIOSPHERE

We study the energy spectra of 0.3-100 MeV protons and find that, at the lower part of the galactic particle spectrum, they are significantly steeper than the J(E) {approx} E spectrum predicted by analytical approximations, such as the force-field model of modulation. We select a series of low-flux periods, and approximate the spectral form by J(E) = AE{sup -{gamma}} + CE{sup {nu}}, where the two terms describe solar/heliospheric and galactic components, respectively. By determining the best fit parameters to energy spectra, correlations are sought with solar activity indices and between the parameters themselves. In the majority of cases, {nu} turns out to be between 1.2 and 1.4, with an average of 1.32 {+-} 0.12, significantly greater than the commonly expected {nu} = 1 predicted by the force-field approximation. In modulation theories {nu} > 1 corresponds to a negative Compton-Getting factor, which poses a challenge. Such an inversion may occur if the radial diffusion coefficient {kappa}{sub rr} < rV (where r is heliocentric distance and V solar wind speed), in which case a large fraction of the 10-100 MeV protons reaching 1 AU would have been cooled down within 1 AU and subsequently convected outward by the solar wind. We also find that the position of the intensity minimum of the proton spectrum dividing the solar and galactic populations shifts toward higher values with increasing solar activity. Correlations obtained with solar activity indicate that the slope of the solar/heliospheric spectrum is practically independent of solar activity. Observations are compared with numerical solutions of the modulation equation adopting simple spherical models. Possible interpretations are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Weak localization with nonlinear bosonic matter waves

We investigate the coherent propagation of dilute atomic Bose-Einstein condensates through irregularly shaped billiard geometries that are attached to uniform incoming and outgoing waveguides. Using the mean-field description based on the nonlinear Gross-Pitaevskii equation, we develop a diagrammatic theory for the self-consistent stationary scattering state of the interacting condensate, which is combined with the semiclassical representation of the single-particle Green function in terms of chaotic classical trajectories within the billiard. This analytical approach predicts a universal dephasing of weak localization in the presence of a small interaction strength between the atoms, which is found to be in good agreement with the numerically computed reflection and transmission probabilities of the propagating condensate. The numerical simulation of this quasi-stationary scattering process indicates that this interaction-induced dephasing mechanism may give rise to a signature of weak antilocalization, which we attribute to the influence of non-universal short-path contributions. - Highlights: Black-Right-Pointing-Pointer Numerical simulation of scattering of Bose-Einstein condensate through billiards. Black-Right-Pointing-Pointer Novel analytical semiclassical theory for nonlinear coherent scattering. Black-Right-Pointing-Pointer Inversion of weak localization due to mean-field interaction within the condensate. Black-Right-Pointing-Pointer Relevance of non-universal short-path contributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mixing dynamics in bubbling fluidized beds

Solids mixing affects thermal and concentration gradients in fluidized bed reactors and is, therefore, critical to their performance. Despite substantial effort over the past decades, understanding of solids mixing continues to be lacking because of technical limitations of diagnostics in large pilot and commercial‐scale reactors. This study is focused on investigating mixing dynamics and their dependence on operating conditions using computational fluid dynamics simulations. Toward this end, fine‐grid 3D simulations are conducted for the bubbling fluidization of three distinct Geldart B particles (1.15 mm LLDPE, 0.50 mm glass, and 0.29 mm alumina) at superficial gas velocities U/Umf = 2–4 in a pilot‐scale 50 cm diameter bed. The Two‐Fluid Model (TFM) is employed to describe the solids motion efficiently while bubbles are detected and tracked using MS3DATA. Detailed statistics of the flow‐field in and around bubbles are computed and used to describe bubble‐induced solids micromixing: solids upflow driven in the nose and wake regions while downflow along the bubble walls. Further, within these regions, the hydrodynamics are dependent only on particle and bubble characteristics, and relatively independent of the global operating conditions. Based on this finding, a predictive mechanistic, analytical model is developed which integrates bubble‐induced micromixing contributions over their size and spatial distributions to describe the gross solids circulation within the fluidized bed. Finally, it is shown that solids mixing is affected adversely in the presence of gas bypass, or throughflow, particularly in the fluidization of heavier particles. This is because of inefficient gas solids contacting as 30–50% of the superficial gas flow escapes with 2–3× shorter residence time through the bed. This is one of the first large‐scale studies where both the gas (bubble) and solids motion, and their interaction, are investigated in detail and the developed framework is useful for predicting solids mixing in large‐scale reactors as well as for analyzing mixing dynamics in complex reactive particulate systems. © 2017 American Institute of Chemical Engineers AIChE J , 63: 4316–4328, 2017

Bakshi, A.↗

Development and Assessment of CTFFuel

This milestone report documents an approach that seeks to develop a standalone fuel solver to simulate steady-state and transient thermo-mechanical responses of light-water reactor fuel rods and assess its temperature predictions against analytical solutions, and a well-known and validated fuel performance code. This manuscript focuses on five major topics: (1) development of the CTFFuel, (2) discussion of the newly implemented fuel deformation models in dynamic gap conductance model of CTFFuel, (3) enabling the dynamic gap conductance model in VERA-CS, (4) verification of the solid conduction solution and comparison of the CTFFuel with respect to the workhorse code, CTF, and (5) benchmark of the thermal conductivity degradation using the burnup dependent fuel thermal conductivity models and benchmark of the fuel centerline temperature using the dynamic gap conductance model with Halden measurements and FRAPCON-4.0 simulation results. CTFFuel temperature predictions are compared to predictions from FRAPCON-4.0 for UO 2 , UO 2 +Gd 2 O 3 and MOX fuels through reactor’s life. The main motivation is to understand discrepancies between CTFFuel and FRAPCON-4.0’s temperature predictions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A three-year dataset supporting research on building energy management and occupancy analytics

Abstract This paper presents the curation of a monitored dataset from an office building constructed in 2015 in Berkeley, California. The dataset includes whole-building and end-use energy consumption, HVAC system operating conditions, indoor and outdoor environmental parameters, as well as occupant counts. The data were collected during a period of three years from more than 300 sensors and meters on two office floors (each 2,325 m 2 ) of the building. A three-step data curation strategy is applied to transform the raw data into research-grade data: (1) cleaning the raw data to detect and adjust the outlier values and fill the data gaps; (2) creating the metadata model of the building systems and data points using the Brick schema; and (3) representing the metadata of the dataset using a semantic JSON schema. This dataset can be used in various applications—building energy benchmarking, load shape analysis, energy prediction, occupancy prediction and analytics, and HVAC controls—to improve the understanding and efficiency of building operations for reducing energy use, energy costs, and carbon emissions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Syntrophic Co-Cultures of Clostridium Organisms to Produce Higher Alcohols & Other C6-C8 Metabolites (Final Report)

The goal of this project was to advance the systems biology understanding and predictive modeling of synthetic & syntrophic Clostridium microbial consortia, focusing on elucidation of metabolic networks and environmental signals in the consortia. The project has direct applicability to lignocellulosic-biomass based production of higher alcohols as advanced biofuels and C6-C8 metabolites, that can be used as chemicals or serve as biofuel precursors. This project studied and optimized three synthetic syntrophic systems. It examined the population dynamics using flow-cytometry, time-lapse microscopy and PCR analysis. 13C-based tracer analyses was used to examine the metabolite exchange between the syntrophic cell populations and the impact of those interactions on the transcriptome of the individual populations. To enhance our analytical and predictive capabilities, genome-scale models (GSMs) for these syntrophies was developed. RNAseq data for these syntrophic coculture systems were acquired to enable a molecular level understanding of the syntrophies aiming to identify the genetic networks of each organism in the co-culture and compare those against the networks of pure cultures.

09 BIOMASS FUELS↗

Catalytic Deoxydehydration of Carbohydrates and Polyols to Chemicals and Fuels

As the world's fossil fuel resources are being depleted and their costs increase, there is an urgent need to discover and develop new processes for the conversion of renewable, biomass resources into fuels and chemical feedstocks. Research and development in this area have been given high priority by both governmental agencies and industry. To increase the energy content and decrease the boiling points of biomass-derived carbohydrates and polyols to the useful liquid range it is necessary to chemically remove water (dehydrate) and, preferably, oxygen (deoxygenate/reduce). The poly-hydroxylic nature of carbohydrates is attractive for their use as functionalized chemical building blocks, but it presents a daunting challenge for their selective conversion to single product chemicals or fuels. The long term, practical objective of this project is to develop catalytic processes for the deoxydehydration (DODH) of biomass-derived carbohydrates and polyols to produce unsaturated alcohols and hydrocarbons of value as chemical feedstocks and fuels; DODH: polyol + reductant --(LMOx catalyst)--> unsaturate + oxidized reductant + H2O. Limited prior studies have established the viability of the DODH process with expensive phosphine reductants and rhenium-catalysts. Initial studies in the PI's laboratory have now demonstrated: 1) the moderately efficient conversion of glycols to olefins by the economical sulfite salts is catalyzed by MeReO3 and Z+ReO4-; 2) effective phosphine-based catalytic DODH of representative glycols to olefins by cheap LMoO2 complexes; and 3) computational studies (with K. Houk, UCLA) have identified several Mo-, W-, and V-oxo complexes that are likely to catalyze glycol DODH. Seeking practically useful DODH reactions of complex polyols and new understanding of the reactivity of polyoxo-metal species with biomass-oxygenates we will employ a two-pronged approach: 1) investigate experimentally the reactivity, both stoichiometric and catalytic, of polyoxo-complexes and practical reductants with representative polyols to establish structure/reactivity relationships and reaction mechanisms; and b) carry out parallel computational studies of these reactions and their mechanisms- both analytical and predictive. Our prioritized action plan is: (1) to optimize the catalytic efficiency, assess the substrate scope/selectivity, and address key mechanistic aspects of Re-catalyzed, sulfite-driven DODH reactions; (2) use the findings from (1), together with computational predictions, to discover new, effective non-precious metal catalysts for sulfite-driven DODH reactions; and 3) to initiate exploratory studies of CO- and H2-driven DODH. Successful execution of this research project will: 1) provide practical chemical processes for the conversion of biomass into useful chemicals and fuels: 2) bring fundamental new understanding of chemical reactions involving metal-oxo catalysts; and 3) provide educational and technical training of future energy scientists.

09 BIOMASS FUELS↗

Modeling Interfacial Fracture in Sierra

This report summarizes computational efforts to model interfacial fracture using cohesive zone models in the SIERRA/SolidMechanics (SIERRA/SM) finite element code. Cohesive surface elements were used to model crack initiation and propagation along predefined paths. Mesh convergence was observed with SIERRA/SM for numerous geometries. As the funding for this project came from the Advanced Simulation and Computing Verification and Validation (ASC V&V) focus area, considerable effort was spent performing verification and validation. Code verification was performed to compare code predictions to analytical solutions for simple three-element simulations as well as a higher-fidelity simulation of a double-cantilever beam. Parameter identification was conducted with Dakota using experimental results on asymmetric double-cantilever beam (ADCB) and end-notched-flexure (ENF) experiments conducted under Campaign-6 funding. Discretization convergence studies were also performed with respect to mesh size and time step and an optimization study was completed for mode II delamination using the ENF geometry. Throughout this verification process, numerous SIERRA/SM bugs were found and reported, all of which have been fixed, leading to over a 10-fold increase in convergence rates. Finally, mixed-mode flexure experiments were performed for validation. One of the unexplained issues encountered was material property variability for ostensibly the same composite material. Since the variability is not fully understood, it is difficult to accurately assess uncertainty when performing predictions.

Brown, Arthur A.↗