Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element simulation”

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

Two-tier model reduction of viscoelastically damped finite element models

Transient simulations of linear viscoelastically damped structures necessitate critical computational resources to directly integrate the full-order finite element model with time-stepping algorithms. Traditional modal reduction techniques are not directly applicable to these systems since viscoelastic materials depend on time and frequency. A more appropriate reduction basis is obtained from the nonlinear, complex eigenvalue problem, whose eigenvectors capture the appropriate kinematics and enable frequency-based mode selection; unfortunately, the computational cost is prohibitive for computing these modes from large-scale engineering models. To address this shortcoming, this work proposes a novel two-tier reduction procedure to reduce the upfront cost of solving the complex, nonlinear eigenvalue problem. The first reduction step reduces the full-order model with real mode shapes linearized about various centering frequencies to capture the kinematics over a full range of viscoelastic material behavior (glassy, rubbery, and glass-transition zones). This tier-one reduction preserves time-temperature superposition and allows the equations to depend parametrically on operating temperature. The second-level reduction then solves the complex, nonlinear eigenmode solutions in the tier-one reduced space about a fixed temperature to further reduce the equations-of-motion. The method is demonstrated on a cantilevered sandwich plate to display its accuracy and efficiency in comparison to full-order model predictions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Mechanical Properties of Niobium Cavities

The mechanical stability of bulk Nb cavity is an important aspect to be considered in relation to cavity material, geometry and treatments. Mechanical properties of Nb are typically obtained from uniaxial tensile tests of small samples. In this contribution we report the results of measurements of the resonant frequency and local strain along the contour of single-cell cavities made of ingot and fine-grain Nb of different purity subjected to increasing uniform differential pressure, up to 6 atm. Measurements have been done on cavities subjected to different heat treatments. Good agreement between finite element analysis simulations and experimental data in the elastic regime was obtained with a single set of values of Young’s modulus and Poisson’s ratio. The experimental results indicate that the yield strength of medium-purity ingot Nb cavities is higher than that of fine-grain, high-purity Nb.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Multidisciplinary Design Of Hot Composite Structures

Unified computer code developed to implement multidisciplinary approach to design and analysis of composite-material structures that must withstand high temperatures. Code modular: includes executive module communicating with and coordinating other modules performing calculations pertaining to traditionally separate disciplines like those of acoustics, structural vibrations, structural loads, and thermal effects. Essential feature, finite-element numerical simulation of relevant physical phenomena according to applicable disciplines. Same finite-element mesh used in thermal, vibrational, and structural analyses; minimizing data-preparation time and eliminating errors incurred in transforming temperatures from one finite-element mesh to another.

Chamis, Christos C.↗
DOE OSTI2016

Interfacing VPSC with finite element codes. Demonstration of irradiation growth simulation in a cladding tube

This Milestone report shows good progress in interfacing VPSC with the FE codes ABAQUS and MOOSE, to perform component-level simulations of irradiation-induced deformation in Zirconium alloys. In this preliminary application, we have performed an irradiation growth simulation in the quarter geometry of a cladding tube. We have benchmarked VPSC-ABAQUS and VPSC-MOOSE predictions with VPSC-SA predictions to verify the accuracy of the VPSCFE interface. Predictions from the FE simulations are in general agreement with VPSC-SA simulations and also with experimental trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Optimization Using Metamodeling in the Context of Integrated Computational Materials Engineering (ICME)

