Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixed 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Compatibility Condition in Theory of Solid Mechanics (Elasticity, Structures, and Design Optimization)

The strain formulation in elasticity and the compatibility condition in structural mechanics have neither been understood nor have they been utilized. This shortcoming prevented the formulation of a direct method to calculate stress. We have researched and understood the compatibility condition for linear problems in elasticity and in finite element analysis. This has lead to the completion of the method of force with stress (or stress resultant) as the primary unknown. The method in elasticity is referred to as the completed Beltrami-Michell formulation (CBMF), and it is the integrated force method (IFM) in structures. The dual integrated force method (IFMD) with displacement as the primary unknown has been formulated. IFM and IFMD produce identical responses. The variational derivation of the CBMF yielded the new boundary compatibility conditions. The CBMF can be used to solve stress, displacement, and mixed boundary value problems. The IFM in structures produced high-fidelity response even with a modest finite element model. The IFM has influenced structural design considerably. A fully utilized design method for strength and stiffness limitation has been developed. The singularity condition in optimization has been identified. The CBMF and IFM tensorial approaches are robust formulations because of simultaneous emphasis on the equilibrium equation and the compatibility condition.

Patnaik, Surya N.↗
DOE OSTI2017

A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields

In this paper, we propose an alternative method to generate samples of a spatially correlated random field with applications to large-scale problems for forward propagation of uncertainty. A classical approach for generating these samples is the Karhunen--Loève (KL) decomposition. However, the KL expansion requires solving a dense eigenvalue problem and is therefore computationally infeasible for large-scale problems. Sampling methods based on stochastic partial differential equations provide a highly scalable way to sample Gaussian fields, but the resulting parametrization is mesh dependent. We propose a multilevel decomposition of the stochastic field to allow for scalable, hierarchical sampling based on solving a mixed finite element formulation of a stochastic reaction-diffusion equation with a random, white noise source function. Lastly, numerical experiments are presented to demonstrate the scalability of the sampling method as well as numerical results of multilevel Monte Carlo simulations for a subsurface porous media flow application using the proposed sampling method.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Eulerian-Lagrangian Simulations of Transonic Flutter Instabilities

This paper presents an overview of recent applications of Eulerian-Lagrangian computational schemes in simulating transonic flutter instabilities. This approach, the fluid-structure system is treated as a single continuum dynamics problem, by switching from an Eulerian to a Lagrangian formulation at the fluid-structure boundary. This computational approach effectively eliminates the phase integration errors associated with previous methods, where the fluid and structure are integrated sequentially using different schemes. The formulation is based on Hamilton's Principle in mixed coordinates, and both finite volume and finite element discretization schemes are considered. Results from numerical simulations of transonic flutter instabilities are presented for isolated wings, thin panels, and turbomachinery blades. The results suggest that the method is capable of reproducing the energy exchange between the fluid and the structure with significantly less error than existing methods. Localized flutter modes and panel flutter modes involving traveling waves can also be simulated effectively with no a priori knowledge of the type of instability involved.

Bendiksen, Oddvar O.↗
NASA NTRS1981

Analysis of cracks emanating from a circular hole in unidirectional fiber reinforced composites, part 2

An analytical method is developed for cracks emanating from a circular hole in an off-axis unidirectional fiber-reinforced composite. The method which is formulated by using conservation laws of elasticity and fundamental relationships in anisotropic fracture mechanics, provides a convenient and accurate means to examine the complicated crack behavior, when used in conjunction with a suitable numerical scheme such as the finite element method. The formulation is eventually reduced to a system of linear algebraic equations of mixed-mode stress intensity factors. Fracture parameters, describing crack-tip deformation and fracture in the composite, are obtained explicitly. Effects of material anisotropy and crack/hole geometry are examined also. Of particular interest are the energy release rates associated with crack extension; their values are evaluated for various cases. Results show that mixed-mode stress intensity factors and energy release rates associated with the cracks emanating from a hole change very appreciably with fiber orientation in the composite. K sub 1 and G increase monotonically with increasing theta; but K sub 2 reaches its maximum at theta = 45 deg, and then decreases gradually as theta increases further.

Wang, S. S.↗
DOE OSTI2023

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗
DOE OSTI2015

Highly non-linear solid core photonic crystal fiber with one nano hole

The numerical study of newly designed solid core photonic crystal fiber (SCPCF) having three hexagonal air hole rings in cladding region and one small nano hole at the center are presented. By using full vectorial finite element method (FV-FEM), we analyses the optical properties like effective area, nonlinearity and confinement loss of the proposed PCF. Results show that the change in core diameter controls the effective area, nonlinearity and confinement loss. A low effective area (3.34 µm{sup 2}), high nonlinearity (36.34 W{sup −1}km{sup −1}) and low confinement loss (0.00106 dB/km) are achieved at the communication wavelength 1.55 µm for the SCPCF having core air hole diameter 0.10 µm, cladding air holes diameter 1.00 µm and pitch 2.50 µm. This type of PCF is very useful in non-linear applications such as supercontinuum generation, four wave mixing, second harmonic generation etc.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Mixed time integration methods for transient thermal analysis of structures, appendix 5

