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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2015

Alignment of graphene oxide nanostructures between microgap electrodes via dielectrophoresis for hydrogen gas sensing applications

Graphene oxide (GO) nanostructures have been aligned between conducting electrodes via dielectrophoresis (DEP) with different electrical configurations. The arrangement of ground with respect to peak-to-peak voltage (V{sub pp}) plays a crucial role in manipulating the GO nanostructures. Grounds on both sides of the V{sub pp} electrode give an excellent linking of GO nanostructures which is explained by scanning electron microscopy and current-voltage characteristics. A finite element method simulation explains the electric field and voltage variation profile during DEP process. The optimized aligned GO nanostructures are used as hydrogen gas sensor with a sensitivity of 6.0% for 800 ppm hydrogen gas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Finite Element Modeling of Spent Nuclear Fuel in Dry Storage under Seismic Excitation

Spent nuclear fuel (SNF) is currently stored at independent spent fuel storage installations (ISFSIs) at nuclear power plant sites across the US. The US Department of Energy (DOE) Spent Fuel and Waste Science and Technology (SFWST) program is currently working to understand the response of SNF during a seismic event. This investigation is part of a long-term effort by DOE to quantify the external mechanical loads on SNF. The SFWST program has recently completed campaigns evaluating SNF behavior during normal conditions of transport and in 30 cm package drop scenarios; both campaigns included a combination of testing and modeling. Likewise, the investigation into seismic response of SNF will include full-size shake table tests of a vertical storage overpack containing instrumented surrogate fuel assemblies, as well as pre- and post-test finite element modeling. This work focuses on the development of finite element models simulating the planned test configuration. The seismic response is calculated at two scales: at the package level, and at the fuel level. The package level model evaluates the behavior of the storage overpack and its contents in response to the ground motions of a representative earthquake. The motion of the package is input to the fuel level model, which consists of a detailed finite element model of a 17x17 or 16x16 pressurized water reactor (PWR) fuel assembly or of a 10x10 boiling water reactor (BWR) fuel assembly. The stress and strain on individual components such as fuel rods and guide tubes are extracted from this model. This work details the SNF response to various input ground motions, considering representative earthquakes in the Central and Eastern US under various site conditions. Also discussed is the development of a “reduced scope” finite element model of the SNF, which consists of a single fuel rod while retaining the dynamic response of the entire fuel assembly. This reduced scope model addresses the significant computational expense of modeling the SNF under long-duration input ground motions and facilitates simulating additional parameter variations of the model beyond the as-tested configuration. The modeling results provide important insight that will help to inform future testing and modeling efforts by the SFWST program to close the knowledge gap related to external loads on SNF.

Kadooka, Kevin↗
DOE OSTI2018

Simulation of Conduction Cooled SRF Cavity

SRF cavities are used heavily in high energy physics research, and simplifying operation of SRF cavities could enable their use in many fields of industry. Liquid helium coolant systems are a notable complexity of operating SRF cavities, and replacing liquid helium with cryocooler conduction cooling is an important research and development interest. As part of this effort, a system consisting of a 650 MHz single cell cavity, a commercial grade pulse tube refrigerator, and an aluminum thermal link was simulated through finite element analysis. The electromagnetic fields inside the cavity, heat flow through the system, and temperature distribution of the thermal link and cavity were calculated. Based on the simulations, the pulse tube refrigerator is capable of supplying adequate cooling power to a cavity with \ce{Nb3Sn} coating operating with an accelerating gradient between 6 and 12 MV/m. The maximum simulated temperature of the cavity was below 7 K, and the oscilla tions in temperature introduced by the pulse tube refrigerator were shown to dampen in amplitude before reaching the cavity surface.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Simulation of Conduction Cooled SRF Cavity

