Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model prediction”

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

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

At least 1,405 records · Page 78

Data-driven models of nonautonomous systems

Nonautonomous dynamical systems are characterized by time-dependent inputs, which complicates the discovery of predictive models describing the spatiotemporal evolution of the state variables of quantities of interest from their temporal snapshots. When dynamic mode decomposition (DMD) is used to infer a linear model, this difficulty manifests itself in the need to approximate the time-dependent Koopman operators. Our approach is to approximate the original nonautonomous system with a modified system derived via a local parameterization of the time-dependent inputs. The modified system comprises a sequence of local parametric systems, which are subsequently approximated by a parametric surrogate model using the DRIPS (dimension reduction and interpolation in parameter space) framework. The offline step of DRIPS relies on DMD to build a linear surrogate model, endowed with reduced-order bases for the observables mapped from training data. The online step interpolates on suitable manifolds to construct a sequence of iterative parametric surrogate models; the target/test parameter points on these manifolds are specified by a local parameterization of the test time-dependent inputs. Here, we use numerical experimentation to demonstrate the robustness of our method and compare its performance with that of deep neural networks.

97 MATHEMATICS AND COMPUTING↗

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗

Assessing the predictive impact of factor fixing with an adaptive uncertainty-based approach .

Despite widespread use of factor fixing in environmental modeling, its effect on model predictions has received little attention and is instead commonly presumed to be negligible. We propose a proof-of-concept adaptive method for systematically investigating the impact of factor fixing. The method uses Global Sensitivity Analysis methods to identify groups of sensitive parameters, then quantifies which groups can be safely fixed at nominal values without exceeding a maximum acceptable error, demonstrated using the 21-dimensional Sobol’ G-function. Furthermore, three error measures are considered for quantities of interest, namely Relative Mean Absolute Error, Pearson Product-Moment Correlation and Relative Variance. Results demonstrate that factor fixing may cause large errors in the model results unexpectedly, when preliminary analysis suggests otherwise, and that the default value selected affects the number of factors to fix. To improve the applicability and methodological development of factor fixing, a new research agenda encompassing five opportunities is discussed for further attention.

54 ENVIRONMENTAL SCIENCES↗

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗

A MPET 2 -mPBPK model for subcutaneous injection of biotherapeutics with different molecular weights: From local scale to whole-body scale

Subcutaneous injection of biotherapeutics has attracted considerable attention in the pharmaceutical industry. However, there is limited understanding of the mechanisms underlying the absorption of drugs with different molecular weights and the delivery of drugs from the injection site to the targeted tissue. Here, we propose the MPET 2 -mPBPK model to address this issue. This multiscale model couples the MPET 2 model, which describes subcutaneous injection at the local tissue scale from a biomechanical view, with a post-injection absorption model at injection site and a minimal physiologically-based pharmacokinetic (mPBPK) model at whole-body scale. Utilizing the principles of tissue biomechanics and fluid dynamics, the local MPET 2 model provides solutions that account for tissue deformation and drug absorption in local blood vessels and initial lymphatic vessels during injection. Additionally, we introduce a model accounting for the molecular weight effect on the absorption by blood vessels, and a nonlinear model accounting for the absorption in lymphatic vessels. The post-injection model predicts drug absorption in local blood vessels and initial lymphatic vessels, which are integrated into the whole-body mPBPK model to describe the pharmacokinetic behaviors of the absorbed drug in the circulatory and lymphatic system. We establish a numerical model which links the biomechanical process of subcutaneous injection at local tissue scale and the pharmacokinetic behaviors of injected biotherapeutics at whole-body scale. With the help of the model, we propose an explicit relationship between the reflection coefficient and the molecular weight and predict the bioavalibility of biotherapeutics with varying molecular weights via subcutaneous injection. The considered drug absorption mechanisms enable us to study the differences in local drug absorption and whole-body drug distribution with varying molecular weights. This model enhances the understanding of drug absorption mechanisms and transport routes in the circulatory system for drugs of different molecular weights, and holds the potential to facilitate the application of computational modeling to drug formulation.

59 BASIC BIOLOGICAL SCIENCES↗

Examining Tropical Convection Features at Storm-Resolving Scales Over the Maritime Continent Region

