SEARCH · Engineering Papers

Results for “finite element methods”

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 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2024

High-Power Polyphase PCB-Type Inductive Coupler for Wireless Electric Vehicle Charging

The cost and manufacturing difficulty of traditionally wound coils for wireless charging of electric vehicles has led to the exploration of alternative manufacturing for scalable deployment at a range of rated powers. A rotating-field polyphase wireless coupler is developed with balanced fully pitched coils which are not possible with conventional Litz wire. The usage of printed circuit board (PCB) technology is proposed for the innovative construction at large dimensions due to reduced cost, increased ease of manufacturability, and potential geometric customization compared to would coils. Design considerations specific to PCBs are detailed for wireless power transfer coils including advantages such as the potential for many turns and transposition with planar interconnections. A design procedure is developed which may be used to investigate the frequency to power to dimension correlation and limits associated with PCB coils. Designs are simulated in 3D electromagnetic finite element analysis (FEA) for the polyphase coil to evaluate the coupling between coils with varying geometric parameters at 100mm and 500mm outer diameters. A discussion is included on methods to increase PCB power ratings through a mixture of layer stacking, connections, and material selection.

Lewis, Donovin D.↗

DOE OSTI2024

Fabrication, Modeling, and Testing of a Prototype Thermal Energy Storage Containment

Increasing penetration of variable renewable energy resources requires the deployment of energy storage at a range of durations. Long-duration energy storage (LDES) technologies will fulfill the need to firm variable renewable energy resource output year round; lithium-ion batteries are uneconomical at these durations. Thermal energy storage (TES) is one promising technology for LDES applications because of its siting flexibility and ease of scaling. Particle-based TES systems use low-cost solid particles that have higher temperature limits than the molten salts used in traditional concentrated solar power systems. A key component in particle-based TES systems is the containment silo for the high-temperature (>1100 degrees C) particles. This study combined experimental testing and computational modeling methods to design and characterize the performance of a particle containment silo for LDES applications. A laboratory-scale silo prototype was built and validated the congruent transient finite element analysis (FEA) model. The performance of a commercial-scale silo was then characterized using the validated model. The commercial-scale model predicted a storage efficiency above 95% after 5 days of storage with a design storage temperature of 1200 degrees C. Insulation material and concrete temperature limits were considered as well. The validation of the methodology means the FEA model can simulate a range of scenarios for future applications. This work supports the development of a promising LDES technology with implications for grid-scale electrical energy storage, but also for thermal energy storage for industrial process heating applications.

clean energy↗

DOE OSTI2021

Method and apparatus for shape-based energy analysis of solids

A computer-readable medium stores instructions including storing a model of a physical structure and defining a mesh for the model. The mesh includes finite elements, each defined by a respective set of edges. The instructions include, for each finite element, identifying a governing differential equation and a set of complementary functions that exactly satisfy the governing differential equation. The instructions include determining an applied physical stimulus for the physical structure. The instructions include generating an energy optimization model that minimizes a difference between internal energy and external energy of the finite elements in response to the applied physical stimulus. The model includes a matrix of respective scalar multipliers for the complementary functions for each finite element. The instructions include transforming the matrix and calculating a physical parameter of interest. The instructions include, in response to the physical parameter not satisfying a design parameter, updating the model of the physical structure.

Spears, Robert E.↗

DOE OSTI2024

Delamination and Buckling Analysis of a Laminated Component in a High-speed Permanent Magnet Motor

High-speed permanent magnet (PM) machines are widely used because of their high-power density and high efficiency. The high rotation speed also inevitably subjects the PMs to high centrifugal load, which might damage them due to their inherent mechanical vulnerability, such as a much lower tensile strength than the compressive strength. To robustly transfer the torque from the magnet to the shaft, the outer diameter of the laminated rotor core is larger than the inner diameter of the rotor frame to ensure tight contact while working at 20,000 rpm. Motor manufacturing requires the shrink-fit method to assemble the rotor frame and rotor core. However, after the shrink-fit assembly, unexpected local delamination and buckling are observed on the rotor core part. Utilizing finite element simulation, we study the internal stress of assembling these two parts at the provided interference. Simulation results indicate the reasons for the delamination and buckling of the rotor core part and provide suggestions for improving the assembly.

Lin, Lianshan↗

DOE OSTI2022

Identification of Distribution of Elastic Constants of SiC/SiC Braided Tubes

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiCf/SiCm) composite tubes due to their superior mechanical properties. There are many parameters in the manufacturing of SiCf/SiCm composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. A sensitivity analysis will be carried out to determine the dominant parameters among the abovementioned ones. We intend to determine the variability in elastic constants of SiCf/SiCm composite tubes considering the variability in dominant parameters. A chosen number of combinations of dominant variables will be sampled through Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in circumferential and longitudinal direction, Poisson's ratio and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used a surrogate to FE analysis. once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The variability can be calculated based on MCS-generated data which will be presented at the conference.

42 ENGINEERING↗

DOE OSTI2020

An Initial Assessment of the Design Margins of Different ASME Section III, Division 5 Design Rules

