Engineering Papers⌕ Search

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.

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

At least 559 records · Page 31

Advances in Integrally Formed Metallic Structures

NASA Langley Research Center (LaRC) is applying flow forming technology to aerospace aluminum alloys for the manufacture of integrally stiffened cylinders that meet aircraft fuselage or launch vehicle specifications. High-strength 2xxx series alloys, in various temper/product form conditions, are being compared with medium-strength, very formable AA 6061 as the benchmark. Screening of candidate materials traditionally relies on empirical forming trials, but more recently on metrics derived from constant nominal strain rate tensile testing. In this work, the established hierarchy of formability metrics is adopted as the baseline; tensile area reduction, work hardening exponent, and modulus of resilience. The suite of metrics can be expanded based on phenomenological and mechanistic considerations of plastic flow instabilities to incorporate anisotropy ratio, rate of work hardening, and dynamic strain aging effects. Separate research thrusts at LaRC are aimed at formulating a protocol for scalable production of fully formed, defect-free components: 1) quantifying the extended set of formability metrics from constant and jump strain rate tensile tests with digital image correlation (DIC); 2) using statistically designed experiments on a vertical spin/flow forming facility (WF VUD-600®) to isolate the effect of convenient process variables, such as axial feed rate and thickness reduction; and 3) predicting the effect of feed ratio (roller translation relative to mandrel rotation speeds) on material response via process modeling using the DEFORM® simulation software. Current results will be discussed in terms of correlating intrinsic material properties with the bulk workability characteristics to successfully fabricate cylinders with complex cross-sections.

Aluminum↗

Evidence of Hexadecapole Deformation in Uranium-238 at the Relativistic Heavy Ion Collider

State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic 238 U + 238 U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding 238 U ions. Here, we show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like 238 U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of 238 U on high-energy collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Lie-algebraic Kähler sigma models with U(1) isotropy

We discuss various questions that emerge in connection with the Lie-algebraic deformation of the sigma model in two dimensions. First, we supersymmetrize the original model endowing it with the minimal and extended supersymmetries. Then we derive the general hypercurrent anomaly in both cases. In the latter case this anomaly is one-loop but is somewhat different from the standard expressions one can find in the literature because the target manifold is nonsymmetric. We also show how to introduce the twisted masses and the term, and study the Bogomol’nyi–Prasad–Sommerfield equation for instantons, in particular the value of the topological charge. Then we demonstrate that the second loop in the function of the nonsupersymmetric Lie-algebraic sigma model is due to an infrared effect. To this end we use a supersymmetric regularization. We also conjecture that the above statement is valid for higher loops too, similar to the parallel phenomenon in four-dimensional super-Yang-Mills. In the second part of the paper we develop a special dimensional reduction—namely, starting from the two-dimensional Lie-algebraic model we arrive at a quasi-exactly solvable quantum-mechanical problem of the Lamé type.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dislocation line tension model to improve understanding of the effects of hydrogen on the deformation of structural materials

Decarbonization efforts highlight hydrogen as an attractive alternative to fossil fuels, but its tendency to embrittle structural metals demands careful consideration when designing hydrogen infrastructure. Moreover, the mechanisms by which hydrogen degrades these materials are still being elucidated. The current work develops new computational tools to quantify the different contributions of hydrogen to the energy barrier of cross-slip, a key deformation mechanism. Novel features are implemented to a line tension model, which include the use of non-singular dislocation interactions, character-dependent dislocation energies and simulations of the constriction configurations. A new molecular dynamics technique is developed to calculate the interaction energy between the partials of a dissociated dislocation via fixing the centers of mass of the regions below and above the Shockley partials and performing time-averaged calculations. Hydrogen is found to impact the stacking fault width of dislocations in different ways depending on their characters: it decreases for dislocations with a character θ > 30°, remains unchanged for θ = 30° and increases for θ < 30°. The latter regime is a newly identified mechanism by which hydrogen inhibits cross-slip. Moreover, formation of nano-hydrides is predicted to occur around screw dislocations for high hydrogen concentrations, a phenomenon previously identified only in dislocations with an edge component. If nano-hydrides develop, their influence extending the equilibrium stacking fault width and increasing both the constriction and cross-slip energy barriers dominate over all other hydrogen contributions. The theory and tools developed will pave the way towards a comprehensive understanding of hydrogen-dislocation interactions in structural metals.

08 HYDROGEN↗

Parametric Finite Element Analysis of Naturally Corroded Steel Specimens Using 3D Surface Laser Scans

Corrosion is considered a uniform thickness reduction design guideline of the maritime industry. However, additionally, the corroded and irregular morphology of the surface affects the steel's load-bearing capacity and its impact on the strength and elongation behaviour of the steel is not yet fully understood. These effects on the local behaviour of steel structures under tensile loading were investigated with tensile tests on naturally corroded steel specimens and nonlinear finite element simulations including the corroded surface morphology with a uniform surface idealation. The models also include the deformed specimen shape. The developed approach led to highly accurate parametric finite element models predicting the ultimate tensile strength and longitudinal position of fracture. The results show that all included aspects are essential for accurate simulations, while solely the maximum available surface resolution was not as decisive.

