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 469 records · Page 26

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

A review of composite polymer-ceramic electrolytes for lithium batteries

All solid-state lithium batteries are garnering attention in both academia and industry. Lithium-ion conductive polymers and lithium-ion conductive ceramics are the two major classes of solid electrolytes that have prevalently been pursued for many years. However, each of them has its own advantages and disadvantages. One approach to overcome the disadvantages and get the best out of each of those materials is a solid composite electrolyte that combines the advantages of inorganic ceramic electrolytes and solid polymer electrolytes. Such composite electrolytes can offer acceptable ionic conductivity, high mechanical strength, and favorable interfacial contact with electrodes, which can greatly improve the electrochemical performance of all-solid-state batteries compared to cells based on a polymer electrolyte alone or a ceramic electrolyte alone. We present in this review the state-of-the-art composite polymer-ceramic electrolytes in view of their electrochemical and physical properties for the applications in lithium batteries. The review mainly encompasses the polymer matrices, various ceramic filler materials, and the polymer/ceramics composite systems. Specifically, the structures, ionic conductivities, electrochemical/chemical stabilities, and fabrications of solid composite electrolytes are discussed in-depth. On the basis of previous work, a perspective on future research directions is highlighted for developing high-performance composite polymer-ceramic electrolytes.

25 ENERGY STORAGE↗

Parametric reduced order models for graded lattice structures

Graded lattice structures, characterized by smoothly varying mechanical properties, hold significant promise for optimizing material distribution in advanced engineering applications. However, accurately modeling these structures poses substantial computational challenges due to the continuous geometric variations within their unit cells. Here, to address these challenges, this paper introduces a novel Efficient Reduced Order Model (EROM) that integrates the Matrix Discrete Empirical Interpolation Method (MDEIM) and Discrete Empirical Interpolation Method (DEIM) with polynomial regression to efficiently manage geometric parametrization in lattice structures. Unlike traditional reduced order models (ROMs) that require extensive precomputed libraries for each geometric configuration, our approach enables continuous geometric variations through a flexible algebraic formulation, significantly reducing computational costs while preserving high accuracy. The method constructs projection matrices for individual unit cells that can be efficiently assembled into global systems, leveraging the repetitive nature of lattice structures. Numerical studies demonstrate that our EROM achieves displacement errors below 1% and von Mises stress prediction errors below 4%, coupled with computational speedups exceeding two orders of magnitude compared to full-order simulations. The proposed method's modularity and scalability make it particularly suitable for design optimization and real-time simulation of functionally graded lattice structures, with applications spanning aerospace to biomedical engineering.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Efficient kinetic thermal inverse modeling for organic material decomposition

The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