Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “element”

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 577 records · Page 32

Peridynamic Multiscale Finite Element Methods

The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the art of local models with the flexibility and accuracy of the nonlocal peridynamic model. In the mixed locality method this coupling occurs across scales, so that the nonlocal model can be used to communicate material heterogeneity at scales inappropriate to local partial differential equation models. Additionally, the computational burden of the weak form of the peridynamic model is reduced dramatically by only requiring that the model be solved on local patches of the simulation domain which may be computed in parallel, taking advantage of the heterogeneous nature of next generation computing platforms. Addition- ally, we present a novel Galerkin framework, the 'Ambulant Galerkin Method', which represents a first step towards a unified mathematical analysis of local and nonlocal multiscale finite element methods, and whose future extension will allow the analysis of multiscale finite element methods that mix models across scales under certain assumptions of the consistency of those models.

42 ENGINEERING↗

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗

Elemental Analysis of Bio-Oils by Inductively Coupled Plasma Optical Emission Spectroscopy (ICP-OES). Laboratory Analytical Procedure (LAP), Issue Date: May 13, 2022

Concentrations of inorganic elements is a key quality metric for bio-oils as certain elements impact upgrading processes and product quality. Unless reduced or removed during production and processing, alkali and alkaline metals native to lignocellulosic biomass can carry over into bio-oils contributing to ash content and degraded catalyst performance during upgrading to hydrocarbon fuels or chemical products. Non-metallic elements such as sulfur and phosphorus can also negatively impact upgrading catalysts and product quality. Inductively coupled plasma optical emission spectroscopy (ICP-OES) can be used to measure inorganic elements of interest in bio-oils. This procedure covers the preparation and analysis of fast pyrolysis (FP) and catalytic fast pyrolysis (CFP) bio-oils. The concentrations of these elements can indicate potential bio-oil quality from the perspective of deoxygenation processes. The implications of elemental composition will depend on process parameters such as upgrading catalyst sensitivities. This Laboratory Analytical Procedure (LAP) covers two methods for quantification of inorganic elements by ICP-OES: Procedure A uses microwave assisted digestion with concentrated nitric acid, and Procedure B is an organic ICP-OES method utilizing a diglyme solvent.

09 BIOMASS FUELS↗

Verification and Regression Testing of a Physically Stabilized Layered Solid Element Formulation in DYNA3D/Paradyn

Recent modeling of filament-wound composite structures drove the need for a new element type within the finite element code DYNA3D/Paradyn. This new layered solid element was implemented and designed to capture the kinematics and constitutive be havior of various lamina layers defined with arbitrary orientations and volume fractions within the element to accurately model a laminated composite material. The new el ement uses a single integration point in each of the lamina layers defined to capture the constitutive behavior, and a novel physical stabilization routine is used to prevent hourglassing while mitigating shear and volumetric locking. This report documents the verification testing and the regression tests added to the DYNA3D/Paradyn software quality assurance test suite to assess the proper implementation of the new layered solid element. There were a total of 27 tests added the DYNA3D test suite, which consist of problems using the layered solid element with a single layer or multiple layers for isotropic and orthotropic material models. These tests are a combination of simple kinematically driven patch tests, beam bending tests, plate bending tests, and more complicated problems. Ultimately, the testing done verifies that the element behaves as expected and is implemented correctly.

42 ENGINEERING↗

Preserving Superconvergence of Spectral Elements for Curved Domains [Slides]

Finite Element Methods (FEM) and Spectral Element Methods (SEM) are crucial for solving partial differential equations (PDEs) on complex geometries. SEM offers superior accuracy due to potential superconvergence for simple domains. Challenges persist for domains with curved boundaries, restricting SEM’s advantages in real-world applications. A proposed solution is the introduction of a novel strategy to enhance accuracy and maintain superconvergence of SEM in curved domains. The strategy includes a mesh-generation procedure with geometrically refined elements near curved boundaries and a post-processing phase using the Adaptive Extended Stencil Finite Element Method (AES-FEM). The method, named AES-FEM post-processed Spectral Element Method (ApSEM), aligns the accuracy of non-tensor-product elements with superconvergent spectral elements.

97 MATHEMATICS AND COMPUTING↗

