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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Assessment of Static Delamination Propagation Capabilities in Commercial Finite Element Codes Using Benchmark Analysis

With capabilities for simulating delamination growth in composite materials becoming available, the need for benchmarking and assessing these capabilities is critical. In this study, benchmark analyses were performed to assess the delamination propagation simulation capabilities of the VCCT implementations in Marc TM and MD NastranTM. Benchmark delamination growth results for Double Cantilever Beam, Single Leg Bending and End Notched Flexure specimens were generated using a numerical approach. This numerical approach was developed previously, and involves comparing results from a series of analyses at different delamination lengths to a single analysis with automatic crack propagation. Specimens were analyzed with three-dimensional and two-dimensional models, and compared with previous analyses using Abaqus . The results demonstrated that the VCCT implementation in Marc TM and MD Nastran(TradeMark) was capable of accurately replicating the benchmark delamination growth results and that the use of the numerical benchmarks offers advantages over benchmarking using experimental and analytical results.

Orifici, Adrian C.↗
NASA NTRS2016

Tangle-Free Finite Element Mesh Motion for Ablation Problems

In numerical simulations involving boundaries that evolve in time, the primary challenge is updating the computational mesh to reflect the physical changes in the domain. In particular, the fundamental objective for any such \mesh motion" scheme is to maintain mesh quality and suppress unphysical geometric anamolies and artifacts. External to a physical process of interest, mesh motion is an added component that determines the specifics of how to move the mesh given certain limited information from the main system. This paper develops a set of boundary conditions designed to eliminate tangling and internal collision within the context of PDE-based mesh motion (linear elasticity). These boundary conditions are developed for two- and three-dimensional meshes. The paper presents detailed algorithms for commonly occuring topological scenarios and explains how to apply them appropriately. Notably, the techniques discussed herein make use of none of the specifics of any particular formulation of mesh motion and thus are more broadly applicable. The two-dimensional algorithms are validated by an extensive verification procedure. Finally, many examples of diverse geometries in both two- and three-dimensions are shown to showcase the capabilities of the tangle-free boundary conditions.

Droba, Justin↗
DOE OSTI2018

High Performance GPU-Based FEM Code for Simulating Welding and Additive Manufacturing

An in-house HPC finite element code with novel computational schemes provides fast and high-fidelity simulation capability for welding processes and allied manufacturing processes such as metal additive manufacturing. Benchmark comparison confirmed that, for the testing cases studied, the newly developed in-house code resulted in increases in computational speed in the ranges of 100x to 1000x depending on the size and complexity of the testing cases, when compared to today's commercial welding process simulation tools. For a reference problem in nuclear reactor repair welding, ORNL's in-house high fidelity code was able to complete the simulation in 3.2 days, which would take estimated 10 to 20 years using today's commercial code. Similar order of magnitudes of reduction in computational times for metal additive manufacturing process were also achieved.

Feng, Zhili↗
DOE OSTI2019

Effect of α/β Forging on Microstructure and Texture Inhomogeneity in a Ti-1023 Forged Disk

The microstructure and texture distribution of a Ti-1023 forged disk were investigated by scanning electron microscopy and synchrotron-based high-energy X-ray diffraction. The finite element method was used to simulate temperature and strain distribution in order to investigate the relationship of the α/β forging process with microstructure and texture distribution. A bimodal microstructure and rolling textures with large inhomogeneity were observed in the disk. A plate-like and a necklace-like morphology and volume fraction variations of the primary α phase were observed in different regions with different forging conditions, such as temperature, deformation, position, and high density of flow lines. Texture sharpness distribution of the β phase was in good agreement with the strain distribution, which suggests deformation may play the most important part in the texture inhomogeneity. A weak cube texture was obtained near the center of the disk, and fiber textures were found near the rim of the disk. The primary a phase also exhibited a transverse texture, which is favored by large deformation.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Crash Testing and Simulation of a Cessna 172 Aircraft: Pitch Down Impact Onto Soft Soil

