Finite-dimensional models for equations of mathematical physics
Finite dimensional models for obtaining finite difference methods with same properties
FINITE DIFFERENCE METHOD↗
SEARCH · Engineering Papers
Results for “Finite”
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.
Asynchronous finite state models of manual control systems
Asynchronous finite-state models of manual control systems in which human operator behaves as finite-state machine
HUMAN PERFORMANCE↗
Symmetric and nonsymmetric buckling of finitely deformed eccentrically stiffened shells of revolution.
Finite difference method stability analysis of deformed eccentrically stiffened shells of revolution, accounting for finite prebuckling rotations
FINITE DIFFERENCE METHOD↗
Symmetric and nonsymmetric buckling of finitely deformed eccentrically stiffened shells of revolution.
Finite difference method stability analysis of deformed eccentrically stiffened shells of revolution, accounting for finite prebuckling rotations
FINITE DIFFERENCE METHOD↗
Large Scale Finite Element Modeling Using Scalable Parallel Processing
An iterative solver for use with finite element codes was developed for the Cray T3D massively parallel processor at the Jet Propulsion Laboratory. Finite element modeling is useful for simulating scattered or radiated electromagnetic fields from complex three-dimensional objects with geometry variations smaller than an electrical wavelength.
finite element modeling parallel processing iterat↗
Azimuthally-dependent Finite Element Solution to the Cylindrical Resonator
The cylindrical cavity resonator loaded with an anisotropic dielectric is analyzed as a two-dimensional problem using a finite element approach that assumes sinusoidal dependence in azimuth. This methodology allows the first finite element treatment of the technically important case of a resonator containing a sapphire element with a cylindrically aligned c axis. Second order trial functions together with quadrilateral elements are adopted in the calculations. The method was validated through comparisons with the analytical solutions for the hollow metal cavity and a coaxial cavity, as well as through measurements on a shielded sapphire resonator.
cylindrical resonator two-dimensional finite eleme↗
Real-time simulation of dynamically deformable finite element models using modal analysis and spectral Lanczos decomposition methods
We present two efficient methods for simulating real-time behavior of a dynamically deformable 3D object modeled by finite element equations. The first method is based on modal analysis, which uses the most significant vibration modes of the object to compute the deformations in real-time for applied forces. The second method uses the spectral Lanczos decomposition to obtain the explicit solutions of the finite element equations.
medical↗
The Discontinuous Galerkin Finite Element Method
The purpose of this report is to present details of the Discontinuous Galerkin (DG) Finite Element Method (DG FEM). First the weighted-residual (WR) form is introduced and then the Galerkin Finite Element (FE) (GFE) and the Petrov-Galerkin FE (PG FE) methods are discussed. The details of the implementation of the DG FEM are presented along with two examples, 2nd order and 4th order differential equations, and the performance of the method is discussed.
Finite element↗
Finite Element Modeling for Compression Strength Prediction of Honeycomb Cores with Geometric Imperfections Measured using X-ray CT Imaging
Aluminum honeycomb cores have been used extensively in composite sandwich panels due to having high bending rigidity while maintaining low density. During manufacturing of the aluminum honeycomb, the thin metal walls are susceptible to imperfections that deviate from an ideal honeycomb shape. Computational tools to quantify imperfections and investigate their effects on quasi-static compression and impact response can inform design criteria to construct safer and lighter launch vehicle structures. This paper presents the results of a study on finite element modeling of metallic honeycomb cores (HCC) with geometric imperfections measured using X-ray Computer Tomography (CT), to predict the compression response. Finite element (FE) models are constructed with measured imperfections with appropriate boundary conditions. It is shown that the average behavior of larger 10x11 cell models can be predicted by sampling single cell models from the larger domain. Amplification of imperfections for cells at the center of sampled specimens are mode switched from the larger out-of-plane imperfection to the mode shape of the adjacent cells resulting in strengthening of the center specimens.
aluminum honeycomb↗
Finite Element Modeling of Extravehicular Mobility Units for Use With Human Body Models – Motivation, Major Challenges, Use Cases and Preliminary Work
Finite element modeling of pressurized spacesuits and implementation with human body models offers many advantages over physical experimentation, but also presents significant challenges. A model of pressurized spacesuit softgoods was developed, integrated with an existing hardgoods model, and fitted to a human body model. Two representative loading scenarios were simulated: dynamic external suit loading and internal occupant-driven loading, to serve as proof-of-concept for the modeling techniques employed. The modeled interactions behaved as intended and illustrate that finite element models of pressurized spacesuits can be used in conjunction with human body models to assess the biomechanical behavior. Further model development and experimental validation are needed.
Finite Element↗
Finite Element Modeling of Extravehicular Mobility Units for Use With Human Body Models – Motivation, Major Challenges, Use Cases and Preliminary Work
Finite element modeling of pressurized spacesuits and implementation with human body models offers many advantages over physical experimentation, but also presents significant challenges. A model of pressurized spacesuit softgoods was developed, integrated with an existing hardgoods model, and fitted to a human body model. Two representative loading scenarios were simulated: dynamic external suit loading and internal occupant-driven loading, to serve as proof-of-concept for the modeling techniques employed. The modeled interactions behaved as intended and illustrate that finite element models of pressurized spacesuits can be used in conjunction with human body models to assess the biomechanical behavior. Further model development and experimental validation are needed.
Finite Element↗
A Full-wave Model for Wave Propagation and Dissipation in the Inner Magnetosphere Using the Finite Element Method
A wide variety of plasma waves play an important role in the energization and loss of particles in the inner magnetosphere. Our ability to understand and model wave-particle interactions in this region requires improved knowledge of the spatial distribution and properties of these waves as well as improved understanding of how the waves depend on changes in solar wind forcing and/or geomagnetic activity. To this end, we have developed a two-dimensional, finite element code that solves the full wave equations in global magnetospheric geometry. The code describes three-dimensional wave structure including mode conversion when ULF, EMIC, and whistler waves are launched in a two-dimensional axisymmetric background plasma with general magnetic field topology. We illustrate the capabilities of the code by examining the role of plasmaspheric plumes on magnetosonic wave propagation; mode conversion at the ion-ion and Alfven resonances resulting from external, solar wind compressions; and wave structure and mode conversion of electromagnetic ion cyclotron waves launched in the equatorial magnetosphere, which propagate along the magnetic field lines toward the ionosphere. We also discuss advantages of the finite element method for resolving resonant structures, and how the model may be adapted to include nonlocal kinetic effects.
43 PARTICLE ACCELERATORS↗
A simple finite element method for non-divergence form elliptic equation
Here, we develop a simple finite element method for solving second order elliptic equations in non-divergence form by combining least squares concept with discontinuous approximations. This simple method has a symmetric and positive definite system and can be easily analyzed and implemented. We could have also used general meshes with polytopal element and hanging node in the method. We prove that our finite element solution approaches to the true solution when the mesh size approaches to zero. Numerical examples are tested that demonstrate the robustness and flexibility of the method.
97 MATHEMATICS AND COMPUTING↗
COMPARISONS OF THE FINITE-ELEMENT-WITH-DISCONTIGUOUS-SUPPORT METHOD TO CONTINUOUS-ENERGY MONTE CARLO FOR PIN-CELL PROBLEMS
The standard multigroup (MG) method for energy discretization of the transport equation can be sensitive to approximations in the weighting spectrum chosen for cross-section averaging. As a result, MG often inaccurately treats important phenomena such as self-shielding variations across a material. From a finite-element viewpoint, MG uses a single fixed basis function (the pre-selected spectrum) within each group, with no mechanism to adapt to local solution behavior. In this work, we introduce the Finite-Element-with-Discontiguous-Support (FEDS) method, whose only approximation with respect to energy is that the angular flux is a linear combination of unknowns multiplied by basis functions. A basis function is non-zero only in the discontiguous set of energy intervals associated with its energy element. Discontiguous energy elements are generalizations of bands and are determined by minimizing a norm of the difference between snapshot spectra and their averages over the energy elements. We begin by presenting the theory of the FEDS method. We then compare to continuous-energy Monte Carlo for one-dimensional slab and two-dimensional pin-cell problem. We find FEDS to be accurate and efficient at producing quantities of interest such as reaction rates and eigenvalues. Results show that FEDS converges at a rate that is approximately first-order in the number of energy elements and that FEDS is less sensitive to weighting spectrum than standard MG.
Cross sections↗
A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS
This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.
36 MATERIALS SCIENCE↗
Extended Finite Element Based Approach in Additive Manufacturing Modeling for Optimizing Highly Complex Manifold in Protonic Ceramic Electrochemical Cells
The objective of this project is to develop a novel numerical approach based on the extended finite element method (XFEM) for modeling moving boundaries of material deposited in additive manufacturing (AM) process with improved accuracy and reduced computational cost. One of the major challenges in simulating AM processes is that the boundaries of the computational domain need to change as material is deposited. Previous approaches have achieved this by activating new finite elements, but this requires a high level of mesh refinement to capture small movement of the boundary. XFEM allows the solution boundary to move independently of the mesh, permitting smooth representation of the evolution of the boundary as material is deposited. As a result, this new method will significantly improve the simulation accuracy while reducing the computational cost compared with existing approaches. Furthermore, this project will bring significant improvements for the design and AM process optimization in an iterative fashion among numerical simuation, parameter optimization, and experimental validation. Specifically, this project will support the protonic ceramic electrochemical cells (PCEC) stack development at INL by providing insights into the density, residual stress, thermal and mechanical properties of the PCEC manifold and interconnect, which is critical in enhancing the system lifetime and reducing the overall cost. Upon the success of proposed development and validation of the proposed approach, simulation will be applied to determine an optimal set of AM process parameters for the PCEC manifold production. This project will consolidate Idaho National Laboratory (INL)’s simulation capabilities provided by the MOOSE framework and the Valhalla AM simulation application, with encouraging expansion to emerging PCEC applications and AM technologies at INL and beyond.
97 MATHEMATICS AND COMPUTING↗
GPU-accelerated 3D neutron diffusion code based on finite difference method
Finite difference method, as a traditional numerical solution to neutron diffusion equation, although considered simpler and more precise than the coarse mesh nodal methods, has a bottle neck to be widely applied caused by the huge memory and unendurable computation time it requires. In recent years, the concept of General-Purpose computation on GPUs has provided us with a powerful computational engine for scientific research. In this study, a GPU-Accelerated multi-group 3D neutron diffusion code based on finite difference method was developed. First, a clean-sheet neutron diffusion code (3DFD-CPU) was written in C++ on the CPU architecture, and later ported to GPUs under NVIDIA's CUDA platform (3DFD-GPU). The IAEA 3D PWR benchmark problem was calculated in the numerical test, where three different codes, including the original CPU-based sequential code, the HYPRE (High Performance Pre-conditioners)-based diffusion code and CITATION, were used as counterpoints to test the efficiency and accuracy of the GPU-based program. The results demonstrate both high efficiency and adequate accuracy of the GPU implementation for neutron diffusion equation. A speedup factor of about 46 times was obtained, using NVIDIA's Geforce GTX470 GPU card against a 2.50 GHz Intel Quad Q9300 CPU processor. Compared with the HYPRE-based code performing in parallel on an 8-core tower server, the speedup of about 2 still could be observed. More encouragingly, without any mathematical acceleration technology, the GPU implementation ran about 5 times faster than CITATION which was speeded up by using the SOR method and Chebyshev extrapolation technique. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Nonlinear acceleration of a continuous finite element discretization of the self-adjoint angular flux form of the transport equation
Nonlinear acceleration of a continuous finite element (CFE) discretization of the transport equation requires a modification of the transport solution in order to achieve local conservation, a condition used in nonlinear acceleration to define the stopping criterion. In this work we implement a coarse-mesh finite difference acceleration for a CFE discretization of the second-order self adjoint angular flux (SAAF) form of the transport equation and use a post processing to enforce local conservation. Numerical results are given for one-group source calculations of one-dimensional slabs. We also give a formal derivation of the boundary conditions for the SAAF. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