Processing MCNP Elemental Edit Outputs

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. A UM geometry is a collection of elements representing a solid geometry. The first step of MCNP UM modeling is using other software packages to create a finite element mesh representation of a solid 3D geometry. Computer-aided design (CAD) or computer-aided manufacturing (CAM) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. MCNP can process a UM model consisting of several different element types including linear tetrahedral or hexahedral elements and calculate quantities of interest such as flux and energy deposition at elements. An MCNP UM simulation provides high-fidelity elemental edit (i.e., tally) outputs, which can be further used in multiphysics calculations. The MCNP UM feature was used for multiphysics simulations where quantities of interest calculated by MCNP are used as inputs for heat transfer calculations in Abaqus. MCNP6.3 can produce two types of elemental edit output (EEOUT) file formats: ASCII and HDF5. An EEOUT file type must be requested on an EMBED card while output type (flux or energy deposition) must be requested on an EMBEE card. We wrote Python3 scripts to extract energy deposition values in an ASCII or HDF5 EEOUT file and compute a heat flux profile for an Abaqus heat transfer calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

How Many Elements Matter?

Some studies of stars' multielement abundance distributions suggest at least 5–7 significant dimensions, but others show that many elemental abundances can be predicted to high accuracy from [Fe/H] and [Mg/Fe] (or [Fe/H] and age) alone. We show that both propositions can be, and are, simultaneously true. We adopt a machine-learning technique known as normalizing flow to reconstruct the probability distribution of Milky Way disk stars in the space of 15 elemental abundances measured by APOGEE. Conditioning on T eff and $\mathrm{log}\,g$ minimizes the differential systematics. After further conditioning on [Fe/H] and [Mg/Fe], the residual scatter for most abundances is σ [X/H] ≲ 0.02 dex, consistent with APOGEE's reported statistical uncertainties of ~0.01–0.015 dex and intrinsic scatter of 0.01–0.02 dex. Despite the small scatter, residual abundances display clear correlations between elements, which we show are too large to be explained by measurement uncertainties or by the finite sampling noise. We must condition on at least seven elements to reduce the correlations to a level consistent with the observational uncertainties. Our results demonstrate that cross-element correlations are a much more sensitive probe of a hidden structure than dispersion, and they can be measured precisely in a large sample even if the star-by-star measurement noise is comparable to the intrinsic scatter. We conclude that many elements have an independent story to tell, even for the mundane disk stars and elements produced by the core-collapse and Type Ia supernovae. The only way to learn these lessons is to measure the abundances directly, and not merely infer them.

79 ASTRONOMY AND ASTROPHYSICS↗

The R-Process Alliance: Abundance Universality among Some Elements at and between the First and Second R-Process Peaks*

Abstract We present new observational benchmarks of rapid neutron-capture process ( r -process) nucleosynthesis for elements at and between the first ( A ∼ 80) and second ( A ∼ 130) peaks. Our analysis is based on archival ultraviolet and optical spectroscopy of eight metal-poor stars with Se ( Z = 34) or Te ( Z = 52) detections, whose r -process enhancement varies by more than a factor of 30 (−0.22 ≤ [Eu/Fe] ≤ +1.32). We calculate ratios among the abundances of Se, Sr through Mo (38 ≤ Z ≤ 42), and Te. These benchmarks may offer a new empirical alternative to the predicted solar system r -process residual pattern. The Te abundances in these stars correlate more closely with the lighter r -process elements than the heavier ones, contradicting and superseding previous findings. The small star-to-star dispersion among the abundances of Se, Sr, Y, Zr, Nb, Mo, and Te (≤0.13 dex, or 26%) matches that observed among the abundances of the lanthanides and third r -process-peak elements. The concept of r -process universality that is recognized among the lanthanide and third-peak elements in r -process-enhanced stars may also apply to Se, Sr, Y, Zr, Nb, Mo, and Te, provided the overall abundances of the lighter r -process elements are scaled independently of the heavier ones. The abundance behavior of the elements Ru through Sn (44 ≤ Z ≤ 50) requires further study. Our results suggest that at least one relatively common source in the early Universe produced a consistent abundance pattern among some elements spanning the first and second r -process peaks.

74 ATOMIC AND MOLECULAR PHYSICS↗

