Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element Method”

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

A layered solid finite element formulation with interlaminar enhanced displacements for the modeling of laminated composite structures

Accurate modeling of layered composite structures often requires the use of detailed finite element models which can sufficiently resolve the kinematics and material behavior within each layer of the composite. However, individually discretizing each material layer into finite elements presents a prohibitive computational expensive given the large number of thin layers comprising some laminated composites. To address these challenges, an 8-node layered solid hexahedral finite element is formulated with the aim of striking an appropriate balance between efficiency and fidelity. The element is discretized into an arbitrary number of distinct material layers, and employs reduced in-plane integration within each layer. The chosen reduced integration scheme is supplemented by a novel physical stabilization approach which includes layerwise enhancements to mitigate various forms of locking phenomena. The proposed framework additionally supports the inclusion of interlaminar enhanced displacements to better represent the kinematics of general layered composite materials. Finally, the described element formulation has been implemented in the ParaDyn finite element code, and its efficacy for modeling laminated composite structures is demonstrated on a variety of verification problems.

42 ENGINEERING↗
DOE OSTI2020

Linear viscoelastic characterization of electrically conductive adhesives used as interconnect in photovoltaic modules

Electrically conductive adhesives (ECAs) are incorporated into recent designs of photovoltaic (PV) modules and replace the traditional metallic solders as interconnects. This transition depicts a significant material change, and a proper understanding of the interconnects' mechanical response has not yet been established. However, such an understanding is necessary to (a) identify the driving forces for module degradation and failure, (b) allow for module design optimization, and (c) enable accurate lifetime predictions. This study summarizes the framework for the mechanical materials characterization and modeling of ECAs for PV applications. Only high-fidelity material models are able to capture the rate and temperature dependency of the ECA interconnect and allow for accurate modeling of the materials response. Furthermore, a linear viscoelastic representation is found to describe the mechanical response of the ECAs sufficiently well. The effects of curing conditions and environmental exposure are investigated, and material models for a variety of ECAs are reported and prepared for the use in numerical simulations. Overall, a finite element simulation of a generic submodel of a shingled cell module is used to highlight the need for high-fidelity material models and demonstrates the error made in the predicted stress states by using less sophisticated models.

14 SOLAR ENERGY↗
DOE OSTI2023

Multiscale numerical investigation of ratchet growth damage effects in PBX 9502

This paper presents results of numerical experiments conducted on the high explosive PBX 9502 to investigate how recently observed grain-scale damage mechanisms of ratchet growth affect uniaxial compression measurements. Simulations are multiscale in the sense of directly resolving grains, pores, cracks, and grain-interfaces based upon scanning electron microscope (SEM) images of damaged and undamaged samples. The combined finite-discrete element method (FDEM) is utilized to resolve both grain-scale microfracture and elastoplastic deformation of solid grains. Pristine (undamaged) and damaged microstructures are compared in simulation of unconfined compression tests of the same material from the literature. Here, the simulation results show the observed microscale mechanisms of damage, specifically microfracture predominantly around and sometimes through grains and crack-associated pore growth, can well-explain the effective degradation of strength and stiffness observed in the laboratory measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Thermomechanical Modeling of High-Temperature Bonded Interface Materials

In an automotive power electronics package, bonded interface materials play a critical role in providing a smooth heat dissipation path from the device to the coolant. However, these materials are susceptible to failure when the package gets exposed to repeated changes in its thermal environment and internal temperature variations. Under these thermal conditions, the coefficient of thermal expansion mismatch between the different component layers results in thermally induced stresses and consequently the delamination of the material from its adjoining surfaces or crack initiation and propagation within the material. Thermomechanical modeling offers a path to study the complex nature of these bonded materials, guides their design and selection, and, in conjunction with experimental results, plays an important role in predicting their lifetime. In this chapter, two modeling strategies - similar in approach and implementation but different in the underlying theory - are discussed in the context of high-temperature bonded materials. In addition, a short review of sintered silver is provided, as it is widely seen as a potentially promising and reliable bonded material for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simple, accurate surrogate models of the elastic response of three-dimensional open truss micro-architectures with applications to multiscale topology design

