Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗
DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

The many-body Wigner Monte Carlo method for time-dependent ab-initio quantum simulations

The aim of ab-initio approaches is the simulation of many-body quantum systems from the first principles of quantum mechanics. These methods are traditionally based on the many-body Schrödinger equation which represents an incredible mathematical challenge. In this paper, we introduce the many-body Wigner Monte Carlo method in the context of distinguishable particles and in the absence of spin-dependent effects. Despite these restrictions, the method has several advantages. First of all, the Wigner formalism is intuitive, as it is based on the concept of a quasi-distribution function. Secondly, the Monte Carlo numerical approach allows scalability on parallel machines that is practically unachievable by means of other techniques based on finite difference or finite element methods. Finally, this method allows time-dependent ab-initio simulations of strongly correlated quantum systems. In order to validate our many-body Wigner Monte Carlo method, as a case study we simulate a relatively simple system consisting of two particles in several different situations. We first start from two non-interacting free Gaussian wave packets. We, then, proceed with the inclusion of an external potential barrier, and we conclude by simulating two entangled (i.e. correlated) particles. The results show how, in the case of negligible spin-dependent effects, the many-body Wigner Monte Carlo method provides an efficient and reliable tool to study the time-dependent evolution of quantum systems composed of distinguishable particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Exploring numerical blow-up phenomena for the Keller–Segel–Navier–Stokes equations

The Keller–Segel–Navier–Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial organism density mass is below 2π there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality. Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfy a priori bounds as well as lower and L 1 (Ω) bounds for the organism and chemoattractant densities. In particular, these latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures. Our findings show that the existence threshold value 2π encountered for the organism density mass may not be optimal and hence it is conjectured that the critical threshold value 4π may be inherited from the fluid-free Keller–Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Implementation and Evaluation of the Virtual Fields Method: Determining Constitutive Model Parameters From Full-Field Deformation Data.

The Virtual Fields Method (VFM) is an inverse method for constitutive model parameter identication that relies on full-eld experimental measurements of displacements. VFM is an alternative to standard approaches that require several experiments of simple geometries to calibrate a constitutive model. VFM is one of several techniques that use full-eld exper- imental data, including Finite Element Method Updating (FEMU) techniques, but VFM is computationally fast, not requiring iterative FEM analyses. This report describes the im- plementation and evaluation of VFM primarily for nite-deformation plasticity constitutive models. VFM was successfully implemented in MATLAB and evaluated using simulated FEM data that included representative experimental noise found in the Digital Image Cor- relation (DIC) optical technique that provides full-eld displacement measurements. VFM was able to identify constitutive model parameters for the BCJ plasticity model even in the presence of simulated DIC noise, demonstrating VFM as a viable alternative inverse method. Further research is required before VFM can be adopted as a standard method for constitu- tive model parameter identication, but this study is a foundation for ongoing research at Sandia for improving constitutive model calibration.

Kramer, Sharlotte Lorraine Bolyard↗
DOE OSTI2018

A mesh-independent method for planar three-dimensional crack growth in finite domains: A mesh-independent method for crack growth in 3D finite domains

The discrete crack mechanics (DCM) method is a dislocation-based crack modeling technique where cracks are constructed using Volterra dislocation loops. The method allows for the natural introduction of displacement discontinuities, avoiding numerically expensive techniques. Mesh dependence in existing computational modeling of crack growth is eliminated by utilizing a superposition procedure. The elastic field of cracks in finite bodies is separated into two parts: the infinite-medium solution of discrete dislocations and an finite element method solution of a correction problem that satisfies external boundary conditions. In the DCM, a crack is represented by a dislocation array with a fixed outer loop determining the crack tip position encompassing additional concentric loops free to expand or contract. Solving for the equilibrium positions of the inner loops gives the crack shape and stress field. The equation of motion governing the crack tip is developed for quasi-static growth problems. Additionally, convergence and accuracy of the DCM method are verified with two- and three-dimensional problems with well-known solutions. Crack growth is simulated under load and displacement (rotation) control. In the latter case, a semicircular surface crack in a bent prismatic beam is shown to change shape as it propagates inward, stopping as the imposed rotation is accommodated.

42 ENGINEERING↗
DOE OSTI2025

Adhesion of Self-Complementary, Sinusoidal Surfaces Fabricated Using Two-Photon Polymerization