Predictive Design Technologies, LLC (PDT) proposed to employ Integrated Computational Materials Engineering (ICME) tools to help the manufacturing industry in the United States regain the competitive advantage in the global economy. ICME uses computational materials science tools within a holistic system in order to accelerate materials development, improve design optimization, and unify design and manufacturing. With the advent of accurate modeling and simulation along with significant increases in high performance computing (HPC) power, virtual design and manufacturing using ICME tools provide the means to reduce product development time and cost by alleviating costly trial-and-error physical design iterations while improving overall quality and manufacturing efficiency. To reduce the computational cost necessary for the large-scale HPC simulations and to make the methodology accessible for small and medium-sized manufacturers (SMMs), metamodels are employed. Metamodels are approximate models (functional relationships between input and output variables) that can reduce the simulation times by one to two orders of magnitude. In Phase I, PDT, partnered with Mississippi State University (MSU), demonstrated the feasibility of the proposed methodology by employing MSU?s internal state variable (ISV) plasticity-damage model with the help of metamodels to optimize the microstructure-process-property-cost for tube manufacturing processes used by Plymouth Tube Company (PTC), which involves complicated temperature and mechanical loading histories. PDT quantified the microstructure-property relationships for PTC?s SAE J525 electric resistance-welded cold drawn low carbon hydraulic 1010 steel tube manufacturing processes at seven different material states and calibrated the ISV plasticity material parameters to fit experimental tensile stress-strain curves. PDT successfully performed large scale finite element (FE) simulations in an HPC environment using the ISV plasticity model in Abaqus FE analyses of the tube forming, sizing, drawing, welding, and normalizing processes. The simulation results coupled with the manufacturing cost data were used to develop prototype metamodeling (quick response) codes which could be used to predict and optimize the microstructure-process-property-cost relationships. The developed ICME metamodeling toolkits are flexible enough to be applied to other manufacturing processes (e.g. forging, forming, casting, extrusion, rolling, stamping, and welding/joining) and metamodeling codes can run on laptop computers. Based on the work completed in Phase I, in Phase II, PDT proposes to continue to refine the ISV model by correlating and incorporating the uncertainties in the microstructure, mechanical testing, and modeling. Following the model refinement, FE analyses will be simulated and will provide even more realistic predictions as they include an appropriate window of uncertainty. Using the HPC output (FE analyses) as input, the quick-response metamodel codes will more accurately predict and optimize the microstructure-process-property-cost relationships. Furthermore, PDT propose to employ the ICME metamodeling toolkits to help develop a new tube product using entirely new high strength steel. The modeling of the high strength steel manufacturing process will replace the costly and time consuming trial-and-error methods that were used in the tubing industry previously. This simulation-based process prototyping will greatly benefit our industrial partners by opening up new market spaces due to new products with greater capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ExaAM: Metal additive manufacturing simulation at the fidelity of the microstructure

Additive manufacturing (AM), or 3D printing, of metals is transforming the fabrication of components, in part by dramatically expanding the design space, allowing optimization of shape and topology. However, although the physical processes involved in AM are similar to those of welding, a field with decades of experimental, modeling, simulation, and characterization experience, qualification of AM parts remains a challenge. The availability of exascale computational systems, particularly when combined with data-driven approaches such as machine learning, enables topology and shape optimization as well as accelerated qualification by providing process-aware, locally accurate microstructure and mechanical property models. We describe the physics components comprising the Exascale Additive Manufacturing simulation environment and report progress using highly resolved melt pool simulations to inform part-scale finite element thermomechanics simulations, drive microstructure evolution, and determine constitutive mechanical property relationships based on those microstructures using polycrystal plasticity. We report on implementation of these components for exascale computing architectures, as well as the multi-stage simulation workflow that provides a unique high-fidelity model of process–structure–property relationships for AM parts. In addition, we discuss verification and validation through collaboration with efforts such as AM-Bench, a set of benchmark test problems under development by a team led by the National Institute of Standards and Technology.

3D printing↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Finite element approximation of the radiative transport equation in a medium with piece-wise constant refractive index

The radiative transport equation can be used as a light transport model in a medium with scattering particles, such as biological tissues. In the radiative transport equation, the refractive index is assumed to be constant within the medium. However, in biomedical media, changes in the refractive index can occur between different tissue types. In this work, light propagation in a medium with piece-wise constant refractive index is considered. Light propagation in each sub-domain with a constant refractive index is modeled using the radiative transport equation and the equations are coupled using boundary conditions describing Fresnel reflection and refraction phenomena on the interfaces between the sub-domains. The resulting coupled system of radiative transport equations is numerically solved using a finite element method. The approach is tested with simulations. The results show that this coupled system describes light propagation accurately through comparison with the Monte Carlo method. It is also shown that neglecting the internal changes of the refractive index can lead to erroneous boundary measurements of scattered light.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Multibody instantly centered moving Lagrangian observer schemes. I - Formulation. II - Application to vehicular simulations

A finite element modeling method employing an instantly centered moving Lagrangian observer is proposed to solve multibody problems involving several rotating components, each with its own rotational history. The present technique allows the steady-state behavior of the problem simulation to be rendered stationary and time independent. In the second part, the method is extended to model the steady and transient response of ground-based automotive-type vehicular systems, including the modeling of vehicular obstruction rollover events.

Padovan, Joseph↗
NASA NTRS2005

Finite Element Analysis of Adaptive-Stiffening and Shape-Control SMA Hybrid Composites