Elastic meta-materials are those whose unique properties come from their micro-architecture, rather than, e.g., from their chemistry. The introduction of such architecture, which is increasingly able to be fabricated due to advances in additive manufacturing, expands the design domain and enables improved design, from the most complex multi-physics design problems to the simple compliance design problem that is our focus. Unfortunately, concurrent design of both the micro-scale and the macroscale is computationally very expensive when the former can vary spatially, particularly in three dimensions. Instead, we provide simple, accurate surrogate models of the homogenized linear elastic response of the isotruss, the octet truss, and the ORC truss based on high-fidelity continuum finite element analyses. These surrogate models are relatively accurate over the full range of relative densities, in contrast to analytical models in the literature, which we show lose accuracy as relative density increases. The surrogate models are also simple to implement, which we demonstrate by modifying Sigmund’s 99-line code to solve a three-dimensional, multiscale compliance design problem with spatially varying relative density. Here, we use this code to generate examples in both two and three dimensions that illustrate the advantage of elastic meta-materials over structures with a single length scale, i.e., those without micro-architectures.

42 ENGINEERING↗
DOE OSTI2018

Optimal Control of a Degenerate PDE for Surface Shape

Controlling the shapes of surfaces provides a novel way to direct self-assembly of colloidal particles on those surfaces and may be useful for material design. This motivates the investigation of an optimal control problem for surface shape in this paper. Specifically, we consider an objective (tracking) functional for surface shape with the prescribed mean curvature equation in graph form as a state constraint. The control variable is the prescribed curvature. We prove existence of an optimal control, and using improved regularity estimates, we show sufficient differentiability to make sense of the first order optimality conditions. This allows us to rigorously compute the gradient of the objective functional for both the continuous and discrete (finite element) formulations of the problem. Numerical results are shown to illustrate the minimizers and optimal controls on different domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Examination of Mean Stress Calculation Approaches in Rock Mechanics

The mean stress, as a fundamental statistical property of a group of stress data, is essential for stress variability characterisation. However, currently in rock mechanics, the mean stress is customarily and erroneously calculated by separately averaging the principal stress magnitudes and orientations. In order to draw the attention of our community to the appropriate approach for stress variability characterisation, here we compare the customary scalar/vector mean with that obtained by the mathematically rigorous tensorial approach—the Euclidean mean. Calculation of mean stress using both a small group of actual in situ stress measurement results and a large group of simulated stress data (obtained using the combined finite– discrete element method, FEMDEM) demonstrates that the two approaches yield different results. Further investigation of these results shows that the scalar/vector approach may yield non-unique and non-orthogonal mean principal stresses, and these may deviate significantly from those of the Euclidean mean. Our calculations and comparisons reveal that the scalar/ vector approach is deficient because it processes the principal stress magnitudes and orientations separately as independent quantities and ignores the connection between them. Conversely, the tensorial approach appropriately averages the tensors that simultaneously carry not only the information of stress magnitudes and orientations, but also the inherent relations between them. Therefore, arbitrarily using scalar/vector mean stress of in situ stress measurements as input in further rock engineering analyses may yield significantly erroneous results. Finally, we advise that stress data should be statistically processed in a tensorial manner using tensors referred to a common Cartesian coordinate system.

58 GEOSCIENCES↗
DOE OSTI2018

Impacts of hydrogeological characteristics on groundwater-level changes induced by earthquakes