Mixed time integration methods for transient thermal analysis of structures are studied. An efficient solution procedure for predicting the thermal behavior of aerospace vehicle structures was developed. A 2D finite element computer program incorporating these methodologies is being implemented. The performance of these mixed time finite element algorithms can then be evaluated employing the proposed example problem.

Liu, W. K.↗
NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗
DOE OSTI2015

Steady-state benchmarks of DK4D: A time-dependent, axisymmetric drift-kinetic equation solver

The DK4D code has been written to solve a set of time-dependent, axisymmetric, finite-Larmor-radius drift-kinetic equations (DKEs) for the non-Maxwellian part of the electron and ion distribution functions using the full, linearized Fokker–Planck–Landau collision operator. The plasma is assumed to be in the low- to finite-collisionality regime, as is found in the cores of modern and future magnetic confinement fusion experiments. Each DKE is formulated such that the perturbed distribution function carries no net density, parallel momentum, or kinetic energy. Rather, these quantities are contained within the background Maxwellians and would be evolved by an appropriate set of extended magnetohydrodynamic (MHD) equations. This formulation allows for straight-forward coupling of DK4D to existing extended MHD time evolution codes. DK4D uses a mix of implicit and explicit temporal representations and finite element and spectral spatial representations. These, along with other computational methods used, are discussed extensively. Steady-state benchmarks are then presented comparing the results of DK4D to expected analytic results at low collisionality, qualitatively, and to the Sauter analytic fits for the neoclassical conductivity and bootstrap current, quantitatively. These benchmarks confirm that DK4D is capable of solving for the correct, gyroaveraged distribution function in stationary magnetic equilibria. Furthermore, the results presented demonstrate how the exact drift-kinetic solution varies with collisionality as a function of the magnetic moment and the poloidal angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

An Optimal Order Nonnested Mixed Multigrid Method for Generalized Stokes Problems

A multigrid algorithm is developed and analyzed for generalized Stokes problems discretized by various nonnested mixed finite elements within a unified framework. It is abstractly proved by an element-independent analysis that the multigrid algorithm converges with an optimal order if there exists a 'good' prolongation operator. A technique to construct a 'good' prolongation operator for nonnested multilevel finite element spaces is proposed. Its basic idea is to introduce a sequence of auxiliary nested multilevel finite element spaces and define a prolongation operator as a composite operator of two single grid level operators. This makes not only the construction of a prolongation operator much easier (the final explicit forms of such prolongation operators are fairly simple), but the verification of the approximate properties for prolongation operators is also simplified. Finally, as an application, the framework and technique is applied to seven typical nonnested mixed finite elements.

Deng, Qingping↗
NASA NTRS1991

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗
DOE OSTI2014

Simulating the spectral response of quantum dot-in-well infrared photodetectors from eight band k.p method

Conduction band energy levels in quantum-dot-in-a-well structures are computed by eight band k.p method (Burt-Foreman Hamiltonian) using finite element software. Optical absorption spectrum due to intersubband transitions is simulated using Fermi golden rule. The use of contact pair boundary condition in strain calculation and criteria for choosing band mixing parameter (E{sub p}) to avoid the spurious solutions are examined in this paper. The simulated intersubband optical absorption spectrum of different structures reported in the literature is in close agreement with the experimentally measured photoconductive absorption region and shows that the method can be used as an effective modeling for quick design of the heterostructures based infrared photodetectors for various wavelengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Development of an Efficient Meso- scale Multi-phase Flow Solver in Nuclear Applications

The proposed research aims at formulating a predictive high-order Lattice Boltzmann Equation for multi-phase flows relevant to nuclear energy related application - namely, saturated and sub-cooled boiling in reactors, and liquid- liquid mixing and extraction for fuel cycle separation. An efficient flow solver will be developed based on the Finite Element based Lattice Boltzmann Method (FE- LBM), accounting for phase-change heat transfer and capable of treating multiple phases over length scales from the submicron to the meter. A thermal LBM will be developed in order to handle adjustable Prandtl number, arbitrary specific heat ratio, a wide range of temperature variations, better numerical stability during liquid-vapor phase change, and full thermo-hydrodynamic consistency. Two-phase FE-LBM will be extended to liquid–liquid–gas multi-phase flows for application to high-fidelity simulations building up from the meso-scale up to the equipment sub-component scale. While several relevant applications exist, the initial applications for demonstration of the efficient methods to be developed as part of this project include numerical investigations of Critical Heat Flux (CHF) phenomena in nuclear reactor fuel bundles, and liquid-liquid mixing and interfacial area generation for liquid-liquid separations. In addition, targeted experiments will be conducted for validation of this advanced multi-phase model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2016