corrosion↗

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗

Implementation of a Tabulated Failure Model Into a Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

polymer matrix composites↗

Wind-tunnel/flight correlation study of aerodynamic characteristics of a large flexible supersonic cruise airplane CXB-70-1). 1: Wind-tunnel tests of a 0.03-scale model at Mach numbers from 0.6 to 2.53

The longitudinal and lateral forces and moments for a 0.03 scale deformed rigid, static force model of the XB-70-1 airplane were determined. Control effectiveness was determined for the elevon in pitch and roll, for the canard, and for the rudders. Component effects of the canard, deflected with tips, variable position canopy, bypass doors, and bleed dump fairing were measured. The effects of small variations in inlet mass flow ratio and small amounts of asymmetric deflection of the wing tips were assessed.

Daugherty, J. C.↗

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗

The New IERS Special Bureau for Loading (SBL)

Currently, the establishment of the International Earth Rotation Service (IERS) Special Bureau for Loading (SBL) is in progress as part of the IERS Global Geophysical Fluids Center (GGFC). The main purpose of the SBL is to provide reliable, consistent model predictions of loading signals that have been thoroughly tested and validated. The products will describe at least the surface deformation, gravity signal and geo-center variations due to the various surface loading processes in reference frames relevant for direct comparison with existing geodetic observing techniques. To achieve these goals, major scientific advances are required with respect to the Earth model, the theory and algorithms used to model deformations of the Earth as well as improvements in the observational data related to surface loading.

vanDam, Tonie↗

Distributed deformation and block rotation in three dimensions

A model is developed that attempts to explain distributed deformation, particularly the relationship between fault slip and rotation of faults and blocks in a 3D stress field. Regions of distributed deformation (such as Southern California) are organized in complex arrays of contemporaneously active block-faulted domains. A 3D block rotation model is proposed here that tracks the orientation of blocks and their bounding faults during rotation. It is shown how only one set of faults, preexisting and rotating in an irrotational strike-slip stress field, can account for the three major phases of deformation observed in the Western Transverse Range domain, Southern California.

Scotti, Oona↗

Finite Element Methods for real-time Haptic Feedback of Soft-Tissue Models in Virtual Reality Simulators

We have applied the linear elastic finite element method to compute haptic force feedback and domain deformations of soft tissue models for use in virtual reality simulators. Our results show that, for virtual object models of high-resolution 3D data (>10,000 nodes), haptic real time computations (>500 Hz) are not currently possible using traditional methods. Current research efforts are focused in the following areas: 1) efficient implementation of fully adaptive multi-resolution methods and 2) multi-resolution methods with specialized basis functions to capture the singularity at the haptic interface (point loading). To achieve real time computations, we propose parallel processing of a Jacobi preconditioned conjugate gradient method applied to a reduced system of equations resulting from surface domain decomposition. This can effectively be achieved using reconfigurable computing systems such as field programmable gate arrays (FPGA), thereby providing a flexible solution that allows for new FPGA implementations as improved algorithms become available. The resulting soft tissue simulation system would meet NASA Virtual Glovebox requirements and, at the same time, provide a generalized simulation engine for any immersive environment application, such as biomedical/surgical procedures or interactive scientific applications.

Frank, Andreas O.↗

Sensitivity Study of Multiscale and Phenomenological Elasto-Viscoplastic Grade 91 Material Models for Component-Scale Response

Many advanced nuclear reactor concepts currently being developed are targeting higher operating temperatures relative to the current fleet of light water nuclear reactors, for efficiency gains and other operational considerations. The design of high temperature structural components with reliable long-term operational performance will depend on material models that accurately capture the inelastic deformation mechanisms active in these environments. In this work, we perform a detailed parameter sensitivity analysis of two unified elasto-viscoplastic Grade 91 material models capable of capturing long term high temperature creep deformation. The first model is a phenomelogical material model from the Nuclear Engineering Material Library (NEML) developed at Argonne National Lab. The NEML model parameters and their uncertainty were fit to a range of Grade 91 experimental data using Bayesian Markov Chain Monte Carlo analysis. The second model is a LAROMance data-driven surrogate material model developed at Los Alamos National Lab. The LAROMance model is fit to a large database of responses produced by a mechanistic crystal plasticity based polycrystal model. Parameters for the LAROMance surrogate material model reflect the pedigree of the Grade 91 microstructure. Both material models have been integrated into the Grizzly code, based on the open-source MOOSE multiphysics simulation framework, to simulate both the progression of aging mechanisms and the effects of that aging on nuclear power plant structures. Grizzly is used analyze a three-dimensional Grade 91 piping system to compare the long-term inelastic response predicted by these two fundamentally different models and assess the sensitivity of the material model input parameters on this quantity of interest.

42 ENGINEERING↗

Solution and sensitivity analysis of nonlinear equations using a hypercomplex-variable Newton-Raphson method