Global Storm Resolving Models (GSRMs) provide a way to understand weather and climate events across scales for better-informed climate impacts. In this work, we apply the recently developed and validated CAM (Community Atmosphere Model)—MPAS (Model for Prediction Across Scales) modeling framework, based on the open-source Community Earth System Model (CESM2), to examine the tropical convection features at the storm resolving scale over the Maritime Continent region at 3 km horizontal spacing. We target two global numerical experiments during the winter season of 2018 for comparison with observation in the region. Here we focus on the investigation of the representations of the convective systems, precipitation statistics, and tropical cyclone behaviors. We found that regional-refined experiments show more accurate precipitation distributions, diurnal cycles, and better agreement with observations for tropical cyclone features in terms of intensity and strength statistics. We expect the exploration of this work will further advance the development and use of the storm-resolving model in precipitation predictions across scales.

54 ENVIRONMENTAL SCIENCES↗

Reconstruction and identification of boosted di-τ systems in a search for Higgs boson pairs using 13 TeV proton-proton collision data in ATLAS

In this paper, a new technique for reconstructing and identifying hadronically decaying τ+τ- pairs with a large Lorentz boost, referred to as the di-τ tagger, is developed and used for the first time in the ATLAS experiment at the Large Hadron Collider. A benchmark di-τ tagging selection is employed in the search for resonant Higgs boson pair production, where one Higgs boson decays into a boosted $b\overline b$ pair and the other into a boosted τ+τ- pair, with two hadronically decaying τ-leptons in the final state. Using 139 fb -1 of proton-proton collision data recorded at a centre-of-mass energy of 13 TeV, the efficiency of the di-τ tagger is determined and the background with quark- or gluon-initiated jets misidentified as di-τ objects is estimated. The search for a heavy, narrow, scalar resonance produced via gluon-gluon fusion and decaying into two Higgs bosons is carried out in the mass range 1–3 TeV using the same dataset. No deviations from the Standard Model predictions are observed, and 95% confidence-level exclusion limits are set on this model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measuring thermomechanical response of large-format printed polymer composite structures via digital image correlation

Large-format additive manufacturing (LFAM) is a branch of additive manufacturing (AM) research with the ability to create large structures typically measuring several meters in scale. LFAM is advantageous for tooling applications, not only because it offers the ability to create complex geometries not easily made using subtractive manufacturing processes, but the cost savings of pelletized feedstock used by these systems result in larger parts printed at faster speeds than traditional AM systems. Fiber reinforced polymer (FRP) is a commonly used feedstock material in LFAM structures because it reduces the distortion experienced during printing. However, FRP introduces highly anisotropic thermomechanical properties and contributes to a nonhomogeneous microstructure that can result in critical distortion of dimensions during tooling. Measuring the global thermomechanical response of LFAM structures requires a more representative method that accounts for not only anisotropic properties but also the nonhomogeneous nature of the final part. This is where traditional techniques to measure thermomechanical response, such as thermomechanical analysis (TMA), fall short as they assume homogeneity. This study evaluated the coefficient of thermal expansion (CTE) of LFAM structures as measured by TMA as compared to a novel digital image correlation oven (DIC Oven) system. The LFAM structures were made from 20 % by weight carbon fiber reinforced acrylonitrile butadiene styrene (CF-ABS). TMA measurements showed significant variations in CTE across a single LFAM bead, confirming the need for a global technique that captures overall thermomechanical response. The CTE values measured using the DIC Oven compared well to average TMA values obtained from localized measurements across the sample. The DIC Oven was also used to quantify the effects of different layer orientations on thermomechanical properties, which cannot be easily captured using TMA. A predictive model was also developed by using localized TMA values across an LFAM bead to predict the overall thermomechanical response of an LFAM structure.

36 MATERIALS SCIENCE↗

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗

Search for non-resonant production of semi-visible jets using Run 2 data in ATLAS

Semi-visible jets, with a significant contribution to the event's missing transverse momentum, can arise in strongly interacting dark sectors. This results in an event topology where one of the jets can be aligned with the direction of the missing transverse momentum. The first search for semi-visible jets produced via a t-channel mediator exchange is presented. The analysis uses proton-proton collisions with an integrated luminosity of 139 fb —1 and a centre-of-mass energy of 13 TeV, collected with the ATLAS detector during the Run 2 of the LHC. No excess over Standard Model predictions is observed. Assuming a coupling strength of unity between the mediator, a Standard Model quark and a dark quark, mediator masses up to 2.7 TeV are excluded at the 95% confidence level. Upper limits on the coupling strength are also derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Resolving dynamic mineral-organic interactions in the rhizosphere by combining in-situ microsensors with plant-soil reactive transport modeling