During the summer of 2015, NASA Langley Research Center conducted three full-scale crash tests of Cessna 172 (C-172) aircraft at the NASA Langley Landing and Impact Research (LandIR) Facility. The first test represented a flare-to-stall emergency or hard landing onto a rigid surface. The second test, which is the focus of this paper, represented a controlled-flight-into-terrain (CFIT) with a nose-down pitch attitude of the aircraft, which impacted onto soft soil. The third test, also conducted onto soil, represented a CFIT with a nose-up pitch attitude of the aircraft, which resulted in a tail strike condition. These three crash tests were performed for the purpose of evaluating the performance of Emergency Locator Transmitters (ELTs) and to generate impact test data for model validation. LS-DYNA finite element models were generated to simulate the three test conditions. This paper describes the model development and presents test-analysis comparisons of acceleration and velocity time-histories, as well as a comparison of the time sequence of events for Test 2 onto soft soil.

Fasanella, Edwin L.↗
DOE OSTI2016

Combined near- and far-field high-energy diffraction microscopy dataset for Ti-7Al tensile specimen elastically loaded in situ

High-energy diffraction microscopy (HEDM) constitutes a suite of combined X-ray characterization methods, which hold the unique advantage of illuminating the microstructure and micromechanical state of a material during concurrent in situ mechanical deformation. The data generated from HEDM experiments provides a heretofore unrealized opportunity to validate meso-scale modeling techniques, such as crystal plasticity finite element modeling (CPFEM), by explicitly testing the accuracy of these models at the length scales where the models predict their response. Combining HEDM methods with in situ loading under known and controlled boundary conditions represents a significant challenge, inspiring the recent development of a new high-precision rotation and axial motion system for simultaneously rotating and axially loading a sample. In this paper, we describe the initial HEDM dataset collected using this hardware on an alpha-titanium alloy (Ti-7Al) under in situ tensile deformation at the Advanced Photon Source, Argonne National Laboratory. We present both near-field HEDM data that maps out the grain morphology and intragranular crystallographic orientations and far-field HEDM data that provides the grain centroid, grain average crystallographic orientation, and grain average elastic strain tensor for each grain. Finally, we provide a finite element mesh that can be utilized to simulate deformation in the volume of this Ti-7Al specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tusas: A fully implicit parallel approach for coupled phase-field equations

In this study, we develop a fully-coupled, fully-implicit approach for phase-field modeling of solidification in metals and alloys. Predictive simulation of solidification in pure metals and metal alloys remains a significant challenge in the field of materials science, as microstructure formation during the solidification process plays a critical role in the properties and performance of the solid material. Our simulation approach consists of a finite element spatial discretization of the fully-coupled nonlinear system of partial differential equations at the microscale, which is treated implicitly in time with a preconditioned Jacobian-free Newton-Krylov method. The approach is algorithmically scalable as well as efficient due to an effective preconditioning strategy based on algebraic multigrid and block factorization. We implement this approach in the open-source Tusas framework, which is a general, flexible tool developed in C++ for solving coupled systems of nonlinear partial differential equations. The performance of our approach is analyzed in terms of algorithmic scalability and efficiency, while the computational performance of Tusas is presented in terms of parallel scalability and efficiency on emerging heterogeneous architectures. We demonstrate that modern algorithms, discretizations, and computational science, and heterogeneous hardware provide a robust route for predictive phase-field simulation of microstructure evolution during additive manufacturing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Predicting the flow stress of Zircaloy-4 under in-reactor accident conditions - 2016-0075

During loss of coolant accident (LOCA) and reactivity-initiated accident (RIA), nuclear fuel rods experience high heating rates that change the microstructure and properties of zirconium cladding materials. The aim was to determine how different fast heating rates affect the kinetics of transformation and whether its consequences on the material yield stresses in the dual-phase region could be predicted. Phase fraction and texture of cold-rolled Zircaloy-4 were thus measured at fast heating rates similar to those experienced during LOCA and RIA by in-situ high-energy synchrotron x-ray diffraction; the yield stress was also measured at different heating rates to the same dual-phase temperature by electrical resistance method. Crystal plasticity finite element model was used to simulate the flow stresses with the inputs of the measured phase fraction and texture. The model is calibrated using the flow stress at the heating rate of 10 deg. Cs{sup -1} and used to predict the flow stress at the heating rate of 50 deg. Cs{sup -1} at the same temperature. A significant shift of the β transus temperature was identified for fast heating rates in the range of 10-100 deg. Cs{sup -1}. There is great similarity in texture at different heating rates. In the dual-phase temperature region, the a texture is almost identical to the starting texture and the β texture is nearly random. Yield stresses of cold-rolled Zircaloy-4 decrease at faster heating rate at the same temperature of 920 deg. C. A good agreement was found between the experimental and predicted flow stress at different heating rates, indicating this methodology has potential to be used to predict flow behavior in a fast transient regime. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2016