SRF cavities are used heavily in high energy physics research, but the operational complexities of SRF cavities prevent their widespread use in industry. One notable complexity of operating an SRF cavity is liquid helium coolant, and replacing liquid helium with cryocooler conduction cooling is an important research and development interest. As part of this effort, a system consisting of a 650 MHz single cell cavity, a commercial grade pulse tube refrigerator, and a thermal conduction link was simulated through finite element analysis. The electromagnetic fields inside the cavity, heat flow through the system, and temperature distribution of the thermal link and cavity were calculated. Based on the simulation, we found that a 2 W @ 4.2 K pulse tube refrigerator is capable of cooling a single cell Nb 3 Sn coated cavity operating at up to 12 MV/m accelerating gradient. The maximum simulated temperature of the cavity was below 7 K, and the thermal conduction link was shown to significantly attenuate the cryocooler thermal oscillations at the location of the cavity

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Evolution of A-Type Macrosegregation in Large Size Steel Ingot After Multistep Forging and Heat Treatment

A-type macrosegregation refers to the channel chemical heterogeneities that can be formed during solidification in large size steel ingots. In this research, a combination of experiment and simulation was used to study the influence of open die forging parameters on the evolution of A-type macrosegregation patterns during a multistep forging of a 40 metric ton (MT) cast, high-strength steel ingot. Macrosegregation patterns were determined experimentally by macroetch along the longitudinal axis of the forged and heat-treated ingot. Mass spectroscopy, on more than 900 samples, was used to determine the chemical composition map of the entire longitudinal sectioned surface. FORGE NxT 1.1 finite element modeling code was used to predict the effect of forging sequences on the morphology evolution of A-type macrosegregation patterns. For this purpose, grain flow variables were defined and implemented in a large scale finite element modeling code to describe oriented grains and A-type segregation patterns. Examination of the A-type macrosegregation showed four to five parallel continuous channels located nearly symmetrical to the axis of the forged ingot. In some regions, the A-type patterns became curved or obtained a wavy form in contrast to their straight shape in the as-cast state. Mass spectrometry analysis of the main alloying elements (C, Mn, Ni, Cr, Mo, Cu, P, and S) revealed that carbon, manganese, and chromium were the most segregated alloying elements in A-type macrosegregation patterns. The observed differences were analyzed using thermodynamic calculations, which indicated that changes in the chemical composition of the liquid metal can affect the primary solidification mode and the segregation intensity of the alloying elements. Finite element modeling simulation results showed very good agreement with the experimental observations, thereby allowing for the quantification of the influence of temperature and deformation on the evolution of the shape of the macrosegregation channels during the open die forging process.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gyrokinetic electromagnetic particle simulations in triangular meshes with C1 finite elements

The triangular mesh-based gyrokinetic scheme enables comprehensive axis-to-edge studies across the entire plasma volume. Our approach employs triangular finite elements with first-derivative continuity (C1), building on previous work to facilitate gyrokinetic simulations. Additionally, we have adopted the mixed variable/pullback scheme for gyrokinetic electromagnetic particle simulations. The filter-free treatment in the poloidal cross-section with triangular meshes introduces unique features and challenges compared to previous treatments using structured meshes. Our implementation has been validated through benchmarks using ITPA-toroidicity-induced Alfvén eigenmode parameters, showing its capability in moderate to small electron skin depth regimes. Additional examinations using experimental parameters confirm its applicability to realistic plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Material Variability Effects on Damage Development Within Composite Adhesively Bonded Joints

Use of composite adhesively bonded joints (ABJ) is of critical importance to the adoption of composite materials in the automotive industry, as ABJ enable lower stress concentrations as compared to conventional mechanically fastened joints. ABJ are better suited for joining composite materials as compared to fastened joints because fastened joints require drilling of holes which locally affect composite material structure. Composite materials and adhesives are subject to unavoidable stochastic local material variations which make different failure scenarios possible. An experimentally tested ABJ configuration is simulated using finite element analysis (FEA). Experimentally, under tension, the joints failed by three major failure modes, with peak loads ranging from 13.0-16.1 kips. Progressive failure analysis tools are used to simulate damage development within each material within the joint. The simulation agreed well with the average experimental peak load. Stochastically occurring adhesive porosity and matrix-fiber micro-disbonding were numerically simulated. The simulations revealed a similar trend as observed experimentally: joints which failed at higher peak loads had lower levels of damage within the face-sheets of the composite panels which were adhesively bonded; these joints which failed at higher peak loads also had greater damage in the doubler of the experimentally tested double lap joint configuration.

Richard Larson↗
DOE OSTI2015

