Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Implementation of Improved Transverse Shear Calculations and Higher Order Laminate Theory Into Strain Rate Dependent Analyses of Polymer Matrix Composites

A numerical procedure has been developed to investigate the nonlinear and strain rate dependent deformation response of polymer matrix composite laminated plates under high strain rate impact loadings. A recently developed strength of materials based micromechanics model, incorporating a set of nonlinear, strain rate dependent constitutive equations for the polymer matrix, is extended to account for the transverse shear effects during impact. Four different assumptions of transverse shear deformation are investigated in order to improve the developed strain rate dependent micromechanics model. The validities of these assumptions are investigated using numerical and theoretical approaches. A method to determine through the thickness strain and transverse Poisson's ratio of the composite is developed. The revised micromechanics model is then implemented into a higher order laminated plate theory which is modified to include the effects of inelastic strains. Parametric studies are conducted to investigate the mechanical response of composite plates under high strain rate loadings. Results show the transverse shear stresses cannot be neglected in the impact problem. A significant level of strain rate dependency and material nonlinearity is found in the deformation response of representative composite specimens.

Zhu, Lin-Fa↗
DOE OSTI2024

Multiscale Modeling of Silicon Carbide Cladding for Nuclear Applications: Thermal Performance Modeling

The complex multiscale and anisotropic nature of silicon carbide (SiC) ceramic matrix composite (CMC) makes it difficult to accurately model its performance in nuclear applications. The existing models for nuclear grade composite SiC do not account for the microstructural features and how these features can affect the thermal and structural behavior of the cladding and its anisotropic properties. In addition to the microstructural features, the properties of individual constituents of the composites and fiber tow architecture determine the bulk properties. Models for determining the relationship between the individual constituents’ properties and the bulk properties of SiC composites for nuclear applications are absent, although empirical relationships exist in the literature. Here, a hierarchical multiscale modeling approach was presented to address this challenge. This modular approach addressed this difficulty by dividing the various aspects of the composite material into separate models at different length scales, with the evaluated property from the lower-length-scale model serving as an input to the higher-length-scale model. The multiscale model considered the properties of various individual constituents of the composite material (fiber, matrix, and interphase), the porosity in the matrix, the fiber volume fraction, the composite architecture, the tow thickness, etc. By considering inhomogeneous and anisotropic contributions intrinsically, our bottom-up multiscale modeling strategy is naturally physics-informed, bridging constitutive law from micromechanics to meso-mechanics and structural mechanics. The effects that these various physical attributes and thermo-physical properties have on the composite’s bulk thermal properties were easily evaluated and demonstrated through the various analyses presented herein. Since silicon carbide fiber-reinforced SiC CMCs are also promising thermal–structural materials with a broad range of high-end technology applications beyond nuclear applications, we envision that the multiscale modeling method we present here may prove helpful in future efforts to develop and construct reinforced CMCs and other advanced composite nuclear materials, such as MAX phase materials, that can service under harsh environments of ultrahigh temperatures, oxidation, corrosion, and/or irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗
NASA NTRS1985

An internal variable assisted formulation of moisture sorption models for polymers and composites

The purpose of the paper is to demonstrate that constitutive equations for moisture sorption and swelling can be obtained with an internal variable method which is rooted in the theory of irreversible thermodynamics. Realism is injected into the model through the concept of free volume. The power of the method is illustrated by the solution of the formidable nonlinear boundary value problem for diffusion with swelling-stress dependent diffusivity.

Hiel, C. C.↗
DOE OSTI2024

Laboratory Upgrade Point Absorber WEC-Sim Model with MoorDyn Moorings

This dataset includes a WEC-Sim and MoorDyn model of the Laboratory Upgrade Point Absorber (LUPA). LUPA is an open-source wave energy converter designed and tested by Oregon State University. The files provided here constitute a stable LUPA configuration with three mooring lines. This model is 1/20 scale, optimized for the O.H. Hinsdale Wave Lab at Oregon State University. This model of LUPA adds MoorDyn functionality for more accurate mooring predictions and uses a more stable, updated version of LUPA's current physical configuration. This model is for WEC-Sim Version 6.0. A recent update of WEC-Sim has changed some functionality of MoorDyn such that this model will not work with WEC-Sim Version 6.1.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

The Ghent Hybrid model in NuWro: a new neutrino single-pion production model in the GeV regime