Predicting the flow stress of Zircaloy-4 under in-reactor accident conditions - 2016-0075

During loss-of-coolant accident (LOCA) and reactivity-initiated accident (RIA), nuclear fuel rods experience high heating rates that change the microstructure and properties of zirconium cladding materials. We sought to determine how different fast heating rates affect the kinetics of transformation and whether their consequences on the material yield stresses in the dual-phase region could be predicted. Phase fraction and texture of cold-rolled Zircaloy-4 were thus measured at fast heating rates similar to those experienced during LOCA and RIA by in situ high-energy synchrotron X-ray diffraction; the yield stress was also measured at different heating rates to the same dual-phase temperature by the electrical resistance method. A crystal plasticity finite-element model was used to simulate the flow stresses with the inputs of the measured phase fraction and texture. The model was calibrated using the flow stress at the heating rate of 10 deg. C s{sup -1} and used to predict the flow stress at the heating rate of 50 deg. C s{sup -1} at the same temperature. A significant shift of the β transus temperature was identified for fast heating rates in the range of 10 to 100 deg. C s{sup -1}. There was great similarity in texture at different heating rates. In the dual-phase temperature region, the α texture was almost identical to the starting texture, and the b texture was nearly random. Yield stresses of cold-rolled Zircaloy-4 decreased at a faster heating rate at the same temperature of 920 deg. C. A good agreement was found between the experimental and predicted flow stress at different heating rates, indicating this methodology can potentially be used to predict flow behavior in a fast transient regime. (authors)

36 MATERIALS SCIENCE↗
NASA NTRS1998

Numerical Modeling of Solidification in Space With MEPHISTO-4

A pre-flight analysis of the directional solidification of BiSn with MEPHISTO-4 is presented. Simplified Bridgman growth under microgravity conditions is simulated using a two dimensional finite element model. This numerical model is a single domain, pseudo-steady state model, and includes the effects of both thermal and solutal convection. The results show that for all orientations of the applied steady state gravity vector, of magnitude 1 micro-g, the directional solidification process remains diffusion controlled. The maximum convective velocity was found to be 4.424 x 10(exp -5) cm/s for the horizontal Bridgman growth configuration. This value is an order of magnitude lower than the growth velocity. The maximum and minimum values or solute concentration in the liquid at the crystal-melt interface were 13.867 at.% and 13.722 at.%, respectively. This gives a radial segregation value of xi = 1.046% at the interface. A secondary objective of this work was to compare the results obtained to those that consider thermal convection only (no solutal convection). It was found that the convective flow patterns in simulations which included solutal convection were significantly different from those which ignored solutal convection. The level of radial segregation predicted by the current simulations is an order of magnitude lower than that found in simulations which ignore solutal convection. The final aim was to investigate the effect of g-jitter on the crystal growth process. A simulation was performed to calculate the system response to a 1 second, 100 micro-g gravity impulse acting normal to the direction of growth. This pulse is consistent with that induced by Orbiter thruster firings. The results obtained indicate that such a gravity pulse causes an increase in the level of radial solute segregation at the interface from the steady state values. The maximum value of solute concentration in the liquid was found to be 13.888 at.%, the minimum value calculated was 13.706 at.%, yielding a radial segregation value of xi = 1.31% at the interface. These values occurred 126 seconds after the pulse terminated. Thus it is anticipated that the process will remain diffusion controlled even when subjected to such g-jitter.

Simpson, James E.↗
DOE OSTI2015

Effect of moisture content and distribution on radiation transport in concrete - 14426

