SEARCH · Engineering Papers

Results for “model deformation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

All sources NASA NTRS DOE OSTI

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↗

NASA NTRS1990

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗

NASA NTRS2003

Small Scale Polygons and the History of Ground Ice on Mars

Recent progress on polygon modeling has focused on the diameter and surface relief that we expect of thermal-contraction polygons in martian permafrost. With this in mind, we developed a finite-element model of thermal-contraction-crack behavior in permafrost in a martian climate. This model was generated from a finite element code by Jay Melosh (called TECTON) originally developed for terrestrial and planetary crustal-deformation studies. We adapted this model to martian permafrost by including time (and temperature) dependent rheologies, boundary conditions, and isotropic thermal-contraction, as well as several small adaptations to a martian environment. We tested our model extensively, including comparison to an analytic solution of pre-fracture stress. We recently published an analysis of two potential sources of water for forming the recent gullies. In this work we first evaluated the potential for near-surface ground ice (in the top meter or so of soil) to melt under conditions of solar heating on sloped surfaces at high obliquity, utilizing both thermal and diffusion-based ground-ice-stability models; our results suggested that the ground ice will sublimate, and the ice table will recede to greater depths before the melting temperature can be reached. An exception can occur only for extremely salt-rich ice, depressing the freezing point.

Mellon, Michael T.↗

NASA NTRS1990

Nonlinear analysis of AS4/PEEK thermoplastic composite laminate using a one parameter plasticity model

A one-parameter plasticity model was shown to adequately describe the orthotropic plastic deformation of AS4/PEEK (APC-2) unidirectional thermoplastic composite. This model was verified further for unidirectional and laminated composite panels with and without a hole. The nonlinear stress-strain relations were measured and compared with those predicted by the finite element analysis using the one-parameter elastic-plastic constitutive model. The results show that the one-parameter orthotropic plasticity model is suitable for the analysis of elastic-plastic deformation of AS4/PEEK composite laminates.

Sun, C. T.↗

NASA NTRS2001

Fully Coupled Micro/Macro Deformation, Damage, and Failure Prediction for SiC/Ti-15-3 Laminates

The deformation, failure, and low cycle fatigue life of SCS-6/Ti-15-3 composites are predicted using a coupled deformation and damage approach in the context of the analytical generalized method of cells (GMC) micromechanics model. The local effects of inelastic deformation, fiber breakage, fiber-matrix interfacial debonding, and fatigue damage are included as sub-models that operate on the micro scale for the individual composite phases. For the laminate analysis, lamination theory is employed as the global or structural scale model, while GMC is embedded to operate on the meso scale to simulate the behavior of the composite material within each laminate layer. While the analysis approach is quite complex and multifaceted, it is shown, through comparison with experimental data, to be quite accurate and realistic while remaining extremely efficient.

Bednarcyk, Brett A.↗

DOE OSTI2025

A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials

Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2023

Precision test of gauge/gravity duality in D0-brane matrix model at low temperature

We test the gauge/gravity duality between the matrix model and type IIA string theory at low temperatures with unprecedented accuracy. To this end, we perform lattice Monte Carlo simulations of the Berenstein-Maldacena-Nastase (BMN) matrix model, which is the one-parameter deformation of the Banks-Fischler-Shenker-Susskind (BFSS) matrix model, taking both the large N and continuum limits. We leverage the fact that sufficiently small flux parameters in the BMN matrix model have a negligible impact on the energy of the system while stabilizing the flat directions so that simulations at smaller N than in the BFSS matrix model are possible. Hence, we can perform a precision measurement of the large N continuum energy at the lowest temperatures to date. The energy is in perfect agreement with supergravity predictions including estimations of α'-corrections from previous simulations. At the lowest temperature where we can simulate efficiently (T = 0.25λ 1/3 , where λ is the ’t Hooft coupling), the difference in energy to the pure supergravity prediction is less than 10%. Furthermore, we can extract the coefficient of the 1/N 4 corrections at a fixed temperature with good accuracy, which was previously unknown.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Non-Synchronous Rotation on Europa Driven by Ocean Currents

It has been suggested that the ice shell of Jupiter's moon Europa may drift non-synchronously due to tidal torques. Here we argue that torques applied by the underlying ocean are also important and can result in non-synchronous rotation. The resulting spin rate can be slightly slower than the synchronous angular rate that would have kept the same point of the ice shell facing Jupiter. We develop an ice shell rotation model, driven by ocean stress calculated using a high-resolution state-of-the-art ocean general circulation model, and take into account the viscoelastic deformation of the ice shell. We use the ice shell model results together with observed limits on the ice shell drift speed to constrain ice shell parameters such as effective viscosity, which is currently uncertain by at least four orders of magnitude. Our results suggest, at best, sluggish ice shell convection. Depending on the relaxation time scale of the ice shell and on the ocean currents, the ice shell may exhibit negligible drift, constant drift, or oscillatory drift superimposed on random fluctuations. The expected rotation rate exceeds ~30 m/yr; future spacecraft observations can be used to test these predictions and yield insight into the properties of the ice shell and underlying ocean.

15 GEOTHERMAL ENERGY↗

NASA NTRS1989

Numerical analysis of supersonic flow through oscillating cascade sections by using a deforming grid

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small-perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small-perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a section of the rotor on a supersonic throughflow compressor designed at NASA Lewis Research Center. Preliminary results indicate that two-dimensional, flat plate analysis predicts conservative flutter boundaries.

Huff, Dennis L.↗

NASA NTRS1989