Here, the classical Newton-Raphson (NR) method for solving nonlinear equations is enhanced in two ways through the use of hypercomplex variables and algebra. In particular, i) the Jacobian is computed in a highly accurate and automated way, and ii) the derivative of the solution to the nonlinear equations is computed with respect to any parameter contained within the system of equations. These advances provide two significant enhancements in that it is straightforward to provide an accurate Jacobian and to construct a reduced order model (ROM) of arbitrary order with respect to any parameter of the system. The ROM can then be used to approximate the solution for other parameter values without requiring additional solutions of the nonlinear equations. Several case studies are presented including 1D and 2D academic examples with fully functioning Python code provided. Additionally, a case of study of the catenary of an elastic cable subject to its own weight and a vertical point load. Derivatives up to 10th order were computed with respect to material, loading, and geometrical parameters. The derivatives were used to generate reduced order models of the cable deformation and reaction forces at its ends with respect to multiple input parameters. Results show that from a single hypercomplex evaluation of the cable under a single vertical point load, it is possible to generate an accurate reduced order model capable of predicting the cable deformation with 1.5 times the load in the opposite direction and with 3.5 times the load in the same direction without resolving the system of equations.

97 MATHEMATICS AND COMPUTING↗

Velocity distributions of ion beams in the plasma sheet boundary layer

The value of measured ion-beam velocity distributions as diagnostic tools to characterize the propagation and acceleration of particle beams in the plasma-sheet boundary layer is assessed. Models based on adiabatic deformation of a flowing Maxwellian, acceleration by field-aligned potentials, and current-sheet acceleration are fitted to observational data, and the results are presented in extensive diagrams and graphs. The data are found to be consistent with models involving field-aligned or cross-tail potential drops, but not with models based solely on the adiabatic deformation of initially isotropic distributions.

Eastman, T. E.↗

Evolution of attached and detached slabs and their associated mantle dynamics

Over the two years of the NASA grant, this project has produced a significant amount of research results related to the plate subduction process and the surface crustal deformation at convergent boundaries (i.e., above subduction zones). While some research objectives are completely accomplished, other research tasks remain active and continue to be investigated at present. A steady state analytic thermal model for subducting slabs was used to examine the torques acting on a descending slab. It is found that gravitational torque vanishes when a slab is dipping either vertically or horizontally, unlike previous studies indicating that the magnitude of gravitational torque decreases as dip angle increases. Subsequently, a new time-dependent, analytic thermal model for a subducting slab was developed. The new model enables us to study transient phenomena associated with plate subduction analytically. On the basis of this model, the nature of slab dip angles was evaluated. Slab dip angles are found to be transient features. As they penetrate into the mantle and increase their lengths, the associated gravitational torque also increases resulting in a downward pulling of the slab to the steeper dip angle. This is especially true once a slab penetrates the olivine-spinel phase boundary at about 400 km depth. However, if the phase transformation does not follow the equilibrium condition, the gravitational torque may have a different behavior. This problem was investigated. Except for fast descending slabs, non-equilibrium phase transformation can only slow down the transient increase of slab dip angles discussed earlier. Its effect is not sufficiently strong to reverse the downward pulling for most of the slabs. However, when slabs subducting at 10 cm/yr or faster, a sufficient amount of metastable olivine can exist beneath 400 km. Because of its low density compared with the surrounding spinel, an upward buoyancy is produced resulting in an upward bending of the slab and possibly an upward rotation of the slab such that smaller dip angles are formed. Seismic studies of the Japanese Slab seem to support this interpretation. The development of oroclinal geometries at convergent boundaries was also examined to study plate obduction which is an important ingredient to the initiation of plate subduction. Although the study suggests that surface features are better modeled by block models, the large scale deformation can be adequately studied by viscous models. Such a model is now under development to complete our original objective to study the initiation of plate subduction. Finally, a three-dimensional, finite element, spherical convective model is developed to study dynamic plate subductions. The model development is now complete and it is being tested to ensure its proper operation. The model is able to generate convection results with a viscosity contrast of about 100. Our research continues to push the viscosity contrast to a level that is appropriate for a subducting slab.

Hsui, Albert T.↗

Volumetric Density-Equalizing Reference Map with Applications

The density-equalizing map, a technique developed for cartogram creation, has been widely applied to data visualization but only for 2D applications. Here, we propose a novel method called the volumetric density-equalizing reference map (VDERM) for computing density-equalizing map for volumetric domains. Given a prescribed density distribution in a volumetric domain in $\mathbb{R}^3$, the proposed method continuously deforms the domain, with different volume elements enlarged or shrunk according to the density distribution. With the aid of the proposed method, medical and sociological data can be visualized via deformations of 3D objects. The method can also be applied to adaptive remeshing and shape modeling. Furthermore, by exploiting the time-dependent nature of the proposed method, applications to shape morphing can be easily achieved. Experimental results are presented to demonstrate the effectiveness of the proposed method.

97 MATHEMATICS AND COMPUTING↗