SEARCH · Engineering Papers

Results for “FEM 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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2020

Investigating Local Buckling and Plastic Collapse in Multifunctional Sandwich Composite Cores using a Unit Cell Method

The primary goal of developing analytical models of composite sandwich cores has been to obtain simple expressions for elastic and failure response to improve the rapid analysis of composite structures. An earlier investigation by the authors was made to verify the accuracy of widely-used analytical models that calculate equivalent elastic moduli of a uniform homogeneous single layer of material in place of the complex honeycomb core properties in sandwich composites. This transformation of elastic properties permits a coarsening of the finite element model (FEM) of the core to a fraction of the computational cost in simulating the original complex core architecture. In addition, the core was assumed constructed of multi-layered isotropic walls to support potential additional functionality beyond simple load-carrying capability.

Erik Saether↗

NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗

DOE OSTI2020

Phase Transformation of Thermally Aged and Neutron Irradiated Duplex Stainless Steel Used in LWRs. Final report

The lifetime of reactor components made of duplex stainless steels can be limited by the embrittlement from thermal aging, neutron irradiation or a synergistic effect between thermal aging and neutron irradiation. There is still a large scientific knowledge gap in understanding the phase evolution in duplex stainless steels upon thermal aging with or without neutron irradiation, the synergistic effect between thermal aging and neutron irradiation, and the impacts of structural evolution on material mechanical responses. In this project, totally, five major tasks were completed through this project: a) the in-situ WAXS tensile tests of thermally aged cast duplex stainless steels, b) the EXAFS investigation of thermally aged cast austenitic stainless steel and welds of austenite sainless steel, c) the characterization of irradiation effect on neutron irradiated CASS using APT, d) characterize the ferrite hardening of neutron irradiated CASS using nanoindentation, and e) the FEM modelling to validate the in-situ WAXS tensile tests and to simulate the stress/strain distribution during deformation for thermally aged duplex stainless steel. In summary, the thermal aging without neutron irradiation can significantly increase the strength of ferrite phase through the spinodal decomposition and G-phase precipitations. The in-situ X-ray experiments show that the hardening effect in ferrite phase can significantly alter the load partition between the ferrite and austenite phases. Atom probe tomography characterizations on samples irradiated to 0.08, 3, 5, 10, 20, and 40 dpa shows that neutron irradiation has dominated the microstructural evolution after 3 dpa, and the spinodal decomposition has saturated yet even at 40 dpa based on the quantified wavelength and amplitude.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Development of a Subchannel Capability for Liquid-Metal Fast Reactors in Pronghorn

This report details the development and demonstration of an entirely new capability in Pronghorn, namely the ability to model liquid-metal fast reactor (LMFR) flow conditions on the engineering scale. We developed two modeling approaches for LMFR that can be used separately or be combined into a hybrid simulation: (1) a modern subchannel capability called Pronghorn-Subchannel for square and hexagonal lattices, and (2) a porous flow capability for LMFR geometries. The report emphasizes the novel aspects of the developed subchannel capability and the interoperability of the subchannel capability, porous flow capability, and multiphysics tools within the multiphysics object oriented simulation environment (MOOSE). Here we demonstrate the ability to: (1) Accurately model subchannel flow in hexagonal lattices; (2) Couple the subchannel flow model to multidimensional finite-element method (FEM) or finite-volume method (FVM) heat conduction models; (3) Model LMFRs using Pronghorn’s porous media FVM approach; (4) Couple porous flow FVM and subchannel models in a single simulation; (5) Explicitly model inter-wrapper flows along with conjugate heat transfer from the intra-element flow; and (6) Demonstrate the numerical robustness of the subchannel algorithm by simulating intra-element flow recirculation in a high-buoyancy, low-flow fuel element.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2019

Evaluation of Impact Energy Attenuators and Composite Material Designs of a UAM VTOL Concept Vehicle

