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 1,405 records · Page 78

An improved three-point bending test method for the investigation of nanosecond laser dicing of ultrathin Si dies with Cu stabilization layer

Highlights: • A selective chemical etching method has been developed for removal of metallization layers on ultrathin Si dies for three-point bending test. • The complex elastoplastic stress distribution in the die is simplified to an elastic stress distribution. • The die sidewall fracture strength can be determined using the standard three-point bending tensile stress equation. • The characteristic fracture strengths of ultrathin Si dies were correlated to the laser-induced microstructures of the die sidewall. - Abstract: The sidewall fracture strength of ultrathin Si dies with frontside and backside metallization layers cannot be determined directly from three-point bending (3PB) test. The complex 3PB elastoplastic stress distribution makes it difficult to determine the die sidewall fracture strength by analytical method or finite element analysis due to various theoretical and practical constraints. In this work, a selective chemical etching method has been developed for complete removal of the frontside and backside metallization layers on the ultrathin die without affecting the critical microstructure and strength of the die sidewall induced by laser dicing. By removing all the metallization layers on the die, the complex 3PB elastoplastic stress distribution is simplified to an elastic stress distribution. This enables a straightforward calculation of the die sidewall fracture strength using the standard 3PB tensile stress equation for a homogeneous linear elastic material with a rectangular cross section. Two chemical etchants were evaluated and found to be suitable but FeCl{sub 3} is preferred because of the much faster etching time. With this improved 3PB test method, the effect of nanosecond laser dicing on ultrathin Si dies with backside Cu layer was investigated. The die frontside and backside characteristic fracture strengths of ultrathin Si dies with varying backside Cu thickness were determined and correlated to TEM analyses of the compositions and microstructures at the fracture origins.

36 MATERIALS SCIENCE↗

Numerical Simulation of Textile Composite Stamping On Double Dome

Stamping is one of the most effective ways to form textile composites in industry for providing high-strength, low-weight and cost-effective products. This paper presents a fully continuum mechanics-based approach for stamping simulation of textile fiber reinforced composites by using finite element (FE) method. A previously developed non-orthogonal constitutive model is used to represent the anisotropic mechanical behavior of textile composites under large deformation during stamping. Simulation are performed on a balanced plain weave composite with 0 deg./90 deg. and {+-}45 deg. as initial yarn orientation over a benchmark double dome device. Simulation results show good agreement with experimental output in terms of a number of parameters selected for comparison.

36 MATERIALS SCIENCE↗

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗

Measurement of Work Hardening Behavior of Pure Titanium Sheet Using A Servo-Controlled Tube Bulge Testing Apparatus

Biaxial stress tests of rolled pure titanium sheet (JIS 1, 0.5 mm thick) have been carried out in order to investigate the anisotropic plastic deformation under biaxial tension. Rolled pure titanium sheet was bent and welded to make tubular specimens. Combined tension-internal pressure was applied to the tubular specimens using the servo-controlled tube bulge testing apparatus developed by one of the authors [Kuwabara, T., Yoshida, K., Narihara, K., Takahashi S., Int. J. Plasticity 21 (1), 101-117 (2002)], so that the strain rate ratio, {epsilon}{sub {phi}}:{epsilon}{theta}, in the axial ({phi}) and circumferential ({theta}) directions of the specimen was controlled to be constant. Contours of plastic work at different levels of plastic strain and stress paths under constant strain rate ratios have been observed in the first quadrant of stress space. It is found that the test material exhibits significant differential work hardening behavior with the increase of plastic work.

36 MATERIALS SCIENCE↗

A Forward Time Stepping Heat Conduction Model for Spot Melt Additive Manufacturing

