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 829 records · Page 46

Modeling and Analysis of DC Microgrids as Stochastic Hybrid Systems

This study proposes a method of predicting the influence of random load behavior on the dynamics of dc microgrids and distribution systems. This is accomplished by combining stochastic load models and deterministic microgrid models. Together, these elements constitute a stochastic hybrid system. The resulting model enables straightforward calculation of dynamic state moments, which are used to assess the probability of desirable operating conditions. Specific consideration is given to systems based on the dual active bridge (DAB) topology. Bounds are derived for the probability of zero voltage switching (ZVS) in DAB converters. A simple example is presented to demonstrate how these bounds may be used to improve ZVS performance as an optimization problem. In conclusion, predictions of state moment dynamics and ZVS probability assessments are verified through comparisons to Monte Carlo simulations.

42 ENGINEERING↗

The Microstructure Aware Plasticity Model: Formulation and Usage Guide

Herein, the formulation, parameter sensitivities, and usage methods for the Microstructure-Aware Plasticity (MAP) model are presented. This document is intend to serve as a reference for the underlying theory that constitutes the MAP model and as a practical guide for analysts and future developers on how aspects of this material model influence generalized mechanical behavior.

42 ENGINEERING↗

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗

Tailoring ductile-phase toughened tungsten hierarchical microstructures for plasma-facing materials

Biological materials, such as bones, nacre, etc., have been found to show attractive and unique combinations of stiffness, strength and fracture toughness. Research has been conducted in various engineering areas to mimic the naturally hierarchical microstructures or nanostructures of these materials to produce materials with optimum stiffness, high strength and fracture toughness for their intended structural applications. In the same objective, this work applies a multiscale microstructural approach recently developed to investigate the deformation and fracture behavior of ductile phase toughened tungsten (W) materials such as tungsten-nickel/iron (W-Ni/Fe) composites that possess lamellar-like and hierarchical “brick-and-mortar” (BAM) microstructures. First, the approach is used to simulate tensile loading of W-Ni/Fe specimens cut out from hot-rolled W-Ni/Fe plates which possess lamellar-like microstructures. The finite element model of the gage section of the specimen consists of a W-Ni/Fe dual-phase microstructural domain in which the constitutive behaviors of W and Ni/Fe phases are described by an elastic-plastic damage model. The predicted material stress-strain response and crack pattern development as a function of loading are then compared to the corresponding experimental results to determine the constitutive model parameters for W and Ni/Fe. Subsequently, the model parameters are used to analyze W-Ni/Fe BAM microstructures that are artificially created to investigate the effects of microstructural features and morphologies on the composite stress-strain response, damage, and fracture patterns. The modeling shows that regular BAM microstructures that experience bridging mechanism combined with crack penetration across W-phase regions exhibit significantly higher strengths than the rather random lamellar-like microstructures. More significantly, however adjusting the brick’s length-to-height ratio, the BAM microstructure can be designed to allow a more distributed damage zone which leads to increased strength, ductility and fracture energy.

36 MATERIALS SCIENCE↗

1D Modeling for the Development of a Digital Twin for Hydropower Systems with Applications for a Testing Rig

Modeling and simulation constitute two important parts in constructing a sensible digital twin to mimic the dynamics of hydropower systems. Because of the physical nature of a hydropower system, the basic modeling should cover the water flow systems from the reservoir to the penstock (inlet water pipes), penstock to hydro turbine, and hydro turbine to generator and from the linkage of the hydropower systems to the grid. This report describes initial attempts to model these dynamic components, including the formulation of the linearized state space model for hydro turbine systems using the well-known six-coefficients linearization method and the formulation of the voltage and power dynamical models of a synchronous generator. To validate the accuracy of the proposed model, data collected from Norwegian University of Science and Technology were used to obtain relevant parameters for the model, and the desired simulation results were obtained. Additionally, neural network modeling and learning were also developed and applied to model the generation torque and water flow rate subsystems to demonstrate a potential learning scheme, which can be used in the development of a digital twin.

13 HYDRO ENERGY↗

Dynamic and Thermal Turbulent Time Scale Modelling for Homogeneous Shear Flows

