Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “difficult”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Krylov subspace recycling for evolving structures

Krylov subspace recycling is a powerful tool when solving a long series of large, sparse linear systems that change only slowly over time. In PDE constrained shape optimization, these series appear naturally, as typically hundreds or thousands of optimization steps are needed with only small changes in the geometry. In this setting, however, applying Krylov subspace recycling can be a difficult task. As the geometry evolves, in general, so does the finite element mesh defined on or representing this geometry, including the numbers of nodes and elements and element connectivity. This is especially the case if re-meshing techniques are used. As a result, the number of algebraic degrees of freedom in the system changes, and in general the linear system matrices resulting from the finite element discretization change size from one optimization step to the next. Changes in the mesh connectivity also lead to structural changes in the matrices. In the case of re-meshing, even if the geometry changes only a little, the corresponding mesh might differ substantially from the previous one. Obviously, this prevents any straightforward mapping of the approximate invariant subspace of the linear system matrix (the focus of recycling in this work) from one optimization step to the next; similar problems arise for other selected subspaces. In this paper, we present an algorithm to map an approximate invariant subspace of the linear system matrix for the previous optimization step to an approximate invariant subspace of the linear system matrix for the current optimization step, for general meshes. This is achieved by exploiting the map from coefficient vectors to finite element functions on the mesh, combined with interpolation or approximation of functions on the finite element mesh. We demonstrate the effectiveness of our approach numerically with several proof of concept studies for a specific meshing technique.

42 ENGINEERING↗
DOE OSTI2022

Learning high-dimensional parametric maps via reduced basis adaptive residual networks

We propose a scalable framework for the learning of high-dimensional parametric maps via adaptively constructed residual network (ResNet) maps between reduced bases of the inputs and outputs. When just few training data are available, it is beneficial to have a compact parametrization in order to ameliorate the ill-posedness of the neural network training problem. By linearly restricting high-dimensional maps to informed reduced bases of the inputs, one can compress high-dimensional maps in a constructive way that can be used to detect appropriate basis ranks, equipped with rigorous error estimates. A scalable neural network learning framework is thus to learn the nonlinear compressed reduced basis mapping. Unlike the reduced basis construction, however, neural network constructions are not guaranteed to reduce errors by adding representation power, making it difficult to achieve good practical performance. Inspired by recent approximation theory that connects ResNets to sequential minimizing flows, we present an adaptive ResNet construction algorithm. This algorithm allows for depth-wise enrichment of the neural network approximation, in a manner that can achieve good practical performance by first training a shallow network and then adapting. We prove universal approximation of the associated neural network class for $L^2_v$ functions on compact sets. Our overall framework allows for constructive means to detect appropriate breadth and depth, and related compact parametrizations of neural networks, significantly reducing the need for architectural hyperparameter tuning. Numerical experiments for parametric PDE problems and a 3D CFD wing design optimization parametric map demonstrate that the proposed methodology can achieve remarkably high accuracy for limited training data, and outperformed other neural network strategies we compared against.

42 ENGINEERING↗
DOE OSTI2023

Computational design of metamaterials with self contact

Inverse homogenization in combination with contact modeling, topology optimization and shape optimization is used to design metamaterials with optimized macroscopic response. The homogenization assumes length scale separation which allows the non-linear macroscopic behavior to be obtained by analyzing a single unit cell in a lattice structure. Self contact in the unit cell, which is modeled using a third medium contact method, is leveraged to obtain a complex homogenized response. The inverse homogenization problem is initially formulated as a topology optimization problem, where the macroscopic stress–strain behavior is tuned to our liking. However, it is well known that boundary phenomena are difficult to model in topology optimization and that interface modeling is crucial to accurately analyze contact. For that reason, the boundary representation of the topology optimized design is extracted and used as initial design in a subsequent shape optimization. The behaviors of our designs are verified by performing rigorous post-processing analyzes using conforming meshes and conventional contact formulations.

42 ENGINEERING↗
DOE OSTI2024

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗
DOE OSTI2025

Boosting efficiency and reducing graph reliance: Basis adaptation integration in Bayesian multi-fidelity networks

The computational cost of high-fidelity numerical models makes outer-loop analysis, which requires repeated interrogation of the model such as uncertainty quantification, computationally demanding. Multi-fidelity methods, which construct a surrogate model using data from an ensemble of models of varying cost and accuracy, can substantially reduce the cost of outer-loop analysis. However, these methods can be difficult to apply when the model ensemble does not admit a clear hierarchy a priori and the correlations between models are low. Consequently, in this paper, we present a multi-fidelity method that leverages dimension reduction to enhance the correlation between models, thereby reducing the amount of data needed to train a surrogate from an unordered ensemble of models. Our method utilizes basis adaptation to build low-dimensional polynomial chaos expansions of each model and employs Multi-fidelity Networks to encode the relationships among models. We show that the resulting method exhibit two notable advantages over its counterpart: (1) enhanced accuracy (both reduced bias and variance); and (2) reduced dependency on the graph structure encoding relationships among models. We demonstrate the approach on an analytical test problem and a challenging finite element model for a spent nuclear fuel. Our method produces a surrogate model that is significantly more accurate than either a single-fidelity surrogate or a multi-fidelity surrogate constructed without basis adaptation.