Multi-element least square HDMR methods and their applications for stochastic multiscale model reduction

Stochastic multiscale modeling has become a necessary approach to quantify uncertainty and characterize multiscale phenomena for many practical problems such as flows in stochastic porous media. The numerical treatment of the stochastic multiscale models can be very challengeable as the existence of complex uncertainty and multiple physical scales in the models. To efficiently take care of the difficulty, we construct a computational reduced model. To this end, we propose a multi-element least square high-dimensional model representation (HDMR) method, through which the random domain is adaptively decomposed into a few subdomains, and a local least square HDMR is constructed in each subdomain. These local HDMRs are represented by a finite number of orthogonal basis functions defined in low-dimensional random spaces. The coefficients in the local HDMRs are determined using least square methods. We paste all the local HDMR approximations together to form a global HDMR approximation. To further reduce computational cost, we present a multi-element reduced least-square HDMR, which improves both efficiency and approximation accuracy in certain conditions. To effectively treat heterogeneity properties and multiscale features in the models, we integrate multiscale finite element methods with multi-element least-square HDMR for stochastic multiscale model reduction. This approach significantly reduces the original model's complexity in both the resolution of the physical space and the high-dimensional stochastic space. We analyze the proposed approach, and provide a set of numerical experiments to demonstrate the performance of the presented model reduction techniques. - Highlights: • Multi-element least square HDMR is proposed to treat stochastic models. • Random domain is adaptively decomposed into some subdomains to obtain adaptive multi-element HDMR. • Least-square reduced HDMR is proposed to enhance computation efficiency and approximation accuracy in certain conditions. • Integrating MsFEM and multi-element least square HDMR can significantly reduce computation complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Design optimization of MR-compatible rotating anode x-ray tubes for stable operation

Purpose: Hybrid x-ray/MR systems can enhance the diagnosis and treatment of endovascular, cardiac, and neurologic disorders by using the complementary advantages of both modalities for image guidance during interventional procedures. Conventional rotating anode x-ray tubes fail near an MR imaging system, since MR fringe fields create eddy currents in the metal rotor which cause a reduction in the rotation speed of the x-ray tube motor. A new x-ray tube motor prototype has been designed and built to be operated close to a magnet. To ensure the stability and safety of the motor operation, dynamic characteristics must be analyzed to identify possible modes of mechanical failure. In this study a 3D finite element method (FEM) model was developed in order to explore possible modifications, and to optimize the motor design. The FEM provides a valuable tool that permits testing and evaluation using numerical simulation instead of building multiple prototypes.Methods: Two experimental approaches were used to measure resonance characteristics: the first obtained the angular speed curves of the x-ray tube motor employing an angle encoder; the second measured the power spectrum using a spectrum analyzer, in which the large amplitude of peaks indicates large vibrations. An estimate of the bearing stiffness is required to generate an accurate FEM model of motor operation. This stiffness depends on both the bearing geometry and adjacent structures (e.g., the number of balls, clearances, preload, etc.) in an assembly, and is therefore unknown. This parameter was set by matching the FEM results to measurements carried out with the anode attached to the motor, and verified by comparing FEM predictions and measurements with the anode removed. The validated FEM model was then used to sweep through design parameters [bearing stiffness (1×10{sup 5}–5×10{sup 7} N/m), shaft diameter (0.372–0.625 in.), rotor diameter (2.4–2.9 in.), and total length of motor (5.66–7.36 in.)] to increase the fundamental frequency past the operating range at 50 Hz.Results: The first large vibration during the prototype motor operation was obtained at 21.64±0.68 Hz in the power spectrum. An abrupt decrease in acceleration occurred at 21.5 Hz due to struggling against the resonance vibrations. A bearing stiffness of 1.2×10{sup 5} N/m in the FEM simulation was used to obtain a critical speed of 21.4 Hz providing 1.1% error. This bearing stiffness value and the 3D model were then confirmed by the experiments with the anode removed, demonstrating an agreement within 6.4% between simulation results and measurements. A calculated first critical frequency (fundamental frequency) of 68.5 Hz was obtained by increasing the bearing stiffness to 1×10{sup 7} N/m and increasing the shaft diameter by 68.0%. Reducing the number of bearings in the design permits decreasing the total length of the motor by 1.7 in., and results in a fundamental frequency of 68.3 Hz in concert with additional changes (shaft diameter of 0.625 in., rotor diameter of 2.4 in., and bearing stiffness of 1×10{sup 6} N/m).Conclusions: An FEM model of the x-ray tube motor has been implemented and experimentally validated. A fundamental frequency above the operational rotation speed can be achieved through modification of multiple design parameters, which allows the motor to operate stably and safely in the MR environment during the repeated acceleration/deceleration cycles required for an interventional procedure. The validated 3D FEM model can now be used to investigate trade-offs between generated torque, maximum speed, and motor inertia to further optimize motor design.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2019