Associations between minerals and organic matter represent one of the most important carbon storage mechanisms in soils. Plant roots are major sources of soil carbon, and resolving the dynamics and dominance of microbial consumption versus mineral sorption of root-derived carbon is critical to understanding soil carbon storage. Here we integrate in-situ rhizosphere microsensor and plant physiological measurements with a 3-D plant-soil reactive transport model to explore the fate of dissolved organic carbon (DOC) in the rhizosphere, particularly its microbial consumption and interaction with Fe oxide minerals. Over several days, a microdialysis probe sampling pore water at the root-soil interface of growing Vicia faba roots in live soil, revealed clear diel patterns of DOC concentration. Daytime DOC spikes coincided with peaks in leaf-level photosynthesis rates and were accompanied by declining redox potential and dissolved oxygen as well as increasing pH in the rhizosphere. Incorporating microsensor data into our modeling framework showed that the measured rapid loss of DOC after each mid-day spike could not be explained by consumption via aerobic respiration, nor via anaerobic respiration dominated by Fe oxide reduction. Rather, in the model, a large fraction of rhizosphere DOC was rapidly immobilized each day by adsorption to Fe oxides. Further, modeled microbial Fe reduction (fueled by DOC) did not mobilize significant organic carbon from Fe oxides during the day. Instead, the model predicted equilibrium desorption of organic carbon from Fe oxides at night. This new mechanistic modeling framework, which couples aboveground plant physiological measurements with non-destructive high-resolution monitoring of rhizosphere processes, has great potential for exploring the dynamics and balance of the various microbial reactions and mineral interactions controlling carbon transformations and storage in soils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning for endoleak detection after endovascular aortic repair

Diagnosis of endoleak following endovascular aortic repair (EVAR) relies on manual review of multi-slice CT angiography (CTA) by physicians which is a tedious and time-consuming process that is susceptible to error. We evaluate the use of a deep neural network for the detection of endoleak on CTA for post-EVAR patients using a novel data efficient training approach. 50 CTAs and 20 CTAs with and without endoleak respectively were identified based on gold standard interpretation by a cardiovascular subspecialty radiologist. The Endoleak Augmentor, a custom designed augmentation method, provided robust training for the machine learning (ML) model. Predicted segmentation maps underwent post-processing to determine the presence of endoleak. The model was tested against 3 blinded general radiologists and 1 blinded subspecialist using a held-out subset (10 positive endoleak CTAs, 10 control CTAs). Model accuracy, precision and recall for endoleak diagnosis were 95%, 90% and 100% relative to reference subspecialist interpretation (AUC = 0.99). Accuracy, precision and recall was 70/70/70% for generalist1, 50/50/90% for generalist2, and 90/83/100% for generalist3. The blinded subspecialist had concordant interpretations for all test cases compared with the reference. In conclusion, our ML-based approach has similar performance for endoleak diagnosis relative to subspecialists and superior performance compared with generalists.

60 APPLIED LIFE SCIENCES↗

Partial-wave decomposition of the Keldysh ionization amplitude

We present an alternative way of calculating the Keldysh amplitude, i.e., the length-gauge form of the ionization amplitude in the strong-field approximation. The amplitude is evaluated exactly by expanding it in Fourier components and partial waves. Comparisons of the semianalytic model predictions with results of ab initio numerical simulations of the time-dependent Schrödinger equation for the interaction of electrons in short-range potentials with intense laser light yield excellent agreement, for wavelengths from the single photon to the multiphoton to the tunneling regime. Specifically, for ionization from initial states with higher angular momentum quantum number, e.g., p states, a significant improvement over predictions based on the popular saddle-point approximation is found. Furthermore, the current model rate allows for interpretation of the strong-field ionization process in terms of multiphoton absorption pathways and angular momentum selection rules.

74 ATOMIC AND MOLECULAR PHYSICS↗

Differentiable Predictive Control with Safety Guarantees: A Control Barrier Function Approach

In this paper, we develop a novel form of differentiable predictive control (DPC) with safety and robustness guarantees. DPC is a form of approximate model predictive control (MPC), wherein the control policy is a neural network that learns a receding horizon, optimal control law. The proposed approach exploits a new form of sampled-data barrier function to enforce safety, while only interrupting the neural network-based controller near the boundary of the safe set. The effectiveness of the proposed approach is demonstrated in simulation.

Shaw Cortez, Wenceslao E.↗