The development of Vertical Take-off and Landing (VTOL) vehicles for the Urban Air Mobility (UAM) markets presents a need for light weight vehicle structures with effective occupant protection capabilities. The National Aeronautics and Space Administration (NASA) has been working to fill that need, recently developing a cadre of concept vehicles to help characterize UAM design feasibility. This paper describes a study, using these concept vehicles, to evaluate the use of advanced composite structure and energy attenuating designs in the UAM vehicle design space. A finite element model (FEM) of a single passenger quadrotor concept vehicle was developed in LS-Dyna® and simulated under nominal and off-nominal vertical impact conditions. A variety of energy attenuating design mechanisms were implemented within this model to quantify their effectiveness in improving occupant safety. The use of carbon composites in both the energy attenuation mechanisms and vehicle structure was evaluated. The results of this study found significant reduction in occupant injury risk with the implementation of energy absorbing composite crush tubes and landing gear within the vehicle design. Additionally the use of a carbon fiber as a structural material was found to provide significant weight reduction while maintaining similar occupant loads to that predicted with an aluminum structure. This work provides a preliminary evaluation of design mechanisms and materials that may be used to optimize occupant protection capabilities within the UAM market.

Putnam, Jacob↗

NASA NTRS2020

Thermal crosstalk measurements and simulations for an X-ray microcalorimeter array

Arrays of high-density microcalorimeters require careful heat sinking in order to minimize the thermal crosstalk between nearby pixels. For the array of microcalorimeters developed for the Athena X-ray Integral Field Unit instrument, which has more than 3000 pixels on a 275 µm pitch, it is essential to address this problem in order to meet the energy-resolution requirements. The instrument’s energy-resolution budget requires that the impact of the thermal crosstalk on the energy resolution be a contribution that, added in quadrature to other energy-resolution contributions, is less than 0.2 eV. This value results in a derived requirement that the ratio between the amplitude of the crosstalk signal to an X-ray pulse (for example at 6 keV) is less than 1 × 10(exp −3) (for the first neighbor), less than 4 × 10(exp −4) (for the diagonal neighbor) and less than 8 × 10(exp −5) (for the second nearest neighbor). We have measured the thermal crosstalk levels between pixels in various geometries and configurations. The results show a crosstalk ratio which is at least a factor of 4 lower than the derived requirement. We also developed a finite element (FEM) 2D thermal model to predict the thermal behavior of large-scale arrays. This model successfully simulates the measured data in terms of pulse amplitude and time constants.

Antoine R. Miniussi↗

DOE OSTI2023

Simulation toolkit for digital material characterization of large image-based microstructures

In this paper, an efficient image-based simulation toolkit for material characterization is presented, which is scalable to work from personal computers to workstations. The effective thermal conductivity, elasticity, and permeability are evaluated employing a computational homogenization framework based on the Finite Element Method (FEM). Two complementary open-source packages are presented: one developed in Python, which can convert digital images into voxel meshes (pyTomoviewer); the other developed in Julia, that can run numerical simulations to compute effective material properties (chpack). Also, a CUDA C version of chpack is provided (chfem_gpu). They were designed to deal with large multi-phase models, so strategies were devised to minimize their memory footprint, while avoiding a high toll on execution time. The voxel-based approach significantly simplifies the FEM meshes and allows efficient matrix-free implementations. In that sense, to handle large linear systems of equations, the element-by-element (EBE) technique is adopted, in conjunction with a low-memory implementation of the Preconditioned Conjugate Gradient (PCG) method. Finally, the code was thoroughly tested on an artificial geometry made of a square array of cylinders, for which analytical solutions exist, as well as on a real micro-tomographic reconstruction of FiberForm TM , a carbon preform commonly used in thermal protection systems.

36 MATERIALS SCIENCE↗

NASA NTRS2014

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗

DOE OSTI2022

Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea

We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.

depressurization↗

NASA NTRS2005

Method for Estimating Thread Strength Reduction of Damaged Parent Holes with Inserts