Changes in groundwater level during earthquakes have been reported worldwide. In this study, field observations of co-seismic groundwater-level changes in wells under different aquifer conditions and sampling intervals due to near-field earthquake events in Taiwan are presented. Sustained changes, usually observed immediately after earthquakes, are found in the confined aquifer. Oscillatory changes due to the dynamic strain triggered by passing earthquake waves can only be recorded by a high-frequency data logger. While co-seismic changes recover rapidly in an unconfined aquifer, they can sustain for months or longer in a confined aquifer. Three monitoring wells with long-term groundwater-level data were examined to understand the association of co-seismic changes with local hydrogeological conditions. The finite element software ABAQUS is used to simulate the pore-pressure changes induced by the displacements due to fault rupture. The calculated co-seismic change in pore pressure is related to the compressibility of the formation. The recovery rate of the change is rapid in the unconfined aquifer due to the hydrostatic condition at the water table, but slow in the confined aquifer due to the less permeable confining layer. Fracturing of the confining layer during earthquakes may enhance the dissipation of pore pressure and induce the discharge of the confined aquifer. The study results indicated that aquifer characteristics play an important role in determining groundwater-level changes during and after earthquakes.

58 GEOSCIENCES↗
DOE OSTI2018

Hydro-mechanically coupled finite-element analysis of the stability of a fractured-rock slope using the equivalent continuum approach: a case study of planned reservoir banks in Blaubeuren, Germany

Construction works associated with the building of reservoirs in mountain areas can damage the stability of adjacent valley slopes. Seepage processes caused by the filling and drawdown operations of reservoirs also affect the stability of the reservoir banks over time. The presented study investigates the stability of a fractured-rock slope subjected to seepage forces in the lower basin of a planned pumped-storage hydropower (PSH) plant in Blaubeuren, Germany. The investigation uses a hydro-mechanically coupled finite-element analyses. For this purpose, an equivalent continuum model is developed by using a representative elementary volume (REV) approach. To determine the minimum required REV size, a large number of discrete fracture networks are generated using Monte Carlo simulations. These analyses give a REV size of 28 × 28 m, which is sufficient to represent the equivalent hydraulic and mechanical properties of the investigated fractured-rock mass. The hydro-mechanically coupled analyses performed using this REV size show that the reservoir operations in the examined PSH plant have negligible effect on the adjacent valley slope.

58 GEOSCIENCES↗
DOE OSTI2018

Experimental Validation of the Earthquake Resistance of a Full-Scale Model of RBMK-1000 Graphite Stack with Post-Repair Characteristics

In providing justification for the operational safety of RBMK-1000 power units with graphite stack distortion, earthquake resistance tests were performed on a full-scale model of RBMK graphite stack under seismic platform conditions. The tests were performed for three mutually perpendicular axes of the model successively for seismic loads corresponding to magnitudes 6, 7 and 8 on the MSK-64 scale. The tests confirmed the earthquake resistance of the graphite stack for prescribed loads. The results were used to verify the FEMGR code (Kurchatov Institute) and the LEGAK-DF finite-element software system (RFYaTs – VNIIEF) for purposes of computational validation of the earthquake resistance of the post-repair restored RMBK-1000 core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Analysis of the Resonance Phenomena in the Operation of the Lower Bureya HPP During an Earthquake

To justify that hydroelectric power plants operating in seismic areas are safe, it is necessary to analyze the resonance phenomena in the generating unit. The analysis was performed using a finite-element model of the powerhouse-foundation-equipment system of the Lower Bureya HPP. The resonance phenomena occurring when the disturbance frequency at the elevation of the stator supports is close to the first natural frequency and speed of the water turbine were studied.

42 ENGINEERING↗
DOE OSTI2018

Study of Dissipative Losses in AC-Biased Mo/Au Bilayer Transition-Edge Sensors

We are developing kilo-pixel arrays of transition-edge sensors (TESs) for the X-ray Integral Field Unit on ESA’s Athena observatory. Previous measurements of AC-biased Mo/Au TESs have highlighted a frequency-dependent loss mechanism that results in broader transitions and worse spectral performance compared to the same devices measured under DC bias. In order to better understand the nature of this loss, we are now studying TES pixels in different geometric configurations. We present measurements on devices of different sizes and with different metal features used for noise mitigation and X-ray absorption. Our results show how the loss mechanism is strongly dependent upon the amount of metal in close proximity to the sensor and can be attributed to induced eddy current coupling to these features. We present a finite element model that successfully reproduces the magnitude and geometry dependence of the losses. Using this model, we present mitigation strategies that should reduce the losses to an acceptable level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Stability Considerations Using a Microscopic Stability Model Applied to a 2G Thin Film Coated Superconductor