Numerical analysis of supersonic flow through oscillating cascade sections by using a deforming grid

Huff, Dennis L.↗

DOE OSTI2021

The Corrected Distortion model for Lagrangian spray simulation of transcritical fuel injection

In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.

42 ENGINEERING↗

DOE OSTI2023

From Well Log to Formation Model: A Novel Laboratory Calibrated Methodology with Demonstration

This work demonstrates how the characterization and modeling of both elastic and creep properties are essential to describe zones in layered rock formations as either low-stress targets for stimulation or high-stress barriers to fracture growth. Prediction of fracture height is critical for designing stimulation operations in oil and gas wells. Ideally, fractures are placed in target zones which will produce hydrocarbons and should not propagate into zones expected to be unproductive or to produce unwanted fluids such as water which in turn must be treated and/or disposed. The essential task in designing stimulation plans is predicting which zones have low horizontal stresses and which will be high-stress barriers to fracture growth. Despite this importance, there are gaps in current knowledge and a complete workflow from laboratory characterization to a finite element model which includes time dependent rock deformation is required. While the research and methodology presented here also have application to CO 2or hydrogen storage, wastewater injection, and geothermal applications, the focus will be on hydrocarbon extraction. This thesis presents the results of a characterization-to-prediction workflow for the Caney shale, which is an emerging hydrocarbon resource in Oklahoma, USA. It begins with an investigation to enable critical evaluation of the Caney zonation into nominally “brittle” and “ductile” zones based on properties observed from well logs. It shows none of the zones are consistently “brittle” or “ductile” mechanical behavior based on the variety of definitions of these terms. However, the nominally ductile zones are weaker and more prone to creep. A laboratory investigation of samples including strength, elastic, and creep properties, is then used in a finite element model of stress evolution. The model includes both elastic deformation and viscoplastic creep. Results predict the least creep-prone layers to have the lowest horizontal stresses, therefore comprising hydraulic fracturing targets. The most creep-prone layers attain a horizontal stress similar to the vertical stress and therefore are predicted to be high stress barriers to hydraulic fracture stimulation. In addition to defining stimulation target intervals, the model shows how as tectonic strain rate increases, there is a transition from creep-dominated stresses to stresses dominated by elasticity.

Benge, Margaret↗

NASA NTRS1992

Prediction of dislocation generation during Bridgman growth of GaAs crystals

Dislocation densities are generated in GaAs single crystals due to the excessive thermal stresses induced by temperature variations during growth. A viscoplastic material model for GaAs, which takes into account the movement and multiplication of dislocations in the plastic deformation, is developed according to Haasen's theory. The dislocation density is expressed as an internal state variable in this dynamic viscoplastic model. The deformation process is a nonlinear function of stress, strain rate, dislocation density and temperature. The dislocation density in the GaAs crystal during vertical Bridgman growth is calculated using a nonlinear finite element model. The dislocation multiplication in GaAs crystals for several temperature fields obtained from thermal modeling of both the GTE GaAs experimental data and artificially designed data are investigated.

Tsai, C. T.↗

NASA NTRS2020

Comparison of Vertical Surface Deformation Estimates Derived From Space-Based Gravimetry, Ground-Based GPS, and Model-Based Hydrologic Loading Over Snow-Dominated Watersheds in the United States

Spatiotemporal variability in Earth's terrestrial water storage (TWS) causes changes in surface deformation. The potential for using ground-based Global Positioning System (GPS) vertical displacement observations for estimating TWS is explored through a comparison of vertical displacements derived from space-based gravimetric retrievals, ground-based GPS, and model-based hydrologic estimates. The study presented here focuses on two snow-dominated basins in the Western United States for the years 2003–2016. Seasonal variations are observed in the vertical displacements derived from all three data sets, and the variation is coherent with the changes in hydrologic loading. Good consistency is observed between any two of the three data sets with gravimetric retrievals and hydrologic model estimates providing the highest level of agreement (i.e., all examined stations with correlation coefficient R > 0.70). Vertical displacements derived from gravimetric retrievals and ground-based GPS yielded R > 0.70 for more than 89% of the stations. In addition, it is found that both GPS-derived and space-based, gravimetry-derived vertical displacements clearly reflected the impact of climate variation (i.e., heavy precipitation during 2010–2011 winter followed by prolonged drought). Vertical displacements derived from the hydrologic model highlighted the relatively large precipitation convergence phase during late 2010 to early 2011 at some stations but not the prolonged drought that followed. The results indicate that ground-based GPS observations of vertical displacement have the capability to capture variations in TWS changes, which can be systematically merged in conjunction with Gravity Recovery and Climate Experiment (GRACE) into a land surface model to improve TWS estimates in a follow-up study.

Gaohong Yin↗

NASA NTRS1987

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗

NASA NTRS1987

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

Huff, Dennis L.↗

DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗

NASA NTRS1992

Hierarchical tapered bar elements undergoing axial deformation

A method is described to model the dynamics of tapered axial bars of various cross sections based on the well-known Craig/Bampton component mode synthesis technique. This element is formed in terms of the static constraint modes and interface restrained normal modes. This is in contrast with the finite elements as implemented in NASTRAN where the interface restrained normal modes are neglected. These normal modes are in terms of Bessel functions. Restoration of a few of these modes leads to higher accuracy with fewer generalized coordinates. The proposed models are hierarchical so that all lower order element matrices are embedded in higher order element matrices. The advantages of this formulation compared to standard NASTRAN truss element formulation are demonstrated through simple numerical examples.

Ganesan, N.↗