During normal assembly and disassembly of bolted-joint components, thread damage and/or deformation may occur. If threads are overloaded, thread damage/deformation can also be anticipated. Typical inspection techniques (e.g. using GO-NO GO gages) may not provide adequate visibility of the extent of thread damage. More detailed inspection techniques have provided actual pitch-diameter profiles of damaged-hardware holes. A method to predict the reduction in thread shear-out capacity of damaged threaded holes has been developed. This method was based on testing and analytical modeling. Test samples were machined to simulate damaged holes in the hardware of interest. Test samples containing pristine parent-holes were also manufactured from the same bar-stock material to provide baseline results for comparison purposes. After the particular parent-hole thread profile was machined into each sample a helical insert was installed into the threaded hole. These samples were tested in a specially designed fixture to determine the maximum load required to shear out the parent threads. It was determined from the pristine-hole samples that, for the specific material tested, each individual thread could resist an average load of 3980 pounds. The shear-out loads of the holes having modified pitch diameters were compared to the ultimate loads of the specimens with pristine holes. An equivalent number of missing helical coil threads was then determined based on the ratio of shear-out loads for each thread configuration. These data were compared with the results from a finite element model (FEM). The model gave insights into the ability of the thread loads to redistribute for both pristine and simulated damage configurations. In this case, it was determined that the overall potential reduction in thread load-carrying capability in the hardware of interest was equal to having up to three fewer threads in the hole that bolt threads could engage. One- half of this potential reduction was due to local pitch-diameter variations and the other half was due to overall pitch-diameter enlargement beyond Class 2 fit. This result was important in that the thread shear capacity for this particular hardware design was the limiting structural capability. The details of the method development, including the supporting testing, data reduction and analytical model results comparison will be discussed hereafter.

Johnson, David L.↗

DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Computational Evaluation of Thermal Response of Open-Cell Foam With Circular Pore (Computional Evaluation of Thermal Response of Open-Cell Foam with Circular Pore)

The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Here, two-dimensional, steady state analyses were performed on unit cells with centered, circular pores using a second order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data.

36 MATERIALS SCIENCE↗

NASA NTRS2014

Calibration of International Space Station (ISS) Node 1 Vibro-Acoustic Model-Report 2

Reported here is the capability of the Energy Finite Element Method (E-FEM) to predict the vibro-acoustic sound fields within the International Space Station (ISS) Node 1 and to compare the results with simulated leak sounds. A series of electronically generated structural ultrasonic noise sources were created in the pressure wall to emulate leak signals at different locations of the Node 1 STA module during its period of storage at Stennis Space Center (SSC). The exact sound source profiles created within the pressure wall at the source were unknown, but were estimated from the closest sensor measurement. The E-FEM method represents a reverberant sound field calculation, and of importance to this application is the requirement to correctly handle the direct field effect of the sound generation. It was also important to be able to compute the sound energy fields in the ultrasonic frequency range. This report demonstrates the capability of this technology as applied to this type of application.

Zhang, Weiguo↗

DOE OSTI2021

Rapid simulations of hyperspectral near-field images of three-dimensional heterogeneous surfaces

The scattering-type scanning near-field optical microscope (s-SNOM) has emerged as a powerful tool for resolving nanoscale inhomogeneities in laterally heterogeneous samples. However, most analytical models used to predict the scattering near-field signals are assuming homogenous landscapes (bulk materials), resulting in inconsistencies when applied to samples with more complex configurations. In this work, we combine the point-dipole model (PDM) to the finite-element method (FEM) to account for the lateral and vertical heterogeneities while keeping the computation time manageable. Full images, spectra, or hyperspectral line profiles can be simulated by calculating the self-consistent dipole radiation demodulated at higher harmonics of the tip oscillation, mimicking real experimental procedures. Using this formalism, we clarify several important yet puzzling experimental observations in near-field images on samples with rich typography and complex material compositions, heterostructures of two-dimensional material flakes, and plasmonic antennas. The developed method serves as a basis for future investigations of nano-systems with nontrivial topography.

36 MATERIALS SCIENCE↗

NASA NTRS2001

SGBEM-FEM Alternating Method for Analyzing 3D Non-Planar Cracks and Their Growth in Structural Components

