Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A direct-adjoint approach for material point model calibration with application to plasticity

Here, this paper proposes a new approach for the calibration of material parameters in local elastoplastic constitutive models. The calibration is posed as a constrained optimization problem, where the constitutive model evolution equations for a single material point serve as constraints. The objective function quantifies the mismatch between the stress predicted by the model and corresponding experimental measurements. To improve calibration efficiency, a novel direct-adjoint approach is presented to compute the Hessian of the objective function, which enables the use of second-order optimization algorithms. Automatic differentiation is used for gradient and Hessian computations. Two numerical examples are employed to validate the Hessian matrices and to demonstrate that the Newton–Raphson algorithm consistently outperforms gradient-based algorithms such as L-BFGS-B.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Random field optimization

Herein we present a new modeling paradigm for optimization that we call random field optimization. Random fields are a powerful modeling abstraction that aims to capture the behavior of random variables that live on infinite-dimensional spaces (e.g., space and time) such as stochastic processes (e.g., time series, Gaussian processes, and Markov processes), random matrices, and random spatial fields. This paradigm involves sophisticated mathematical objects (e.g., stochastic differential equations and space-time kernel functions) and has been widely used in neuroscience, geoscience, physics, civil engineering, and computer graphics. Despite of this, however, random fields have seen limited use in optimization; specifically, existing optimization paradigms that involve uncertainty (e.g., stochastic programming and robust optimization) mostly focus on the use of finite random variables. This trend is rapidly changing with the advent of statistical optimization (e.g., Bayesian optimization) and multi-scale optimization (e.g., integration of molecular sciences and process engineering). Our work extends a recently-proposed abstraction for infinite-dimensional optimization problems by capturing more general uncertainty representations. Moreover, we discuss solution paradigms for this new class of problems based on finite transformations and sampling, and identify open questions and challenges.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Matrix-free subcell residual distribution for Bernstein finite elements: Monolithic limiting

This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. Finally, the convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Extracting single fiber transverse and shear moduli from off-axis misalignment fiber tensile testing

Small diameter (<100 μm) fibers (e.g. carbon fibers, Kevlar, and fiberglass) and wires (e.g. ultrafine copper and aluminum wires) are frequently used in many different engineering applications, such as for light weighting structures, electromagnetic shielding for aircraft/infrastructure/EVs, vibration damping, biological sensors, aerospace electrical devices, and electric windings just to name a few. Due to the manufacturing process, the fibers and wires are pulled and stretched to produce a preferential alignment. Therefore, thin fibers and wires typically display different properties along the length of the fiber as opposed to their cross section and many fibers/wires are considered transversely isotropic. The axial properties of fibers/wires can be ascertained via tensile testing of single-filaments or fiber tows, but the radial properties require much more effort to measure. Knowing these properties is important for the accurate prediction of micromechanical models and manipulation of fibers during micromanufacturing. In this paper, a new technique was developed to determine the transverse/shear moduli and strength of a material by conducting tensile tests of the material at increasing misalignment angles from the tensile axis. Due to the transversely isotropic nature of the material, the transverse/shear moduli and strength influence the experimental results recorded by the test machine to different degrees based on the amount of misalignment in the test setup. An equation was derived to determine the influence of each of the material properties based on the misalignment angle by manipulating the stiffness matrix for transversely isotropic materials using the transformation matrices. Then, curve fitted coefficients were used to identify the material properties. Here, the strengths were similarly determined by curve fitting an off-axis Tsai-Hill failure criteria to determine the influence of transverse, shear, and tensile strengths based on the complex loading condition provided by the off-axis tensile test. Zoltek Panex 35 carbon fibers were used to demonstrate this new technique and the determined properties were then compared to those obtained from nanoindentation and from literature. Fracture surfaces provide insight into the different failure mechanisms at various misalignment angles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quest for environmentally sustainable materials: A case for animal-based fillers and fibers in polymeric biocomposites

This review explores the potential of animal-based fillers and fibers as eco-friendly alternatives to conventional synthetic ones. Examining materials such as wool, silk, feather, hair and beak, the review elucidates their chemistry, structure, properties and sources, emphasizing biodegradability and renewability. It also discusses the compatibility of these materials with polymer matrices and their mechanical, acoustic and thermal performances. Further, the review critically analyzes environmental and ethical implications, presenting challenges and opportunities in the emerging field. By addressing ecological and performance aspects, it contributes to global efforts in fostering sustainability in materials science. Future research to address gaps and enhance the design, manufacture and application of animal-based reinforcements in various industries are clearly outlined at the end of the review.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Recent advancements of plant-based natural fiber–reinforced composites and their applications

