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 217 records · Page 12

Automated design of an additive manufactured compact broadband antenna for plasma reflectometry

Broadband antennas operating in the gigahertz frequency range are regularly used for plasma reflectometry diagnostics. Due to a lack of space and unique diagnostic constraints, these antennas are often custom in design and frequency range. Recent advances in additive manufacturing of high temperature copper alloys allow for expanded freedom in design of these diagnostic antennas. In this work, a heuristic simulated annealing algorithm is used alongside 3-D finite element simulation to automate the design of a double ridged rectangular horn antenna for a reflectometry diagnostic on the DIII-D tokamak. Optimization of antenna performance given the design constraints results in a compact broadband (6-20+ GHz) antenna design. Measured transmission from the additively manufactured antenna matches simulation within reasonable error, and experimental plasma electron density profiles from the DIII-D high-field side scrape-off layer are shown.

Additive manufacturing↗

Phase Shift Control of a Three-Phase Inverter for Balanced Secondary Currents in Misaligned Three-Phase Inductive Power Transfer Systems

This paper presents a strategy for generating fundamental reference voltage commands for a double-sided LCC tuned three-phase wireless power transfer (WPT) system to achieve balanced output currents when the receiver and transmitter are not perfectly aligned. Formulas for translating the normalized fundamental current commands into phase-shift commands for two of the phase legs in a two-level three-phase inverter are derived. The methodology is validated through circuit simulations performed using coupler models derived from finite-element simulations. Simulation results for a 1.5 k W system design are included, which indicate near perfect correction of imbalance in the secondary currents for a typical misalignment case.

Pries, Jason↗

Analytical Simulations of Energy-Absorbing Impact Spheres for a Mars Sample Return Earth Entry Vehicle

Nonlinear dynamic finite element simulations were performed to aid in the design of an energy-absorbing impact sphere for a passive Earth Entry Vehicle (EEV) that is a possible architecture for the Mars Sample Return (MSR) mission. The MSR EEV concept uses an entry capsule and energy-absorbing impact sphere designed to contain and limit the acceleration of collected samples during Earth impact without a parachute. The spherical shaped impact sphere is composed of solid hexagonal and pentagonal foam-filled cells with hybrid composite, graphite-epoxy/Kevlar cell walls. Collected Martian samples will fit inside a smaller spherical sample container at the center of the EEV's cellular structure. Comparisons were made of analytical results obtained using MSC.Dytran with test results obtained from impact tests performed at NASA Langley Research Center for impact velocities from 30 to 40 m/s. Acceleration, velocity, and deformation results compared well with the test results. The correlated finite element model was then used for simulations of various off-nominal impact scenarios. Off-nominal simulations at an impact velocity of 40 m/s included a rotated cellular structure impact onto a flat surface, a cellular structure impact onto an angled surface, and a cellular structure impact onto the corner of a step.

Billings, Marcus Dwight↗

Crash Simulation of a Boeing 737 Fuselage Section Vertical Drop Test

A 30-ft/s vertical drop test of a fuselage section of a Boeing 737 aircraft was conducted in October of 1999 at the FAA Technical Center in Atlantic City, NJ. This test was performed to evaluate the structural integrity of a conformable auxiliary fuel tank mounted beneath the floor and to determine its effect on the impact response of the airframe structure and the occupants. The test data were used to compare with a finite element simulation of the fuselage structure and to gain a better understanding of the impact physics through analytical/experimental correlation. To perform this simulation, a full-scale 3-dimensional finite element model of the fuselage section was developed using the explicit, nonlinear transient-dynamic finite element code, MSC.Dytran. The emphasis of the simulation was to predict the structural deformation and floor-level acceleration responses obtained from the drop test of the B737 fuselage section with the auxiliary fuel tank.

Fasanella, Edwin L.↗

Impact Test and Simulation of Energy Absorbing Concepts for Earth Entry Vehicles

Nonlinear dynamic finite element simulations have been performed to aid in the design of an energy absorbing concept for a highly reliable passive Earth Entry Vehicle (EEV) that will directly impact the Earth without a parachute. EEV's are designed to return materials from asteroids, comets, or planets for laboratory analysis on Earth. The EEV concept uses an energy absorbing cellular structure designed to contain and limit the acceleration of space exploration samples during Earth impact. The spherical shaped cellular structure is composed of solid hexagonal and pentagonal foam-filled cells with hybrid graphite- epoxy/Kevlar cell walls. Space samples fit inside a smaller sphere at the center of the EEV's cellular structure. Comparisons of analytical predictions using MSC,Dytran with test results obtained from impact tests performed at NASA Langley Research Center were made for three impact velocities ranging from 32 to 40 m/s. Acceleration and deformation results compared well with the test results. These finite element models will be useful for parametric studies of off-nominal impact conditions.

Billings, Marcus D.↗

Crushing behaviors of unidirectional carbon fiber reinforced plastic composites under dynamic bending and axial crushing loading

