Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic model”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Dissociated vacancies and screw dislocations in MgO and UO 2 : atomistic modeling and linear elasticity analysis

Understanding the effect of dislocations on the mass transport in ionic ceramics is important for understanding the behavior of these materials in a variety of contexts. In particular, the dissociated nature of vacancies at screw dislocations, or more generally, at a wide range of low-angle twist grain-boundaries, has ramifications for the mechanism of defect migration and thus mass transport at these microstructural features. In this paper, a systematic study of the dissociated vacancies at screw dislocations in MgO is carried out. The important role of stress migration in the atomistic modeling study is identified. Another aspect of the current work is a rigorous treatment of the linear elasticity model. As a result, good agreement between the atomistic modeling results and the linear elasticity model is obtained. Furthermore, we demonstrate that the proposed vacancy dissociation mechanism can also be extended to more complicated ionic ceramics such as UO 2 , highlighting the generality of the mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A 3D Orthotropic Strain-Rate Dependent Elastic Damage Material Model.

A three dimensional orthotropic elastic constitutive model with continuum damage and cohesive based fracture is implemented for a general polymer matrix composite lamina. The formulation assumes the possibility of distributed (continuum) damage followed b y localized damage. The current damage activation functions are simply partially interactive quadratic strain criteria . However, the code structure allows for changes in the functions without extraordinary effort. The material model formulation, implementation, characterization and use cases are presented.

English, Shawn Allen↗
DOE OSTI2021

An optimization-based strategy for peridynamic-FEM coupling and for the prescription of nonlocal boundary conditions

We develop and analyze an optimization-based method for the coupling of a static peridynamic (PD) model and a static classical elasticity model. The approach formulates the coupling as a control problem in which the states are the solutions of the PD and classical equations, the objective is to minimize their mismatch on an overlap of the PD and classical domains, and the controls are virtual volume constraints and boundary conditions applied at the local-nonlocal interface. Our numerical tests performed on three-dimensional geometries illustrate the consistency and accuracy of our method, its numerical convergence, and its applicability to realistic engineering geometries. We demonstrate the coupling strategy as a means to reduce computational expense by confining the nonlocal model to a subdomain of interest, and as a means to transmit local (e.g., traction) boundary conditions applied at a surface to a nonlocal model in the bulk of the domain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Joint Electricity and Natural Gas Transmission Planning with Endogenous Market Feedbacks

Recent trends in gas-fired power plant installation has increased the connections between the electric power and natural gas industries. Despite these dependencies, both industries must meet commercial, political, operational and technical requirements that often force the industries to plan and operate in isolation. As a result, undesired situations may arise, such as those experienced by both systems during the winter of 2013/2014 in the northeastern United States. In this work, we consider the technical challenges and present a Combined Electricity and Gas Expansion (CEGE) planning model. The CEGE model minimizes the cost of meeting gas and electricity demand during high stress conditions and introduces an elasticity model for analysis of gas-price volatility caused by congestion. The underlying mathematical formulation considers recent advances in convex approximations to make the problem computationally tractable when applied to large-scale CEGE instances. We conduct an indepth analysis on integrated test systems that include the New England area.

03 NATURAL GAS↗
DOE OSTI2018

Stressor-layer-induced elastic strain sharing in SrTiO 3 complex oxide sheets

A precisely selected elastic strain can be introduced in submicron-thick single-crystal SrTiO 3 sheets using a silicon nitride stressor layer. A conformal stressor layer deposited using plasma-enhanced chemical vapor deposition produces an elastic strain in the sheet consistent with the magnitude of the nitride residual stress. Synchrotron x-ray nanodiffraction reveals that the strain introduced in the SrTiO 3 sheets is on the order of 10 -4 , matching the predictions of an elastic model. Using this approach to elastic strain sharing in complex oxides allows the strain to be selected within a wide and continuous range of values, an effect not achievable in heteroepitaxy on rigid substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio calculations of third-order elastic coefficients