Neutrino-induced single-pion production constitutes an essential interaction channel in modern neutrino oscillation experiments, with its products building up a significant fraction of the observable hadronic final states. Frameworks of oscillation analyses strongly rely on Monte Carlo neutrino event generators, which provide theoretical predictions of neutrino interactions on nuclear targets. Thus, it is crucial to integrate state-of-the-art single-pion production models with Monte Carlo simulations to prepare for the upcoming systematics-dominated landscape of neutrino measurements. In this work, we present the implementation of the Ghent Hybrid model for neutrino-induced single-pion production in the NuWro Monte Carlo event generator. The interaction dynamics includes coherently-added contributions from nucleon resonances and a non-resonant background, merged into the pythia branching predictions in the deep-inelastic regime, as instrumented by NuWro. This neutrino-nucleon interaction model is fully incorporated into the nuclear framework of the generator, allowing it to account for the influence of both initial- and final-state nuclear medium effects. We compare the predictions of this integrated implementation with recent pion production data from accelerator-based neutrino experiments. The results of the novel model show improved agreement of the generator predictions with the data and point to the significance of the refined treatment of the description of pion-production processes beyond the ∆ region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The Structural Basis of T4 Phage Lysis Control: DNA as the Signal for Lysis Inhibition

Optimal phage propagation depends on the regulation of the lysis of the infected host cell. In T4 phage infection, lysis occurs when the holin protein (T) forms lesions in the host membrane. However, the lethal function of T can be blocked by an antiholin (RI) during lysis inhibition (LIN). LIN sets if the infected cell undergoes superinfection, then the lysis is delayed until host/phage ratio becomes more favorable for the release of progeny. It has been thought that a signal derived from the superinfection is required to activate RI. Here we report structures that suggest a radically different model in which RI binds to T irrespective of superinfection, causing it to accumulate in a membrane as heterotetrameric 2RI–2T complex. Moreover, we show the complex binds non-specifically to DNA, suggesting that the gDNA from the superinfecting phage serves as the LIN signal and that stabilization of the complex by DNA binding is what defines LIN. Finally, we show that soluble domain of free RI crystallizes in a domain-swapped homotetramer, which likely works as a sink for RI molecules released from the RI–T complex to ensure efficient lysis. These results constitute the first structural basis and a new model not only for the historic LIN phenomenon but also for the temporal regulation of phage lysis in general.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

An asymmetric electron-scattering photosphere around optical tidal disruption events

We report a star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events have now been identified in optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three tidal disruption events. The continuum polarization appears independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code, we simulate the polarization degree as a function of disk mass and viewing angle and we show that the predicted levels are compatible with the observations for extended reprocessing envelopes of ~1,000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for tidal disruption event models, and opens an important new line of exploration in the study of tidal disruption events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2019

Comparison of Space Radiation GCR Models to AMS Heavy Ion Data

Galactic cosmic rays (GCR) are a constant source of radiation that constitutes one of the major hazards during deep space exploration missions for both astronauts and hardware. In this work, GCR models commonly used by the space radiation protection community are compared with recently published high-precision, high- resolution measurements of cosmic ray lithium, beryllium, boron, carbon, nitrogen, and oxygen fluxes along with their ratios (Li/B, Li/C, Li/O, Be/B, Be/C, Be/O, B/C, B/O, C/O, N/B, N/O) from the Alpha Magnetic Spectrometer (AMS). All of the models were developed and calibrated prior to the publication of this AMS data, therefore this is an opportunity to validate the models against an independent data set. This paper is a compliment to the previously published comparison of GCR models with AMS hydrogen, helium, and the boron-to-carbon ratio.

Whitman, Kathryn↗
NASA NTRS1975

Model of Birkeland currents associated with an auroral arc

The ionospheric configuration of the vertical Birkeland currents associated with an auroral arc is calculated from a model based on the assumption that the auroral electrons constitute a net current into the ionosphere. It is further assumed that the return vertical current out of the ionosphere is impeded along the field lines which connect the arc to the magnetosphere. The magnitude, location, and structure of the return Birkeland current calculated from the model agree well with measurements of these parameters by a rocket-borne experiment.

Vondrak, R. R.↗
DOE OSTI2024

Low-energy radiative backgrounds in CCD-based dark-matter detectors