Demands for reducing energy consumption and environmental impacts are the major driving factors for the development of natural fiber–reinforced composites (NFRCs) in many sectors. Compared with synthesized fiber, natural fiber provides several advantages in terms of biodegradability, light weight, low price, life-cycle superiority, and satisfactory mechanical properties. However, the inherent features of plant-based natural fibers have presented challenges to the development and application of NFRCs, such as variable fiber quality, limited mechanical properties, water absorption, low thermal stability, incompatibility with hydrophobic matrices, and propensity to agglomeration. Substantial research has recently been conducted to address these challenges for improved performance of NFRCs and their applications. This article reviews the recent advancements of plant-based NFRCs, focusing on strategies and breakthroughs in enhancing the NFRCs’ performance, including fiber modification, fiber hybridization, lignocellulosic fillers incorporation, conventional processing techniques, additive manufacturing (3D printing), and new fiber source exploration. In conclusion, the sustainability of plant-based NFRCs using life-cycle assessment and the burgeoning applications of NFRCs with emphasis on the automotive industry are also discussed.

3D reinforcement↗
DOE OSTI2024

Understanding interfacial crystallization dynamics on carbon fiber reinforced polypropylene composite manufacturing

Reinforcing polymers with discontinuous fibers improves mechanical properties, such as strength and stiffness, and in some cases achieve isotropic properties, rendering them suitable for various engineering applications. Matrix materials are generally highly engineered thermosets (e.g. crosslinked epoxies), bonded to the fiber periphery by proprietary surface and sizing chemistries. Semicrystalline thermoplastic matrices are less utilized due to poor fiber-matrix bonding resulting in inefficient interfacial load-transfer in reinforced composites. However, flexibility with melt-processing or molding conditions can be leveraged to promote non-covalent interfacial bonding between matrix and fiber via crystallization of the matrix onto fiber surface. In the present study, we utilize a co-mingle chopped carbon and isotactic polypropylene fibers to form isotropic composites, tailoring interfacial immobilized matrix or interphase morphology to optimize performance through precise control of thermal processing/molding windows. Calorimetry and optical microscopy were employed to investigate the impact of carbon fiber at various volume fractions (10, 20, and 30 %) on isotactic polypropylene crystallization and mechanical performance. Variations in mechanical properties correspond to the structural evolution of the interfacial region and are correlated to underlying microstructural attributes using wide-angle X-ray scattering, thermal analysis, and low-field nuclear magnetic resonance spectroscopy. These results provide a practical framework for the manufacturing of thermoplastic matrix composites. In conclusion, the results presented provide a guide for the strategic optimization of interphase design, showcasing tailorable tensile strengths which outperform any isotactic polypropylene carbon fiber composites previously reported in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crashworthiness of recycled carbon fiber composite sinusoidal structures at dynamic rates

Fiber reinforced polymer composites are finding applications in the automotive space for structurally critical applications including crash management. The use of recycled carbon fiber reinforcement can greatly reduce carbon emissions. In this study, recycled carbon fiber composites were manufactured into a self-supporting geometry using three matrices: polyphenylene sulfide (PPS), acrylonitrile butadiene styrene (ABS), and a structural epoxy, and subsequently crushed between flat platens at dynamic rates ranging from 4.6 to 9.1 ms -1 at temperatures in a range of -40°C to 80°C. Load-displacement data was used to evaluate their specific energy absorption, crush efficiency, and steady-state crush stress. Further, the energy absorption of the ABS composites was strongly sensitive to temperature, while the PPS composites exhibited strong crush efficiency dependence with both rate and temperature. The epoxy composites exhibited stable crush behavior at dynamic crush rates but exhibited a dramatic reduction in crush efficiency relative to quasi-static tests. The results of this study indicate that recycled fiber composites can achieve very high energy absorption levels (50–80kJkg -1 at room temperature) that make them an excellent alternative to more expensive and less environmentally friendly continuous virgin fiber laminates.

42 ENGINEERING↗
DOE OSTI2020

Influence of mineralogical and chemical compositions on alkali-silica-reaction of Tennessee limestones

The paper aims to identify potential correlations between mineralogical and chemical compositions of AAR reactive dolomitic limestone aggregates from Tennessee and their expansion on ASTM C1293 for quick evaluation of reactivity. The chemical compositions of 26 limestones, namely the silica content (SiO 2 ), showed a strong correlation with the expansion in the concrete prism test. The early and highly expansive limestones are characterized by a high clay/ mica content (>5%). The reactivity of Tennessee limestones is mainly caused by crypto/microcrystalline quartz, chert, and finely disseminated silica in carbonate and clay/mica matrices, and it can be reasonably evaluated using X-ray diffraction and chemical analysis coupled with petrographic examinations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Biotechnological solutions for critical mineral recovery from unconventional feedstocks

