Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element problems”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A priori mesh quality metrics for three-dimensional hybrid grids

Use of general hybrid grids to attain complex-geometry field simulations poses a challenge on estimation of their quality. Apart from the typical problems of non-uniformity and non-orthogonality, the change in element topology is an extra issue to address. The present work derives and evaluates an a priori mesh quality indicator for structured, unstructured, as well as hybrid grids consisting of hexahedra, prisms, tetrahedra, and pyramids. Emphasis is placed on deriving a direct relation between the quality measure and mesh distortion. The work is based on use of the Finite Volume discretization for evaluation of first order spatial derivatives. The analytic form of the truncation error is derived and applied to elementary types of mesh distortion including typical hybrid grid interfaces. The corresponding analytic expressions provide relations between the truncation error and the degree of stretching, skewness, shearing, torsion, expansion, as well as the type of grid interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Performance Portability of the Aeras Atmosphere Model to Next Generation Architectures using Kokkos

The subject of this report is the performance portability of the Aeras global atmosphere dynamical core (implemented within the Albany multi-physics code) to new and emerging architecture machines using the Kokkos library and programming model. We describe the process of refactoring the finite element assembly process for the 3D hydrostatic model in Aeras and highlight common issues associated with development on GPU architectures. After giving detailed build and execute instructions for Aeras with MPI, OpenMP and CUDA on the Shannon cluster at Sandia National Laboratories and the Titan supercomputer at Oak Ridge National Laboratory, we evaluate the per- formance of the code on a canonical test case known as the baroclinic instability problem. We show a speedup of up to 4 times on 8 OpenMP threads, but we were unable to achieve a speedup on the GPU due to memory constraints. We conclude by providing methods for improving the performance of the code for future optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Application of parasound data for sediment study on methane seep site at Simeulue basin

The Parasound data presents sea depth and sub-bottom profiler. In terms of geological terminology, parasound data represents significant recent surface sedimentary structures that valuable for the selection of subsequent sampling site such as sampling at methane seep site. Therefore, Parasound is used to detailing methane seep at surface sediment following seismic data interpretation. In this study, parasound is used to focus observe area especially for sediment study on methane seep site. The Parasound systems works both as narrow beam sounder use high frequency and as sediment echosounder use low frequency. Parasound acquisition applies parametric effect. It produces additional frequency by nonlinear acoustic interaction of finite amplitude waves. Parasound transducers have 128 elements on 1 m2 and need transmission power up to 70 kW. The results of this study are discovered large seep carbonate with porous surface which means there are gas expulsions passing through that rock.

58 GEOSCIENCES↗
DOE OSTI2015

Failure of wooden sandwich beam reinforced with glass/epoxy faces

The mechanical properties and the failure of wooden beam strengthened with two faces from glass/epoxy composite and a wooden beam without strengthening was studied. Stresses and deflections on both beams, which are imposed in three point bending loading. On the idealized geometry of the specimens with detailed nonlinear orthotropic analysis was performed with a finite elements program. The failure study of the wooden beams was performed, applying the criterion of Tsai-Hill. The shear strength of the adhesive was taken into account. All the specimens were tested with three point bending loading and the experimental results were compared to those of the theoretical approach with the finite elements analysis. Comparing the results, the advantage of strengthened wooden beam against the simple wooden beam becomes obvious. Theoretical predictions were in good agreement with experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Inverse point source location with the Helmholtz equation on a bounded domain

The problem of recovering acoustic sources, more specifically monopoles, from point-wise measurements of the corresponding acoustic pressure at a limited number of frequencies is addressed. To this purpose, a family of sparse optimization problems in measure space in combination with the Helmholtz equation on a bounded domain is considered in this study. A weighted norm with unbounded weight near the observation points is incorporated into the formulation. Optimality conditions and conditions for recovery in the small noise case are discussed, which motivates concrete choices of the weight. The numerical realization is based on an accelerated conditional gradient method in measure space and a finite element discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Characterization of Flame Cut Heavy Steel: Modeling of Temperature History and Residual Stress Formation