Third-order elasticity (TOE) theory predicts strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the natural, or zero pressure, reference state. This definition is inconvenient for predictions of SOECs under significant initial stresses. Therefore, when TOE theory is necessary to study the strain dependence of elasticity, the seismological community has resorted to an empirical version of the theory. This study reviews the thermodynamic definition of the third-order elastic tensor and proposes using an “effective” third-order elastic tensor. An explicit expression for the effective third-order elastic tensor is given and verified. Additionally, we extend the ab initio approach to calculate third-order elastic tensors under finite pressure and apply it to two cubic systems, namely, NaCl and MgO. As applications and validations, we evaluate (a) strain-induced changes in SOECs and (b) pressure derivatives of SOECs based on ab initio calculations. Good agreement between third-order elasticity-based predictions and numerically calculated values confirms the validity of our theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Residual Stress and Opening-Mode Fracture Analysis of Multilayered Structures Subjected to Thermal Loading

Temperature change after formation commonly results in thermal residual stress in multilayered structures due to the different thermal and mechanical properties of each layer. In this paper, a three-dimensional (3D) elastic model is developed to study the residual stress and opening-mode fractures (OMFs) in a multilayered structure consisting of arbitrary number of layers under temperature change. The general solution of displacement field in the multilayered structure is derived by solving the elastic boundary value problem. In order to verify the proposed model, the elastic field in the advanced polymeric solar reflectors that consist of four layers is solved by applying the present model and compared with the finite-element (FE) simulation. In addition, parametric studies are conducted to investigate the effect of the thickness ratio between each layer on the accuracy of the developed model. As a result, based on the obtained elastic field, the fracture energy release rate (ERR) in the surface layer of the advanced polymeric reflector is obtained and used to study the fracture initiation, infilling, and saturation successfully.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Implementation of the Hybrid-Elastic-Plastic Material Model in DYNA3D

The Hybrid–Elastic–Plastic (HEP) material model was developed by the U.S. Army Engineer Research and Development Center (ERDC) and L3/Titan Research and Technology (Akers, et al., 1995; Zimmerman et al., 1987) to simulate geological materials. The present effort implements HEP into LLNL’s explicit finite element code DYNA3D (Lin and Zywicz, 2010) [version 11.2 – revision 1.1921] as material model 72. The HEP material model includes a library of 26 sets of pre–fit material parameters and a mechanism to read HEP material parameters from auxiliary files. The library and external file coding were included in the DYNA3D implementation.

42 ENGINEERING↗
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.↗
NASA NTRS2009

Solar Thermal Vacuum Test of Deployable Astromesh Reflector

On September 10, 2008, a 36-hour Solar Thermal Vacuum Test of a 5m deployable mesh reflector was completed in JPL's 25' Space Simulator by the Advanced Deployable Structures Group at JPL. The testing was performed under NASA's Innovative Partnership Program (IPP) as a risk reduction effort for two JPL Decadal Survey Missions: DESDynI and SMAP. The 5.0 m aperture Astromesh reflector was provided by Northrop Grumman Aerospace Systems (NGAS) Astro Aerospace, our IPP industry partner. The testing utilized a state-of-the-art photogrammetry system to measure deformation of the reflector under LN2 cold soak, 0.25 Earth sun, 0.5 sun and 1.0 sun. An intricate network of thermocouples (approximately 200 in total) was used to determine the localized temperature across the mesh as well as on the perimeter truss of the reflector. Half of the reflector was in a fixed shadow to maximize thermal gradients. A mobility system was built for remotely actuating the cryo-vacuum capable photogrammetry camera around the circumference of the Solar Simulator. Photogrammetric resolution of 0.025 mm RMS (0.001") was achieved over the entire 5 meter aperture for each test case. The data will be used for thermo-elastic model correlation and validation, which will benefit the planned Earth Science Missions.

Photogrammetry↗
DOE OSTI2014

Development of an Offshore Direct-Drive Wind Turbine Model by Using a Flexible Multibody Simulation (Poster)

