A hybridized formulation for the weak Galerkin mixed finite element method
Not Available
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “finite element problems”
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.
An a posteriori error estimator for the weak Galerkin least-squares finite-element method
Not Available
97 MATHEMATICS AND COMPUTING↗
A velocity/stress mixed stabilized nodal finite element for elastodynamics: Analysis and computations with strongly and weakly enforced boundary conditions
A new nodal mixed finite element is proposed for the simulation of linear elastodynamics and wave propagation problems in time domain. Our method is based on equal-order interpolation discrete spaces for both the velocity (or displacement) and stress (or strain) tensor variables. The mixed form is derived using either the velocity/stress or velocity/strain pair of unknowns, the latter being instrumental in extensions of the method to nonlinear mechanics. The proposed approach works equally well on hexahedral or tetrahedral grids and, for this reason, it is suitable for time-domain engineering applications in complex geometry. The peculiarity of the proposed approach is the use of the rate form of the stress update equation, which yields a set of governing equations with the structure of a non-dissipative space/time Friedrichs’ system. We complement standard traction boundary conditions for the stress with strongly and weakly enforced boundary conditions for the velocity (or displacement). Weakly enforced boundary conditions are particularly suitable when considering complex geometrical shapes, because they do not require dedicated data structures for the imposition of the boundary degrees of freedom, but, rather, they utilize the structure of the variational formulation. Here we also show how the framework of weakly enforced boundary conditions can be used to develop variational forms for multi-domain simulations of heterogeneous media. A complete analysis including stability and convergence proofs is included, in the case of a space–time variational approach. A series of computational tests are used to demonstrate and verify the performance of the proposed approach.
42 ENGINEERING↗
Carbon nanotube biconvex microcavities
Developing highly efficient microcavities with predictive narrow-band resonance frequencies using the least amount of material will allow the applications in nonlinear photonic devices. We have developed a microcavity array that comprised multi-walled carbon nanotubes (MWCNT) organized in a biconvex pattern. The finite element model allowed designing microcavity arrays with predictive transmission properties and assessing the effects of the microarray geometry. The microcavity array demonstrated negative index and produced high Q factors. 2–3 μm tall MWCNTs were patterned as biconvex microcavities, which were separated by 10 μm in an array. The microcavity was iridescent and had optical control over the diffracted elliptical patterns with a far-field pattern, whose properties were predicted by the model. It is anticipated that the MWCNT biconvex microcavities will have implications for the development of highly efficient lenses, metamaterial antennas, and photonic circuits.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Three-point bending analysis of doubly clamped silicon nanowire beams; Young's modulus, initial stress, and crystal orientation
A non-linear model is introduced describing the force-deflection relation of doubly clamped beams, including initial stress. Several approximations for the exact model are developed and compared, revealing the importance of considering the initial stress during 3-point bending measurements analysis. A novel approximation is found to be better than others, and both the exact model and this approximation are in perfect agreement with finite element simulations. A brief experimental example of silicon nanowires is presented in which the Young's modulus, the initial stress, and the crystallographic growth orientation are extracted by 3-point bending analysis.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
FY18 Progress on Improvements and V&V of MC 2 -3 and PROTEUS
The PROTEUS code is a high-fidelity capable deterministic neutron transport code based on unstructured finite element meshes, which solves the steady-state and transient neutron transport problem using the 2nd-order discrete ordinate method (SN2ND), the method of characteristics (MOC), and the advanced nodal transport method (NODAL). The PROTEUS system in a broader sense consists of cross section generation tools, mesh generation tools, and output post processing tools as well as neutron transport solvers.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Stochastic Design Optimization of Microstructural Features Using Linear Programming for Robust Design
Microstructure design can have a substantial effect on the performance of critical components in numerous aerospace applications. However, the stochastic nature of metallic microstructures leads to deviations in material properties from the design point and alters the performance of these critical components. In this paper, a novel stochastic linear programming (LP) methodology is developed for microstructure design accounting for uncertainties in desired properties. The metallic microstructure is represented using a finite element discretized form of the orientation distribution function (ODF). The inverse LP problem solves the mean values and covariance matrix of the ODFs to maximize the mean values of a property, given the statistical constraints on other properties. The highlight is an analytical uncertainty quantification model via a Gaussian distribution to model propagation of microstructural uncertainties to the properties. Examples illustrate maximization of the yield strength and magnetostriction of a galfenol alloy when constrained by uncertainties in a set of stiffness constants.
42 ENGINEERING↗
Transient MOC with frequency transform and DSA on unstructured mesh
We present an implementation of the transient method of characteristics (MOC) with isotropic time derivatives, accelerated by diffusion synthetic acceleration (DSA). The fully implicit frequency transform method is used to solve the transient problem with analytic precursor integration. The code works on meshes composed of almost any of the commonly used non-curvilinear finite element types, and can handle the deformation of geometry in time-dependent transport calculations. We present results of a continuous Fourier analysis for the transient multigroup DSA problem, and representative benchmarking results are presented for the C5G7-TD benchmark in 2D showing reasonable performance and agreement compared to other codes. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
A stable potential-based time-domain method for wideband elec- tromagnetic analysis
In previous research, the frequency-domain A-ϕ formulation has been validated using the finite element method for electromagnetic simulations of low-frequency and multi- scale problems, demonstrating excellent numerical accuracy, good matrix condition, and high computational efficiency. Time- domain simulations provide significant advantages for modeling wideband problems and are crucial for multiphysics applications. In this paper, the frequency-domain A-ϕ formulation is extended to the time domain. The central difference scheme is employed for temporal discretization to ensure both accuracy and stability. A numerical example is presented to demonstrate the capability of the proposed time-domain method in wideband electromagnetic analysis.
Mekonnen, Minyichil↗
A stochastic approach to uncertainty in the equations of MHD kinematics
The magnetohydrodynamic (MHD) kinematics model describes the electromagnetic behavior of an electrically conducting fluid when its hydrodynamic properties are assumed to be known. In particular, the MHD kinematics equations can be used to simulate the magnetic field induced by a given velocity field. While prescribing the velocity field leads to a simpler model than the fully coupled MHD system, this may introduce some epistemic uncertainty into the model. If the velocity of a physical system is not known with certainty, the magnetic field obtained from the model may not be reflective of the magnetic field seen in experiments. Additionally, uncertainty in physical parameters such as the magnetic resistivity may affect the reliability of predictions obtained from this model. By modeling the velocity and the resistivity as random variables in the MHD kinematics model, we seek to quantify the effects of uncertainty in these fields on the induced magnetic field. We develop stochastic expressions for these quantities and investigate their impact within a finite element discretization of the kinematics equations. We obtain mean and variance data through Monte Carlo simulation for several test problems. Toward this end, we develop and test an efficient block preconditioner for the linear systems arising from the discretized equations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Verification Problems for Parameterized Load Curves in DYNA3D
Two verification problems are presented which test a set of newly implemented load curve features in the explicit finite element program DYNA3D. Each test consists of a single element, with the velocity of each node specified via a different load curve type. The first test covers initialization and run-time execution of these new features, while the second test covers redefinition of load curves upon restart. The proposed problems may be easily extended to test additional load curve features.
42 ENGINEERING↗
The Target-Matrix Optimization Paradigm for High-Order Meshes
Here, we describe a framework for controlling and improving the quality of high-order finite element meshes based on extensions of the Target-Matrix Optimization Paradigm (TMOP) of [P. Knupp, Eng. Comput., 28 (2012), pp. 419--429]. This approach allows high-order applications to have a very precise control over local mesh quality, while still improving the mesh globally. We address the adaption of various TMOP components to the settings of general isoparametric element mappings, including the mesh quality metric in 2D and 3D, the selection of sample points and the solution of the resulting mesh optimization problem. We also investigate additional practical concerns, such as tangential relaxation and restricting the deviation from the original mesh. The benefits of the new high-order TMOP algorithms are illustrated on a number of test problems and examples from a high-order arbitrary Lagrangian--Eulerian (ALE) application [BLAST: High-order curvilinear finite elements for shock hydrodynamics, http://www.llnl.gov/CASC/blast]. Our implementation is freely available in an open-source library form [MFEM: Modular parallel finite element methods library, http://mfem.org].
97 MATHEMATICS AND COMPUTING↗
Analysis of a Computational Framework for Bayesian Inverse Problems: Ensemble Kalman Updates and MAP Estimators under Mesh Refinement
This paper analyzes a popular computational framework to solve infinite-dimensional Bayesian inverse problems, discretizing the prior and the forward model in a finite-dimensional weighted inner product space. We demonstrate the benefit of working on a weighted space by establishing operator-norm bounds for finite element and graph-based discretizations of Matérn-type priors and deconvolution forward models. For linear-Gaussian inverse problems, we develop a general theory to characterize the error in the approximation to the posterior. We also embed the computational framework into ensemble Kalman methods and MAP estimators for nonlinear inverse problems. Furthermore, our operator-norm bounds for prior discretizations guarantee the scalability and accuracy of these algorithms under mesh refinement.
Bayesian inverse problem↗
A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines
We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.
97 MATHEMATICS AND COMPUTING↗
Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment
We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.
97 MATHEMATICS AND COMPUTING↗
Boundary coefficient determination for an eddy current problem based on the potential field method
We study a recovery problem for an unknown boundary coefficient relating to one material characteristic in an eddy current field. The field equations are represented in terms of the potential field method (T-ψ method) and can be solved numerically by the nodal finite element method. We introduce a measurement as an additional condition and prove the existence and uniqueness of the weak solution. Further, we present an iteration algorithm for the recovery problem and validate its efficiency by two numerical experiments.
97 MATHEMATICS AND COMPUTING↗
Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces
Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Shell Element Material Model Verification Problems for DYNA3D
A suite of shell element material model verification tests were developed for the explicit finite element program DYNA3D. Each test consists of a single quadrilateral shell element, whose nodal displacements and rotations have been prescribed such that the deformation of the material is fully controlled throughout the analysis. At one of the element’s integration points, the numerically computed stress, through-thickness strain, and equivalent plastic strain are compared against closed-form reference solutions that are derived to be consistent with the theory of finite deformation elasto-plasticity. A collection of 31 feature-specific verification tests are proposed, each of which may be used to verify all models which share in common the tested feature of interest. For each material model, a subset of these tests are selected and used to verify the correct behavior of all features of the model. The collective suite of tests cover 26 of the 27 currently available material models in DYNA3D. In the course of developing the proposed test suite, 18 implementation bugs were identified and have subsequently been resolved.
42 ENGINEERING↗