Microscale, pick-and-place assembly is a non-lithographic assembly method poised to impact diverse fields including flexible electronics, microfluidics and robotics. However, a major technological challenge is the need to deterministically control adhesion between parts. Here, switchable adhesion involving 3D-printed, self-complementary surfaces is demonstrated. Mechanical properties of metasurfaces pressed against flat, rigid substrates are modeled using finite element methods. A series of flat slabs and metastructured slabs with 2D sinusoidal surfaces are printed using two-photon polymerization (2PP) of a shape-memory resin. The surface frequency of featured slabs was varied between $3.\bar3$ mm −1 and $26.\bar6$ mm −1 with similar amplitudes. Adhesion between printed metasurfaces and glass and between printed, self-complementary metasurfaces is studied above and below the cured resin’s glass transition temperature (∼45 °C). Simple heating of adhering surfaces to above 60 °C lowers adhesion, and compression of surfaces while above the glass transition temperature followed by cooling to room temperature elevates adhesion. The nominal adhesive strength between printed, self-complementary surfaces, as determined by the maximum observable pull-off stress, exceeds 3 MPa. Further tailoring complementary surfaces for adhesion control may facilitate microscale disassembly for recovery of components or precious metals.

amorphous materials↗
DOE OSTI2016

OpenACC directive-based GPU acceleration of an implicit reconstructed discontinuous Galerkin method for compressible flows on 3D unstructured grids

In this study, we use a combination of modeling techniques to describe the relationship between fracture radius that might be accomplished in a hypothetical enhanced geothermal system (EGS) and drilling distance required to create and access those fractures. We use a combination of commonly applied analytical solutions for heat transport in parallel fractures and 3D finite-element method models of more realistic heat extraction geometries. For a conceptual model involving multiple parallel fractures developed perpendicular to an inclined or horizontal borehole, calculations demonstrate that EGS will likely require very large fractures, of greater than 300 m radius, to keep interfracture drilling distances to ~10 km or less. As drilling distances are generally inversely proportional to the square of fracture radius, drilling costs quickly escalate as the fracture radius decreases. It is important to know, however, whether fracture spacing will be dictated by thermal or mechanical considerations, as the relationship between drilling distance and number of fractures is quite different in each case. Information about the likelihood of hydraulically creating very large fractures comes primarily from petroleum recovery industry data describing hydraulic fractures in shale. Those data suggest that fractures with radii on the order of several hundred meters may, indeed, be possible. The results of this study demonstrate that relatively simple calculations can be used to estimate primary design constraints on a system, particularly regarding the relationship between generated fracture radius and the total length of drilling needed in the fracture creation zone. Comparison of the numerical simulations of more realistic geometries than addressed in the analytical solutions suggest that simple proportionalities can readily be derived to relate a particular flow field.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

CU-BENs: A structural modeling finite element library

The present work discusses capabilities within the finite element library CU-BENs. CU-BENs focuses on applying the finite element method to structural mechanics problems encountered within the context of inverse problems and partitioned fluid–structure interaction; thus element formulations are primarily of a structural type — truss, frame, and triangular discrete Kirchhoff theory shells. CU-BENs defaults to the skyline sparse storage scheme for the system matrix, but also supports other storage schemes when using external libraries such as LAPACK and UMFPACK. CU-BENs includes built-in nonlinear solution strategies, such as the Newton Raphson method and modified spherical arc length method, that are available within static analyses as well within the context of transient dynamic analyses involving a generalized-α implementation of the Newmark implicit time integration scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Investigating the Interaction Between Circulating Tumor Cells and Local Hydrodynamics via Experiment and Simulations

Introduction: The biological and mechanical properties of circulating tumor cells (CTCs) in combination with the hemodynamics affect the preference of metastatic sites in the vasculature. Despite the extensive literature on the effects of biological properties on cell adhesion, the effects of hydrodynamic forces on primary attachment remains an active area of research. Here, using simulations in conjunction with experimentation, we provide new insight into the interplay of CTCs dynamics and local hydrodynamics. Methods: A flow experiment of CTC attachment was performed within a bioprinted, double branching endothelialized vessel. Simulations of fluid flow and CTC transport in the reconstructed and idealized bifurcated vessel were respectively performed by HARVEY, our in-house massively parallel computational fluid dynamics solver. HARVEY is based on the lattice Boltzmann and finite element methods to model the fluid and cells dynamics. The immersed boundary method is employed for resolving the fluid–structure interaction. Results: CTC attachment was quantified experimentally at all regions of the complex vessel. The results demonstrate a clear preference for CTCs to attach at the branch points. To elucidate the effect of the vessel topology on the location of attachment, a fluid-only simulation was performed assessing the differences in the hydrodynamics along the vessel. CTC transport in idealized bifurcated vessels was subsequently studied to examine the effects of cell deformability on the local hydrodynamics patterns and, thus, the preference of attachment sites. Conclusions: The current work provides evidence on the correlation of the hydrodynamics forces arising from the vessel topology and CTC properties on the attachment regions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A variational method for the sheath potential of hypersonic leading edges with space-charge limitations