This report develops a method for assessing the design margin of the ASME Section III, Division 5 Class A design rules, focusing on a definition of margin as the ratio between the actual, expected component life and the ASME design life. Full inelastic finite element simulations with a complete damage model capturing the available creep, creep-fatigue, fatigue, and tension test failure data produce the expected component lives, which can then be compared to corresponding ASME design calculations. The focus of this particular work is on Alloy 617 and the ASME creep-fatigue design rules, but the approach is general and could be applied to other design limits and to other materials. A margin assessment of a representative Alloy 617 component shows the ASME creep-fatigue design rules applied to this component have a margin of 10 on design life. The report describes several challenges in completing a more comprehensive margin assessment of the ASME rules, in particular challenges in developing probabilistic assessment methods for high temperature structural components. The report details these challenge and discusses the immediate possibilities of the deterministic margin assessment method developed here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗

DOE OSTI2020

Initial development and verification of a primary load design method based on elastic-perfectly plastic analysis

This report describes a new primary load design method to supplement the current design-by-elastic-analysis rules in Section III, Division 5, Subsection HB, Subpart B of the ASME Boiler and Pressure Vessel Code, covering the design and construction of Class A high temperature nuclear structural components. The main objectives for the new design method are to provide a procedure that can be applied to components with complicated geometries and to simplify the primary load design process by providing rules more compatible with modern finite element analysis. The report includes a complete set of design rules, presented as a draft ASME Code Case, as well as a commentary on the rules and a set of verification problems. The verification design problems demonstrate the new primary load design rules produce safe, efficient components when compared to the current Division 5 approach, while greatly simplifying the design analysis process.

42 ENGINEERING↗

DOE OSTI2021

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗

DOE OSTI2022

Mach Conference, 2022.

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiC/SiC) composite tubes due to their superior mechanical properties. The SiC/SiC tubes are known to remain stable even after prolonged exposure to radiation. However, there are many variables in the manufacturing of SiC/SiC composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. The intention is to identify the distribution of elastic constants of SiC/SiC tubes which would help model the performance of the novel cladding material. To do so, a sensitivity analysis will be carried out to determine the dominant variables that influence the elastic constants of SiC/SiC tubes. A chosen number of combinations of dominant variables will be sampled through the Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in the circumferential and longitudinal direction, Poisson's ratio, and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used as a surrogate to FE analysis. Once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The uncertainty can be quantified by calculating the coefficient of variation (CV) based on MCS-generated data which will be presented at the conference.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Bipolar Membranes with an Electrospun 3D Junction

Freestanding bipolar membranes (BPMs) with an entirely new and transformative morphology were fabricated and characterized. The key disruptive design element was a water-splitting/water-generation junction layer of finite thickness composed of interwoven and interlocking electrospun anion-exchange polymer (AEP) fibers and cation-exchange polymer (CEP) fibers. Methods were developed to sandwich the 3D junction layer between two pre-formed dense AEP and CEP films., with the closure of all interfiber voids, where the processing steps are easily incorporated into a commercial roll-to-roll membrane manufacturing scheme. A series of membranes were made using different polymers and junction layers to identify the optimum morphology and composition for water-splitting and water-generation applications. For water splitting, the 3D junction BPMs worked remarkably well with operating current densities at/above 1.0 A/cm2 (10X greater than commercial BPMs) at a transmembrane voltage drop of only 1.1 V. In water generation mode, the 3D junction BPMs operated at 0.5 A/cm2, a world-record current density. The membranes were found to exhibit outstanding durability and can be manufactured at scale for low cost. The high operating current densities and stable long-term operation are due to the 3D junction layer design, where there is a high interfacial area for water splitting or water generation reactions and where the interlocking fibers prevent delamination of the outer films. These new BPMs are ideal candidate materials for both existing water-splitting electrodialysis separations and new electrochemical processes, such as large-scale direct air capture of CO 2 , reactors for CO 2 reduction, self-hydrating fuel cells, and redox flow batteries.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2023

A hybrid deep neural operator/finite element method for ice-sheet modeling

One of the most challenging and consequential problems in climate modeling is to provide probabilistic projections of sea level rise. A large part of the uncertainty of sea level projections is due to uncertainty in ice sheet dynamics. At the moment, accurate quantification of the uncertainty is hindered by the cost of ice sheet computational models. In this work we develop a hybrid approach to approximate existing ice sheet models at a fraction of their cost. Our approach consists of replacing the finite element model for the momentum equations for the ice velocity, the most expensive part of an ice sheet model, with a Deep Operator Network, while we retain a classic finite element discretization for the evolution of the ice thickness. We show that the resulting hybrid model is very accurate and it is an order of magnitude faster than the traditional finite element model. Further, a distinctive feature of the proposed model, compared to other neural network approaches, is that it can handle high-dimensional parameter spaces (parameter fields) such as the basal friction at the bed of the glacier and can therefore be used for generating samples for uncertainty quantification. Further, we study the impact of hyper-parameters, number of unknowns and correlation length of the parameter distribution on the training and accuracy of the Deep Operator Network on a synthetic ice sheet model. We then target the evolution of the Humboldt glacier in Greenland and show that our hybrid model can provide accurate statistics of the glacier mass loss and can be effectively used to accelerate the quantification of uncertainty.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Quantifying leakage fields at ionic grain boundaries using off-axis electron holography