42 ENGINEERING↗
DOE OSTI2025

A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. Here, the two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Automatic differentiation↗
DOE OSTI2026

Instabilities and phase transitions in architected metamaterials: a gradient-enhanced continuum approach

Architected metamaterials such as foams and lattices exhibit a wide range of properties governed by microstructural instabilities and emerging phase transitions. Their macroscopic response–including energy dissipation during impact, large recoverable deformations, morphing between configurations, and auxetic behavior–remains difficult to capture with conventional continuum models, which often rely on discrete approaches that limit scalability. In this work, we propose a nonlocal continuum formulation that captures both stable and unstable responses of elastic architected metamaterials. The framework extends anisotropic hyperelasticity by introducing nonlocal variables and internal length scales reflective of microstructural features. Local polyconvex free-energy models are systematically augmented with two families of non-(poly)convex energies, enabling both metastable and bistable responses. Implementation in a finite element framework enables solution using a hybrid monolithic–staggered strategy. Simulations capture densification fronts, forward and reverse transitions, hysteresis loops, imperfection sensitivity, and globally coordinated auxetic modes. Overall, this framework provides a robust foundation for accelerated modeling of instability-driven phenomena in architected metamaterials, while enabling extensions to anisotropic, dissipative, and active systems as well as integration with data-driven and machine learning approaches.

42 ENGINEERING↗
DOE OSTI2023

Interlaboratory Study: Testing Reproducibility of Solid Biofuels Component Identification Using Reflected Light Microscopy

Considering global market trends and concerns about climate change and sustainability, increased biomass use for energy is expected to continue. As more diverse materials are being utilized to manufacture solid biomass fuels, it is critical to implement quality assessment methods to analyze these fuels thoroughly. One such method is reflected light microscopy (RLM), which has the potential to complement and enhance current standard testing, leading to improving fuel quality assessment and, ultimately, preventing avoidable air pollution. An interlaboratory study (ILS) was conducted to test the reproducibility of biomass fuels component identification using a reflected light microscopy technique. The exercise was conducted on thirty photomicrographs showing biomass and various undesired components (like plastics or mineral matter), which were purposely added (by the ILS organizers) to contaminate wood pellets and charcoal-based grilling fuels. Forty-six participants had various levels of difficulty identifying the marked components, and as a result, the percentage of correct answers ranged from 52.2 to 94.4%. Among the most difficult components to distinguish were petroleum products and inorganic matter. Various reasons led to the misidentification, including insufficient morphological descriptions of the components provided to participants, ambiguities of the nomenclature, limitations of the analytical and exercise method, and insufficient experience of the participants. Overall, the results indicate that RLM has the potential to enhance the quality assessment of biomass fuels. However, they also demonstrate that the petrographic classification used in this exercise requires further refinement before it can be standardized. While a new simplified classification of solid biomass fuels components was created as an outcome of this study, future research is necessary to refine the nomenclature, develop a microscopic morphological description of the components, and verify the accuracy of component identification with a follow-up ILS.

biomass↗
DOE OSTI2023

High-gradient magnetic separation of colloidal uranium oxide particles from soil components in aqueous suspensions

The separation of uranium oxide (UO 2 ) particles from soil-surrogate particles in aqueous suspensions was achieved using filtration enhanced by a magnetic field. Enhanced attraction of paramagnetic UO 2 colloids to a ferromagnetic stainless-steel filter placed in a strong magnetic field arises because of the positive magnetic susceptibility of the particles and the high-gradient field generated near ferromagnetic fibers. Enhanced uptake of smaller particles over larger ones occurs through Brownian motion that promotes the collision of particles with the ferromagnetic fibers of the filter. Hence, this work focused on UO 2 particles in the colloidal size range. Experiments used a water-cooled electromagnet and an array of permanent magnets. Chemical analysis showed that the magnetic field increased the capture efficiency of uranium particles from a range of 27–53% with the magnet off up to 98% with the magnet on after a single pass of the suspension through the filter. Further, the recovery of the UO 2 particles from the filter, however, was more difficult to achieve. Small amounts of UO 2 , together with significant amounts of background SiO 2 particles, were removed from the filter during a first flush with the magnetic field on. A much larger recovery of UO 2 was not observed until a second out-of-field flush was performed, which also released some SiO 2 . The degree to which particle separation was enhanced through the use of multi-stage filtration compared to single pass-through filtration was also examined. A design was suggested that could be used to optimize the separation efficiency for a continuous process.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A statistical model to predict ignition probability