A new turbulence model, based upon dynamic and thermal turbulent time scale transport equations, is developed and applied to homogeneous shear flows with constant velocity and temperature gradients. The new model comprises transport equations for k, the turbulent kinetic energy; tau, the dynamic time scale; k(sub theta), the fluctuating temperature variance; and tau(sub theta), the thermal time scale. It offers conceptually parallel modeling of the dynamic and thermal turbulence at the two equation level, and eliminates the customary prescription of an empirical turbulent Prandtl number, Pr(sub t), thus permitting a more generalized prediction capability for turbulent heat transfer in complex flows and geometries. The new model also incorporates constitutive relations, based upon invariant theory, that allow the effects of nonequilibrium to modify the primary coefficients for the turbulent shear stress and heat flux. Predictions of the new model, along with those from two other similar models, are compared with experimental data for decaying homogeneous dynamic and thermal turbulence, homogeneous turbulence with constant temperature gradient, and homogeneous turbulence with constant temperature gradient and constant velocity gradient. The new model offers improvement in agreement with the data for most cases considered in this work, although it was no better than the other models for several cases where all the models performed poorly.

Schwab, John R.↗

Measuring and modeling the backscattering cross section of a leaf

Leaves are a significant feature of any vegetation canopy, and for remote sensing purposes it is important to develop an effective model for predicting the scattering from a leaf. From measurements of the X band backscattering cross section of a coleus leaf in varying stages of dryness, it is shown that a uniform resistive sheet constitutes such a model for a planar leaf. The scattering is determined by the (complex) resistivity which is, in turn, entirely specified by the gravimetric moisture content of the leaf. Using an available asymptotic expression for the scattering from a rectangular resistive plate which includes, as a special case, a metallic plate whose resistivity is zero, the computed backscattering cross sections for both principal polarizations are found to be in excellent agreement with data measured for rectangular sections of leaves with different moisture contents. If the resistivity is sufficiently large, the asymptotic expressions do not differ significantly from the physical optics ones, and for naturally shaped leaves as well as rectangular sections, the physical optics approximation in conjunction with the resistive sheet model faithfully reproduces the dominant feataures of the scattering patterns under all moisture conditions.

Senior, T. B. A.↗

Two-Dimensional Quantum Model of a Nanotransistor

A mathematical model, and software to implement the model, have been devised to enable numerical simulation of the transport of electric charge in, and the resulting electrical performance characteristics of, a nanotransistor [in particular, a metal oxide/semiconductor field-effect transistor (MOSFET) having a channel length of the order of tens of nanometers] in which the overall device geometry, including the doping profiles and the injection of charge from the source, gate, and drain contacts, are approximated as being two-dimensional. The model and software constitute a computational framework for quantitatively exploring such device-physics issues as those of source-drain and gate leakage currents, drain-induced barrier lowering, and threshold voltage shift due to quantization. The model and software can also be used as means of studying the accuracy of quantum corrections to other semiclassical models.

Govindan, T. R.↗

New Developments in the CompDam Fiber Kinking Model for the Interaction

This paper describes two new additions to the CompDam fiber kinking model: 1) a failure index for fiber kinking to facilitate interaction between kinking and matrix cracking and 2) a strain-based criterion for fiber rupture to account for fiber breakage before or after fiber kinking. The model formulation is described with the aid of the graphical Considère construction, which is used to determine the conditions for fiber kinking. Verification studies are conducted with single element models to show the behavior of the model in light of the new developments and to evaluate the effect of numerical parameters of the model. The model is used to analyze open hole compression (OHC) specimens with soft, quasi-isotropic, and hard layups. Analysis results are shown for a conventional continuum damage mechanics model using a trilinear constitutive law for fiber damage and for the proposed model to highlight the differences between the two. Comparison with test data highlights the capabilities and remaining gaps in the quantitative prediction of the damage and failure response. While both models can predict many features of the OHC response, accurate prediction of kinking onset in the OHC specimens remains elusive. The results suggest that it is necessary to predict the combined effects of compression and shear loading as well as the structural effects (e.g., ply thickness, strain gradient effect) accurately to match the experimental observations of fiber kinking onset in the three OHC laminates.

