Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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 667 records · Page 37

Coupling High-Throughput Experiments and Regression Algorithms to Optimize PGM-Free ORR Electrocatalyst Synthesis

Over the past decades, significant improvement has been achieved in the performance of platinum group metal-free (PGM-free) materials as an alternative to Pt-based electrocatalysts for oxygen reduction reaction (ORR). However, further progress in ORR activity requires evaluation of precursors and synthesis approaches. In response to this challenge, we generated a first of its kind experimental data set of 36 samples using high-throughput synthesis and activity measurements. Several control parameters (e.g., Fe precursor identity, the precursor content, and pyrolysis temperature) were varied. We then developed several state-of-the-art machine learning (ML) based regression models to predict ORR activity, dependent on selected synthesis variables. Through an iterative algorithm, higher prediction accuracy (smaller root-mean-square error) was achieved. We identified that gradient boosting regression (GBR) and support vector regression (SVR), among several methods, work best for this data set. Aided by our ML-based surrogate models, we decided to alter catalyst synthesis conditions, which resulted in a 36% increase in measured ORR activity in comparison to the maximum ORR mass activity value of 21.9 A/g catalyst in the original data set. Overall, this combined experiment and machine learning approach represents a promising path forward toward developing highly efficient next-generation ORR electrocatalysts and, more generally, functional materials.

25 ENERGY STORAGE↗

Resolving experimental biases in the interpretation of diffusion experiments with a user-friendly numerical reactive transport approach

The reactive transport code CrunchClay was used to derive effective diffusion coefficients (D e ), clay porosities (ε), and adsorption distribution coefficients (K D ) from through-diffusion data while considering accurately the influence of unavoidable experimental biases on the estimation of these diffusion parameters. These effects include the presence of filters holding the solid sample in place, the variations in concentration gradients across the diffusion cell due to sampling events, the impact of tubing/dead volumes on the estimation of diffusive fluxes and sample porosity, and the effects of O-ring-filter setups on the delivery of solutions to the clay packing. Doing so, the direct modeling of the measurements of (radio)tracer concentrations in reservoirs is more accurate than that of data converted directly into diffusive fluxes. While the above-mentioned effects have already been described individually in the literature, a consistent modeling approach addressing all these issues at the same time has never been described nor made easily available to the community. A graphical user interface, CrunchEase, was created, which supports the user by automating the creation of input files, the running of simulations, and the extraction and comparison of data and simulation results. While a classical model considering an effective diffusion coefficient, a porosity and a solid/solution distribution coefficient (D e –ε–K D ) may be implemented in any reactive transport code, the development of CrunchEase makes it easy to apply by experimentalists without a background in reactive transport modeling. CrunchEase makes it also possible to transition more easily from a D e –ε–K D modeling approach to a state-of-the-art process-based understanding modeling approach using the full capabilities of CrunchClay, which include surface complexation modeling and a multi-porosity description of the clay packing with charged diffuse layers.

58 GEOSCIENCES↗

Ultraefficient resistance switching between charge ordered phases in 1 T -TaS 2 with a single picosecond electrical pulse

Progress in high-performance computing demands significant advances in memory technology. Among novel memory technologies that promise efficient device operation on a sub-ns timescale, resistance switching between charge ordered phases of 1T-TaS 2 has shown to be potentially useful for development of high-speed, energy efficient nonvolatile memory devices. Measurement of the electrical operation of such devices in the picosecond regime is technically challenging and hitherto still largely unexplored. Here, we use an optoelectronic “laboratory-on-a-chip” experiment for measurement of ultrafast memory switching, enabling accurate measurement of electrical switching parameters with 100 fs temporal resolution. Photoexcitation and electro-optic sampling on a (Cd,Mn)Te substrate are used to generate and, subsequently, measure electrical pulse propagation with intra-band excitation and sub-gap probing, respectively. We demonstrate high contrast nonvolatile resistance switching from high to low resistance states of a 1T-TaS 2 device using single sub-2ps electrical pulses. Using detailed modeling, we find that the switching energy density per unit area is exceptionally small, E A =9.4 fJ/μm 2 . Finally, the speed and energy efficiency of an electronic “write” process place the 1T-TaS 2 devices into a category of their own among new generation nonvolatile memory devices.