The reach of sub-GeV dark-matter detectors is at present severely affected by low-energy events from various origins. We present the theoretical methods to compute the single- and few-electron events that arise from secondary radiation emitted by high-energy particles as they pass through detector materials and perform a detailed simulation to quantify them at (Skipper) CCD-based experiments, focusing on the SENSEI data collected at Fermilab near the MINOS cavern. The simulations account for the generation of secondaries from Cherenkov and luminescent recombination radiation; photo-absorption in the bulk, backside layer, pitch adapter, and epoxy; the photon reflection and refraction at interfaces; thin-film interference; the roughness of the interfaces; the dynamics of charges produced in the highly doped CCD-backside-layers; and the partial charge collection on the CCD backside. We consider several systematic uncertainties, notably those stemming from the backside modeling, which we estimate with a “fiducial” and an “extreme” charge-diffusion model, with the former model being preferred due to better agreement with partial-charge collection data. We find that Cherenkov photons constitute about 30% of the observed single-electron events for both diffusion models; radiative recombination contributes negligibly to the event rate for the fiducial model, although it can dominate over Cherenkov for the extreme model. We also estimate the fraction of 2-electron events that arise from 1-electron event coincidences in the same pixel, finding that the entire 2-electron rate can be explained by coincidences of radiative events and spurious charge. Accounting for both radiative and non-radiative backgrounds, we project the sensitivity of future Skipper-CCD-based experiments to different dark-matter models. For light-mediator models with dark-matter masses of 1, 5, and 10 MeV, we find that future experiments with 10-kg-year exposures and successful background mitigation could have a sensitivity that is larger by 9, 3, and 2 orders of magnitude, respectively, when compared to an experiment without background improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2012

Computational Implementation of a Thermodynamically Based Work Potential Model For Progressive Microdamage and Transverse Cracking in Fiber-Reinforced Laminates

A continuum-level, dual internal state variable, thermodynamically based, work potential model, Schapery Theory, is used capture the effects of two matrix damage mechanisms in a fiber-reinforced laminated composite: microdamage and transverse cracking. Matrix microdamage accrues primarily in the form of shear microcracks between the fibers of the composite. Whereas, larger transverse matrix cracks typically span the thickness of a lamina and run parallel to the fibers. Schapery Theory uses the energy potential required to advance structural changes, associated with the damage mechanisms, to govern damage growth through a set of internal state variables. These state variables are used to quantify the stiffness degradation resulting from damage growth. The transverse and shear stiffness of the lamina are related to the internal state variables through a set of measurable damage functions. Additionally, the damage variables for a given strain state can be calculated from a set of evolution equations. These evolution equations and damage functions are implemented into the finite element method and used to govern the constitutive response of the material points in the model. Additionally, an axial failure criterion is included in the model. The response of a center-notched, buffer strip-stiffened panel subjected to uniaxial tension is investigated and results are compared to experiment.

Pineda, Evan J.↗
DOE OSTI2024

Implications of a Pervasive Climate Model Bias for Low‐Cloud Feedback

Abstract How low clouds respond to warming constitutes a key uncertainty for climate projections. Here we observationally constrain low‐cloud feedback through a controlling factor analysis based on ridge regression. We find a moderately positive global low‐cloud feedback (0.45 W , 90% range 0.18–0.72 W ), about twice the mean value (0.22 W ) of 16 models from the Coupled Model Intercomparison Project. We link this discrepancy to a pervasive model mean‐state bias: models underestimate the low‐cloud response to warming because (a) they systematically underestimate present‐day tropical marine low‐cloud amount, and (b) the low‐cloud sensitivity to warming is proportional to this present‐day low‐cloud amount. Our results hence highlight the importance of reducing model biases in both the mean state of clouds and their sensitivity to environmental factors for accurate climate change projections.

58 GEOSCIENCES↗
NASA NTRS2002

Finite-Strain Fractional-Order Viscoelastic (FOV) Material Models and Numerical Methods for Solving Them

Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.

Freed, Alan D.↗
DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on ultralight axion dark matter through galaxy cluster number counts