A Multi-Physics Study on High-Specific Power Li-O2 Batteries for Electric Aircraft

Commercialization of lithium-air batteries faces many challenges, such as electrolyte decomposition, short cycle life, low energy and power density, etc. However, commercialization of Li-O2 batteries for aeronautics is much more challenging due to additional safety constraints on cyclability and performance (high specific power and specific energy). For this presentation, we will discuss inter-related aspects of physics-based modeling of a pack: cell and battery model calibration. In addition, we will evaluate and present optimal battery designs for high discharge current density, high discharge time, and low battery mass using simulation-based optimization.The Finite Element Model (FEM) used to simulate a Li-O2 cell is based on the work of Bevara [1]. The different aspects of the model are based on: porous electrode theory and concentrated electrolyte theory; quantum tunneling model for the resistance of conformal layer of discharge product (Li2O2) [1]; Butler-Volmer kinetics for electrochemical reaction; Fick's diffusion for oxygen transport; and an oxygen dissolution model is applied at the air/electrolyte interface [2]. The electrolyte properties such as ion conductivity, ion diffusion, oxygen diffusion, and mass density of the electrolyte were taken from Molecular Dynamics (MD) simulations [3]; while the other model parameters, which includes mass of cell components, were calibrated to match experiments at high discharge current densities. The cell mass includes the anode, cathode, separator, electrolyte, and other components (such as current collector). This calibrated model is used to perform parametric studies on cathode thickness, porosity, tortuosity, carbon particle size, electrolyte transport and material properties, partial pressure of oxygen, discharge time, and discharge current density to study optimal designs for high specific power and energy. References:1. Bevara, V. & Andrei, P. (2014), J. Electrochem. Soc. 161 (14), A2068-A2079.2.Mehta, M. & Andrei, P. (2015), J. Power Sources. 286, 299-308.3.Liyana-Arachchi, T.; Haskins, J.; Burke, C.; Diederichsen, K.; McCloskey, B.; & Lawson, J. (2018), J. Phys. Chem. B. 122 (36), 8548 - 8559.4.Choi, W.; Kikumoto, H.; Choudhary, R. & Ooka, R. (2018), Applied Energy, 209, 306-321.

Mehta, Mohit↗
NASA NTRS2002

Comparisons of Kinematics and Dynamics Simulation Software Tools

Kinematic and dynamic analyses for moving bodies are essential to system engineers and designers in the process of design and validations. 3D visualization and motion simulation plus finite element analysis (FEA) give engineers a better way to present ideas and results. Marshall Space Flight Center (MSFC) system engineering researchers are currently using IGRIP from DELMIA Inc. as a kinematic simulation tool for discrete bodies motion simulations. Although IGRIP is an excellent tool for kinematic simulation with some dynamic analysis capabilities in robotic control, explorations of other alternatives with more powerful dynamic analysis and FEA capabilities are necessary. Kinematics analysis will only examine the displacement, velocity, and acceleration of the mechanism without considering effects from masses of components. With dynamic analysis and FEA, effects such as the forces or torques at the joint due to mass and inertia of components can be identified. With keen market competition, ALGOR Mechanical Event Simulation (MES), MSC visualNastran 4D, Unigraphics Motion+, and Pro/MECHANICA were chosen for explorations. In this study, comparisons between software tools were presented in terms of following categories: graphical user interface (GUI), import capability, tutorial availability, ease of use, kinematic simulation capability, dynamic simulation capability, FEA capability, graphical output, technical support, and cost. Propulsion Test Article (PTA) with Fastrac engine model exported from IGRIP and an office chair mechanism were used as examples for simulations.

