Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

The shifted boundary method for embedded domain computations. Part I: Poisson and Stokes problems

We propose a new finite element method for embedded domain computations, which falls in the category of surrogate/approximate boundary algorithms. The key feature of the proposed approach is the idea of shifting the location where boundary conditions are applied from the true to the surrogate boundary, and to appropriately modify the shifted boundary conditions, enforced weakly, in order to preserve optimal convergence rates of the numerical solution. This process yields a method which, in our view, is simple, efficient, and also robust, since it is not affected by the small-cut-cell problem. Although general in nature, here we apply this new concept to the Poisson and Stokes problems. We present in particular the full analysis of stability and convergence for the case of the Poisson operator, and numerical tests for both the Poisson and Stokes equations, for geometries of progressively higher complexity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Simulation of the RPI Reactor Critical Facility Core Using High-Fidelity Neutronics Code PROTEUS

The PROTEUS code is a high-fidelity three-dimensional (3D) deterministic neutron transport code developed by ANL under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The code contains the SN2ND (second order discrete ordinates method) and two MOC transport solvers, MOCFE and MOCEX (first order methods of characteristics), based on unstructured finite element meshes which allow users to model complex or unconventional geometry reactor problems. Among the two MOC solvers which normally perform better than SN for heterogeneous geometry problems, MOCEX is based on a rigorous formulation using the two-dimensional (2D) MOC radially and the discontinuous Galerkin finite element method axially to overcome the drawbacks of MOCFE with the full 3D MOC representation which requires enormous memory and computational efforts. The Reactor Critical Facility (RCF), an open tank research reactor at the Rensselaer Polytechnic Institute (RPI), became operational in 1956 under the American Locomotive company and has been operating using low-enriched uranium fuel pins, which makes it the only university research facility to use fuel rods similar to operating commercial light water reactors (LWRs). This paper presents modeling and simulation results of the RCF using PROTEUS for code verification and validation which also supports an ongoing NEUP project in which a set of benchmark validation experiments for multiphysics coupling is being developed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020: Part 3

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: Task2c: Explicit modeling of pebble transient temperature response. In this simulation, the 400 MWth Pebble-Bed Modular Reactor (PBMR) design, PBMR- 400, experiences a 20-second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling.; Extended scope part 1: An assessment of the computational efficiency of the Discontinuous Finite Element Method (DFEM) heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the Finite Element Method (FEM) solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems.; Extended scope part 2: The GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables.; Extended scope part 3: Demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat-up transient and used in the GapHeatTransferInterface model. A full-domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench-marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1991

Parallel Methods on Large-Scale Structural Analysis and Physics Applications; Symposium, Hampton, VA, Feb. 5, 6, 1991, Selected Papers

Recent advances in parallel methods and algorithms integrated into large-scale codes are presented. Consideration is given to problem decomposition (substructuring), efficient matrix solution algorithms for shared memory architectures, dynamic and transient analysis algorithms for shared memory architectures, and algorithms for distributed and massively parallel architectures. Particular attention is given to partitioning of unstructured problems for parallel processing, parallel-vector computation for linear-structural analysis and nonlinear unconstraint optimization problems, a parallel-vector equation solver for unsymmetric matrices on supercomputers, parallel nonlinear finite element dynamic response, multigrid algorithms for solving structural mechanics problems on supercomputers, structural analysis on massively parallel computers, explicit finite element methods with contact-impact on SIMD computers, and the impact of mapping and sparsity on parallelized finite element method modules.

Storaasli, Olaf O.↗
DOE OSTI2019

Experimental and 3D Micromechanical Analysis of Stress–Strain Behavior and Damage Initiation in Dual-Phase Steels