CYRANO3 the EDF fuel code performance: Global overview and recent developments on MOX fuel

CYRANO3 is the thermal mechanical industrial code developed and used by EDF to simulate the nuclear fuel rod performance in pile, and during the disposal storage period as well. This code has already been used successfully for the last ten years by the EDF engineering teams in justifying the overall reloading options covering various fuels: UO{sub 2}, UO{sub 2} + gadolinium, MOX, UO{sub 2} plus additives, and various cladding as well proposed by the nuclear fuel suppliers. Its particularity is related to a construction around a standard one-dimensional thermal and mechanical finite elements kernel with an appropriate architecture to plug in any new functionalities or models. Since the beginning it is developed as a module able integrating within a core calculation chain. Other constrain is to always guaranty a running time consistent with the calculation of a core within a reasonable time. The objective is to provide the EDF designers with a convivial and qualified code, with robust algorithms and able to quickly evolve with computer performances. Algorithms are able to account for adiabatic transients and are robust enough to avoid any undesirable overflows. The permanent challenge is indeed to provide with the more descriptive simulations not degrading to much the running time. Developments are obviously performed within the strict Quality Insurance rules. The code is then associated to a large documentation base, including model specification, realization, validation and test procedures. Maintenance is managed through a system of anomaly and improvement forms. Computers are still evolving quickly. Therefore, after an architecture adaptation, the CYRANO3 code is now integrated in the PLEIADES platform as a PWR industrial application. The platform is an exchange space allowing sharing functionalities, models and properties. In this paper, we will give a global overview of the capabilities of CYRANO3 code to simulate the thermomechanical behavior of fuel rods, and we will focus on recent developments realized on MOX fuel. In CYRANO3 code, heterogeneous models (neutronic model, micromechanical behavior law and fission gas model) are developed for MOX fuels and integrated in order to account for differential evolutions between phases under irradiation. The methodology and these models are presented in this paper. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

BISON: A Flexible Code for Advanced Simulation of the Performance of Multiple Nuclear Fuel Forms

BISON is a nuclear fuel performance application built using the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library. One of its major goals is to have a great amount of flexibility in how it is used, including in the types of fuel it can analyze, the geometry of the fuel being modeled, the modeling approach employed, and the dimensionality and size of the models. Fuel forms that can be modeled include standard light water reactor fuel, emerging light water reactor fuels, tri-structural isotropic fuel particles, and metallic fuels. BISON is a platform for research in nuclear fuel performance modeling while simultaneously serving as a tool for the analysis of nuclear fuel designs. Recent research in BISON includes techniques such as the extended finite element method for fuel cracking, exploration of high-burnup light water reactor fuel behavior, swelling behavior of metallic fuels, and central void formation in mixed-oxide fuel. BISON includes integrated documentation for each of its capabilities, follows rigorous software quality assurance procedures, and has a growing set of rigorous verification and validation tests.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

FY 2016 Status Report: Documentation of All CIRFT Data including Hydride Reorientation Tests (Draft M2)

The first portion of this report provides a detailed description of fiscal year (FY) 2015 test result corrections and analysis updates based on FY 2016 updates to the Cyclic Integrated Reversible-Bending Fatigue Tester (CIRFT) program methodology, which is used to evaluate the vibration integrity of spent nuclear fuel (SNF) under normal conditions of transport (NCT). The CIRFT consists of a U-frame test setup and a real-time curvature measurement method. The three-component U-frame setup of the CIRFT has two rigid arms and linkages connecting to a universal testing machine. The curvature SNF rod bending is obtained through a three-point deflection measurement method. Three linear variable differential transformers (LVDTs) are clamped to the side connecting plates of the U-frame and used to capture deformation of the rod. The second portion of this report provides the latest CIRFT data, including data for the hydride reorientation test. The variations in fatigue life are provided in terms of moment, equivalent stress, curvature, and equivalent strain for the tested SNFs. The equivalent stress plot collapsed the data points from all of the SNF samples into a single zone. A detailed examination revealed that, at the same stress level, fatigue lives display a descending order as follows: H. B. Robinson Nuclear Power Station (HBR), LMK, and mixed uranium-plutonium oxide (MOX). Just looking at the strain, LMK fuel has a slightly longer fatigue life than HBR fuel, but the difference is subtle. The third portion of this report provides finite element analysis (FEA) dynamic deformation simulation of SNF assemblies . In a horizontal layout under NCT, the fuel assembly’s skeleton, which is formed by guide tubes and spacer grids, is the primary load bearing apparatus carrying and transferring vibration loads within an SNF assembly. These vibration loads include interaction forces between the SNF assembly and the canister basket walls. Therefore, the integrity of the guide tubes and spacer grids critically affects the vibration intensity of the fuel assembly during transport and must be considered when developing the multipurpose purpose canister (MPC) design for safe SNF transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