A numerical, finite element simulation of local temperature fields, critical current and transport current distributions is applied to a 2G thin film coated high-temperature superconductor. The focus is on simulation of quenches originating from transport current locally exceeding critical current density. As in previous reports, a statistical treatment of superconductor parameters (critical temperature, current and magnetic field) is applied for this analysis; it shall take into account uncertainties possibly arising from shortcomings in conductor manufacture or during measurements of their properties. The results of the calculations are used as quasi-boundary (driving) conditions for a subsequent transient microscopic numerical stability analysis. Traditionally, stability analysis usually applies analytic stability functions that incorporate conventional, phonon-related timescales, t and disturbances. The question is whether decay of electron pairs and subsequent relaxation of the excited state to a new dynamic equilibrium carrier concentration, the “electron aspect” of the stability problem, under the same disturbances, proceeds on another timescale, t{sup ′}. Is this timescale identical to the traditional (phonon) timescale, t? If not, how large are the differences? The recently reported microscopic stability model, now applied to a thin film conductor, is consulted to find an answer to these questions. A time limit, t{sub Quench}, can be extracted from the simulations, as the time of immediate onset of a quench. This is a new approach in stability considerations since, conventionally, a temperature limit, T < T{sub Crit}, is set as stability criterion. If electrical operation or cooling conditions cannot immediately respond to this challenge, the thin film superconductor, if this time limit is exceeded, will hardly be able to return to zero loss current transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Investigation on Effect of Magnetic Field Dependency Coefficient of Critical Current Density on the AC Magnetizing Loss in HTS Tapes Exposed to External Field

High-temperature superconducting (HTS) technology is the most promising, advanced, and beneficial technological developments of the last two decades in terms of fabricating more efficient, compact, and reliable electrical power apparatus. The second-generation HTS wires are now quite competitive with traditional conductors from the cost point of view. In addition, they are of high interest in electrical machine manufacturing due to their unique merits such as high current density, high magnetic field, and low loss. The most significant issue for any superconducting devices which work with alternating current or magnetic field is the AC loss of the HTS tapes. Since HTS tapes are always exposed to an external magnetic field in electrical machine applications, it is vital to consider the effect of this external field on their AC magnetizing (ACMG) loss. It is well-known that external magnetic field decreases the current carrying capability of the wire. In real life, critical current density of HTS tapes is dependent to magnetic field but this dependency sometime was not considered in modeling stage in order to make it faster and simple. In this paper, the electromagnetic properties of HTS tapes for use in superconducting electric machines have been investigated with a particular focus on the ACMG loss of YBCO tapes. For this purpose, a two-dimensional finite element model has been implemented based on the H formulation in order to study the effect of dependency of critical current density to magnetic field on the variation of ACMG loss. The modeling results will be helpful for electrical machine designers to have a better understanding about this effect and have a more accurate estimate of ACMG loss in design stage.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Calculation of AC Magnetizing Loss of ReBCO Superconducting Tapes Subjected to Applied Distorted Magnetic Fields