In the recent decade, dual-phase steel has been used extensively in the industry due to its high strength and formability. In this paper, deformation pattern and mechanical behavior of DP steel are predicted using the 2D and 3D representative volume element (RVE) and finite elements method. Also, the obtained results were compared with the experimental SEM images. The experiments are performed in three stages to obtain the stress–strain curve and failure pattern in the specimens with the different state of stress distribution at the test sections. In the 3D state of stress, the state of applied stress to the RVE is obtained from macro-modeling of the test specimens and this model was analyzed using the finite element method. Furthermore, the effects of volume fraction and mesh size in 3D simulation on the stress–strain behavior are investigated. Metallography and SEM images are also used to access the failure mechanism at the microscale. It is shown that the martensite particles remain almost in their original form and their deformations are very small for the specimens with the thickness-to-width ratio of 0.8 (3D state of stress). This phenomenon is completely different from those obtained for smaller thickness-to-width ratio that the martensite particles were pulled and its length was increased. On the other hand, SEM images show that voids are located in the phase boundaries. According to the experimental and numerical results, although the main factor in the formation and growth of 2D voids is shear stress, in the 3D state of stress, the voids grow in spherical form due to the existence of large hydrostatic stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deep Learning-Assisted Real-Time Forward Modeling of Electromagnetic Logging in Complex Formations

Higher dimensional (i.e., 2-D and 3-D) modeling is indispensable to correctly evaluate the responses of electromagnetic logging tools in complex formation environments. However, limited by the high computational cost of rigorous modeling such as finite difference method and the finite element method, the real-time applications in the well logging industry primarily rely on the 1-D forward solver, which would result in erroneous formation evaluation for complex scenarios. As a result, aiming at realizing fast modeling for electromagnetic logging tools in complex formations, this paper proposes a general framework assisted by deep neural networks (DNNs). The framework consists of three modules: earth model classification, parameter extraction, and surrogate construction. Separate DNNs are trained and tested for different modules. The accuracy and efficiency of the DNN assisted fast modeling are validated by several experiments. Here, this study finds that the fast modeling assisted by DNNs is able to calculate the tool responses and reconstruct the subsurface formations in real-time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Sierra/RKPM: Current State and Applications

The Reproducing Kernel Particle Method (RKPM), a meshfree method, has been implemented in Sandia's Sierra/Solid Mechanics in a collaboration between Sandia and the University of California San Diego's Center for Extreme Events Research (UCSD/CEER). Meshfree methods, like RKPM, have an advantage over mesh-based methods, like the Finite Element Method (FEM), in applications where achieving or maintaining a quality mesh becomes burdensome or impractical. For example, using FEM for problems with very large deformations will result in poor element Jacobians which causes problems with the parametric mapping. RKPM constructs the approximation functions in the physical domain, circumventing the parametric mapping issue. Also, reconstructing the approximation functions at very large deformations is straight-forward. RKPM has an advantage over traditional meshfree methods such as Smoothed-Particle Hydrodynamics (SPH) due to its ability to reproduce linear or higher-order functions exactly. This removes the tensile instabilities that are present in SPH and allows preservation of angular momentum. The point of this memo is to explore the capabilities and limitations of the current implementation by testing it on three different applications: 1) a quasi-static ductile shearing problem 2) a dynamic concrete panel perforation problem and 3) a set of dynamic metal panel perforation problems. In summary, areas where RKPM appears to be a promising alternative to current methods have been identified. Also, outstanding inefficiencies and issues (bugs) with code are noted, ways to improve the capabilities using material from literature are mentioned and areas deserving of new research are highlighted.

42 ENGINEERING↗
DOE OSTI2019

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1981

The stability of pseudospectral-Chebyshev methods

The pseudospectral-Chebyshev methods are shown to be convergent in variable coefficient problems and, in some cases, hyperbolic problems. The analysis demonstrates that the rate of convergence is greater for finite difference methods or the finite element method. For a single first-order hyperbolic equation, the method is seen as remaining stable even when the coefficient changes sign, although in this case it is specified that care must be taken to have adequate spatial resolution. It is noted that this fact, combined with the fact that collocation methods are easy to apply in the nonlinear case, shows that the pseudospectral method is in general preferable to the Galerkin or Tau methods.

Gottlieb, D.↗
NASA NTRS1984

Megamechanics research Consortium

The mechanical behavior of large space lattice systems is investigated. Continuum modeling methods such as the Load Correction Method (LCM) and finite element operations aid in solving structural analysis problems. Static analysis of frames and trusses, dynamic analysis of trusses, analysis of thermal and mechanical loads on frames and trusses, and load corrections for three dimensional complex geometries is possible by using the LCM. The finite element method and the LCM are used to extend Saint-Venant's principle via the Boussinesq problem into the regime of large repetitive lattice systems.