Andrew C. Bergan↗

Comparison of continuum and particle simulations of expanding rarefied flows

Comparisons of Navier-Stokes solutions and particle simulations for a simple two-dimensional model problem at a succession of altitudes are performed in order to assess the importance of rarefaction effects on the base flow region. In addition, an attempt is made to include 'Burnett-type' extensions to the Navier-Stokes constitutive relations. The model geometry consists of a simple blunted wedge with a 0.425 meter nose radius, a 70 deg cone half angle, a 1.7 meter base length, and a rounded shoulder. The working gas is monatomic with a molecular weight and viscosity similar to air and was chosen to focus the study on the continuum and particle methodologies rather than the implementation of thermo-chemical modeling. Three cases are investigated, all at Mach 29, with densities corresponding to altitudes of 92 km, 99 km, and 105 km. At the lowest altitude, Navier-Stokes solutions agree well with particle simulations. At the higher altitudes, the Navier-Stokes equations become less accurate. In particular, the Navier-Stokes equations and particle method predict substantially different flow turning angle in the wake near the after body. Attempts to achieve steady continuum solutions including 'Burnett-type' terms failed. Further research is required to determine whether the boundary conditions, the equations themselves, or other unknown causes led to this failure.

Lumpkin, Forrest E., III↗

Rate Dependent Deformation and Strength Analysis of Polymer Matrix Composites

A research program is being undertaken to develop rate dependent deformation and failure models for the analysis of polymer matrix composite materials. In previous work in this program, strain-rate dependent inelastic constitutive equations used to analyze polymers have been implemented into a mechanics of materials based composite micromechanics method. In the current work, modifications to the micromechanics model have been implemented to improve the calculation of the effective inelastic strain. Additionally, modifications to the polymer constitutive model are discussed in which pressure dependence is incorporated into the equations in order to improve the calculation of constituent and composite shear stresses. The Hashin failure criterion is implemented into the analysis method to allow for the calculation of ply level failure stresses. The deformation response and failure stresses for two representative uniaxial polymer matrix composites, IM7/977-2 and AS4-PEEK, are predicted for varying strain rates and fiber orientations. The predicted results compare favorably to experimentally obtained values.

Goldberg, Robert K.↗

The star-formation history in the last 10 billion years from CIB cross-correlations

ABSTRACT The cosmic infrared background (CIB) traces the emission of star-forming galaxies throughout all cosmic epochs. Breaking down the contribution from galaxies at different redshifts to the observed CIB maps would allow us to probe the history of star formation. In this paper, we cross-correlate maps of the CIB with galaxy samples covering the range z ≲ 2 to measure the bias-weighted star-formation rate (SFR) density 〈bρSFR〉 as a function of time in a model independent way. This quantity is complementary to direct measurements of the SFR density ρSFR, giving a higher weight to more massive haloes, and thus provides additional information to constrain the physical properties of star formation. Using cross-correlations of the CIB with galaxies from the DESI Legacy Survey and the extended Baryon Oscillation Spectroscopic Survey, we obtain high signal-to-noise ratio measurements of 〈bρSFR〉, which we then use to place constraints on halo-based models of the star-formation history. We fit halo-based SFR models to our data and compare the recovered ρSFR with direct measurements of this quantity. We find a qualitatively good agreement between both independent data sets, although the details depend on the specific halo model assumed. This constitutes a useful robustness test for the physical interpretation of the CIB, and reinforces the role of CIB maps as valuable astrophysical probes of the large-scale structure. We report our measurements of 〈bρSFR〉 as well as a thorough account of their statistical uncertainties, which can be used to constrain star-formation models in combination with other data.

Astronomy & Astrophysics↗

graphenv: a Python library for reinforcement learning on graph search spaces