Occlusion-Perturbed Deep Learning for Probabilistic Solar Forecasting via Sky Images

Solar forecasting is shifting to the probabilistic paradigm due to the inherent uncertainty within the solar resource. Input uncertainty quantification is one of the widely used and best-performing ways to model solar uncertainty. However, compared to other sources of inputs, such as numerical weather prediction models, pure sky image-based probabilistic solar forecasting lags behind. In this research, an occlusion-perturbed convolutional neural network, named the PSolarNet, is developed. The PSolarNet provides very short-term deterministic forecasts, forecast scenarios, and probabilistic forecasts of the global horizontal irradiance from sky image sequences. Case studies based on 6 years of open-source data show that the developed PSolarNet is able to generate accurate 10-minute ahead deterministic forecasts with a 5.62% normalized root mean square error, realistic and diverse forecast scenarios with a 0.966 average correlation with the actual time series, and reliable and sharp probabilistic forecasts with a 2.77% normalized continuous ranked probability score.

Bayesian model averaging↗

Learned implicit representations of aerosol chemistry and physics for enhancing the predictability of water cycle extreme events

Focal Area(s): Focal area 2: Predictive modeling through the use of AI-derived model components; and Focal area 3: physics-guided AI. Science Challenge: Comprehensive models of many geophysical processes require a large number of state variables, disqualifying them for inclusion in Earth System Models (ESMs) for the foreseeable future. This white paper proposes the paradigm shift from the use of state variables representing explicit physical quantities to the use of machine-learning to create and integrate a compact, implicit representation of a physical system with a smaller number of state variables. As an initial application, we propose to create parsimonious machine-learned surrogate models of gas- and aerosol-phase chemistry and physics with the purpose of improving the accuracy of the representation of cloud formation and cloud-aerosol interaction in the E3SM model. This type of improved representation of cloud microphysics is critical for enhancing the predictability of water cycle extremes.

54 ENVIRONMENTAL SCIENCES↗

Low-Cost, High-Performance Carbon Fiber for Compressed Natural Gas Storage Tanks (Final Technical Report – Down Select Report)

The aim of this project is to reduce the cost of Type IV, carbon fiber (CF) composite overwrap compressed gas storage tanks by reducing the cost of CF and CF composites. The project team worked to reduce the cost of CF by exploring and testing opportunities for a low-cost alternative precursor material for CF production to supplant market-dominant and costly polyacrylonitrile (PAN). Concurrently, the team aimed to reduce the cost of the tanks at the composite level by improving the interfacial adhesion between the fibers and the matrix resin through the incorporation of low-cost nanoparticles recycled from waste materials, which would reduce the volume of costly CF required to achieve the same tank performance. At the end of the first year, the project team selected mesophase pitch as the primary precursor candidate from a field of materials based on the superior mechanical performance and cost-saving potential. During the second year, the team produced CFs derived from mesophase pitch achieving an average tensile strength of 365.6 ksi and average tensile modulus of 40.74 Msi. Facility availability for spinning and converting these fibers at greater scale has hindered scale-up demonstration, but the team has identified opportunities to conduct this work in the near term. Cost modeling shows that these mesophase pitch-derived CFs can be up to 40% less expensive than PAN-derived CFs due to the lower cost of the feedstock material, higher throughput, greater conversion yield, and lower cost spinning method and compared to PAN. Additionally, the team has demonstrated at lab-scale that nanoparticle coating CFs can significantly increase the interfacial shear strength and load transfer efficiency of CFs in a matrix. Single filament pull-out testing showed a 27% average increase in max interfacial shear strength due to this coating. A continuous method of applying these coatings to a tow of CF has been developed for scale-up. 26 m tows of coated CFs were produced using this system and formed into composite ring samples for ASTM ring burst testing. Issues with the testing protocol have limited assessment of these results. A prototype Type IV tank was designed to meet ANSI HGV2 standards, and the design criteria set out by DOE, using the CF properties developed by the team paired with a proprietary resin matrix, a polyamide liner, and aluminum end bosses. The tank weighs 153.1 kg and with a total capacity of 5.8 kg H2 (5.6 kg usable), which yields a gravimetric capacity of 1.17 kWh/kg. Cost modeling predicts that the tank will have a projected cost of $15.73/kWh. Tank performance modeling does not include considerations for fiber-matrix load transfer efficiency improvements offered by nanoparticle coating method.

08 HYDROGEN↗