Solidification dynamics are crucial for determining microstructure development in additively manufactured parts. Multiphysics models based on finite element or finite volume methods may help gain insight for complicated phenomena such as fluid flow, keyholing, and porosity but are too computationally expensive to use for simulating actual builds. Recent analytic and semi-analytic solutions for moving heat sources in a semi-infinite three-dimensional space provide a way to accurately estimate the solidification conditions for entire builds. The downside to these methods is that, unlike finite element or finite volume methods, they cannot use the temperature distribution of the previous timesteps to march the solution forward in time. This paper provides the mathematical formulation and implementation of a forward time stepping (FTS) approach to an existing semi-analytic solution. The speed and accuracy of the two methods are then compared for various scan patterns. The result is that, for spot melts, the forward time-stepping model provides improvements in both speed and accuracy. This is especially true for longer simulations, where the simulation can be orders of magnitude faster. Finally, the longest simulation analyzed in this paper was roughly 30× faster when using the forward time-stepping model versus the straightforward implementation of the semi-analytic solution.

36 MATERIALS SCIENCE↗

Collaborative Proposal: Transforming How Climate System Models are Used: A Global, Multi-Resolution Approach

Despite the great interest in regional modeling for both weather and climate applications, regional modeling is not yet at the stage that it can be used routinely and effectively for climate modeling of the ocean. The overarching goal of this project is to transform how climate models are used by developing and implementing a robust, efficient, and accurate global approach to regional ocean modeling. To achieve this goal, we will use theoretical and computational means to resolve several basic modeling and algorithmic issues.

54 ENVIRONMENTAL SCIENCES↗

Characterization of high order spatial discretizations and lumping techniques for discontinuous finite element SN transport

We examine several possible methods of mass matrix lumping for discontinuous finite element discrete ordinates transport using a Lagrange interpolatory polynomial trial space. Though positive outflow angular flux is guaranteed with traditional mass matrix lumping in a purely absorbing 1-D slab cell for the linear discontinuous approximation, we show that when used with higher degree interpolatory polynomial trial spaces, traditional lumping does yield strictly positive outflows and does not increase in accuracy with an increase in trial space polynomial degree. As an alternative, we examine methods which are 'self-lumping'. Self-lumping methods yield diagonal mass matrices by using numerical quadrature restricted to the Lagrange interpolatory points. Using equally-spaced interpolatory points, self-lumping is achieved through the use of closed Newton-Cotes formulas, resulting in strictly positive outflows in pure absorbers for odd power polynomials in 1-D slab geometry. By changing interpolatory points from the traditional equally-spaced points to the quadrature points of the Gauss-Legendre or Lobatto-Gauss-Legendre quadratures, it is possible to generate solution representations with a diagonal mass matrix and a strictly positive outflow for any degree polynomial solution representation in a pure absorber medium in 1-D slab geometry. Further, there is no inherent limit to local truncation error order of accuracy when using interpolatory points that correspond to the quadrature points of high order accuracy numerical quadrature schemes. (authors)

97 MATHEMATICS AND COMPUTING↗

Modelling of Superplastic Forming of AZ31 Magnesium Alloy

In this study the constitutive equation of the superplastic AZ31 magnesium-based alloy is modelled by the power law relationship between the stress, the strain and the strain-rate and an accurate procedure for determining the constants of the material is presented. Moreover, the problem of optimizing the pressure-time load curve of a free forming process is investigated and resolved by means of a pressure jump forming process. The experimental tests, carried out to support the finite-element modelling, have shown good agreement between the numerical results and the experimental data.

36 MATERIALS SCIENCE↗

A model of ultrasonic testing for cracks near a non-planar surface

Based on the methodology of the 2D models presented at QNDE 2008 and QNDE 2009, this paper presents an extension to 3D of a model of ultrasonic testing for cracks near a non-planar back surface. The elastic wave scattering in a thick-walled component by an interior rectangular crack located near a back surface of arbitrary geometry is considered. The 3D wave scattering problem is solved using a hybrid boundary integral equation method (BIEM): a displacement boundary element method (BEM) for the back surface displacements is combined with an analytical technique for the hypersingular traction integral equation for the crack-opening displacements. The action of ultrasonic transducers in transmission and reception is taken into account in the model, and a few numerical examples illustrating the influence of the back surface geometry are given.

42 ENGINEERING↗

Development on electromagnetic impedance function modeling and its estimation

