Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum mechanics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

On the Fidelity of Computational Models for the Flow of Milled Loblolly Pine: A Benchmark Study on Continuum-Mechanics Models and Discrete-Particle Models

The upstream of bioenergy industry has suffered from unreliable operations of granular biomass feedstocks in handling equipment. Computational modeling, including continuum-mechanics models and discrete-particle models, offers insightful understandings and predictive capabilities on the flow of milled biomass and can assist equipment design and optimization. This paper presents a benchmark study on the fidelity of the continuum and discrete modeling approaches for predicting granular biomass flow. We first introduce the constitutive law of the continuum-mechanics model and the contact law of the coarse-grained discrete-particle model, with model parameters calibrated against laboratory characterization tests of the milled loblolly pine. Three classical granular material flow systems (i.e., a lab-scale rotating drum, a pilot-scale hopper, and a full-scale inclined plane) are then simulated using the two models with the same initial and boundary conditions as the physical experiments. The close agreement of the numerical predictions with the experimental measurements on the hopper mass flow rate, the hopper critical outlet width, the material stopping thickness on the inclined plane, and the dynamic angle of repose, clearly indicates that the two methods can capture the critical flow behavior of granular biomass. The qualitative comparison shows that the continuum-mechanics model outperforms in parameterization of materials and wall friction, and large-scale systems, while the discrete-particle model is more preferred for discontinuous flow systems at smaller scales. Industry stakeholders can use these findings as guidance for choosing appropriate numerical tools to model biomass material flow in part of the optimization of material handling equipment in biorefineries.

Jin, Wencheng↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2025

Green's function fast multipole method for continuum mechanics (SM-FMM)

Solid Mechanics Fast Multipole Method based on elastic Green's function and accelerated using FFTs The code calculates the mechanical fields (stress and strain) for a heterogeneous elasto-plastic material unit cell under quasi-static conditions (no dynamic effects). The fields are calculated using a Green's function method, where the strain field is given by a discrete convolution of Green's operator with an auxiliary stress field. The convolution is calculated using the fast multipole method.

Zecevic, Miroslav↗
DOE OSTI2022

Time-dependent density functional theory of charge, energy, and spin transport, and quantum continuum mechanics (Final Report)

The overarching goal of this project, which unfolded over a period of seventeen years, including two no-cost extensions in the last two years, was to extend and deepen our ability to use time-dependent density functional theory and related theoretical methods for the description of the transport of charge, spin, or energy in quantum electronic systems at the nanoscale. The output of our efforts is contained in 52 papers, listed below, which we have published since 2006. There have also been many presentations at international conferences, workshops, colloquia and seminars. This body of work has been carried out with the essential input of 8 postdoctoral students (who were supported by the grant) and several external senior collaborators and their postdocs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fabric-enriched continuum breakage mechanics (F-CBM )

A fabric-enriched continuum breakage mechanics (F-CBM) framework is developed to examine the relation between elastic and inelastic anisotropy in granular materials subjected to high pressure. For this purpose, emphasis is given to the influence of the anisotropic fabric on the plastic cap regime – that is, a macroscopic signature of the inelasticity of granular materials widely dependent on the particle crushability. Anisotropic behaviour prior to yielding is introduced through a symmetric second-order fabric tensor embedded in the expression of the elastic energy potential. It is shown that, thanks to its energy-based formulation, the F-CBM framework provides a platform to predict the rotation and distortion of the yield cap of granular materials as an outcome of anisotropic strain energy storage prior to grain crushing. Parametric analyses indicate that both the degree of anisotropy and the orientation of the fabric have an impact on the energy release resulting from grain crushing and contribute to the resulting stress–strain behaviour. It is shown that, despite the use of very few parameters, the proposed framework can accurately predict the anisotropic yielding and inelastic stress–strain response solely on the basis of the directional properties inferred from measurements of elastic anisotropy.

Engineering↗
DOE OSTI2022

Application of the Continuum Damage Mechanics Wilshire-Cano-Stewart (WCS) Model