42 ENGINEERING↗

HOD-dependent systematics in Emission Line Galaxies for the DESI 2024 BAO analysis

The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we estimate the level of systematics induced in the DESI BAO analysis due the assumed Halo Occupation Distribution (HOD) model for the Emission Line Galaxy (ELG) tracer. We make use of mock galaxy catalogs constructed by fitting various HOD models to early DESI data, namely the One-Percent survey data. Our analysis includes typical HOD models for the ELG tracer used in the literature as well as extensions to the baseline models. Among the extensions, we consider various recipes for galactic conformity and assembly bias. We use 25 AbacusSummit simulations under the ΛCDM cosmology for each HOD model and perform independent analyses in Fourier space and in configuration space. To recover the BAO signal from our mocks we perform BAO reconstruction and apply the control variates technique to reduce sample variance noise. Our BAO analyses can recover the isotropic BAO parameter α iso within 0.1% and the Alcock Paczynski parameter α AP within 0.3%. Overall, we find that the systematic error due to the HOD dependence is below 0.17%, with the Fourier space analysis being more robust against the HOD systematics. We conclude that our analysis pipeline is robust enough against the HOD systematics for the ELG tracer in the DESI 2024 BAO analysis, for the assumptions made.

79 ASTRONOMY AND ASTROPHYSICS↗

Association of optical, ultraviolet, and soft X-ray excess emissions in AGNs

Soft X-ray emission (0.5–2.0 keV) plays a pivotal role in regulating the optical and ultraviolet (UV) emission in the active galactic nuclei (AGNs). We collected a sample of 1811 AGNs from the SDSS database and obtained various parameters of Balmer lines, optical continuum, Mg ii line & UV continuum and studied their dependencies on soft X-ray luminosity. Based on the linear regression analysis, we found that FWHM$_{\rm {Mg\,\,\small {II}}}$ ∝ FWHM$_{\text{H}\beta }^{0.554}$ suggesting that UV emission is arising from a region relatively outside the broad-line region (BLR) associated to the Hβ emission and found a strong correlation between optical and UV luminosities (L$_{\rm {Mg\,\,\small {II}}}$ ∝ $L_{\rm {H}\beta }^{0.822}$). It was noticed that the dependency of optical continuum luminosities on soft excess changes with the redshift (L X ∝ L$^{0.596}_{5100\, \mathring{\rm A}}$ for z < 0.5 and L X ∝ L$^{0.429}_{5100\, \mathring{\rm A}}$ for z > 0.5). The full width at half-maximum components of Hβ and Mg II core components were found to be virialized and is not affected by the soft excess emission whereas the wings of Mg II display a dependency. We estimated a relation viz. L X ∝L$^{0.520}_{3000\, \mathring{\rm A}}$ FWHM$^{0.525}_{\rm {Mg\,\,\small {II}}}$ and found to be well in agreement with a proposed physical scenario. All the derived relations were used to understand the intermodulating association of the BLR and disc in the AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep drilling in the time domain with DECam – II: characterizing the light curves of candidates in the extragalactic fields

ABSTRACT In this second paper on the DECam deep-drilling field (DDF) program, we release 2020 optical $gri$-band light curves for transients and variables in the extragalactic COSMOS and ELAIS fields based on time series observations with a 3-d cadence from semester 2021A through 2023A. In order to demonstrate the wide variety of time domain events detected by the program and encourage others to use the data set, we characterize the sample by presenting a brief analysis of the light-curve parameters such as time span, amplitude, and peak brightness. We also present preliminary light-curve categorizations, and identify potential stellar variables, active galactic nuclei, tidal disruption events, supernovae (SNe) (such as Type Ia, Type IIP, superluminous, and gravitationally lensed SNe), and fast transients. Where relevant, the number of identified transients is compared to the predictions of the original proposal. We also discuss the challenges of analyzing DDF data in the context of the upcoming Vera C. Rubin Observatory and its Legacy Survey of Space and Time, which will include DDFs. Images from the Dark Energy Camera DDF program are available without proprietary period and the light curves presented in this work are publicly available for analysis.