Criticality Emergency Planning (CEP) requires a detailed knowledge of the Facility design and construction in order to undertake radiation transport studies in support of criticality consequence calculations. Such calculations are used at AWE to determine mitigated doses, CEP exclusion zones and in support of Criticality Warning Systems (CWS) detector head placement studies. At AWE, such calculations have traditionally been performed by assuming that structural concrete has a fixed and evenly distributed moisture profile. In order to validate this assumption, work has been undertaken using the finite element analysis tool Femmasse to simulate moisture movement within the concrete and calculate moisture distributions in a 65 cm thick concrete wall as a function of aging and specific environmental conditions. This work identified that considerable variation in the moisture profile may arise, with an equilibrium condition being reached after circa 20 years, where the overall moisture content is reduced to a minimum. Furthermore, the calculations identified that the overall moisture content used in extant CEP assessments was considerably lower than that estimated in this work. In order to determine the effect of the moisture distribution profiles, the Monte Carlo transport code MCNP5 was used to study the transmission of radiation through aged concrete and thus assess its performance as a radiation shield in the event of a criticality. Various moisture distributions were assessed based on the findings of the Femmasse calculations and were compared with that assessed using traditional assumptions. Results were obtained for neutron and photon doses. The paper provides an overview of the findings and makes recommendations on the suitability of using evenly distributed moisture profiles in concrete in support of mitigated doses, CEP and CWS detector head placement studies. (authors)

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

A Workflow to Rapidly Interrogate Multiscale Model Simulation Results Across Multiple Length Scales

Many tools can be used to visualize field and state variables for a single scale analysis so that the influence of relevant mechanisms can be evaluated. Finite element software is often utilized to simulate a unit cell of a material and visualize results at that scale. Material properties can be homogenized from individual constituents and local deformation, damage, and failure mechanisms can be evaluated within the unit cell due to globally applied boundary conditions. Such solutions can produce satisfactory results if a user is only interested in analyzing a single scale. But materials in general contain features across multiple length scales, and assumptions must be made when attempting to account for lower length scale phenomena within a higher length scale model. Multiscale modeling is an attractive means to model materials because detailed material responses can be tracked across multiple disparate length scales while reducing the amount of required assumptions. However, as the complexity of these models increases, a large amount of data can be produced, and data traceability can become increasingly more difficult. Field and state variables, which are naturally dependent on spatial position, may themselves be calculated from one or more lower length scale unit cell models each with their own appropriate field and state variables. The NASA Multiscale Analysis Tool (NASMAT) is one software that can be used to perform a multiscale analysis efficiently and output requested data at all length scales in the analysis. A companion open-source Python software, NASMAT PrePost, can be used to visualize NASMAT model results and rapidly interrogate multiscale data across multiple length scales. This presentation will demonstrate some of the key features of NASMAT PrePost on two multiscale problems by quickly displaying and demonstrating connectivity among multiscale results from large datasets.

Python↗
NASA NTRSDate not supplied

Study of Aircraft Structural Response and Occupant Loading During a Ditching Event Utilizing LS-DYNA Simulation

In this study, the capability of finite element model (FEM) analysis to assess aircraft structural response and occupant injury risk during aircraft water ditching conditions was evaluated. A FEM of the Fokker F28 aircraft, previously validated against full-scale crash test data, was simulated in a realistic water impact environment utilizing recorded conditions from the ditching of US Airways Flight 1549 into the Hudson River. The FEM of the F28 aircraft, was developed from a NASTRAN loads model provided by the vehicle manufacturer, calibrated using component fuselage drop-tests, and validated through test-analysis comparisons from a full-scale crash test of the F28 onto soil. Simulation results are presented in several categories: vehicle kinematics, structural accelerations, damage, occupant responses, and injury risk predictions. Results were consistent with the general outcomes of Flight 1549, indicating that the developed vehicle FEM, validated against land impact test data, may provide an effective tool for evaluating water ditching events.

Jacob Putnam↗
NASA NTRS1989

A NASTRAN/TREETOPS solution to a flexible, multi-body dynamics and controls problem on a UNIX workstation

Demands for nonlinear time history simulations of large, flexible multibody dynamic systems has created a need for efficient interfaces between finite-element modeling programs and time-history simulations. One such interface, TREEFLX, an interface between NASTRAN and TREETOPS, a nonlinear dynamics and controls time history simulation for multibody structures, is presented and demonstrated via example using the proposed Space Station Mobile Remote Manipulator System (MRMS). The ability to run all three programs (NASTRAN, TREEFLX and TREETOPS), in addition to other programs used for controller design and model reduction (such as DMATLAB and TREESEL, both described), under a UNIX Workstation environment demonstrates the flexibility engineers now have in designing, developing and testing control systems for dynamically complex systems.

Benavente, Javier E.↗
DOE OSTI2015

Improving the Accuracy of PCMI Simulations with More Realistic Geometry and Material Models