The dynamic crushing characteristics of unidirectional carbon fiber reinforced plastic composites under two loading types, dynamic three-point bending and axial crushing, are investigated by experiment and finite element simulation. Hat section samples with two different layup orientations are tested at various impact velocities to investigate its effects. The experiment results show that delamination plays a critical role in dynamic bending deformation. Different layup orientations lead to similar crushing bending behaviors, but remarkable variations in delamination pattern evolutions. In contrast, layup orientation exerts more significant effects on failure modes and energy absorption in axial crushing. A finite element model with multiple layers of thick shell element and cohesive element is established, and simulations are performed for both dynamic three-point bending and axial crushing. The model can capture the crushing behaviors in the two loading types reasonably well. The damage evolution and energy absorption through various mechanisms at different loading conditions are discussed.

42 ENGINEERING↗

Verification and Validation of a Coordinate Transformation Method in Axisymmetric Transient Magnetics.

We present a verification and validation analysis of a coordinate-transformation-based numerical solution method for the two-dimensional axisymmetric magnetic diffusion equation, implemented in the finite-element simulation code ALEGRA. The transformation, suggested by Melissen and Simkin, yields an equation set perfectly suited for linear finite elements and for problems with large jumps in material conductivity near the axis. The verification analysis examines transient magnetic diffusion in a rod or wire in a very low conductivity background by first deriving an approximate analytic solution using perturbation theory. This approach for generating a reference solution is shown to be not fully satisfactory. A specialized approach for manufacturing an exact solution is then used to demonstrate second-order convergence under spatial refinement and tem- poral refinement. For this new implementation, a significant improvement relative to previously available formulations is observed. Benefits in accuracy for computed current density and Joule heating are also demonstrated. The validation analysis examines the circuit-driven explosion of a copper wire using resistive magnetohydrodynamics modeling, in comparison to experimental tests. The new implementation matches the accuracy of the existing formulation, with both formulations capturing the experimental burst time and action to within approximately 2%.

97 MATHEMATICS AND COMPUTING↗

Experimental evidence of high-frequency complete elastic bandgap in pillar-based phononic slabs

We present strong experimental evidence for the existence of a complete phononic bandgap, for Lamb waves, in the high frequency regime (i.e., 800 MHz) for a pillar-based phononic crystal (PnC) membrane with a triangular lattice of gold pillars on top. The membrane is composed of an aluminum nitride film stacked on thin molybdenum and silicon layers. Experimental characterization shows a large attenuation of at least 20 dB in the three major crystallographic directions of the PnC lattice in the frequency range of 760 MHz–820 MHz, which is in agreement with our finite element simulations of the PnC bandgap. The results of experiments are analyzed and the physics behind the attenuation in different spectral windows is explained methodically by assessing the type of Bloch modes and the in-plane symmetry of the displacement profile.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗

Deformational Features and Microstructure Evolution of Copper Fabricated by a Single Pass of the Elliptical Cross-Section Spiral Equal-Channel Extrusion (ECSEE) Process

The objectives of this work are to study the deformational feature, textures, microstructures, and dislocation configurations of ultrafine-grained copper processed by the process of elliptical cross-section spiral equal-channel extrusion (ECSEE). The deformation patterns of simple shear and pure shear in the ECSEE process were evaluated with the analytical method of geometric strain. The influence of the main technical parameters of ECSEE die on the effective strain distribution on the surface of ECSEE-fabricated samples was examined by the finite element simulation. The high friction factor could improve the effective strain accumulation of material deformation. Moreover, the pure copper sample fabricated by ECSEE ion shows a strong rotated cube shear texture. The refining mechanism of the dislocation deformation is dominant in copper processed by a single pass of ECSEE. The inhomogeneity of the micro-hardness distribution on the longitudinal section of the ECSEE-fabricated sample is consistent with the strain and microstructure distribution features.

36 MATERIALS SCIENCE↗

Processing Parameter Effects on Residual Stress and Mechanical Properties of Selective Laser Melted Ti6Al4V

Selective laser melting (SLM) process is characterized by large temperature gradients resulting in high levels of residual stress within the additively manufactured metallic structure. SLM-processed Ti6Al4V yields a martensitic microstructure due to the rapid solidification and results in a ductility generally lower than a hot working equivalent. Post-process heat treatments can be applied to SLM components to remove in-built residual stress and improve ductility. Residual stress buildup and the mechanical properties of SLM parts can be controlled by varying the SLM process parameters. This investigation studies the effect of layer thickness on residual stress and mechanical properties of SLM Ti6Al4V parts. This is the first-of-its kind study on the effect of varying power and exposure in conjunction with keeping the energy density constant on residual stress and mechanical properties of SLM Ti6Al4V components. It was found that decreasing power and increasing exposure for the same energy density lowered the residual stress and improved the % elongation of SLM Ti6Al4V parts. Increasing layer thickness resulted in lowering the residual stress at the detriment of mechanical properties. The study is based on detailed experimental analysis along with finite element simulation of the process using ABAQUS to understand the underlying physics of the process.