Graham, Melissa L. (ORCID:0000000291543136)↗

Gamma-ray bursts, supernovae Ia, and baryon acoustic oscillations: A binned cosmological analysis

Abstract Cosmological probes at any redshift are necessary to reconstruct consistently the cosmic history. Studying properly the tension on the Hubble constant, H0, obtained by supernovae type Ia (SNe Ia) and the Planck measurements of the cosmic microwave background radiation would require complete samples of distance indicators at any epoch. Gamma-ray bursts (GRBs) are necessary for the aforementioned task because of their huge luminosity that allows us to extend the cosmic ladder to very high redshifts. However, using GRBs alone as standard candles is challenging, because their luminosity varies widely. To this end, we choose a reliable correlation for GRBs with a very small intrinsic scatter: the so-called fundamental plane correlation for GRB afterglows corrected for selection biases and redshift evolution. We choose a well defined sample: the platinum sample, composed of 50 long GRBs. To further constrain the cosmological parameters, we use baryon acoustic oscillations (BAOs) given their reliability as standard rulers. Thus, we have applied GRBs, SNe Ia, and BAOs in a binned analysis in redshifts so that the GRB contribution is fully included in the last redshift bin, which reaches z = 5. We use the fundamental plane correlation (also known as the 3D Dainotti relation), together with SNe Ia and BAOs, to constrain H0 and the density matter today, ΩM. This methodology allows us to assess the role of GRBs combined with SNe Ia and BAOs. We have obtained results for H0 and ΩM using GRBs+SNe Ia+BAOs with better precision than SNe Ia alone for every bin, thus confirming the beneficial role of BAOs and GRBs added together. In addition, consistent results between GRBs+SNe Ia+BAOs are obtained when compared with SNe Ia+BAOs, showing the importance of GRBs since the distance ladder is extended up to z = 5 with a similar precision obtained with other probes without including GRBs.

Dainotti, Maria Giovanna↗

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Observation of Zero-Field Transverse Resistance in AlO x / SrTi O 3 Interface Devices

Domain walls in AlO x =SrTiO 3 (AlO x =STO) interface devices at low temperatures give a rise to a new signature in the electrical transport of two-dimensional carrier gases formed at the surfaces or interfaces of STO-based heterostructures: a finite transverse resistance observed in Hall bars in zero external magnetic field. This transverse resistance depends on the local domain wall configuration and hence changes with temperature, gate voltage, thermal cycling, and position along the sample and can even change sign as a function of these parameters. The transverse resistance is observed below ≃70 K but grows and changes significantly below ≃40 K, the temperature at which the domain walls become increasingly polar. Surprisingly, the transverse resistance is much larger in (111) oriented heterostructures in comparison to (001) oriented heterostructures. Measurements of the capacitance between the conducting interface and an electrode applied to the substrate, which reflect the dielectric constant of the STO, indicate that this difference may be related to the greater variation of the temperature-dependent dielectric constant with electric field when the electric field is applied in the [111] direction. Further, the finite transverse resistance can be explained inhomogeneous current flow due to the preferential transport of current along domain walls that are askew to the nominal direction of the injected current.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Strong and Weak C P Tests in Sequential Decays of Polarized Σ 0 Hyperons