Today the Electromagnetic methods such as magnetotellurics (MT) and controlled sources audio MT (CSAMT) is used in a broad variety of applications. Its usefulness in poor seismic areas and its negligible environmental impact are integral parts of effective exploration at minimum cost. As exploration was forced into more difficult areas, the importance of MT and CSAMT, in conjunction with other techniques, has tended to grow continuously. However, there are obviously important and difficult problems remaining to be solved concerning our ability to collect process and interpret MT as well as CSAMT in complex 3D structural environments. This talk aim at reviewing and discussing the recent development on MT as well as CSAMT impedance functions modeling, and also some improvements on estimation procedures for the corresponding impedance functions. In MT impedance modeling, research efforts focus on developing numerical method for computing the impedance functions of three dimensionally (3-D) earth resistivity models. On that reason, 3-D finite elements numerical modeling for the impedances is developed based on edge element method. Whereas, in the CSAMT case, the efforts were focused to accomplish the non-plane wave problem in the corresponding impedance functions. Concerning estimation of MT and CSAMT impedance functions, researches were focused on improving quality of the estimates. On that objective, non-linear regression approach based on the robust M-estimators and the Hilbert transform operating on the causal transfer functions, were used to dealing with outliers (abnormal data) which are frequently superimposed on a normal ambient MT as well as CSAMT noise fields. As validated, the proposed MT impedance modeling method gives acceptable results for standard three dimensional resistivity models. Whilst, the full solution based modeling that accommodate the non-plane wave effect for CSAMT impedances is applied for all measurement zones, including near-, transition-as well as the far-field zones, and consequently the plane wave correction is no longer needed for the impedances. In the resulting robust impedance estimates, outlier contamination is removed and the self consistency between the real and imaginary parts of the impedance estimates is guaranteed. Using synthetic and real MT data, it is shown that the proposed robust estimation methods always yield impedance estimates which are better than the conventional least square (LS) estimation, even under condition of severe noise contamination. A recent development on the constrained robust CSAMT impedance estimation is also discussed. By using synthetic CSAMT data it is demonstrated that the proposed methods can produce usable CSAMT transfer functions for all measurement zones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Physics-based shape matching for intraoperative image guidance

Purpose: Soft-tissue deformations can severely degrade the validity of preoperative planning data during computer assisted interventions. Intraoperative imaging such as stereo endoscopic, time-of-flight or, laser range scanner data can be used to compensate these movements. In this context, the intraoperative surface has to be matched to the preoperative model. The shape matching is especially challenging in the intraoperative setting due to noisy sensor data, only partially visible surfaces, ambiguous shape descriptors, and real-time requirements. Methods: A novel physics-based shape matching (PBSM) approach to register intraoperatively acquired surface meshes to preoperative planning data is proposed. The key idea of the method is to describe the nonrigid registration process as an electrostatic–elastic problem, where an elastic body (preoperative model) that is electrically charged slides into an oppositely charged rigid shape (intraoperative surface). It is shown that the corresponding energy functional can be efficiently solved using the finite element (FE) method. It is also demonstrated how PBSM can be combined with rigid registration schemes for robust nonrigid registration of arbitrarily aligned surfaces. Furthermore, it is shown how the approach can be combined with landmark based methods and outline its application to image guidance in laparoscopic interventions. Results: A profound analysis of the PBSM scheme based on in silico and phantom data is presented. Simulation studies on several liver models show that the approach is robust to the initial rigid registration and to parameter variations. The studies also reveal that the method achieves submillimeter registration accuracy (mean error between 0.32 and 0.46 mm). An unoptimized, single core implementation of the approach achieves near real-time performance (2 TPS, 7–19 s total registration time). It outperforms established methods in terms of speed and accuracy. Furthermore, it is shown that the method is able to accurately match partial surfaces. Finally, a phantom experiment demonstrates how the method can be combined with stereo endoscopic imaging to provide nonrigid registration during laparoscopic interventions. Conclusions: The PBSM approach for surface matching is fast, robust, and accurate. As the technique is based on a preoperative volumetric FE model, it naturally recovers the position of volumetric structures (e.g., tumors and vessels). It cannot only be used to recover soft-tissue deformations from intraoperative surface models but can also be combined with landmark data from volumetric imaging. In addition to applications in laparoscopic surgery, the method might prove useful in other areas that require soft-tissue registration from sparse intraoperative sensor data (e.g., radiation therapy)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Variational multiscale element-free Galerkin method combined with the moving Kriging interpolation for solving some partial differential equations with discontinuous solutions

