Search NASA⌕ Search

SEARCH · Search NASA

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 OSTI2021

Overall equilibrium in the coupling of peridynamics and classical continuum mechanics

Coupling peridynamics based computational tools with those using classical continuum mechanics can be very beneficial, because it can provide a means to generate a computational method that combines the efficiency of classical continuum mechanics with the capability to simulate crack propagation, typical of peridynamics. This paper presents an overlooked issue in this type of coupled computational methods: the lack of overall equilibrium. This can be the case even if the coupling strategy satisfies the usual numerical tests involving rigid body motions as well as uniform and linear strain distributions. We focus our investigation on the lack of overall equilibrium in an approach to couple peridynamics and classical continuum mechanics recently proposed by the authors. In our examples, the magnitude of the out-of-balance forces is a fraction of a percent of the applied forces, but it cannot be assumed to be a numerical round-off error. We show analytically and numerically that the main reason for the existence of out-of-balance forces is a lack of balance between the local and nonlocal tractions at the coupling interface. This usually results from the presence of high-order derivatives of displacements in the coupling zone.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2020

Numerical Considerations in the Modeling of a High Explosive Cylinder Experiment Using an ALE Continuum Mechanics Code

A series of experiments involving the detonation of PBX 9501 encased in a copper cylinder are modeled with the objective of evaluating a proposed set of phenomenological parameters for the Wescott–Stewart–Davis reactive burn model. The numerical analysis is conducted using the Los Alamos continuum mechanics code FLAG. Numerical considerations pertaining to various aspects of modeling the experiments using FLAG are discussed. It is shown that use of the proposed set of phenomenological parameters results in predictions of free-surface velocity that match empirically measured velocities reasonably well.

36 MATERIALS SCIENCE↗
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 OSTI2020

A Review of Computational Models for the Flow of Milled Biomass Part II: Continuum-Mechanics Models

The design of efficient material-handling systems for milled lignocellulosic biomass is challenging due to their complex particle morphologies and frictional interactions. Computational modeling, including the discrete element method (DEM) and continuum-based finite-element/volume methods, may offer scientific insight and predictive capabilities for the flow of milled biomass in hoppers and feeders. Herein, this article (Part II) presents a review of current state-of-the-art continuum models for the flow of milled biomass, whereas DEM models are reviewed in a companion article (Part I). Advances of numerical methods to solve the global governing equations are discussed first, followed by a comprehensive review of constitutive models for granular materials, including Drucker–Prager, hypoplastic, Cambridge-type, inertial-rheology, and nonlocal granular fluidity models. Specifically, we provide in-depth discussion on the suitability of those models for milled lignocellulosic biomass materials in terms of nonlinear elasticity, dependence of flow strength on pressure, density and shear rate, and compaction (dilation) associated with hardening (softening). Furthermore, our study shows that, despite the recent advances in continuum granular flow modeling, the most suitable constitutive models still need further development to account for material parametrization, multiflow regimes, and multiscale behavior before they can be reliably used to optimize the design and operation of biomass handling systems.

09 BIOMASS FUELS↗
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 OSTI2020

A continuum damage mechanics (CDM) based Wilshire model for creep deformation, damage, and rupture prediction

In this report, a new constitutive model is derived that combines the Wilshire equations with continuum damage mechanics (CDM) to enable the long-term prediction of creep deformation, damage, and rupture. The Wilshire equations have been demonstrated to accurately extrapolate the stress-rupture, minimum-creep-strain-rate, and time-to-strain of various alloys across decade of life. Recently, the time-to-creep-strain equation has been rearranged to predict creep deformation curves; however, the resulting equation is difficult to calibrate and does not track damage. The CDM-based Sinh constitutive model accurately predicts creep deformation, damage, and stress-rupture for a variety of alloys; however, it lacks an explicit description of stress and temperature co-dependence which limits interpolation and extrapolation ability. In this study, the Wilshire and the CDM-based Sinh equations are combined to create the “WCS” model. The WCS model combines the best features of each model while eliminating their deficiencies. Experimental data for alloy P91 is gathered and the WCS model is calibrated to the data. The WCS model accurately predicts the stress rupture, minimum-creep-strain-rate, creep deformation, and damage of alloy P91 across the experimental stress and temperature range as well as consistency in parametric simulations.

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↗