Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strain gradient modeling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
NASA NTRS1995

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗
DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Direct numerical simulation of a three-dimensional turbulent boundary layer

The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a decrease in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. Terms in the Reynolds stress transport equations were computed. The balances indicate that the decrease in turbulent kinetic energy is a result of a decrease in turbulence production, along with an increase in turbulent dissipation. Intuitive reasoning and current turbulence models would predict an increase in kinetic energy along with increases in production and dissipation rates as a result of increased mean-flow strain rate. Later in the evolution of the flow, both turbulence production and dissipation increase.

Moin, P.↗
NASA NTRS1987

Direct numerical simulation of buoyantly driven turbulence

Numerical simulations of homogeneous turbulence subject to buoyant forcing were performed. The presence of a mean temperature gradient combined with a gravitational field results in a forcing term in the momentum equations. The development of the turbulence was studied and compared to the decay of similar fields in the absence of gravity. In the buoyantly driven field, the vorticity is preferentially aligned with the intermediate eigenvector of the strain-rate tensor and the local temperature gradient is more likely to be aligned with the most compressive eigenvector. These relationships are qualitatively similar to those observed in previous shear flow results studied by Ashurst (1987). A tensor diffusivity model for passive scalar transport developed from shear flow results in Rogers, Moin, and Reynolds (1986) also predicts this buoyant scalar transport, indicating that the relationship between the scalar flux and the Reynolds stress is similar in both flows.

Ashurst, William T.↗
NASA NTRS1989

More accurate predictions with transonic Navier-Stokes methods through improved turbulence modeling

Significant improvements in predictive accuracies for off-design conditions are achievable through better turbulence modeling; and, without necessarily adding any significant complication to the numerics. One well established fact about turbulence is it is slow to respond to changes in the mean strain field. With the 'equilibrium' algebraic turbulence models no attempt is made to model this characteristic and as a consequence these turbulence models exaggerate the turbulent boundary layer's ability to produce turbulent Reynolds shear stresses in regions of adverse pressure gradient. As a consequence, too little momentum loss within the boundary layer is predicted in the region of the shock wave and along the aft part of the airfoil where the surface pressure undergoes further increases. Recently, a 'nonequilibrium' algebraic turbulence model was formulated which attempts to capture this important characteristic of turbulence. This 'nonequilibrium' algebraic model employs an ordinary differential equation to model the slow response of the turbulence to changes in local flow conditions. In its original form, there was some question as to whether this 'nonequilibrium' model performed as well as the 'equilibrium' models for weak interaction cases. However, this turbulence model has since been further improved wherein it now appears that this turbulence model performs at least as well as the 'equilibrium' models for weak interaction cases and for strong interaction cases represents a very significant improvement. The performance of this turbulence model relative to popular 'equilibrium' models is illustrated for three airfoil test cases of the 1987 AIAA Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 'nonequilibrium' turbulence model is currently being applied to wing flows for which similar improvements in predictive accuracy are being realized.

Johnson, Dennis A.↗
DOE OSTI2023

Campylobacter jejuni uses energy taxis and a dehydrogenase enzyme for l -fucose chemotaxis

ABSTRACT Campylobacter jejuni is a leading cause of bacterial diarrhea worldwide, and infection in infants is associated with growth stunting in low- and middle-income countries. Approximately half of all C. jejuni isolates are asaccharolytic, while the remainder encode enzymes capable of l -fucose catabolism. Our previous study suggested that breastfed infants originating from sub-Saharan Africa and South Asia were colonized less frequently with l -fucose-metabolizing C. jejuni isolates compared to asaccharolytic strains. We also showed that FucX is an l -fucose dehydrogenase that binds to its NADP + cofactor even in the absence of l -fucose and that this enzyme is sufficient for l -fucose chemotaxis when introduced into strains lacking the metabolic pathway. This study indicates that FucX may influence l -fucose chemotaxis by increasing cellular NADPH/NADP + ratios, which subsequently affect the energy taxis components, CetABC that respond to shifting gradients of electron acceptors and donors. Furthermore, C. jejuni CetAB and CetAC can complement Escherichia coli aerotaxis. Our data support a model in which the shift in NADP + levels, impacted by the FucX dehydrogenase, affects chemotaxis in response to l -fucose. This swimming behavior can be phenocopied with a homologous l -fucose dehydrogenase from Burkholderia multivorans . Taken together, our work provides a possible explanation for why l -fucose-metabolizing C. jejuni isolates swim away from intestinal epithelial cells toward free fucose in the lumen where they are subsequently cleared by breastfed infants. IMPORTANCE In this study, we identify a separate role for the Campylobacter jejuni l -fucose dehydrogenase in l -fucose chemotaxis and demonstrate that this mechanism is not only limited to C. jejuni but is also present in Burkholderia multivorans . We now hypothesize that l -fucose energy taxis may contribute to the reduction of l -fucose-metabolizing strains of C. jejuni from the gastrointestinal tract of breastfed infants, selecting for isolates with increased colonization potential.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Extended FFT-based micromechanical formulation to consider general non-periodic boundary conditions

Here, this paper presents a new approach for applying non-periodic boundary conditions in the context of FFT-based methods to solve micromechanical problems in heterogeneous solids. The domain of the original problem is extended to satisfy the periodicity requirements at the boundary of the extended domain. The velocity constraint on the boundary of the original domain is replaced by a corresponding constraint on the velocity gradient in the extended volume, and a two-level augmented Lagrangian method is used to enforce the constraint. The proposed method is implemented as an extension of the large-strain elasto-viscoplastic FFT-based (LS-EVPFFT) model of Zecevic et al. (2022). The proposed method is verified in the cases of fully imposed velocity boundary conditions and mixed velocity/traction-free boundary conditions. The accuracy and convergence of the method are studied next, followed by applications to bending and indentation of polycrystals that illustrate the extended capabilities of the proposed formulation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Response tailoring of elasto-plastic trusses