The electrical properties of interfaces in semiconductors and ionic conductors are immensely important in a wide range of applications. Electron holography is ideally suited for the direct measurement of the electrostatic potential of such interfaces. A key challenge with this approach is the contribution of the leakage field from the sample to the observed electron phase shift. This leakage field cannot be a priori independently determined and can cause an overestimation of the phase shift. In this work, we use finite element simulations to compute the three-dimensional electrostatic potential in the vicinity of an interface associated with a given interfacial charge density distribution. We then evaluate the predicted phase shift and demonstrate that the leakage field strongly affects the recovery of the projected interface potential. From the difference between the true potential and uncorrected, recovered potential, we propose a method to correct for this effect. We then demonstrate the application of this methodology to the analysis of experimental off-axis electron holography data acquired from the grain boundaries in lightly doped ceria.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Electromagnetic Modeling of Lattice Structures in Additively Manufactured Electric Machines

Recent advancements in metal additive manufacturing (AM) show great potential to revolutionize the design and manufacturing of electromagnetic components used in the field of electrical engineering. Lattice structures directly printed by AM processes typically offer better structural performance with reduced weight, such as high stiffness, surface area, elongation, energy absorption, and porosity, than the solid counterpart. This paper aims to study the electromagnetic modeling of lattice structures used in additively manufactured magnetic cores or windings for electric machines. Three dimensional (3D) electromagnetic finite element (FE) analysis with high performance computing (HPC) shows the highest fidelity in predicting the electromagnetic performance of designs with lattice structures by preserving complex geometry details. FE-based homogenization methods have also been explored to potentially speed up concept design. A case study based on an additively manufactured axial-flux permanent magnet machine with a Hilbert pattern stator validates the discussed electromagnetic modeling approaches.

Chung, Shi-Uk [ANSYS, Inc.]↗

DOE OSTI2024

Study of a High-speed Outer Rotor Burst Containing Enclosure

The increasing demand for high-speed electric machines in many applications pushes the development of high-power density electric motors. Outer rotor motors, which can be designed with a larger airgap diameter than inner rotor motors for the same overall diameter and, therefore, provide higher torque, distinguish themselves in space-constrained but high torque applications. However, increasing the motor speed and diameter results in significant centrifugal loads for the outer rotor motor due to its higher airgap diameter. Additional attention must be paid at the design stage to consider the extreme situation of possible mechanical failure, such as the sudden burst of the high-speed rotating parts in the motor, and measures must be taken to prevent their damage to other system components in advance. Through the finite element dynamic impact analysis, we studied the damage caused by broken parts on the inner wall of the electric motor housing when the motor rotates at 20,000 RPM. The proposed method is expected to provide reasonable recommendations of the enclosure material and its thickness required to protect other components in the same powertrain system when the high-power density motor fails unexpectedly.

Lin, Lianshan↗

DOE OSTI2023

Stellarator Simplification using Permanent Magnets (PM4Stell)

This Report describes the design and prototyping of an array of rare-Earth permanent magnets to form a stellarator. This effort was motivated by the hypothesis that the usage of permanent magnets, rather than electromagnetic coils with complex geometry, could reduce the cost of stellarator construction and thereby make increase the feasibility of the stellarator as a technology for a fusion-based power plant. In this project, we have developed novel methods for specifying the positions, shapes, and polarizations of the magnets in the array, and have developed designs for mounting structures and tooling for assembly. We have also performed detailed finite-element modeling to qualify the accuracy of the magnetic field produced by the magnet array as designed, and to confirm that the structure can withstand the forces between the magnets. We have also developed techniques for measuring the magnetic field produced by the array once constructed, as well as for correcting errors in the field arising from misalignments and offsets within the tolerances for mounting and fabrication. Finally, we have constructed a tabletop prototype of a section of the array to qualify the concept for assembling and mounting magnets within the array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Stellarator Simplification using Permanent Magnets (PM4Stell)

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Extent of interlocking and metallurgical bonding in friction riveting of aluminum alloy to steel

In this study, the joining of 6061-T6 aluminum alloy and DP590 steel using a M42 steel rivet via friction riveting technique is investigated. The surface morphology and microstructure characterization reveal the formation of an anchor zone that imparts mechanical interlock as well as the formation of metallurgical bonds at the interface of aluminum and steel. A combination of interlocking and bonding results in the achievement of a high load-carrying capacity of 5.7 kN during lap shear testing at room temperature. A finite element-based computational model was developed which accurately predicted the lap shear response of the joint. The model revealed that the metallurgical bond formed during fric-riveting adds 39% peak load strength to the joint. An extensive microstructural investigation, post-lap-shear fractography, and the modeling results, together provided insights on the joint failure mechanism. In conclusion, this study highlights that friction riveting is a promising method for aluminum-to-steel dissimilar joining, which is important for lighweighing automotive vehicles for energy efficiency.

36 MATERIALS SCIENCE↗