Yeu-Sheng Paul Shiue↗
NASA NTRS2022

Laser Patterning Assisted Devitrification and Domain Engineering of Amorphous and Nanocrystalline Alloys

We have investigated laser-patterning and devitrification as a novel method to realize emergent properties in modified regions at the ribbon surface consisting of periodic localized laser spots thermally treated under rapid heating and cooling conditions. Fe-based amorphous ribbons were annealed for systematically varied laser power. Extremely high heating and cooling rates on the order of 10(exp 8)-10(exp 10)K/s have been estimated by finite element analysis simulations. Observations suggest surface melting followed by rapid solidification fast enough to quench and form an amorphous structure upon cooling. Diffusion of heat occurs from the laser irradiated surface region so that a finite volume of material surrounding this region rises in temperature above the crystallization temperature relevant for conventional isothermal annealing experiments. The underlying mechanism of laser annealing and consequences of heat transients involved are discussed in terms of impacts upon micro/nanostructure, residual stresses, and magnetic domain structure surrounding the laser irradiated region. The study illustrates the potential to exploit spatially optimized phase transformations in a scalable manufacturing process of amorphous and nanocrystalline alloys to locally access otherwise inaccessible extreme heating and cooling rates.

Amorphous ribbons↗
NASA NTRS1993

A class of hybrid finite element methods for electromagnetics: A review

Integral equation methods have generally been the workhorse for antenna and scattering computations. In the case of antennas, they continue to be the prominent computational approach, but for scattering applications the requirement for large-scale computations has turned researchers' attention to near neighbor methods such as the finite element method, which has low O(N) storage requirements and is readily adaptable in modeling complex geometrical features and material inhomogeneities. In this paper, we review three hybrid finite element methods for simulating composite scatterers, conformal microstrip antennas, and finite periodic arrays. Specifically, we discuss the finite element method and its application to electromagnetic problems when combined with the boundary integral, absorbing boundary conditions, and artificial absorbers for terminating the mesh. Particular attention is given to large-scale simulations, methods, and solvers for achieving low memory requirements and code performance on parallel computing architectures.

Volakis, J. L.↗
DOE OSTI2016

Theoretical study on the top- and enclosed-contacted single-layer MoS{sub 2} piezotronic transistors

Recently, the piezotronic effect has been observed in two-dimensional single-layer MoS{sub 2} materials, which have potential applications in force and pressure triggered or controlled electronic devices, sensors, and human-machine interfaces. However, classical theory faces the difficulty in explaining the mechanism of the piezotronic effect for the top- and enclosed-contacted MoS{sub 2} transistors, since the piezoelectric charges are assumed to exist only at the edge of the MoS{sub 2} flake that is far from the electronic transport pathway. In the present study, we identify the piezoelectric charges at the MoS{sub 2}/metal-MoS{sub 2} interface by employing both the density functional theory and finite element method simulations. This interface is on the transport pathway of both top- and enclosed-contacted MoS{sub 2} transistors, thus it is capable of controlling their transport properties. This study deepens the understanding of piezotronic effect and provides guidance for the design of two-dimensional piezotronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Mechanical properties of niobium radio-frequency cavities

Radio-frequency cavities made of bulk niobium are one of the components used in modern particle accelerators. The mechanical stability is an important aspect of cavity design, which typically relies on finite-element analysis simulations using material properties from tensile tests on sample. This contribution presents the results of strain and resonant frequency measurements as a function of a uniform pressure up to 722 kPa, applied to single-cell niobium cavities with different crystallographic structure, purity and treatments. In addition, burst tests of high-purity multi-cell cavities with different crystallographic structure have been conducted up to the tensile strength of the material. Finite-element analysis of the single-cell cavity geometry is in good agreement with the observed behavior in the elastic regime assuming a Young's modulus value of 88.5 GPa and a Poisson's ratio of 0.4, regardless of crystallographic structure, purity or treatment. However, the measured yield strength and tensile strength depend on crystallographic structure, material purity and treatment. In particular, the results from this study show that the mechanical properties of niobium cavities with large crystals are comparable to those of cavities made of fine-grain niobium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