Abstract In this work we tailor the response of trusses loaded beyond their yield limit. The truss structures are modeled using finite strain theory and rate-independent elasto-plasticity. We design trusses with a tailored mechanical response that is between “elastic” and “elastic-ideal-plastic” subject to the volume constraint. The design updates are generated by the gradient-based Method of Moving Asymptotes (MMA) solver and the sensitivities of the response functions are computed using a path-dependent adjoint sensitivity analysis. The computations are performed in Matlab.

42 ENGINEERING↗
NASA NTRS1992

Stochastic modeling of turbulent reacting flows

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

Fox, R. O.↗
NASA NTRS1995

Aluminum U-groove weld enhancement based on experimental stress analyses

Though butt-welds are among the most preferred joining methods in aerostructures because of their sealing and assembly integrity and general elastic performance; their inelastic mechanics are generally the least understood. This study investigated experimental strain distributions across a thick aluminum U-grooved weld and identified two weld process considerations for improving the multipass weld strength. The extreme thermal expansion and contraction gradient of the fusion heat input across the tab thickness between the grooves produce severe peaking, which induces bending moment under uniaxial loading. The filler strain hardening decreased with increasing filler pass sequence. These combined effects reduce the weld strength, and a depeaking index model was developed to select filler pass thicknesses, pass numbers, and sequences to improve the welding process results over the current normal weld schedule.

Verderaime, V.↗
NASA NTRS1986

Testing experience with unheated stain-gage balances in the NTF

A series of cryogenic (cryo) cycles was conducted in the cryo chamber at the National Transonic Facility (NTF) in order to identify the cause of apparent strain shifts in axial force with temperature for the Pathfinder I model and to minimize their effects. The results of the investigation indicated that the major cause of axial force end point shifts and thermal hysteresis loops was the thickness of the Teflon insulation on the instrumentation wires crossing the balance. By reducing the thickness of the insulation and the total number and size of the wires, apparent strain values were achieved for the model with instrumentation wires which were nearly identical to those for the model without wires. Because of the special design features used, the balance output was very accurate and repeatable over the entire NTF temperature range, even with balance thermal gradients as large as 64 F and transient conditions as large as 3 F/minute.

Jacobs, Peter F.↗
NASA NTRS1990

Thermomechanical deformation testing and modeling in the presence of metallurgical instabilities

A number of viscoplastic constitutive models were developed to describe deformation behavior under complex combinations of thermal and mechanical loading. Questions remain, however, regarding the validity of procedures used to characterize these models for specific structural alloys. One area of concern is that the majority of experimental data available for this purpose are determined under isothermal conditions. This experimental study is aimed at determining whether viscoplastic constitutive theories characterized using an isothermal data base can adequately model material response under the complex thermomechanical loading conditions typical of power generation service. The approach adopted was to conduct a series of carefully controlled thermomechanical experiments on a nickel-based superalloy, Hastelloy Alloy X. Previous investigations had shown that this material experiences metallurgical instabilities leading to complex hardening behavior, termed dynamic strain aging. Investigating this phenomenon under full thermomechanical conditions leads to a number of challenging experimental difficulties which up to the present work were unresolved. To correct this situation, a number of advances were made in thermomechanical testing techniques. Advanced methods for dynamic temperature gradient control, phasing control and thermal strain compensation were developed and incorporated into real time test control software. These advances allowed the thermomechanical data to be analyzed with minimal experimental uncertainty. The thermomechanical results were evaluated on both a phenomenological and microstructural basis. Phenomenological results revealed that the thermomechanical hardening trends were not bounded by those displayed under isothermal conditions. For the case of Hastelloy Alloy X (and similar dynamic strain aging materials), this strongly suggests that some form of thermomechanical testing is necessary when characterizing a thermoviscoplastic deformation model. Transmission electron microscopy was used to study the microstructural physics, and analyze the unique phenomenological behavior.

Castelli, Michael G.↗
DOE OSTI2025

Emergent topological polarization textures in relaxor ferroelectrics

Relaxor ferroelectrics underpin high-performance actuators and sensors, yet the nature of polar heterogeneities driving their broadband dielectric response remains debated. Using a unified, multimodal structural refinement framework— simultaneously fitting complementary X-ray and neutron total scattering, X-ray absorption spectra, and diffuse scattering—we reconstruct 3D mesoscale polarization maps in the classic relaxor system PbMg 1/3 Nb 2/3 O 3 –PbTiO 3 . We uncover self-organized swirling polarization textures with half-skyrmion (meron) vortices, challenging models of independent polar nanoregions. These textures, characterized by smooth changes in the polarization direction, originate from overlapping volumes in which the projections of locally correlated polarization vectors onto each volume’s long axis share the same sign. Vortex cores correlate strongly with local charge and strain gradients imposed by compositional heterogeneities. In this work, our results suggest that chemical disorder, acting via depolarizing and strain fields, stabilizes topological vortex textures of the polarization field, offering a route for engineering new dielectric and ferroelectric functionalities.

Eremenko, Maksim [National Institute of Standards ↗
NASA NTRSDate not supplied

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↗
DOE OSTI2021

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