Ignition capability is a critical design constraint for aeronautical gas turbines. However the current trend toward overall lean burn is detrimental to the engine ignition and relight and the ignition system must be adapted to ensure a fast and reliable light-round in all circumstances. As ignition is a stochastic phenomenon, the optimization of an ignition system requires to build ignition probability maps, which is difficult and costly with either experiment or numerical simulation as both require many tests. Here, this work proposes a model to predict the ignition probability map, knowing only flow statistics in non-reacting conditions, i.e., with only one test. The originality of the model is to construct statistics of the flame kernel trajectory, which are then combined with local flow indicators to evaluate the ignition probability at the considered sparking location. Application to a swirled burner operated in premixed, non-premixed and spray combustion modes illustrates the model concepts and demonstrates its ability to recover the experimental ignition map with good accuracy.

42 ENGINEERING↗
DOE OSTI2021

Data reduction considerations for spherical R-32(CH 2 F 2 )-air flame experiments

The burning velocities of mixtures of refrigerant R-32 (CH 2 F 2 ) with air over a range of equivalence ratios are studied via shadowgraph images of spherically expanding flames (SEFs) in a large, optically accessible spherical chamber at constant pressure. Numerical simulations of the 1-D, unsteady, spherical flames incorporating an optically thin radiation model and detailed kinetics accurately predict experimental results for a range of equivalence ratios. For these low burning velocity flames, the effects of stretch and radiation occur simultaneously and make extraction of the unstretched burning velocity from the experimental data difficult. Different data reduction approaches are shown to have large effects on the burning velocities inferred from the experiments and simulations. A new flame radius tracking approach for the experimental images is shown to provide improved agreement of the burned gas velocity variation with stretch predicted by the simulations and helps to compensate for mild flame distortion due to buoyancy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Fabrication, oxidation, and combustion of nanoscale magnesium diboride and tetraboride

The difficult ignition and low combustion efficiency of boron particles decrease the performance of boron-loaded, fuel-rich propellants for solid fuel ramjets and ducted rockets. One approach to solving this problem involves the use of magnesium diboride (MgB 2 ), which ignites easier than boron. Magnesium tetraboride (MgB 4 ) offers greater energy density owing to its higher boron content. However, the effect of B/Mg ratio on the ignition and combustion is unknown. Additionally, while nanoscale MgB₂ particles and quasi-2D structures are promising energetic additives, the oxidation and combustion properties of nanoscale MgB₄ have not been explored. To address these knowledge gaps, the present work included synthesis and high-energy ball milling of MgB 2 and MgB 4 powders, thermogravimetric analysis (TGA) of their oxidation, and combustion experiments with thin layers of the obtained powders. Comparison of two synthesis routes (a solid-state reaction in a tube furnace and combustion synthesis) has shown that the former is the superior method for producing magnesium borides. TGA has revealed that oxidation of both MgB 2 and MgB 4 results in a high conversion into the oxides (88–91 %), far exceeding the low conversion of boron (62.5 %). MgB 4 begins to oxidize rapidly at a much lower temperature (∼900 °C) than MgB 2 (∼1200 °C). The burning rates of milled MgB 2 and MgB 4 are about eight and five times, respectively, faster than that of submicron boron. Magnesium borides exhibit a stable, sustained boron flame, needed for high combustion efficiency, whereas physical Mg/B mixtures undergo Mg-driven "flash" combustion.

Boron↗
DOE OSTI2020

Three-dimensional phase field sintering simulations accounting for the rigid-body motion of individual grains

Sintering is a widely used powder processing technique in industrial applications. During sintering, atoms migrate to decrease the energy of the system via two main mechanisms: coarsening and densification, both of which lead to significant morphological variation of the sintered microstructure. When simulating sintering dynamics, the phase-field method has been broadly utilized because of its convenience in tracking morphology evolution. When a large number of grains is involved, it is common to use the same order parameter to describe multiple grains that are not in direct contact with one another (in order to reduce the computational memory demands). However, with this treatment it is difficult to handle the rigid-body motion of individual grains during densification. In this work, an implementation scheme is introduced to overcome the challenge of calculating individual particle motion based on existing equations. It uses a grouping algorithm and sets a cutoff radius on each grain for calculating the particle velocity during densification. This method allows for the incorporation of the densification mechanism, which has been commonly ignored in previous work, into phase-field sintering models in three-dimensional simulations with a large number of particles/grains. Moreover, through combination with the smoothed boundary method, material properties of sintered microstructures, such as the effective diffusivity and Young’s modulus, can be calculated during the sintering processes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An automated cluster surface scanning method for exploring reaction paths on metal-cluster surfaces