In this study, the applications of the continuum damage mechanics-based Wilshire-Cano-Stewart (WCS) model are explored to predict rupture time, minimum-creep-strain-rate (MCSR), damage, damage evolution, and creep deformation. Increase knowledge in manufacturing methods has pushed the limit of material science and the development of new materials. Conventional testing is required to qualify materials against creep which according to the ASME B&PV III code, 10,000+ hours of experiments are necessary for each heat before materials are put into service. This process is costly and not feasible for new materials. As an alternative, models have been employed to predict creep behaviors and reduce the amount of time necessary for material qualification. Many models have been developed to predict distinct creep behaviors and the question of which model is best remains. Amongst current models, the WCS model has emerge with the ability to predict multiple behaviors using an explicit analytical approach with the ability to predict long-term creep. In this study the novel continuum damage mechanics WCS is employed in multiple applications. The goals of the study are (a) to discuss and determine the framework of the WCS model and validated it mathematically and using parametric simulations, (b) applied the model to accelerated creep data to show the capabilities of the model with non-conventional data, and (c) applied a novel numerical method, the datum temperature method (DTM) to show the model extrapolations and interpolations capabilities with limited and reduce data sets. To accomplish these goals, data is gathered for alloy P91 and Inconel 718 to develop and post-audit validate the model. The benefits of using the WCS model is that it provides an explicit stress and temperature dependency ideal for extrapolations, the coupled equations are suitable for finite element analysis (FEA) implementation, and it follows an explicit calibration approach. The model also proves that it can be applied to accelerated testing data and using the DTM.

Cano, Jaime A↗
DOE OSTI2024

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Continuum Damage Mechanics Modeling of High-Temperature Flaw Propagation: Application to Creep Crack Growth In 316H Standardized Specimens and Nuclear Reactor Components

Predicting creep crack growth (CCG) of flaws found during operation in high-temperature alloy components is essential for assessing the remaining lifetime of those components. While defect assessment procedures are available for this purpose in design codes, these are limited in their range of applicability. This study assesses the application of a local damage-based finite-element methodology as a more general technique for the prediction of CCG at high temperatures on a variety of structural configurations. Numerical results for stainless steel 316H, which are validated against experimental data, show the promise of this approach. This integration of continuum damage mechanics (CDM) based methodologies, together with adequate inelastic models, into assessment procedures can therefore inform the characterization of CCG under complex operating conditions, while avoiding excessive conservatism. This article shows that such modeling frameworks can be calibrated to experimental data and used to demonstrate that the degree of tri-axiality ahead of a growing creep crack affects its rate of growth. The framework is also successfully employed in characterizing CCG in a realistic reactor pressure vessel geometry under an arbitrary loading condition. These results are particularly relevant to the nuclear power industry for defect assessment and inspections as part of codified practices of structural components with flaws in high-temperature reactors.

42 ENGINEERING↗
DOE OSTI2023

Multiphysics Modeling Framework for Soft PVC Gel Sensors with Experimental Comparisons

Polyvinyl chloride (PVC) gels have recently been found to exhibit mechanoelectrical transduction or sensing capabilities under compressive loading applications. This phenomenon is not wholly understood but has been characterized as an adsorption-like phenomena under varying amounts and types of plasticizers. A different polymer lattice structure has also been tested, thermoplastic polyurethane, which showed similar sensing characteristics. This study examines mechanical and electrical properties of these gel sensors and proposes a mathematical framework of the underlying mechanisms of mechanoelectrical transduction. COMSOL Multiphysics is used to show solid mechanics characteristics, electrostatic properties, and transport of interstitial plasticizer under compressive loading applications. The solid mechanics takes a continuum mechanics approach and includes a highly compressive Storakers material model for compressive loading applications. The electrostatics and transport properties include charge conservation and a Langmuir adsorption migration model with variable diffusion properties based on plasticizer properties. Results show both plasticizer concentration gradient as well as expected voltage response under varying amounts and types of plasticizers. Experimental work is also completed to show agreeance with the modeling results.

42 ENGINEERING↗
DOE OSTI2021

A numerical study of reversible plasticity using continuum dislocation mechanics