Shape memory alloy hybrid composites with adaptive-stiffening or morphing functions are simulated using finite element analysis. The composite structure is a laminated fiber-polymer composite beam with embedded SMA ribbons at various positions with respect to the neutral axis of the beam. Adaptive stiffening or morphing is activated via selective resistance heating of the SMA ribbons or uniform thermal loads on the beam. The thermomechanical behavior of these composites was simulated in ABAQUS using user-defined SMA elements. The examples demonstrate the usefulness of the methods for the design and simulation of SMA hybrid composites. Keywords: shape memory alloys, Nitinol, ABAQUS, finite element analysis, post-buckling control, shape control, deflection control, adaptive stiffening, morphing, constitutive modeling, user element

Gao, Xiujie↗
NASA NTRS1999

Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms

This paper presents a final report on Wavelet and Multiresolution Analysis for Finite Element Networking Paradigms. The focus of this research is to derive and implement: 1) Wavelet based methodologies for the compression, transmission, decoding, and visualization of three dimensional finite element geometry and simulation data in a network environment; 2) methodologies for interactive algorithm monitoring and tracking in computational mechanics; and 3) Methodologies for interactive algorithm steering for the acceleration of large scale finite element simulations. Also included in this report are appendices describing the derivation of wavelet based Particle Image Velocity algorithms and reduced order input-output models for nonlinear systems by utilizing wavelet approximations.

Kurdila, Andrew J.↗
NASA NTRSDate not supplied

Modeling and Simulation Credibility Assessment of Finite Element Models for Fall From Heights Injury Scenarios

Finite element (FE) computational models of the whole human body are being used to investigate the injury modalities and mechanisms arising from suited astronauts’ ExtraVehicular Activities (EVAs). These FE whole body models – simplified Elemance and THUMS pedestrian versions – have undergone a series of verification and validation (V&V) assessments for applications in automotive, domestic, and contact sports arenas. However, the FE models’ Modeling and Simulation (M&S) credibility assessments for simulating suited astronaut injuries due to fall from heights in altered gravity environments are limited. The current study evaluates the M&S credibility of simplified Elemance and THUMS pedestrian FE models based on eight credibility factors – data pedigree, input pedigree, code verification, solution verification, conceptual validation, referent validation, results uncertainty, and results robustness – as per NASA-STD-7009A [1]. These credibility factors were assigned ordinal scores ranging between 0 to 4, with a score of 0 in the case of insufficient evidence and 4 when all data for the model is known, used appropriately, and traceable to the Real-World System (RWS) with acceptable accuracy, precision, and uncertainty. The credibility factor sufficiency thresholds were assigned upon guidance from NASA EVA subject matter experts. The factor scores were ascertained through the evaluation of M&S credibility assessment procedures described in the published literature and FE model reports. The M&S credibility scores for the Elemance and THUMS FE models vary between 1 and 3 for data pedigree, input pedigree, and referent validation credibility factors, and between 0 and 1 for code verification, solution verification, conceptual validation, results uncertainty and results robustness. The results show that these FE models are credible for simulating fall from heights injury scenarios in specific kinematic and kinetic ranges that correlate with highway and contact sports-related external conditions. Multiple credibility elevation strategies are prescribed to improve the Elemance and THUMS models’ credibility for use in the EVA fall from height injury scenarios.

R A Perkins↗
DOE OSTI2022

Uncertainty quantification in elastic constants of SiC f /SiC m tubular composites using global sensitivity analysis

Silicon carbide fiber and silicon carbide matrix (SiC f /SiC m ) tubes produced through the chemical vapor infiltration process have become a candidate cladding material in nuclear applications. The performance of this composite is influenced by many variables such as braiding angle, porosity, material properties, etc., which vary over a range of values due to the inherent fluctuations in the manufacturing process. In this study, the variability in elastic constants of SiC f /SiC m composite has been quantified through multiscale finite element (FE) simulations, variable screening, and high-fidelity surrogate modeling. The key variables dominantly affecting the elastic constants of SiC f /SiC m tubes were identified using global sensitivity analysis. A surrogate to the high-fidelity FE-based model was used in Monte Carlo simulations to generate a hundred thousand samples from which the uncertainty in elastic constants was assessed. It turned out that the coefficient of variation was less than 10%.

Materials Science↗
DOE OSTI2020

Tailored Graphene Micropatterns by Wafer-Scale Direct Transfer for Flexible Chemical Sensor Platform