Modern wind turbines are complex, highly-coupled systems. The dynamic interaction between various components is especially pronounced for multi-megawatt wind turbines. As a result, design process is generally split in several phases. First step consists of creating a global aero-elastic model that includes essential dynamics of structural components using the minimum-possible number of degrees of freedom (d.o.f.). The most important simplifications concern drivetrain and rotor-nacelle assembly (RNA). This approach has been shown valid for several wind turbine configurations. Nevertheless, with increasing size of wind turbines, any simplified design approach must be validated. The present work deals with the comparison and validation of the two modeling approaches for directdrive offshore wind turbines. ARNA/drivetrain model idealized as collection of lumped masses and springs is compared to a detailed Finite Element Method (FEM) based model. The comparison between models focuses on dynamic loads concerning drivetrain system. The comparison is performed in several operational conditions in order to explore the range of validity of the simplified model. Finally, the paper proposes a numerical-based workflow to assess the validity of simplified models of RNA/drivetrain in an aero-elastic global WT model.

17 WIND ENERGY↗
NASA NTRS1973

Mathematical modeling of spinning elastic bodies for modal analysis.

The problem of modal analysis of an elastic appendage on a rotating base is examined to establish the relative advantages of various mathematical models of elastic structures and to extract general inferences concerning the magnitude and character of the influence of spin on the natural frequencies and mode shapes of rotating structures. In realization of the first objective, it is concluded that except for a small class of very special cases the elastic continuum model is devoid of useful results, while for constant nominal spin rate the distributed-mass finite-element model is quite generally tractable, since in the latter case the governing equations are always linear, constant-coefficient, ordinary differential equations. Although with both of these alternatives the details of the formulation generally obscure the essence of the problem and permit very little engineering insight to be gained without extensive computation, this difficulty is not encountered when dealing with simple concentrated mass models.

Likins, P. W.↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗
NASA NTRS1977

Nonlinear viscoelasticity and relaxation phenomena of polymer solids

In the light of a three-chain model of statistical network theories of rubberlike elastic models, it is assumed that the free energy function of incompressible viscoelastic polymer solids is a separable, symmetric function of the principle stretch ratios and the hidden thermodynamic coordinates along the same directions. This assumption leads to a characterization of those viscoelastic polymer solids which exhibit the property of factorizability between the time and strain functions.

Peng, S. T. J.↗
NASA NTRSDate not supplied

Enhancements, Verification, and VMS Integration of VTOL Concept Vehicle Simulation Models

Advanced Rotorcraft Technology (ART) and the NASA Ames Aeromechanics branch have jointly developed FLIGHTLAB simulation models for Advanced Air Mobility (AAM) VTOL concept vehicles. The overarching purpose of the simulation model development is to establish a set of well defined reference vehicles for FLIGHTLAB users and the rotorcraft community. The ongoing research effort and enhancement of these AAM simulation models to fulfill the role of quality reference vehicles is this paper’s focus. The content of this paper expands on the established characteristics of these AAM models in three primary areas. First, enhancement of the lift+cruise and tiltwing models with elastic airframe properties is discussed. The process of setting up the elastic airframe model in FLIGHTLAB, as well as the impacts on flight characteristics are explained. The introduction of the elastic airframe modeling allows these models to be used in flight dynamics, loads, and vibration analysis of the configuration designs. Next, linear model generation from the enhanced simulation model is covered. Confirming the validity of the linearized models is of importance, as these linear models are utilized for flight control design and tuning for these experimental configurations. For the final focus, the progress towards implementation of these models into the NASA Ames Vertical Motion Simulator (VMS) is described. This task seeks to demonstrate the procedures of integrating a FLIGHTLAB flight simulation in the VMS environment, test fully integrated simulation with communication between flight dynamics, control, and propulsion models, and explore the essential aspects of simulation model integration in a full flight simulator environment. This includes I/O definition, initialization, trim, flying, etc. By expanding the capabilities of the AAM simulation models, they continue to develop as valuable and approachable modeling references.

VMS↗
DOE OSTI2016

Elastically-homogeneous lattice models of damage in geomaterials