Metal-cluster surfaces present a wide variety of unique coordination environments. This complexity makes it difficult to manually probe the surface reactivity of such clusters. Here, we present a simple and automated method to systematically discover reaction pathways on cluster surfaces, based on the automated cluster surface scanning (ACSS) technique for mapping out potential energy surfaces. We showcase our method on 55-atom icosahedral Cu and Ag clusters, where we determine the activation energies of four elementary steps common in heterogeneous catalysis – hydrogen recombination (H* + H* → H 2 * + *), oxygen recombination (O* + O* → O 2 * + *), water formation (OH* + H* → H 2 O* + *), and CO oxidation (CO* + O* → CO 2 * + *) – with density functional theory calculations (DFT-PBE + D3). We show that the ACSS method requires significantly less human effort than the established manually performed climbing-image nudged elastic band (MP + CI-NEB) technique and locates transition states with comparable accuracy (root-mean-squared error of 0.10 eV) and similar computational cost. Rigorous sampling of the potential energy surface with the ACSS method allows one to locate all lowest-energy reaction pathways obtained via the MP+CI-NEB approach, as well as alternative pathways that one may have missed with the MP+CI-NEB approach due to the many possible pathways available on these clusters. The accuracy and efficiency afforded by the ACSS method could enable high-throughput exploration of the diverse reactivity of metal clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deep-freeze graph training for latent learning

Scientific and engineering advances are primarily driven by multi-tier conceptual constructs and conditional theoretical frameworks. The theories allow predictions of hypothetical system responses, given a set of approximate conditions (ranges of applicability) imposed on latent parameters that cannot be measured directly. Learning to estimate the latent variables (Latent Learning) helps to pinpoint the anticipated range-edge anomalies and improves the confidence in interpretation, interpolation and extrapolation of limited experimental data. Due to high dimensionality and extreme non-linearity of the materials science problems, very large datasets are typically required for conventional data-driven model development. The vital experimental data collection, particularly on microstructural phases, is very challenging, which makes it difficult to compile a high-quality database. Incorporation of the domain knowledge into the computational graph structure, initialization and optimization processes presents a viable mechanism for developing accurate models, with limited datasets. Furthermore, this study successfully utilized the approach to build the Deep Freeze Graph (DeepFreG) by mapping known causality relationships and by digitizing empirical domain knowledge for Latent Learning (LL), with specific applications in materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using machine-learning to understand complex microstructural effects on the mechanical behavior of Ti-6Al-4V alloys

Structural materials properties are highly dependent on their microstructure. Their microstructure is in turn affected by multiple fabrication and thermo-mechanical treatment parameters, all of which conform a highly-dimensional parametric space with often hidden correlations that are difficult to extract by experimentation alone. This is particularly true for alloys of the dual-phase Ti-6Al-4V family, with their greatly complex and rich microstructures, which combine several intrinsic length scales associated with multiple grain and subgrain structures, grains with different crystal lattices (α and β phases), and complex chemistry. In this paper we use a comprehensive set of machine learning techniques to develop predictive tools relating the yield strength and hardening rate of these alloys to a set of input parameters covering extensive ranges. The data generator is a finite-element crystal plasticity model for polycrystal deformation that takes into account slip anisotropy and employs standard dislocation evolution models for the α and β phases of Ti-based alloys. Our dataset includes over two thousand independent simulations and is used to train the machine learning models, which are then used to establish correlations between microstructural parameters and the alloys’ mechanical response. Our results point to the most influential parameters affecting yield strength and hardening rate, information that can then be used to guide experimental synthesis and characterization efforts to save time and resources.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding structure-processing relationships in metal additive manufacturing via featurization of microstructural images

Understanding and predicting accurate property-structure-processing relationships for additively manufactured components is important for both forward and inverse design of robust, reliable parts and assemblies. While direct mapping of process parameters to properties is sometimes plausible, it is often rendered difficult due to poor microstructural control. Exploring the direct relationship between processing conditions and microstructural features can thus provide significant physical insights and aid the overall design process. Here, in this study, we develop an automated high-throughput framework to simulate an uncertainty-aware additive manufacturing (AM) process, characterize microstructural images, and extract meaningful features/descriptors. A kinetic Monte Carlo (KMC) based model of the AM process is used to simulate microstructural evolution for a diverse set of experimentally relevant processing conditions. We perform a parametric study to explore the relationship between microstructural features and processing conditions. Our results indicate that a many-to-one mapping can exist between processing conditions and typical descriptors; therefore, multiple descriptors are thus necessary to unambiguously represent microstructural images. Our work provides crucial quantitative and qualitative in-formation that would aid in the selection of features for microstructural images. Featurized microstructures could then be utilized to build data-driven models for predictive control of microstructures and thereby properties of additively manufactured components.

36 MATERIALS SCIENCE↗