2D materials, such as graphene, exhibit great potential as functional materials for numerous novel applications due to their excellent properties. The grafting of conventional micropatterning techniques on new types of electronic devices is required to fully utilize the unique nature of graphene. However, the conventional lithography and polymer-supported transfer methods often induce the contamination and damage of the graphene surface due to polymer residues and harsh wet-transfer conditions. In this work, a novel strategy to obtain micropatterned graphene on polymer substrates using a direct curing process is demonstrated. Employing this method, entirely flexible, transparent, well-defined self-activated graphene sensor arrays, capable of gas discrimination without external heating, are fabricated on 4 in. wafer-scale substrates. Finite element method simulations show the potential of this patterning technique to maximize the performance of the sensor devices when the active channels of the 2D material are suspended and nanoscaled. This study contributes considerably to the development of flexible functional electronic devices based on 2D materials.

2D materials, Graphene, Microscale patterning, Che↗
DOE OSTI2016

The cross-aisle seismic performance of storage rack base connections

Steel storage racks used in retail stores and warehouses are seismically designed as moment resisting frames in the down-aisle direction, and braced frames in the cross-aisle direction. While their down-aisle response is relatively well understood, there is little understanding of their cross-aisle response, especially as it pertains to the desired mode of inelastic deformation and associated design methods. In this paper, results are presented from six full scale tests on braced frames representing storage racks in the cross-aisle direction. These tests investigate the base plate thickness and dimensions, and the upright (column) cross section. The experiments indicate that inelastic deformation in the base plate provides stable hysteretic response with significant ductility and energy dissipation. Ductile tearing is also observed in welds connecting the base plate to the upright. However, it does not appear to negatively influence the hysteretic response. The tests are complemented by Finite Element (FE) simulations of the base connections. These simulations provide insights into internal force distributions within the connections. Based on these insights, analytical equations are proposed for characterizing the backbone curve of the hysteretic response, for use in displacement based design methods. It is determined that the current approach for characterizing design forces in the anchors is unconservative, since it does not incorporate the effects of strain hardening or the membrane action as the base plate undergoes large deformations. A new approach which incorporates these phenomena is presented, and determined to be significantly more accurate. Finally, limitations of the study are outlined and directions for future work are identified.

36 MATERIALS SCIENCE↗
DOE OSTI2015

The Young's modulus of high-aspect-ratio carbon/carbon nanotube composite microcantilevers by experimental and modeling validation

This paper reports the Young's modulus of a carbon nanotube (CNT)-reinforced carbon/CNT (C/CNT) composite microcantilevers measured by laser Doppler vibrometer and validated by finite element method. Also, the microfabrication process of the high-aspect-ratio C/CNT microcantilever arrays based on silicon micromolding and pyrolysis is presented in detail. With the in-plane natural resonant frequencies of the microcantilevers measured by a laser Doppler vibrometer, a single degree of freedom (SDoF) model based on Euler-Bernoulli (E-B) beam theory is used to calculate the Young's modulus of this composite. To figure out whether this SDoF model can be applied to these composite microcantilevers, the finite element (FE) simulation of these microcantilevers was performed. The Young's modulus of C/CNT composite microcantilevers fabricated by the pyrolysis process at 600 °C is 9391 MPa, and a good agreement between the results from experiments and FE simulation is obtained.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

Implementation of the full viscoresistive magnetohydrodynamic equations in a nonlinear finite element code

Numerical simulations form an indispensable tool to understand the behavior of a hot plasma that is created inside a tokamak for providing nuclear fusion energy. Various aspects of tokamak plasmas have been successfully studied through the reduced magnetohydrodynamic (MHD) model. The need for more complete modeling through the full MHD equations is addressed here. Our computational method is presented along with measures against possible problems regarding pollution, stability, and regularity. The problem of ensuring continuity of solutions in the center of a polar grid is addressed in the context of a finite element discretization of the full MHD equations. A rigorous and generally applicable solution is proposed here. Useful analytical test cases are devised to verify the correct implementation of the momentum and induction equation, the hyperdiffusive terms, and the accuracy with which highly anisotropic diffusion can be simulated. A striking observation is that highly anisotropic diffusion can be treated with the same order of accuracy as isotropic diffusion, even on non-aligned grids, as long as these grids are generated with sufficient care. This property is shown to be associated with our use of a magnetic vector potential to describe the magnetic field. Several well-known instabilities are simulated to demonstrate the capabilities of the new method. The linear growth rate of an internal kink mode and a tearing mode are benchmarked against the results of a linear MHD code. The evolution of a tearing mode and the resulting magnetic islands are simulated well into the nonlinear regime. The results are compared with predictions from the reduced MHD model. Finally, a simulation of a ballooning mode illustrates the possibility to use our method as an ideal MHD method without the need to add any physical dissipation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