In recent years, there are fast-increasing concerns on the utilizations of superconducting rotating electrical machines in different application areas, such as ship propulsion systems, aircraft drivers, and wind turbine generators, since these machines exhibit the merits of high current density, compact design, high power density, light weight, high torque density as well as high efficiency. One of the main limitations in front of the vast use of superconducting tapes in the fabrication of electrical machineries is AC magnetizing loss when tapes are exposed to an external magnetic field, which can decrease the critical current density of wires, as well. In the literature, most of the research works have been done on calculation of the AC magnetizing loss under a pure external magnetic field, while in reality, magnetic flux lines in AC electrical machines are usually distorted with harmonics because of different reasons such as distorted leakage flux, distributed coils of a winding in several slots, cogging fields, mechanical faults, etc. Since these distorted fields contain harmonics, then in this paper, the AC magnetizing loss of superconducting tapes has been electromagnetically modeled and calculated when they are subjected to nonsinusoidally distorted external magnetic fields. The magnetic field dependency of critical current density has been considered in a proposed finite element model. The results have shown that the AC magnetizing loss increases significantly under a distorted applied field compared with a sinusoidal one. In addition, the loss increase depends on the harmonic content which would increases drastically with total harmonic distortion of the applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Shallow Embedded Strain Measurements and Analysis for a NbTi Superconducting Sextupole Coil During Cooling, Excitation, and Quench

The present study deals with strain measurements and their numerical analysis for a niobium titanium wire wounded superconducting sextupole coil in a multi-field operation environment, including cooling down, excitation, and quench. Four low-temperature resistance strain gauges with designed compensation circuits were affixed directly on the surface of the superconducting wires but within glass fiber cloth. Before the sextuple coil’s test was carried out, for the purpose of reliability, based on the designed compensation circuit, some basic experiments were carried out during cooling down, excitation, and quench. The experimental results showed that the shallow embedded low-temperature resistance strain gauges were still highly reproducible with an error of 1% for a variable magnetic field and an error of 5% for a variable thermal field by executing various modes. It shows multi-stage strain monitoring is feasible for the superconducting sextupole. In addition, a numerical multiphysics model for the sextuple coil is also built from Maxwell’s equations and the equilibrium equations for mechanical and thermal deformation. A self-tested elastic modulus of NbTi/Cu stack with epoxy, which is cut from a NbTi/Cu sextupole prototype, at low temperature is used in this simulation. The numerical results have reasonable agreement with the experimental observations and show the multiphysics model with self-tested elastic modulus can be applied to analysis strain behaviors under complicated fields before superconducting magnet assembly. The present results will also provide some strain-measuring technology and multiphysics numerical analysis on an LTS/HTS superconducting magnet developed at the Institute of Modern Physics, Chinese Academy of Sciences (IMP, CAS).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Effect Laws of Different Factors on Levitation Characteristics of High-T{sub c} Superconducting Maglev System with Numerical Solutions

The interaction between the high-temperature superconducting (HTS) bulk and permanent magnet guideway (PMG) has been a topic among the existing HTS magnetic levitation systems. Up to now, the theories about the electromagnetic and force characteristics have been relatively consummated. Here, using the finite element model by COMSOL Multiphysics 5.3a, we investigated the impact of the intrinsic and external factors on the levitation characteristics of HTS maglev systems, the application laws are concluded. Factors including the critical current density and geometric dimensions of the bulk superconductor, as well as the moving velocity, are discussed. The effects of the multiple back and forth movements are also analyzed. By comparing calculation results of the distributions of the induced current and the levitation force of the bulk superconductor, we can clearly confirm the effect laws of various factors on the electromagnetic characteristics and levitation force. These observations support the physical parameters that are difficult to be tested by experiments. And they also have a significant contribution to the experimental condition selections and optimization of the high-temperature superconducting magnetic levitation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Remanent Magnetic Flux Distribution in Superconducting-Ferromagnetic Layered Heterostructures

Stacks of superconducting tape, made of second-generation coated conductors, present an interesting alternative to bulk superconductors. The layered geometry allows incorporation of other functional, namely ferromagnetic, materials. The ferromagnetic material changes the field configuration within the superconductor and can be used to shape or divert the flux density from its interior. This can be used to optimise the remanent flux density after magnetisation. FEM simulations are employed to test the performance of a variety of superconductor/ferromagnet structures and explain the role of the ferromagnetic material in the composite structure. Cases where trapped field is much higher and lower than remanent magnetisation of the ferromagnetic material are considered, showing that performance can be enhanced or maintained by using ferromagnetic material while reducing the usage of costly superconducting tape.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