In this contribution, an elasto-viscoplastic fast Fourier transform-based (EVPFFT) numerical implementation of the Mesoscale Field Dislocation Mechanics (MFDM) formulation, called MFDM-EVPFFT, is applied to study the reversible plastic behavior of periodic two-phase crystalline composites with an elasto-viscoplastic plastic matrix and a purely elastic second phase. Additionally, periodic laminate microstructures of this kind with different periods (i.e. sizes) are considered to examine the size dependence of the Bauschinger effect and hardening during cyclic loading. Comparisons with classic composite effects obtained with conventional crystal plasticity are discussed. Specifically, the MFDM-EVPFFT results shed light on the hardening mechanisms due to piling-up/unpiling-up of geometrically-necessary dislocations (GND) during reverse loading.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction of nonlinear specific heat during single crystal HMX phase transition

Here, we develop a thermodynamically consistent chemo-thermo-mechanical model for the β → δ phase transition of energetic HMX crystals. In contrast to previous models, which either considered specific heat to be a constant or utilized a calibrated function, this model provides novel expressions for the specific heats at constant volume and constant elastic strains derived directly from continuum mechanics. In addition, the model provides a novel prediction for the critical temperature at which the chemical heating rate achieves its extremum for Arrhenius kinetics. The numerical solution predicts highly nonlinear specific heat behavior including order of magnitude changes.

42 ENGINEERING↗
DOE OSTI2024

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗
DOE OSTI2022

Statistical Estimation of Strain Using Spatial Correlation Functions

Ex-situ estimation of strains from deformed micrographs is not possible as there are no persistent features which can be tracked. Two point spatial statistics enable the rigorous quantification of spatial patterns in heterogeneous media. In this paper, we propose a novel method for estimating strains directly from dissimilar micrographs using a continuum mechanics approach. Rather than operating directly on images from sequential frames, as is done in digital image correlation, we operate on different microstructure realizations. This is made possible by comparing the spatial autocorrelation maps of deformed and undeformed micrographs rather than direct comparison of images. Additionally, a Bayesian framework is proposed for quantifying uncertainty. We first illustrate the efficacy of this method on speckle pattern images from digital image correlation experiments. Then, we demonstrate that the method is capable of operating on dissimilar micrographs using deformed synthetic binary microstructures. Finally, we present a case study on polycrystalline additively manufactured 316L deformed via tension. The proposed method works well and we discuss implications and limitations of the presented work.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The role of deviatoric stress and dislocations on the α to ω phase transformation in Ti

Under extreme conditions, α-Ti becomes unstable and transforms either into β-Ti at high temperature or into ω-Ti at high pressure. In what concerns the α to ω phase transformation (PT), there has been a wide range of experimentally reported transition pressures from approximately 2 to 15 GPa at room temperature. Deviatoric stresses and internal defects are often assumed to be the root cause of this variation. Here, in this study, these postulates are revisited using both continuum mechanics and molecular dynamics (MD) simulations. First, a simple continuum model, assuming linear elasticity and isotropic plasticity, is developed to describe the effects of applied stress and dislocations on the stability of an ω nucleus in an infinite α domain. Second, a new MD simulation method is developed to generate an ω nucleus in the α domain utilizing the displacement field identified from the topological analysis. Results from MD simulations show that despite the fact that phase diagrams typically delineate the limits between two phases in terms of only P and T, deviatoric stress promotes the α to ω phase transformation by reducing the critical radius above which an ω nucleus is stable. Furthermore, the required deviatoric stress to nucleate and stabilize a nanoscale ω nucleus is likely emanating from the internal stress of defects such as dislocations. The MD-informed micromechanics models are used to identify favorable configurations where dislocations help favor the α to ω transformation. These configurations show that the interaction with a basal or prismatic dislocation reduces the critical radius of a ω nucleus by about 10 or 16 %, respectively. In addition, prismatic edge dislocations are found to promote the growth of ω nucleus when interacting with the ($\bar{1}$$\bar{1}$20) α //(0001) ω interface. Importantly, a simple model of the arrival of dislocation at an ω nucleus suggests that PT does not necessarily require a pile-up to be present but could alternatively be mediated by a constant rapid flow of dislocations.

36 MATERIALS SCIENCE↗