Broome, T. H., Jr.↗
NASA NTRS1989

An advanced probabilistic structural analysis method for implicit performance functions

In probabilistic structural analysis, the performance or response functions usually are implicitly defined and must be solved by numerical analysis methods such as finite element methods. In such cases, the most commonly used probabilistic analysis tool is the mean-based, second-moment method which provides only the first two statistical moments. This paper presents a generalized advanced mean value (AMV) method which is capable of establishing the distributions to provide additional information for reliability design. The method requires slightly more computations than the second-moment method but is highly efficient relative to the other alternative methods. In particular, the examples show that the AMV method can be used to solve problems involving non-monotonic functions that result in truncated distributions.

Wu, Y.-T.↗
NASA NTRS1991

A new procedure for calculating contact stresses in gear teeth

A numerical procedure for evaluating and monitoring contact stresses in meshing gear teeth is discussed. The procedure is intended to extend the range of applicability and to improve the accuracy of gear contact stress analysis. The procedure is based upon fundamental solution from the theory of elasticity. It is an iterative numerical procedure. The method is believed to have distinct advantages over the classical Hertz method, the finite-element method, and over existing approaches with the boundary element method. Unlike many classical contact stress analyses, friction effects and sliding are included. Slipping and sticking in the contact region are studied. Several examples are discussed. The results are in agreement with classical results. Applications are presented for spur gears.

Somprakit, Paisan↗
NASA NTRS1991

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches; the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(sup eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(sup eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗
NASA NTRS1990

Monitoring of contact stresses in advanced propulsion systems

This paper discusses a new and innovative numerical procedure for evaluating and monitoring contact stresses in the mating surfaces of propulsion systems. The procedure is intended to extend the range of applicability and the accuracy of contact stress analyses. Originally proposed by Burton Paul, et. al., (1981) the procedure is based upon fundamental solutions from the theory of elasticity. It is an iterative numerical procedure. It is believed that the method has distinct advantages over the classical Hertz method, the finite-element method and over existing approaches with the boundary element method. Several examples are discussed. Favorable comparisons are obtained with classical results. Applications with mating element of propulsion systems are discussed.

Somprakit, Paisan↗
NASA NTRS1992

Modeling of crack bridging in a unidirectional metal matrix composite

The effective fatigue crack driving force and crack opening profiles were determined analytically for fatigue tested unidirectional composite specimens exhibiting fiber bridging. The crack closure pressure due to bridging was modeled using two approaches: the fiber pressure model and the shear lag model. For both closure models, the Bueckner weight function method and the finite element method were used to calculate crack opening displacements and the crack driving force. The predicted near crack tip opening profile agreed well with the experimentally measured profiles for single edge notch SCS-6/Ti-15-3 metal matrix composite specimens. The numerically determined effective crack driving force, Delta K(eff), was calculated using both models to correlate the measure crack growth rate in the composite. The calculated Delta K(eff) from both models accounted for the crack bridging by showing a good agreement between the measured fatigue crack growth rates of the bridged composite and that of unreinforced, unbridged titanium matrix alloy specimens.

Ghosn, Louis J.↗
NASA NTRS1993

Contact stress analysis of spiral bevel gears using nonlinear finite element static analysis

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Bibel, G. D.↗
NASA NTRS1990

Theoretical and numerical difficulties in 3-D vector potential methods in finite element magnetostatic computations

This paper describes the results of application of three well known 3D magnetic vector potential (MVP) based finite element formulations for computation of magnetostatic fields in electrical devices. The three methods were identically applied to three practical examples, the first of which contains only one medium (free space), while the second and third examples contained a mix of free space and iron. The first of these methods is based on the unconstrained curl-curl of the MVP, while the second and third methods are predicated upon constraining the divergence of the MVP 10 zero (Coulomb's Gauge). It was found that the two latter methods cease to give useful and meaningful results when the global solution region contains a mix of media of high and low permeabilities. Furthermore, it was found that their results do not achieve the intended zero constraint on the divergence of the MVP.

Demerdash, N. A.↗