36 MATERIALS SCIENCE↗

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗

Developing parallel GeoFEST(P) using the PYRAMID AMR library

The PYRAMID parallel unstructured adaptive mesh refinement (AMR) library has been coupled with the GeoFEST geophysical finite element simulation tool to support parallel active tectonics simulations. Specifically, we have demonstrated modeling of coseismic and postseismic surface displacement due to a simulated Earthquake for the Landers system of interacting faults in Southern California. The new software demonstrated a 25-times resolution improvement and a 4-times reduction in time to solution over the sequential baseline milestone case. Simulations on workstations using a few tens of thousands of stress displacement finite elements can now be expanded to multiple millions of elements with greater than 98% scaled efficiency on various parallel platforms over many hundreds of processors. Our most recent work has demonstrated that we can dynamically adapt the computational grid as stress grows on a fault. In this paper, we will describe the major issues and challenges associated with coupling these two programs to create GeoFEST(P). Performance and visualization results will also be described.

adaptive mesh refinement (AMR)↗

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↗

High quality factor, fully switchable terahertz superconducting metasurface

We present a complementary THz metasurface realised with Niobium thin film which displays a quality factor Q = 54 and a fully switchable behaviour as a function of the temperature. The switching behaviour and the high quality factor are due to a careful design of the metasurface aimed at maximising the ohmic losses when the Nb is above the critical temperature and minimising the radiative coupling. The superconductor allows the operation of the cavity with high Q and the use of inductive elements with a high aspect ratio. Comparison with three dimensional finite element simulations highlights the crucial role of the inductive elements and of the kinetic inductance of the Cooper pairs in achieving the high quality factor and the high field enhancement.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Three-point bending analysis of doubly clamped silicon nanowire beams; Young's modulus, initial stress, and crystal orientation

A non-linear model is introduced describing the force-deflection relation of doubly clamped beams, including initial stress. Several approximations for the exact model are developed and compared, revealing the importance of considering the initial stress during 3-point bending measurements analysis. A novel approximation is found to be better than others, and both the exact model and this approximation are in perfect agreement with finite element simulations. A brief experimental example of silicon nanowires is presented in which the Young's modulus, the initial stress, and the crystallographic growth orientation are extracted by 3-point bending analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffuse interface approach to modeling crystal plasticity with accommodation of grain boundary sliding

In diffuse-interface (phase field) models of polycrystalline solids, the grain boundaries naturally possess a narrow region of finite thickness. A Ginzburg-Landau type kinetic equation for plastic deformation of a solid is derived that combines crystal plasticity and J 2 plasticity. Within such a framework for a polycrystal the grain interior can be defined by crystal plasticity, and the grain boundary region can be assigned J 2 type plasticity. The two are connected smoothly. This enables the modeling framework to accommodate grain boundary sliding (GBS), an important deformation mechanism for creep. The relaxation of elastic modulus of a polycrystalline solid and the power-law creep (both accommodated by GBS) are reported on by 2D and 3D simulations, respectively, and the results are compared to the 2D finite-element simulations of Ghahremani and Crossman-Ashby. A strong grain shape dependence and orientation dependence (for non-equiaxed grains) of the GBS effect are predicted by this model.

36 MATERIALS SCIENCE↗

Influence of Welding Strength Matching Coefficient and Cold Stretching on Welding Residual Stress in Austenitic Stainless Steel

Austenitic stainless steel is widely used in pressure vessels for the storage and transportation of liquid gases such as liquid nitrogen, liquid oxygen, and liquid hydrogen. Cryogenic pressure vessel manufacturing uses cold stretching technology, which relies heavily on welding joint performance, to construct lightweight and thin-walled vessels. Residual stress from welding is a primary factor in cases of austenitic stainless steel pressure vessel failure. In this paper, on the basis of Visual Environment 10.0 finite element simulation technology, the residual stress resulting from different welding strength matching coefficients (0.8, 1, 1.2, 1.4) for two S30408 plates welded with three-pass butt welds is calculated according to thermal elastoplastic theory. In addition, the stress field was calculated under a loading as high as 410 MPa and after the load was released. Path 1 was set to analyze stress along the welding line, and path 2 was set to analyze stress normal to the welding line. The welding strength matching coefficient strongly affected both the longitudinal residual stress (center of path 1) and the transverse residual stress (both ends of path 1) after the welding was completed. However, the coefficient had little effect on the longitudinal and transverse residual stress of path 2. Under the loading of 410 MPa, the longitudinal and transverse stress decreased and the stress distribution, with different welding strength matching coefficients, was less diverse. After the load was released, longitudinal and transverse stress distribution for both path 1 and path 2 decreased to a low level. Cold stretching could reduce the effect of residual stress to various degrees. Transverse strain along the stretching direction was also taken into consideration. The experimental results validated the reliability of the partial simulation.

36 MATERIALS SCIENCE↗