The J / ψ , ψ ( 3686 ) → Σ 0 Σ ¯ 0 processes and subsequent decays are studied using the world’s largest J / ψ and ψ ( 3686 ) data samples collected with the BESIII detector. The parity-violating decay parameters of the decays Σ 0 → Λ γ and Σ ¯ 0 → Λ ¯ γ , α Σ 0 = − 0.0017 ± 0.0021 ± 0.0018 and α ¯ Σ 0 = 0.0021 ± 0.0020 ± 0.0022 , are measured for the first time. The strong C P symmetry is tested in the decays of the Σ 0 hyperons for the first time by measuring the asymmetry A C P Σ = α Σ 0 + α ¯ Σ 0 = ( 0.4 ± 2.9 ± 1.3 ) × 10 − 3 . The weak C P test is performed in the subsequent decays of their daughter particles Λ and Λ ¯ . Also for the first time, the transverse polarizations of the Σ 0 hyperons in J / ψ and ψ ( 3686 ) decays are observed with opposite directions, and the ratios between the S -wave and D -wave contributions of the J / ψ , ψ ( 3686 ) → Σ 0 Σ ¯ 0 decays are obtained. These results are crucial to understand the decay dynamics of the charmonium states and the production mechanism of the Σ 0 − Σ ¯ 0 pairs. Published by the American Physical Society 2024

Physics↗

Magnetism in metastable and annealed compositionally complex alloys

Compositionally complex materials (CCMs) present a potential paradigm shift in the design of magnetic materials. These alloys exhibit long-range structural order coupled with limited or no chemical order. As a result, extreme local environments exist with a large variations in the magnetic energy terms, which can manifest large changes in the magnetic behavior. In the current work, the magnetic properties of (Cr, Mn, Fe, Ni) alloys are presented. These materials were prepared by room-temperature combinatorial sputtering, resulting in a range of compositions with a single bcc structural phase and no chemical ordering. The combinatorial growth technique allows CCMs to be prepared outside of their thermodynamically stable phase, enabling the exploration of otherwise inaccessible order. The mixed ferromagnetic and antiferromagnetic interactions in these alloys causes frustrated magnetic behavior, which results in an extremely low coercivity (<1mT), which increases rapidly at 50 K. At low temperatures, the coercivity achieves values of nearly 500 mT, which is comparable to some high-anisotropy magnetic materials. Further, commensurate with the divergent coercivity is an atypical drop in the temperature dependent magnetization. These effects are explained by a mixed magnetic phase model, consisting of ferro-, antiferro-, and frustrated magnetic regions, and are rationalized by simulations. A machine-learning algorithm is employed to visualize the parameter space and inform the development of subsequent compositions. Annealing the samples at 600 °C orders the sample, more-than doubling the Curie temperature and increasing the saturation magnetization by as much as 5×. Simultaneously, the large coercivities are suppressed, resulting in magnetic behavior that is largely temperature independent over a range of 350 K. The ability to transform from a hard magnet to a soft magnet over a narrow temperature range makes these materials promising for heat-assisted recording technologies.

36 MATERIALS SCIENCE↗

Autoregressive neural quantum states of Fermi Hubbard models

Neural quantum states (NQSs) have emerged as a powerful ansatz for variational quantum Monte Carlo studies of strongly correlated systems. Here, we apply recurrent neural networks (RNNs) and autoregressive transformer neural networks to the Fermi-Hubbard and the (non-Hermitian) Hatano-Nelson-Hubbard models in one and two dimensions. In both cases, we observe that the convergence of the RNN ansatz is challenged when increasing the interaction strength. We present a physically motivated and easy-to-implement strategy for improving the optimization, namely, by ramping of the model parameters. Furthermore, we investigate the advantages and disadvantages of the autoregressive sampling property of both network architectures. For the Hatano-Nelson-Hubbard model, we identify convergence issues that stem from the autoregressive sampling scheme in combination with the non-Hermitian nature of the model. Our findings provide insights into the challenges of the NQS approach and make the first step towards exploring strongly correlated electrons using this ansatz. Published by the American Physical Society 2025