Heavy steel plates are used in demanding applications that require both high strength and hardness. An important step in the production of such components is cutting the plates with a cost-effective thermal cutting method such as flame cutting. Flame cutting is performed with a controlled flame and oxygen jet, which burns the steel and forms a cutting edge. However, the thermal cutting of heavy steel plates causes several problems. A heat-affected zone (HAZ) is generated at the cut edge due to the steep temperature gradient. Consequently, volume changes, hardness variations, and microstructural changes occur in the HAZ. In addition, residual stresses are formed at the cut edge during the process. In the worst case, unsuitable flame cutting practices generate cracks at the cut edge. The flame cutting of thick steel plate was modeled using the commercial finite element software ABAQUS. The results of modeling were verified by X-ray diffraction-based residual stress measurements and microstructural analysis. The model provides several outcomes, such as obtaining more information related to the formation of residual stresses and the temperature history during the flame cutting process. In addition, an extensive series of flame cut samples was designed with the assistance of the model.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗
DOE OSTI2015

Local time–space mesh refinement for simulation of elastic wave propagation in multi-scale media

This paper presents an original approach to local time–space grid refinement for the numerical simulation of wave propagation in models with localized clusters of micro-heterogeneities. The main features of the algorithm are –the application of temporal and spatial refinement on two different surfaces; –the use of the embedded-stencil technique for the refinement of grid step with respect to time; –the use of the Fast Fourier Transform (FFT)-based interpolation to couple variables for spatial mesh refinement. The latter makes it possible to perform filtration of high spatial frequencies, which provides stability in the proposed finite-difference schemes. In the present work, the technique is implemented for the finite-difference simulation of seismic wave propagation and the interaction of such waves with fluid-filled fractures and cavities of carbonate reservoirs. However, this approach is easy to adapt and/or combine with other numerical techniques, such as finite elements, discontinuous Galerkin method, or finite volumes used for approximation of various types of linear and nonlinear hyperbolic equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Efficient solvers for hybridized three-field mixed finite element coupled poromechanics

We consider a mixed hybrid finite element formulation for coupled poromechanics. A stabilization strategy based on a macro-element approach is advanced to eliminate the spurious pressure modes appearing in undrained/incompressible conditions. The efficient solution of the stabilized mixed hybrid block system is addressed by developing a class of block triangular preconditioners based on a Schur-complement approximation strategy. Robustness, computational efficiency and scalability of the proposed approach are theoretically discussed and tested using challenging benchmark problems on massively parallel architectures.

42 ENGINEERING↗
DOE OSTI2021

Spent Nuclear Fuel Mechanical Loads in the General Package Drop Scenario

The US Department of Energy Spent Fuel and Waste Science and Technology program is performing research to determine the mechanical loading conditions applied to spent nuclear fuel (SNF) during normal conditions of transportation to inform mechanical tests of SNF and close an important knowledge gap related to the practical disposition of SNF in the US. A recent multi-national collaborative test campaign measured SNF assembly impact response to the 30 cm horizontal package drop scenario, which is a common regulatory basis test of SNF package design. Researchers at Pacific Northwest National Laboratory (PNNL) are using the test data to validate explicit finite element models to calculate the mechanical loads and structural response of spent nuclear fuel assemblies in the as-tested 30 cm horizontal package drop scenario. Once the as-tested package drop model is validated, the next step is to apply the model to the general 30 cm drop scenario, which includes all impact angles, all fuel assembly types, and all burnup conditions. PNNL is developing a damage model that will incorporate the results of finite element parametric studies to establish trends in SNF mechanical loading to various input parameters, like impact orientation and burnup. The damage model will have the capability to estimate mechanical loads for any single set of input parameters, but it will first be used to describe the upper bounds of potential SNF loading in the 30 cm drop scenario. This paper describes PNNL’s progress toward developing the general solution to the problem of spent nuclear fuel mechanical loads in the 30 cm package drop scenario and it describes the next steps in closing the knowledge gap.