This paper involves the development of the auxiliary stress approach for producing elastically-homogeneous lattice models of damage in geomaterials. The lattice models are based on random, three-dimensional assemblages of rigid-body-spring elements. Unlike conventional lattice or particle models, the elastic constants of a material (e.g., Young’s modulus and Poisson’s ratio) are represented properly in both global and local senses, without any need for calibration. The proposed approach is demonstrated and validated through analyses of homogeneous and heterogeneous systems under uni- and tri-axial loading conditions. Comparisons are made with analytical solutions and finite element results. Thereafter, the model is used to simulate a series of standard laboratory tests: (a) split-cylinder tests, and (b) uniaxial compressive tests of sedimentary rocks at the Horonobe Underground Research Laboratory in Hokkaido, Japan. Model inputs are based on physical quantities measured in the experiments. Finally, the simulation results agree well with the experimental results in terms of pre-peak stress-strain/displacement responses, strength measurements, and failure patterns.

3D irregular lattice↗
DOE OSTI2017

Transient Simulation of the Multi-SERTTA Experiment with MAMMOTH

This work details the MAMMOTH reactor physics simulations of the Static Environment Rodlet Transient Test Apparatus (SERTTA) conducted at Idaho National Laboratory in FY-2017. TREAT static-environment experiment vehicles are being developed to enable transient testing of Pressurized Water Reactor (PWR) type fuel specimens, including fuel concepts with enhanced accident tolerance (Accident Tolerant Fuels, ATF). The MAMMOTH simulations include point reactor kinetics as well as spatial dynamics for a temperature-limited transient. The strongly coupled multi-physics solutions of the neutron flux and temperature fields are second order accurate both in the spatial and temporal domains. MAMMOTH produces pellet stack powers that are within 1.5% of the Monte Carlo reference solutions. Some discrepancies between the MCNP model used in the design of the flux collars and the Serpent/MAMMOTH models lead to higher power and energy deposition values in Multi-SERTTA unit 1. The TREAT core results compare well with the safety case computed with point reactor kinetics in RELAP5-3D. The reactor period is 44 msec, which corresponds to a reactivity insertion of 2.685% delta k/k$. The peak core power in the spatial dynamics simulation is 431 MW, which the point kinetics model over-predicts by 12%. The pulse width at half the maximum power is 0.177 sec. Subtle transient effects are apparent at the beginning insertion in the experimental samples due to the control rod removal. Additional difference due to transient effects are observed in the sample powers and enthalpy. The time dependence of the power coupling factor (PCF) is calculated for the various fuel stacks of the Multi-SERTTA vehicle. Sample temperatures in excess of 3100 K, the melting point UO$_2$, are computed with the adiabatic heat transfer model. The planned shaped-transient might introduce additional effects that cannot be predicted with PRK models. Future modeling will be focused on the shaped-transient by improving the control rod models in MAMMOTH and adding the BISON thermo-elastic models and thermal-fluids heat transfer.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1986

Viscoelastic deformation near active plate boundaries

Model deformations near the active plate boundaries of Western North America using space-based geodetic measurements as constraints are discussed. The first six months of this project were spent gaining familarity with space-based measurements, accessing the Crustal Dynamics Data Information Computer, and building time independent deformation models. The initial goal was to see how well the simplest elastic models can reproduce very long base interferometry (VLBI) baseline data. From the Crustal Dynamics Data Information Service, a total of 18 VLBI baselines are available which have been surveyed on four or more occasions. These data were fed into weighted and unweighted inversions to obtain baseline closure rates. Four of the better quality lines are illustrated. The deformation model assumes that the observed baseline rates result from a combination of rigid plate tectonic motions plus a component resulting from elastic strain build up due to a failure of the plate boundary to slip at the full plate tectonic rate. The elastic deformation resulting from the locked plate boundary is meant to portray interseismic strain accumulation. During and shortly after a large interplate earthquake, these strains are largely released, and points near the fault which were previously retarded suddenly catch up to the positions predicted by rigid plate models. Researchers judge the quality of fit by the sum squares of weighted residuals, termed total variance. The observed baseline closures have a total variance of 99 (cm/y)squared. When the RM2 velocities are assumed to model the data, the total variance increases to 154 (cm/y)squared.

Ward, S. N.↗