Ibarra-García-Padilla, Eduardo (ORCID:000000019165↗

A science-driven approach to optimize the design for a biological small-angle neutron scattering instrument

Biological small-angle neutron scattering (SANS) instruments facilitate critical analysis of the structure and dynamics of complex biological systems. However, with the growth of experimental demands and the advances in optical systems design, a new neutron optical concept is necessary to overcome the limitations of current instruments. This work presents an approach to include experimental objectives ( i.e. the science to be supported by a specific neutron scattering instrument) in the optimization of the neutron optical concept. The approach for a proposed SANS instrument at the Second Target Station of the Spallation Neutron Source at Oak Ridge National Laboratory, USA, is presented here. Further, the instrument is simulated with the McStas software package. The optimization process is driven by an evolutionary algorithm using McStas output data, which are processed to calculate an objective function designed to quantify the expected performance of the simulated neutron optical configuration for the intended purpose. Each McStas simulation covers the complete instrument, from source to detector, including realistic sample scattering functions. This approach effectively navigates a high-dimensional parameter space that is otherwise intractable; it allows the design of next-generation SANS instruments to address specific scientific cases and has the potential to increase instrument performance compared with traditional design approaches.

47 OTHER INSTRUMENTATION↗

Improved joint X-ray and neutron refinement procedure in Phenix

Neutron diffraction is one of the three crystallographic techniques (X-ray, neutron and electron diffraction) used to determine the atomic structures of molecules. Its particular strengths derive from the fact that H (and D) atoms are strong neutron scatterers, meaning that their positions, and thus protonation states, can be derived from crystallographic maps. However, because of technical limitations and experimental obstacles, the quality of neutron diffraction data is typically much poorer (completeness, resolution and signal to noise) than that of X-ray diffraction data for the same sample. Further, refinement is more complex as it usually requires additional parameters to describe the H (and D) atoms. The increase in the number of parameters may be mitigated by using the `riding hydrogen' refinement strategy, in which the positions of H atoms without a rotational degree of freedom are inferred from their neighboring heavy atoms. However, this does not address the issues related to poor data quality. Therefore, neutron structure determination often relies on the presence of an X-ray data set for joint X-ray and neutron (XN) refinement. In this approach, the X-ray data serve to compensate for the deficiencies of the neutron diffraction data by refining one model simultaneously against the X-ray and neutron data sets. To be applicable, it is assumed that both data sets are highly isomorphous, and preferably collected from the same crystals and at the same temperature. However, the approach has a number of limitations that are discussed in this work by comparing four separately re-refined neutron models. To address the limitations, a new method for joint XN refinement is introduced that optimizes two different models against the different data sets. This approach is tested using neutron models and data deposited in the Protein Data Bank. The efficacy of refining models with H atoms as riding or as individual atoms is also investigated.

36 MATERIALS SCIENCE↗

ARC-DPA 1-MeV Neutron Fluence Equivalence Model for GaAs

Here, the legacy data used in the production of the American Society for Testing and Materials (ASTM) E722 standard for gallium arsenide (GaAs) were fit with the athermal recombination corrected-displacements per atom (ARC-DPA) model for the development of a novel 1-MeV neutron equivalent fluence metric. Improving on previous work, the coefficients for the ARC-DPA model were optimized using a conjugate gradient method, and the uncertainties in the observed damage in different benchmark neutron fields were sampled to quantify their contribution to the damage metric. The optimized parameters of −0.474 ± 0.03 for b arcdpa and 0.016 ± 0.002 for c arcdpa provided good agreement and indicated that the solution is not sensitive to the corresponding fluence uncertainties. Compared to the current 1-MeV neutron fluence equivalence standard for displacement damage in GaAs, our results indicate that for a thermal reactor neutron environment, displacement damage has been underestimated by around 25% and that the saturation value for the damage efficiency had been higher than previously modeled.

1-MeV equivalent fluence↗

hIPPYlib-MUQ: A Bayesian Inference Software Framework for Integration of Data with Complex Predictive Models under Uncertainty

Bayesian inference provides a systematic framework for integration of data with mathematical models to quantify the uncertainty in the solution of the inverse problem. However, the solution of Bayesian inverse problems governed by complex forward models described by partial differential equations (PDEs) remains prohibitive with black-box Markov chain Monte Carlo (MCMC) methods. We present hIPPYlib-MUQ, an extensible and scalable software framework that contains implementations of state-of-the art algorithms aimed to overcome the challenges of high-dimensional, PDE-constrained Bayesian inverse problems. These algorithms accelerate MCMC sampling by exploiting the geometry and intrinsic low-dimensionality of parameter space via derivative information and low rank approximation. The software integrates two complementary open-source software packages, hIPPYlib and MUQ. hIPPYlib solves PDE-constrained inverse problems using automatically-generated adjoint-based derivatives, but it lacks full Bayesian capabilities. MUQ provides a spectrum of powerful Bayesian inversion models and algorithms, but expects forward models to come equipped with gradients and Hessians to permit large-scale solution. By combining these two complementary libraries, we created a robust, scalable, and efficient software framework that realizes the benefits of each and allows us to tackle complex large-scale Bayesian inverse problems across a broad spectrum of scientific and engineering disciplines. To illustrate the capabilities of hIPPYlib-MUQ, we present a comparison of a number of MCMC methods available in the integrated software on several high-dimensional Bayesian inverse problems. These include problems characterized by both linear and nonlinear PDEs, various noise models, and different parameter dimensions. The results demonstrate that large (~ 50×) speedups over conventional black box and gradient-based MCMC algorithms can be obtained by exploiting Hessian information (from the log-posterior), underscoring the power of the integrated hIPPYlib-MUQ framework.

97 MATHEMATICS AND COMPUTING↗

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part II: Idealized Markov Chain Monte Carlo Experiments

Observationally informed development of a new framework for bulk rain microphysics, the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS; described in Part I of this study), is demonstrated. This scheme’s development is motivated by large uncertainties in cloud and weather simulations associated with approximations and assumptions in existing microphysics schemes. Here, a proof-of-concept study is introduced using a Markov chain Monte Carlo sampling algorithm with BOSS to probabilistically estimate microphysical process rates and parameters directly from a set of synthetically generated rain observations. The framework utilized is an idealized steady-state one-dimensional column rainshaft model with specified column-top rain properties and a fixed thermodynamical profile. Different configurations of BOSS—flexibility being a key feature of this approach—are constrained via synthetic observations generated from a traditional three-moment bulk microphysics scheme. The ability to retrieve correct parameter values when the true parameter values are known is illustrated. For cases when there is no set of true parameter values, the accuracy of configurations of BOSS that have different levels of complexity is compared. It is found that addition of the sixth moment as a prognostic variable improves prediction of the third moment (proportional to bulk rain mass) and rain rate. In contrast, increasing process rate formulation complexity by adding more power terms has little benefit—a result that is explained using further-idealized experiments. BOSS rainshaft simulations are shown to well estimate the true process rates from constraint by bulk rain observations, with the further benefit of rigorously quantified uncertainty of these estimates.

58 GEOSCIENCES↗

Identification of Distribution of Elastic Constants of SiC/SiC Braided Tubes

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiCf/SiCm) composite tubes due to their superior mechanical properties. There are many parameters in the manufacturing of SiCf/SiCm composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. A sensitivity analysis will be carried out to determine the dominant parameters among the abovementioned ones. We intend to determine the variability in elastic constants of SiCf/SiCm composite tubes considering the variability in dominant parameters. A chosen number of combinations of dominant variables will be sampled through Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in circumferential and longitudinal direction, Poisson's ratio and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used a surrogate to FE analysis. once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The variability can be calculated based on MCS-generated data which will be presented at the conference.

42 ENGINEERING↗