Klymyshyn, Nicholas A.↗
DOE OSTI2019

An a posteriori residual-based spatial discretization error estimator for S N neutron transport

The “residual source estimator”, an implicit residual-based spatial discretization error estimator for S N neutron transport methods, is introduced. The estimator uses an approximation of the residual as a fixed source for an analogous transport problem whose solution is the error estimate, using the same transport operator as the original numerical solution. The complete derivation of the residual source estimator, the approximation of the residual, and higher-order derivatives is included. The performance of the residual source estimator for 0-order DGFEM is tested against reference solutions. For problems with true solutions discontinuous in the first derivative across the singular characteristics, the residual source estimator is generally accurate, precise, and computationally affordable compared to similar estimators. It suffers in problems where solutions are discontinuous, particularly when the discontinuity inflicts significant inaccuracies on the discrete solution; however, compared to other estimators, which suffer similarly, it remains a competitive option.

42 ENGINEERING↗
DOE OSTI2018

Parameter covariance and non-uniqueness in material model calibration using the Virtual Fields Method

Traditionally, material identification is performed using global load and displacement data from simple boundary-value problems such as uni-axial tensile and simple shear tests. More recently, however, inverse techniques such as the Virtual Fields Method (VFM) that capitalize on heterogeneous, full-field deformation data have gained popularity. In this work, we have written a VFM code in a finite-deformation framework for calibration of a viscoplastic (i.e. strain-rate dependent) material model for 304L stainless steel. Using simulated experimental data generated via finite-element analysis (FEA), we verified our VFM code and compared the identified parameters with the reference parameters input into the FEA. The identified material model parameters had surprisingly large error compared to the reference parameters, which was traced to parameter covariance and the existence of many essentially equivalent parameter sets. This parameter non-uniqueness and its implications for FEA predictions is discussed in detail. Lastly, we present two strategies to reduce parameter covariance – reduced parametrization of the material model and increased richness of the calibration data – which allow for the recovery of a unique solution.

304L stainless steel↗
DOE OSTI2022

Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma

The Naval Research Laboratory ‘Mag Noh problem’, described in this paper, is a self-similar magnetized implosion flow, which contains a fast magnetohydrodynamic (MHD) outward propagating shock of constant velocity. We generalize the classic Noh ( OSTI Tech. Rep. 577058, 1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$ . Our self-similar solutions are indeed realized when we solve the initial value MHD problem with the finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field and a shock wave.

Mechanics↗
DOE OSTI2017

A Variational Nodal Approach to 2D/1D Pin Resolved Neutron Transport for Pressurized Water Reactors

A two-dimensional/one-dimensional (2D/1D) variational nodal approach is presented for pressurized water reactor core calculations without fuel-moderator homogenization. A 2D/1D approximation to the within-group neutron transport equation is derived and converted to an even-parity form. The corresponding nodal functional is presented and discretized to obtain response matrix equations. Within the nodes, finite elements in the x-y plane and orthogonal functions in z are used to approximate the spatial flux distribution. On the radial interfaces, orthogonal polynomials are employed; on the axial interfaces, piecewise constants corresponding to the finite elements eliminate the interface homogenization that has been a challenge for method of characteristics (MOC)-based 2D/1D approximations. The angular discretization utilizes an even-parity integral method within the nodes, and low-order spherical harmonics (P N ) on the axial interfaces. The x-y surfaces are treated with high-order P N combined with quasi-reflected interface conditions. Furthermore, the method is applied to the C5G7 benchmark problems and compared to Monte Carlo reference calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

A novel pillar indentation splitting test for measuring fracture toughness of thin ceramic coatings

Fracture toughness is an important material property that plays a role in determining the in-service mechanical performance and adhesion of thin ceramic films. Unfortunately, measuring thin film fracture toughness is affected by influences from the substrate and the large residual stresses that can exist in the films. In this paper, we explore a promising new technique that potentially overcomes these problems based on nanoindentation testing of micro-pillars produced by focused ion beam milling of the films. By making the pillar diameter approximately equal to its length, the residual stress in the pillar’s upper portion is almost fully relaxed, and when indented with a sharp Berkovich indenter, the pillars fracture by splitting at reproducible loads that are readily quantified by a sudden displacement excursion in the load displacement behavior. Cohesive finite element simulations are used to analyze and develop, for a given material, a simple relation between the critical load at failure, pillar radius, and fracture toughness. The main novel aspect of this work is that neither crack geometries nor crack sizes need to be measured post test. Furthermore, the residual stress can be measured at the same time with toughness, by comparing the indentation results from the stress-free pillars and the as-deposited film. The method is tested on three different hard coatings formed by physical vapor deposition: titanium nitride, chromium nitride, and a CrAlN/Si 3 N 4 nanocomposite. Results compare well to independently measured values of fracture toughness for the three brittle films. The technique offers several benefits over existing methods.

42 ENGINEERING↗
DOE OSTI2025

A Misalignment Tolerant Magnetic Coupler Design and Implementation for Wireless Charging of Unmanned Aerial Vehicles (UAVs)

Wireless power transfer technology is getting more attention these days as a part of expanded battery charging infrastructure and network at a wide range of power levels. Charging Unmanned aerial vehicles (UAVs) is one wireless charging application, and charging their batteries can be challenging because they have short flight times and are frequently charged. For wireless charging of UAVs, the misalignment performance of the magnetic coupler pads is critical because landing on a charging pad at a perfectly aligned position is difficult and can be time-consuming. Different coil shapes and system structures have been proposed in the literature to overcome this challenge. A new design with a high misalignment tolerance is needed; therefore, to fill this gap, a wireless power transfer system has been designed for a 250 W UAV battery charging system with a new magnetic coupler design. A 3D model of the proposed design was created, and finite element simulations were performed. To prove the proposed design, the magnetic coupler prototypes were developed, and experimental results are presented herein. The misalignment tolerance of the proposed system was observed for various rotational positions of the secondary coil, from 0° to 360°. The results demonstrate that DC–DC efficiency varies between 80% and 85% within full 360° rotation of the receiver coil. This result enables the proposed system to address the misalignment problem for wireless charging of UAVs.

UAVs↗
DOE OSTI2016

Crashworthiness Models for Automotive Batteries - Report on Project 2088-A031-15 for DOT/NHTSA

Safety is a key element of any device designed to store energy, in particular electrochemical batteries, which convert energy of chemical reactions to electrical energy. Safety considerations are especially important when applied to large automotive batteries designed for propulsion of electric vehicles (EV). The high amount of energy stored in EV battery packs translates to higher probability of fire in case of severe deformation of battery compartment due to automotive crash or impact caused by road debris. While such demand for safety has resulted in heavier protection of battery enclosure, the mechanisms leading to internal short circuit due to deformation of the battery are not well understood even on the level of a single electrochemical cell. Moreover, not all internal shorts result in thermal runaway, and thus a criterion for catastrophic failure needs to be developed. This report summarizes the effort to pinpoint the critical deformation necessary to trigger a short via experimental study on large format automotive Li-ion cells in a rigid spherical indentation configuration. Cases of single cells and cell stacks undergoing indentation were investigated. Mechanical properties of cell components were determined via experimental testing and served as input for constitutive models of Finite Element (FE) analysis. The ability of the model to predict the behavior of cell(s) under spherical indentation and to predict failure leading to internal short circuit was validated against experiments. The necessity of resolving pairs of negative and positive electrodes in the FE formulation is clearly demonstrated by comparing layer-resolved simulations with simulations involving batteries with homogenized material properties. Finally, a coupled solution of electrochemical-electrical-thermal (EET) problem on a Nissan Leaf battery module was demonstrated towards the goal of extending the simulations to module level.

25 ENERGY STORAGE↗
DOE OSTI2017

Development of a Reduced-Order Three-Dimensional Flow Model for Thermal Mixing and Stratification Simulation during Reactor Transients

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the multi-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

Hu, Rui↗