Assessment of fuel performance codes is an endeavour that requires re-evaluation as new simulation capabilities become available. BISON, the finite element-based fuel performance code in development at Idaho National Laboratory, has been used to simulate a variety of experiments. Calculations from these simulations are compared to experimental measurements for the purpose of evaluating the progress of code development. One such experiment, from the Risø 3 Fission Gas Project, provides measurements of fission gas release and clad diameter. A subset of the Risø3 Fission Gas Project, which is referred to here, consisted of a base-irradiation period of a full-length fuel rod in the Quad-Cities-1 BWR. The base irradiation test was followed by a “bump test”, in which a sub-section of the original rod was subjected to a strongly non-uniform axial power profile. The bump test took place in Risø test reactor DR3 in a water-cooled HP1 rig under BWR conditions. This experiment was chosen because it allows for an evaluation of several aspects of BISON, including fully coupled thermo-mechanics, contact, and several nonlinear material models. A summary of simulation results and comparison to the experiment data has been published previously, however, this paper shows new results with recently added features to BISON. These new features include representing the pellets as discrete entities rather than a smeared column, including primary creep in the clad creep model, and the coupling of fuel gaseous swelling to fission gas release. Comparing to previous results, the addition of theses features has narrowed the gap between BISON calculations and experimental measurements for final clad diameter and fission gas release.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Coupled nonlinear elasticity, plastic slip, twinning, and phase transformation in single crystal titanium for plate impact loading

High purity single crystal titanium (Ti) under shockwave loading in plate impact experiment is modelled and simulated, with help of the finite element. A thermodynamically consistent system of equations is formulated in the frame of large deformation to consider material anisotropy, rate dependence, and multi-physics including nonlinear elasticity, dislocation based plastic slip, deformation twinning, and phase transformation. A novel kinematics is proposed to consider phase transformation in the α twin variants and also to consider the dislocation based slip in all of components of parent material, primary twins and the ω phase. The thermodynamically consistent driving forces for plastic slip, twinning, and phase transformation are developed based upon the second law of thermodynamics. For nonlinear elasticity, the stiffness matrix is dependent on the volume fractions of all components, and the volumetric part of Cauchy stress is obtained from the nonlinear equation of state. Dislocation based plastic slip in multi-variant and multiphase heterogeneous materials is developed and interactions of dislocations among all slip modes are taken into account. A mechanism for dislocation density evolution during twinning and phase transformation is proposed. The shock loading along the [0001] and [10$\bar{1}1$ directions of single crystal high purity Ti is investigated computationally. Very good correspondence between simulation and experiment is obtained, which includes pole figures, volume fractions of components, free surface velocity, peak pressure, and phase transformation pressure. Multiple experimental phenomena are interpreted based upon the progression of dislocation slip, deformation twinning, and phase transformation. In compression with the [10$\bar{1}1$ crystal, a higher volume fraction in the primary twins but a lower secondary twin volume fraction in the [0001] crystal in experiment was observed. Finally, the higher propensity for phase transformation occurs in the [0001] crystal is reproduced. In addition to material texture, distributions of temperature, stresses, and plastic strains dependent on the impact loading directions are revealed.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Frequency dependence of the acoustic field generated from a spherical cavity transducer with open ends

Resolution of high intensity focused ultrasound (HIFU) focusing is limited by the wave diffraction. We have developed a spherical cavity transducer with two open ends to improve the focusing precision without sacrificing the acoustic intensity (App Phys Lett 2013; 102: 204102). This work aims to theoretically and experimentally investigate the frequency dependence of the acoustic field generated from the spherical cavity transducer with two open ends. The device emits high intensity ultrasound at the frequency ranging from 420 to 470 kHz, and the acoustic field is measured by a fiber optic probe hydrophone. The measured results shows that the spherical cavity transducer provides high acoustic intensity for HIFU treatment only in its resonant modes, and a series of resonant frequencies can be choosen. Furthermore, a finite element model is developed to discuss the frequency dependence of the acoustic field. The numerical simulations coincide well with the measured results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Radiation and viscous dissipation effect on square porous annulus

The present study is carried out to investigate the effect of radiation and viscous dissipation in a square porous annulus subjected to outside hot T{sub h} and inside cold T{sub c} temperature. The square annulus has a hollow section of dimension D×D at the interior of annulus. The flow is assumed to obey Darcy law. The governing equations are non-dimensionalised and solved with the help of finite element method. Results are discussed with respect to viscous dissipation parameter, radiation parameter and size of the hollow section of annulus.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