Chemical Doppelgangers in GALAH DR3: The Distinguishing Power of Neutron-capture Elements among Milky Way Disk Stars

The observed chemical diversity of Milky Way stars places important constraints on Galactic chemical evolution and the mixing processes that operate within the interstellar medium. Recent works have found that the chemical diversity of disk stars is low. For example, the Apache Point Observatory Galactic Evolution Experiment (APOGEE) "chemical doppelganger rate," or the rate at which random pairs of field stars appear as chemically similar as stars born together, is high, and the chemical distributions of APOGEE stars in some Galactic populations are well-described by two-dimensional models. However, limited attention has been paid to the heavy elements (Z > 30) in this context. In this work, we probe the potential for neutron-capture elements to enhance the chemical diversity of stars by determining their effect on the chemical doppelganger rate. We measure the doppelganger rate in GALactic Archaeology with HERMES DR3, with abundances rederived using The Cannon, and find that considering the neutron-capture elements decreases the doppelganger rate from ~2.2% to 0.4%, nearly a factor of 6, for stars with –0.1 < [Fe/H] < 0.1. While chemical similarity correlates with similarity in age and dynamics, including neutron-capture elements does not appear to select stars that are more similar in these characteristics. Our results highlight that the neutron-capture elements contain information that is distinct from that of the lighter elements and thus add at least one dimension to Milky Way abundance space. This work illustrates the importance of considering the neutron-capture elements when chemically characterizing stars and motivates ongoing work to improve their atomic data and measurements in spectroscopic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Stepwise chemical abrasion–isotope dilution–thermal ionization mass spectrometry with trace element analysis of microfractured Hadean zircon

Abstract. The Hadean Jack Hills zircons represent the oldest known terrestrial material, providing a unique and truly direct record of Hadean Earth history. This zircon population has been extensively studied via high-spatial-resolution high-throughput in situ isotopic and elemental analysis techniques, such as secondary ionization mass spectrometry (SIMS), but not by comparatively destructive, high-temporal-precision (<0.05 % two-sigma) thermal ionization mass spectrometry (TIMS). In order to better understand the lead loss and alteration history of terrestrial Hadean zircons, we conduct stepwise chemical abrasion–isotope dilution–thermal ionization mass spectrometry with trace element analysis (CA-ID-TIMS-TEA) on manually microfractured Hadean Jack Hills zircon fragments previously dated by SIMS. We conducted three successive HF leaching steps on each individual zircon fragment, followed by column chromatography to isolate U–Pb and trace element fractions. Following isotopic and elemental analysis, the result is an independent age and trace element composition for each leachate of each zircon fragment. We observe ~50 Myr of age heterogeneity in concordant residues from a single zircon grain, along with a protracted history of post-Hadean Pb loss with at least two modes circa ~0 and 2–4 Ga. Meanwhile, stepwise leachate trace element chemistry reveals enrichments of light rare earth elements, uranium, thorium, and radiogenic lead in early leached domains relative to the zircon residue. In addition to confirming the efficacy of the LREE-I alteration index and providing new insight into the mechanism of chemical abrasion, the interpretation and reconciliation of these results suggest that Pb loss is largely driven by low-temperature aqueous recrystallization and that regional thermal events may act to halt – not initiate – Pb loss from metamict domains in the Hadean Jack Hills zircons.

Brenhin Keller, C. (ORCID:0000000174009428)↗

Application of finite-element method in the computation of temperature with emphasis on radiative exchanges.

Analyses pertaining to the solution of heat transfer problems in combined modes based on the finite-element method are presented. Two elements - a triangular element employing two spatial variables and a multi-faceted bar element employing one spatial variable - with nonlinear radiation on the boundaries are detailed. The radiative effects considered on the diffuse-gray surface elements include: (1) directional radiant fluxes from distant sources, (2) radiative exchanges including surfaces of prescribed temperatures, and (3) radiative exchanges including elements whose temperatures are not known a priori. The nonlinear part of the boundary condition is treated in two different ways: (1) a consistent linearization method, and (2) a direct energy distribution method. Applications of these elements together with solution algorithms to three sample problems exhibit solution capability and obtainable accuracy.

Lee, H.-P.↗