Ultralight axions are hypothetical scalar particles that influence the evolution of large-scale structures of the Universe. Depending on their mass, they can potentially be part of the dark matter component of the Universe as candidates commonly referred to as fuzzy dark matter. While strong constraints have been established for pure fuzzy dark matter models, the more general scenario where ultralight axions constitute only a fraction of the dark matter has been limited to only a few observational probes. In this work, we use the galaxy cluster number counts obtained from the first All-Sky Survey (eRASS1) of the SRG/eROSITA mission together with gravitational weak lensing data from the Dark Energy Survey, the Kilo-Degree Survey, and the Hyper Suprime-Cam to constrain the fraction of ultralight axions in the mass range 10 −32 eV to 10 −24 eV. We put upper bounds on the ultralight axion relic density Ω a in independent logarithmic axion mass bins by performing a full cosmological parameter inference. We find an exclusion region in the intermediate ultralight axion mass regime with the tightest bounds reported so far in the mass bins around m a = 10 −27 eV with Ω a < 0.0035 and m a = 10 −26 eV with Ω a < 0.0079; both are at a 95% confidence level. When combined with cosmic microwave background probes, these bounds are tightened to Ω a < 0.0030 in the m a = 10 −27 eV mass bin and Ω a < 0.0058 in the m a = 10 −26 eV mass bin, with both at a 95% confidence level. This is the first time that constraints on ultralight axions have been obtained using the growth of structure measured by galaxy cluster number counts. These results pave the way for large surveys, which can be utilized to obtain tight constraints on the mass and relic density of ultralight axions with better theoretical modeling of the abundance of halos.

cosmological parameters↗
DOE OSTI2023

Oncogenic Ras and ΔNp63α cooperate to recruit immunosuppressive polymorphonuclear myeloid-derived suppressor cells in a mouse model of squamous cancer pathogenesis

Introduction Amplification of human chromosome 3q26-29, which encodes oncoprotein ΔNp63 among other isoforms of the p63 family, is a feature common to squamous cell carcinomas (SCCs) of multiple tissue origins. Along with overexpression of ΔNp63, activation of the protooncogene, RAS , whether by overexpression or oncogenic mutation, is frequently observed in many cancers. In this study, analysis of transcriptome data from The Cancer Genome Atlas (TCGA) demonstrated that expression of TP63 mRNA , particularly ΔNp63 isoforms, and HRAS are significantly elevated in advanced squamous cell carcinomas of the head and neck (HNSCCs), suggesting pathological significance. However, how co-overexpressed ΔNp63 and HRAS affect the immunosuppressive tumor microenvironment (TME) is incompletely understood. Methods Here, we established and characterized an immune competent mouse model using primary keratinocytes with retroviral-mediated overexpression of ΔNp63α and constitutively activated HRAS (v-ras Ha G12R) to evaluate the role of these oncogenes in the immune TME. Results In this model, orthotopic grafting of wildtype syngeneic keratinocytes expressing both v-ras Ha and elevated levels of ΔNp63α consistently yield carcinomas in syngeneic hosts, while cells expressing v-ras Ha alone yield predominantly papillomas. We found that polymorphonuclear (PMN) myeloid cells, experimentally validated to be immunosuppressive and thus representing myeloid-derived suppressor cells (PMN-MDSCs), were significantly recruited into the TME of carcinomas arising early following orthotopic grafting of ΔNp63α/v-ras Ha -expressing keratinocytes. ΔNp63α/v-ras Ha -driven carcinomas expressed higher levels of chemokines implicated in recruitment of MDSCs compared to v-ras Ha -initiated tumors, providing a heretofore undescribed link between ΔNp63α/HRAS-driven carcinomas and the development of an immunosuppressive TME. Conclusion These results support the utilization of a genetic carcinogenesis model harboring specific genomic drivers of malignancy to study mechanisms underlying the development of local immunosuppression.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A mesoscale 3D model of irradiated concrete informed via a 2.5 U-Net semantic segmentation

The concrete biological shield in light-water reactors is exposed to neutron and gamma irradiation, which deteriorates the concrete’s mechanical properties in the long term. To assess the irradiation-induced damage, predictive mechanical models are developed and used in parallel with the characterization of irradiated concrete samples. Realistic 3D simulation domains can drastically improve a model’s prediction. In this work, we utilized x-ray computed tomography (XCT) data of a concrete specimen to reconstruct its 3D microstructure. The XCT data shows low contrast between the concrete’s aggregates and cement paste, resulting in poor image segmentation when using traditional unsupervised techniques. To address this issue, we developed and trained a 2.5D U-Net model on only 24 pre-labeled XCT layers to segment 651 layers of the XCT data. The overall F1-score of the model is approximately 96%. Then, we created a 3D finite element (FE) mesh based on the stack of segmented images. The FE model contains radiation-induced expansion, damage, and creep. The constitutive equations are adapted to each phase (aggregates and cement paste). Here, we simulated the effects of neutron irradiation in the concrete specimen as well as the specimen’s mechanical response to uniaxial compression. Finally, model validation was performed using experimental data on similar concrete specimens in the literature.

2.5D U-Net↗