Many important and challenging problems in combinatorial optimization (CO) can be expressed as graph search problems, in which graph vertices represent full or partial solutions and edges represent decisions that connect them. Graph structure not only introduces strong relational inductive biases for learning (Battaglia et al., 2018) - in this context, by providing a way to explicitly model the value of transitioning (along edges) between one search state (vertex) and the next - but lends itself to problems both with and without clearly defined algebraic structure. For example, classic CO problems on graphs such as the Traveling Salesman Problem (TSP) can be expressed as either pure graph search or integer programs. Other problems, however, such as molecular optimization, do no have concise algebraic formulations and yet are readily implemented as a graph search (V. et al., 2022; Zhou et al., 2019). Such "model-free" problems constitute a large fraction of modern reinforcement learning (RL) research owing to the fact that it is often much easier to write a forward simulation that expresses all of the state transitions and rewards, than to write down the precise mathematical expression of the full optimization problem. In the case of molecular optimization, for example, one can use domain knowledge alongside existing software libraries to model the effect of adding a single bond or atom to an existing but incomplete molecule, and let the RL algorithm build a model of how good a given decision is by "experiencing" the simulated environment many times through. In contrast, a model-based mathematical formulation that fully expresses all the chemical and physical constraints is intractable. In recent years, RL has emerged as an effective paradigm for optimizing searches over graphs and led to state-of-the-art heuristics for games like Go and chess, as well as for classical CO problems such as the TSP. This combination of graph search and RL, while powerful, requires non-trivial software to execute, especially when combining advanced state representations such as Graph Neural Networks (GNN) with scalable RL algorithms.

97 MATHEMATICS AND COMPUTING↗

Modeling Modern Methane Emissions from Natural Wetlands: Model Description and Results - 1

Methane is an important greenhouse gas which contributes about 22 percent to the present greenhouse effect. Natural wetlands currently constitute the biggest methane source and were the major source in preindustrial times. Wetland emissions depend highly on the climate, i.e., on soil temperature and water table. To investigate the response of methane emissions from natural wetlands to climate variations, a process-based model that derives methane emissions from natural wetlands as a function of soil temperature, water table, and net primary productivity is used. For its application on the global scale, global data sets for all model parameters are generated. In addition, a simple hydrologic model is developed in order to simulate the position of the water table in wetlands. The hydrologic model is tested against data from different wetland sites, and the sensitivity of the hydrologic model to changes in precipitation is examined. The global methane­ hydrology model constitutes a tool to study temporal and spatial variations in methane emissions from natural wetlands. The model is applied using high-frequency atmospheric forcing fields from European Center for Medium-range Weather Forecasts (ECMWF) re-analyses of the period from 1982 to 1993. We calculate global annual methane emissions from wetlands to be 260 teragrams per year. Twenty-five percent of these methane emissions originate from wetlands north of 30 degrees North Latitude. Only 60 percent of the produced methane is emitted, while the rest is re-oxidized. A comparison of zonal integrals of simulated global wetland emissions and results obtained by an inverse modeling approach shows good agreement. In a test with data from two wetlands the seasonality of simulated and observed methane emissions agrees well.

Biogeochemical processes↗

Physics augmented machine learning discovery of composition-dependent constitutive laws for 3D printed digital materials

Multi-material 3D printing, particularly through polymer jetting, enables the fabrication of digital materials by mixing distinct photopolymers at the micron scale within a single build to create a composite with tunable mechanical properties. Here, this work presents an integrated experimental and computational investigation into the composition-dependent mechanical behavior of 3D printed digital materials. We experimentally characterize five formulations, combining soft and rigid UV-cured polymers under uniaxial tension and torsion across three strain and twist rates. The results reveal nonlinear and rate-dependent responses that strongly depend on composition. To model this behavior, we develop a physics-augmented neural network (PANN) that combines a partially input convex neural network (pICNN) for learning the composition-dependent hyperelastic strain energy function with a quasi-linear viscoelastic (QLV) formulation for time-dependent response. The pICNN ensures convexity with respect to strain invariants while allowing non-convex dependence on composition. To enhance interpretability, we apply $L_0$ sparsification. For the time-dependent response, we introduce a multilayer perceptron (MLP) to predict viscoelastic relaxation parameters from composition. The proposed model accurately captures the nonlinear, rate-dependent behavior of 3D printed digital materials in both uniaxial tension and torsion, achieving high predictive accuracy for interpolated material compositions. This approach provides a scalable framework for automated, composition-aware constitutive model discovery for multi-material 3D printing.

Constitutive modeling↗