Atmophilic elements in lunar soils

Ten volatile elements are more abundant in soil 14163 than in 14259 by factors of 1.4 or greater, too high to be explained by the 15% higher KREEP content of 14163. The enrichments suggest that these elements are easily volatilized by energetic events generated by solar particles or micrometeorites and as a result move from site to site through the lunar atmosphere. The elements Zn, C, and Cd are correlated with excess Ar-40 in Apollo 15 soils. Since concentrations also increase in the direction of the Apennine Front, the correlation chiefly indicates that the Front regolith is enriched in these elements relative to the Apollo 15 regolith away from the Front. The absence of known lunar rock types showing Zn and Cd concentrations as high as those in the Apollo 15 Front regolith or the Apollo 14 or 16 highlands regoliths and the relatively uniform concentrations of both elements in these soils despite differences in major element composition suggest that large amounts of these elements were released as volatiles during early lunar magmatic activity.

Chou, C.-L.↗

A correlation study of finite-element modeling for vibrations of composite material panels

A correlation study has been made to evaluate the effectiveness of finite-element modeling for vibrations of composite material panels. The three types of composite material panels studied were: (1) symmetrically laminated cross-ply panels, (2) unsymmetrically laminated cross-ply panels, and (3) a 1/6-scale sandwich wing model. Correlations of finite-element predictions for frequencies and vibration modes were made with Rayleigh-Ritz analyses and experimental data. For the symmetrically laminated panels, the finite-element representation predicted vibration modes very accurately. For the unsymmetrically laminated panels it was found that predictions using uncoupled finite-element plate bending elements did not give good agreement with experimental data. The constant-chord, swept-wing model was effectively analyzed by finite-element presentations using either sandwich elements or ribs and spars.

Thornton, E. A.↗

TAP 1: A Finite Element Program for Steady-State Thermal Analysis of Convectively Cooled Structures

The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature dependent thermal parameters is performed using the Newton-Raphson iteration method. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. A companion plotting program for displaying the finite element model and predicted temperature distributions is presented. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

Finite-element methodology for thermal analysis of convectively cooled structures

A finite-element method for steady-state thermal analysis of convectively cooled structures is presented. The method is based on representing the coolant passages by finite elements with fluid bulk temperature nodes and fluid/structure interface modes. Four finite elements are described: two basic elements (mass transport and surface convection) and two integrated elements for applications to discrete tube and plate-fin convectively cooled structural configurations. Comparative finite-element and lumped-parameter thermal analyses of several convectively cooled structures demonstrate the practicality of utilizing finite-element methodology for thermal analysis of realistic convectively cooled structures.

Source record↗

A direct element resequencing procedure

Element by element frontal solution algorithms are utilized in many of the existing finite element codes. The overall computational efficiency of this type of procedure is directly related to the element data input sequence. Thus, it is important to have a pre-processor which will resequence these data so as to reduce the element wavefronts to be encountered in the solution algorithm. A direct element resequencing algorithm is detailed for reducing element wavefronts. It also generates computational by products that can be utilized in pre-front calculations and in various post-processors. Sample problems are presented and compared with other algorithms.

Akin, J. E.↗

TAP 2: A finite element program for thermal analysis of convectively cooled structures

A finite element computer program (TAP 2) for steady-state and transient thermal analyses of convectively cooled structures is presented. The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature-dependent thermal parameters is performed using the Newton-Raphson iteration method. Transient analyses are performed using an implicit Crank-Nicolson time integration scheme with consistent or lumped capacitance matrices as an option. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

A finite element approach for shells of revolution with a local deviation

A finite element model that is suitable for the static analysis of shells of revolution with arbitrary local deviations is presented. The model employs three types of elements: rotational, general, and transitional shell elements. The rotational shell elements are used in the region where the shell is axisymmetric. The general shell element are used in the local region of the deviation. The transitional shell elements connect these two distinctively different types of elements and make it possible to combine them in a single analysis. The form of the global stiffness matrix resulting when different forms of nodal degrees of freedom are combined is illustrated. The coupling of harmonic degrees of freedom due to the locally nonaxisymmetric geometry was studied. The use of a substructuring technique and separate partial harmonic analysis is recommended.

Han, K. J.↗