Secure and sustainable metal recovery from unconventional feedstocks is needed to meet the mineral demands of energy, defense, and electronic technologies. Here, we highlight the potential to leverage nature’s ability to extract and differentiate metal ions in biotechnologies that could become the next generation of mining and refining. We describe bulk and trace processes and then discuss the advances and opportunities of two key bioprocesses: microbially mediated solubilization of metal ions from solid matrices (termed ‘bioleaching’) and bio-based separation of solubilized ions via selective adsorption to proteins. Both biotechnologies have advantages such as reduced energy input for leaching low-grade feedstocks and reduced organic solvent demand for separating ions with similar physiochemical properties but require more development for industrial scale recovery from unconventional feedstocks. Innovation in biological science and engineering may bring timely solutions to key challenges toward recovering critical minerals from unconventional feedstocks.

organic↗
DOE OSTI2021

Calculating elements of matrix functions using divided differences

In this work, we introduce a method for calculating individual elements of matrix functions. Our technique makes use of a novel series expansion for the action of matrix functions on basis vectors that is memory efficient even for very large matrices. We showcase our approach by calculating the matrix elements of the exponential of a transverse-field Ising model and evaluating quantum transition amplitudes for large many-body Hamiltonians of sizes up to 2 64 x 2 64 on a single workstation. We also discuss the application of the method to matrix inverses. We relate and compare our method to the state-of-the-art and demonstrate its advantages. We also discuss practical applications of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Accelerating multigrid with streaming chiral SVD for Wilson fermions in lattice QCD

A modification to the setup algorithm for the multigrid preconditioner of Wilson fermions in lattice QCD is presented. A larger basis of test vectors than that used in regular multigrid is calculated by the smoother and truncated by singular value decomposition on the chiral components of the test vectors. The truncated basis is used to form the prolongation and restriction matrices of the multigrid hierarchy. This modification of the setup method is demonstrated to increase the convergence of linear solvers on an anisotropic lattice with m π ≈ 239 MeV from the Hadron Spectrum Collaboration and an isotropic lattice with m π ≈ 220 MeV from the MILC Collaboration. The lattice volume dependence of the method is also examined. Increasing the number of test vectors improves speedup up to a point, but storing these vectors becomes impossible in limited memory resources such as GPUs. To address storage cost, we implement a streaming singular value decomposition of the basis of test vectors on the chiral components and demonstrate a decrease in the number of fine level iterations by a factor of 1.7 for m q ≈ m crit

Iterative methods↗
DOE OSTI2022

Evaluating the performance of random forest and iterative random forest based methods when applied to gene expression data

Gene-to-gene networks, such as Gene Regulatory Networks (GRN) and Predictive Expression Networks (PEN) capture relationships between genes and are beneficial for use in downstream biological analyses. There exists multiple network inference tools to produce these gene-to-gene networks from matrices of gene expression data. Random Forest-Leave One Out Prediction (RF-LOOP) is a method that has been shown to be efficient at producing these gene-to-gene networks, frequently known as GEne Network Inference with Ensemble of trees (GENIE3). Random Forest can be replaced in this process by iterative Random Forest (iRF), which performs variable selection and boosting. Here we validate that iterative Random Forest-Leave One Out Prediction (iRF-LOOP) produces higher quality networks than GENIE3 (RF-LOOP). We use both synthetic and empirical networks from the Dialogue for Reverse Engineering Assessment and Methods (DREAM) Challenges by Sage Bionetworks, as well as two additional empirical networks created from Arabidopsis thaliana and Populus trichocarpa expression data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Polymer electrolytes based on protic ionic liquids with perfluorinated anions for safe lithium-ion batteries

Here, the quest for safe and high-performance polymer electrolytes in lithium-ion batteries (LIBs) has led researchers to explore protic ionic liquids (PILs) as potential candidates to be entrapped in polymer matrices. In this context, we present an investigation into solid polymeric systems based on poly(methyl methacrylate) (PMMA) as a host for PILs, featuring 1,8-diazabicyclo-[5,4,0]-undec-7-ene (DBU) cation paired with three different anions: bis(trifluoromethanesulfonyl)imide (TFSI – ), trifluoromethanesulfonate (TFO – ), and (trifluoromethanesulfonyl-nonafluorobutylsulfonyl)imide (IM14 – ). Additionally, we explore the lithium-doped IM14-gel-like system to broaden our understanding of these intriguing materials. Through comprehensive thermal analysis, solid-state NMR, and diffusion NMR techniques, we delve into the interactions and structural features of these binary and ternary polymeric systems. Our investigation reveals unique dynamics and ion interactions within the PMMA matrix, shedding light on the potential of these materials for advanced energy storage technologies. Particularly, we highlight the distinctive features of DBUH-IM14 and its specific interaction with the polymeric matrix and the lithium ions, underscoring its significance in advancing safer and more efficient energy storage devices.

25 ENERGY STORAGE↗