The highly accurate and efficient Symmetric Galerkin Boundary Element Method (SGBEM), a Finite Element Method (FEM)-based alternating method, is proposed for analyzing three-dimensional non-planar cracks and their growth. The cracks are modeled using the symmetric Galerkin boundary element method as a distribution of displacement discontinuities, simulating an infinite medium. The finite element method only analyzes the stress for the uncracked body. The solution for the cracked structural component is determined by an iteration procedure. This process alternates between an FEM solution for the uncracked body and the SGBEM solution for a crack in an infinite body. Numerical analysis, and the Java code used, evaluate stress intensity factors and model fatigue crack growth. Examples of non-planar cracks in infinite media and planar cracks in finite bodies, as well as growth under fatigue, show the accuracy of the method.

Nikishkov, G. P.↗

DOE OSTI2024

First considerations on the supporting structures of FCC-ee booster and collider in the arc regions

In 2022, the FCC Feasibility Study management mandated a working group to analyse the best configuration of the FCC-ee tunnel in the arc regions, in view of the construction of a mock-up of the arc half-cell. One of the main and most challenging goals of the study, named FCC-ee Arc Half-Cell Mock-up Project, was to perform a preliminary investigation on the principles of supporting the Short-Straight Sections and dipoles of the half-cells, both for the booster and for the collider machines. This is an important input needed for the choice of the best configuration of the relative placement of the booster with respect to the collider. The structural stiffness, mass and stability of the supporting structures must be optimized to minimize the vibrations transmitted/transferred to the magnetic system of the accelerators by elements such as pumps, water cooling system, beam thermomechanical stresses, powering elements, etc. To perform the study, tools such as CAD software, FEM and analytical techniques were employed. This paper summarizes the preliminary design concepts and the results of the simulations performed.

Accelerator Subsystems and Technologies↗

NASA NTRS2014

Elastic Model Transitions: a Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes both quadratically constrained least squares (LSQI) and Direct Shape Mapping (DSM) algorithms to determine physical displacements. This approach is applicable to the simulation of the elastic behavior of launch vehicles and other structures that utilize multiple LTI finite element model (FEM) derived mode sets that are propagated throughout time. The time invariant nature of the elastic data for discrete segments of the launch vehicle trajectory presents a problem of how to properly transition between models while preserving motion across the transition. In addition, energy may vary between flex models when using a truncated mode set. The LSQI-DSM algorithm can accommodate significant changes in energy between FEM models and carries elastic motion across FEM model transitions. Compared with previous approaches, the LSQI-DSM algorithm shows improvements ranging from a significant reduction to a complete removal of transients across FEM model transitions as well as maintaining elastic motion from the prior state.

Jurenko, Robert J.↗

DOE OSTI2021

Modeling and analyzing parasitic parameters in high frequency converters

This research focuses on electromagnetic interference (EMI) / electromagnetic compatibility (EMC) design and analysis in power electronics systems. To limit the EMI under the standards, different methods and strategies are investigated. Parasitic parameters of high frequency (HF) transformer are analyzed using a novel analytical method, finite element method (FEM), and experimental measurements for different structures and windings arrangements. Also, the magnetic field, electric field, electric displacement, and electric potential distribution are simulated and analyzed. Moreover, a high voltage system is considered and analyzed to improve the EMC. The EMI propagation paths are analyzed. The EMI noise level of the system is obtained and compared to the IEC61800-3 standard. To improve the EMC, the parasitic parameters of the transformer, as the main path of EMI circulation, are analyzed and optimized to block the propagation. Furthermore, the geometry structure of the HF transformer is optimized to lower the parasitics in the system. Three pareto-optimal techniques are investigated for the optimization. The models and results are verified by 3D-FEM and experimental results for several given scenarios. Furthermore, the EMC modeling and conducted EMI analysis are developed for a system including an AC-DC-DC power supply (rectifier and dual active bridge (DAB) converter). Moreover, the common mode (CM) EMI noise propagation through the system is discussed and the noise sources and effect of components on the noise are analyzed. Additionally, the CM impedance of different parts of the system and the noise levels are discussed. Finally, EMI attenuation techniques were applied to the system.

24 POWER TRANSMISSION AND DISTRIBUTION↗