Solving partial differential equations with discontinuous solutions is an important challenging problem in numerical analysis. To this end, there are some methods such as finite volume method, discontinuous Galerkin approach and particle technique that are able to solve these problems. In the current paper, the moving Kriging element-free Galerkin method has been combined with the variational multiscale algorithm to obtain acceptable and high-resolution solutions. For testing this technique, we select some PDEs with discontinuous solution such as Burgers’, Sod’s shock tube, advection–reaction–diffusion, Kuramoto–Sivashinsky, Boussinesq and shallow water equations. First, we obtain a time-discrete scheme by approximating time derivative via finite difference technique. Then we introduce the moving Kriging interpolation and also obtain their shape functions. We use the element-free Galerkin method for approximating the spatial derivatives. This method uses a weak form of the considered equation that is similar to the finite element method with the difference that in the classical element-free Galerkin method test and trial functions are moving least squares approximation (MLS) shape functions. Since the shape functions of moving least squares (MLS) approximation do not have Kronecker delta property, we cannot implement the essential boundary condition, directly. Thus, we employ the shape functions of moving Kriging interpolation and radial point interpolation technique which have the mentioned property. Also, in the element-free Galekin method, we do not use any triangular, quadrangular or other type of meshes. The element-free Galerkin method is a global method while finite element method is a local one. This technique employs a background mesh for integration which makes it different from the truly mesh procedures. Several test problems are solved and numerical simulations are reported which confirm the efficiency of the proposed schemes.

97 MATHEMATICS AND COMPUTING↗

The evolution history of superalloy powders during hot consolidation and plastic deformation

Highlights: • Element segregations in a powder metallurgy superalloy at different states are identified. • The surface energy drives precipitation of oxides, carbides, and γ' particles on powder surface during hot consolidation. • Prior particle boundary structure is distorted and more uniformly refined grains are obtained after hot extrusion. - Abstract: Powder metallurgy superalloys are widely used to produce critical rotating parts in turbine engines. However, the issues like what is the driving force of precipitation on prior particle boundary (PPB) during consolidation and how the PPB structure evolves during plastic deformation still remain to be resolved. In this study, with the characterization of the microstructure evolution from powder, consolidation, to plastic deformation, the segregation of alloying elements on the particle boundaries at elevated temperature was analyzed by equilibrium thermodynamic calculations, and the distortion of PPB structure under large plastic deformation was quantitively simulated by finite element (FE) methods. The results indicated that surface energy serves as the main driving force of PPB precipitation and hot extrusion could distort the PPB structure and yield more uniformly refined grains.

36 MATERIALS SCIENCE↗

SCEPTRE 1.4 Quick Start Guide

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general purpose C++ code for solving the Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates. Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Party Libraries (TPL) to be available, and example scripts for building these TPL’s are provided. The TPL’s needed by SCEPTRE are Trilinos, boost, and netcdf. SCEPTRE uses an autoconf build system, and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance analysis of higher mode spoof surface plasmon polariton for terahertz sensing

We investigated the spoof surface plasmon polaritons (SSPPs) on 1D grooved metal surface for terahertz sensing of refractive index of the filling analyte through a prism-coupling attenuated total reflection setup. From the dispersion relation analysis and the finite element method-based simulation, we revealed that the dispersion curve of SSPP got suppressed as the filling refractive index increased, which cause the coupling resonance frequency redshifting in the reflection spectrum. The simulated results for testing various refractive indexes demonstrated that the incident angle of terahertz radiation has a great effect on the performance of sensing. Smaller incident angle will result in a higher sensitive sensing with a narrower detection range. In the meanwhile, the higher order mode SSPP-based sensing has a higher sensitivity with a narrower detection range. The maximum sensitivity is 2.57 THz/RIU for the second-order mode sensing at 45° internal incident angle. The proposed SSPP-based method has great potential for high sensitive terahertz sensing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