Electron transpiration cooling for the leading edges (LE) of hypersonic aircraft utilizes thermionic emission; however, space-charge effects limit the electron emission rate, potentially diminishing the efficiency of this cooling mechanism. We develop a variational weak form of the Poisson equation that describes the sheath potential and then numerically solve it using the finite element method. This formulation has two main benefits: (1) the space-charge limit condition can be incorporated as a constraint and (2) it allows for the analysis of three-dimensional geometries with complex boundary conditions. We demonstrate that the current emitted from the surface of an LE is generally a small fraction of the Child–Langmuir limit due to space charge. We then propose several methods to enhance the emitted current from the surface and to boost the cooling effect of thermionic emission. These include increasing the plasma density, applying a negative surface potential, and using fringe fields under suitable geometric conditions. For a LaB6 emitting LE, the total emitted current is shown to be minimal and independent of the temperature of a surface with floating potential. However, when a negative potential is applied and the surface is heated, the emitted current follows the Richardson–Dushman relationship up to a critical temperature, beyond which it remains constant. At an applied surface potential of −5 V, the critical temperature is around 1700 K.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Non-perturbative Approach to Computing Seismic Normal Modes in Rotating Planets

In this work, a continuous Galerkin method based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial filtering eigensolver. The relevant elastic-gravitational system of equations, including the Coriolis force, is subjected to a mixed finite-element method, while self-gravitation is accounted for with the fast multipole method. Our discretization utilizes fully unstructured tetrahedral meshes for both solid and fluid regions. The relevant eigenvalue problem is solved by a combination of several highly parallel and computationally efficient methods. We validate our three-dimensional results in the non-rotating case using analytical results for constant elastic balls, as well as numerical results for an isotropic Earth model from standard “radial” algorithms. We also validate the computations in the rotating case, but only in the slowly-rotating regime where perturbation theory applies, because no other independent algorithms are available in the general case. The algorithm and code are used to compute the point spectra of eigenfrequencies in several Earth and Mars models studying the effects of heterogeneity on a large range of scales.

58 GEOSCIENCES↗
DOE OSTI2024

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗
DOE OSTI2013

Full Core Reactor Analysis: Running Denovo on Jaguar

Fully-consistent, full-core, 3D, deterministic neutron transport simulations using the orthogonal mesh code Denovo were run on the massively parallel computing architecture Jaguar XT5. Using energy and spatial parallelization schemes, Denovo was able to efficiently scale to more than 160k processors. Cell-homogenized cross sections were used with step-characteristics, linear-discontinuous finite element, and trilinear-discontinuous finite element spatial methods. It was determined that using the finite element methods gave considerably more accurate eigenvalue solutions for large-aspect ratio meshes than using step-characteristics.

Denovo↗
DOE OSTI2023

Investigation of desorption of hydrogen gas from polymer matrix using thermal desorption analysis and finite element modeling

Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A New Method to Measure Crack Extension in Nuclear Graphite Based on Digital Image Correlation

Graphite components, used as moderators, reflectors, and core-support structures in a High-Temperature Gas-Cooled Reactor, play an important role in the safety of the reactor. Specifically, they provide channels for the fuel elements, control rods, and coolant flow. Fracture is the main failure mode for graphite, and breaching of the above channels by crack extension will seriously threaten the safety of a reactor. In this paper, a new method based on digital image correlation (DIC) is introduced for measuring crack extension in brittle materials. Cross-correlation of the displacements measured by DIC with a step function was employed to identify the advancing crack tip in a graphite beam specimen under three-point bending. The load-crack extension curve, which is required for analyzing the R-curve and tension softening behaviors, was obtained for this material. Furthermore, a sensitivity analysis of the threshold value employed for the cross-correlation parameter in the crack identification process was conducted. Finally, the results were verified using the finite element method.

Lai, Shigang↗