Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERPOLATION”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermodynamics of monoclinic and tetragonal hafnium dioxide (HfO 2 ) at ambient pressure

Accurate and precise thermodynamic models of oxide compounds and their phases are important for calculating the phase stability of oxide materials. We develop and use a coupled quantum mechanical and molecular dynamics approach to create thermodynamic models of hafnia (HfO2) polymorphs from 0 K to 3000 K at ambient pressure. The approach is based on the quasi-harmonic approximation below the Debye temperature and on ab-initio molecular dynamics calculations above the Debye temperature to predict constant pressure heat capacities (C-p). A Bayesian model provides interpolated values between these regimes. As a case study, we develop thermodynamic models of monoclinic and tetragonal HfO2 polymorphs. The predicted heat capacities are in excellent agreement with experiment, and the predicted temperature of the monoclinic to tetragonal phase transition (2173 K) is in good agreement with the experimental value (2078K). These results provide a comprehensive and accurate thermodynamic model of the monoclinic and tetragonal phases of hafnia on a broad range of temperatures and can serve as input to CALPHAD assessment of multi-component hafnia-based phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2019

An improved multislope MUSCL scheme for solving shallow water equations on unstructured grids

In this paper we describe an improved vector manipulation multislope monotone upstream-centred scheme for conservation laws (MUSCL) reconstruction for solving the shallow water equations on unstructured grids. This improved MUSCL reconstruction method includes a bigger stencil for the interpolation and saves time for determining the geometric relations compared to the original vector manipulation method, so it is computationally more efficient and straightforward to implement. Four examples involving an analytical solution, laboratory experiments and field-scale measurements are used to test the performance of the proposed scheme. It has been proven that the proposed scheme can provide comparable accuracy and higher efficiency compared to the original vector manipulation method. With the increasing of the number of cells, the advantage of the proposed scheme becomes more apparent.

42 ENGINEERING↗
DOE OSTI2021

Empirical estimation of densities in NaCl-KCl-UCl 3 and NaCl-KCl-YCl 3 molten salts using Redlich-Kister expansion

Densities of molten KCl-NaCl-UCl 3 and KCl-NaCl-YCl 3 ternary systems have been estimated using a multidimensional Redlich-Kister model. Temperature and composition dependent Redlich-Kister functions have been used to generate binary interaction parameters in the outlined ternary salt systems. These binary interactions have been used in the extrapolation to ternary system densities. The results of the density extrapolations by Muggianu interpolation scheme provide agreement within 2–3% for the NaCl-KCl-YCl 3 liquids and 11% in NaCl-KCl-UCl 3 liquids compared to the available experimental data. Modeling NaCl-KCl-UCl 3 molten phase density with a ternary interaction parameter improved the agreement within 4%. Thermophysical modeling used in this study has shown promising results for use in other material properties, such as viscosity, thermal conductivity, and heat capacity of the molten salts. Lastly, the outlined modeling method applied in these specific molten salt ternaries can be used for quaternary or higher multicomponent molten salt systems.

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

High pressure experimental investigation of clinopyroxene dissolution in a K-basaltic melt

Dissolution of clinopyroxene (cpx) in a K-basaltic melt from the Campi Flegrei Volcanic District (Italy) has been investigated through dissolution and dissolution-crystallization experiments at pressure of 0.8 GPa, superliquidus temperature of 1350 °C, and dwell times between 0.5 and 1 h. The obtained dissolution rates range from 7.9·10 –6 cm s –1 to 6.1·10 –6 cm s –1 as a function of dwell time. In the dissolution-crystallization experiment (1300 °C; 0.8 GPa; 2 h), the formation of overgrowth rims accompanied by new cpx crystals suggests that the injection of recharge magmas in basaltic reservoir may lead to inverse or oscillatory zonation. The interaction between cpx and K-basaltic melt at ~1570 °C was studied by in situ radiography using synchrotron radiation combined with the Paris-Edinburgh press. Here, this resulted in cpx resorption to occur depending on the temperature conditions with respect to the liquidus temperature of the cpx (T cpxL ). The calculated cpx dissolution rates are ~5·10 –3 cm s –1 at T ≤ T cpxL and ~ 3·10 –2 cm s –1 at T ≥ T cpxL . The role of crystal dissolution in the estimation of magma residence times has been also tested for a natural magmatic system by interpolating the dissolution rates (~10 –5 –10 –6 cm s –1 ) with the textural data of cpx phenocrysts from the Agnano-Monte Spina pyroclastic deposit at Campi Flegrei caldera (Campanian region, Italy). Results from calculations indicate that the time required for partial or complete resorption of phenocrysts varies from ~0.5 to ~40 h, and that the effect of crystal dissolution may be relevant to estimate magma residence times whether significant dissolution occurs during magma mixing processes.

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

A velocity/stress mixed stabilized nodal finite element for elastodynamics: Analysis and computations with strongly and weakly enforced boundary conditions

A new nodal mixed finite element is proposed for the simulation of linear elastodynamics and wave propagation problems in time domain. Our method is based on equal-order interpolation discrete spaces for both the velocity (or displacement) and stress (or strain) tensor variables. The mixed form is derived using either the velocity/stress or velocity/strain pair of unknowns, the latter being instrumental in extensions of the method to nonlinear mechanics. The proposed approach works equally well on hexahedral or tetrahedral grids and, for this reason, it is suitable for time-domain engineering applications in complex geometry. The peculiarity of the proposed approach is the use of the rate form of the stress update equation, which yields a set of governing equations with the structure of a non-dissipative space/time Friedrichs’ system. We complement standard traction boundary conditions for the stress with strongly and weakly enforced boundary conditions for the velocity (or displacement). Weakly enforced boundary conditions are particularly suitable when considering complex geometrical shapes, because they do not require dedicated data structures for the imposition of the boundary degrees of freedom, but, rather, they utilize the structure of the variational formulation. Here we also show how the framework of weakly enforced boundary conditions can be used to develop variational forms for multi-domain simulations of heterogeneous media. A complete analysis including stability and convergence proofs is included, in the case of a space–time variational approach. A series of computational tests are used to demonstrate and verify the performance of the proposed approach.

42 ENGINEERING↗
DOE OSTI2019

A configurational force for adaptive re-meshing of gradient-enhanced poromechanics problems with history-dependent variables

Here, we introduce a mesh-adaption framework that employs a multi-physical configurational force and Lie algebra to capture multiphysical responses of fluid-infiltrating geological materials while maintaining the efficiency of the computational models. To resolve sharp gradients of both displacement and pore pressure, we introduce an energy-estimate-free re-meshing criterion by extending the configurational force theory to consider the energy dissipation due to the fluid diffusion and the gradient-dependent plastic flow. To establish new equilibria after remeshing, the local tensorial history-dependent variables at the integration points are first decomposed into spectral forms. Then, the principal values and directions are projected onto smooth fields interpolated by the basis function of the finite element space via the Lie-algebra mapping. Our numerical results indicate that this Lie algebra operator in general leads to a new trial state closer to the equilibrium than the ones obtained from the tensor component mapping approach. A new configurational force for dissipative fluid-infiltrating porous materials that exhibit gradient-dependent plastic flow is introduced such that the remeshing may accommodate the need to resolve the sharp pressure gradient as well as the strain localization. The predicted responses are found to be not influenced by the mesh size due to the micromorphic regularization, while the adaptive meshing enables us to capture the width of deformation bands without the necessity of employing fine mesh everywhere in the domain.

42 ENGINEERING↗
DOE OSTI2020

Simultaneous material, shape and topology optimization

Using three design fields we develop an optimization environment that can simultaneously optimize material, shape and topology. We use the implicit representation of the boundaries with level-set functions that define the shape and topology. Differentiable R-functions allow us to combine these shapes and topology descriptions with Boolean operations. Additionally, we incorporate design dependent-stiffness materials with another design field. Notably, this framework accommodates design dependent loads, has the ability to introduce holes, and ensures the satisfaction of optimality criteria. It builds upon the fictitious domain, ersatz material, material interpolation and level-set methods. Additionally, it also borrows from parameterized density-based topology optimization methods. Since analytical sensitivities can be computed, we use efficient nonlinear programming algorithms to update the design instead of the Hamilton–Jacobi’s scheme of level-set methods. We illustrate the features of our framework by designing a cantilever beam with octet truss microlattice, a dam with design-dependent loads, and a composite clevis plate.

42 ENGINEERING↗
DOE OSTI2021

Anisotropic elastoplastic response of double-porosity media

Here we present a continuum framework to simulate fluid flow through anisotropic elastoplastic media with double porosity. Two effective stress measures σ' and δ" emerge from the thermodynamic formulation, which are energy-conjugate to the elastic and plastic components of strain, respectively. Both effective stress measures can be expressed as a combination of the total Cauchy stress σ and the average pore pressure $\bar{p}$ in the two pore scales. In the effective stress for elasticity, $\bar{p}$ is scaled with a rank-2 Biot tensor, whereas the effective stress for plasticity follows the Terzaghi form in which $\bar{p}$ is scaled by the Kronecker delta. The Biot tensor and storage coefficients are derived as functions of elasticity parameters and porosities. A mixed finite element formulation is introduced to discretize the domain and solve initial boundary-value problems. A stabilization scheme is employed on equal-order interpolation for both displacement and pressure fields throughout the entire range of drainage responses. Numerical simulations reproduce the hydromechanical response of Opalinus shale in one-dimensional consolidation tests throughout the range of primary and secondary consolidation under different external loads. Numerical simulations of the consolidation of a rectangular domain subjected to a strip load demonstrate the efficacy of the proposed stabilization scheme, as well as illustrate the impacts of stress history, mass transfer, and different pore systems on the hydromechanical response.

42 ENGINEERING↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
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 OSTI2023

Modular machine learning-based elastoplasticity: Generalization in the context of limited data

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

42 ENGINEERING↗
DOE OSTI2023

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗
DOE OSTI2023

Reliable extrapolation of deep neural operators informed by physics or sparse observations

Deep neural operators can learn nonlinear mappings between infinite-dimensional function spaces via deep neural networks. As promising surrogate solvers of partial differential equations (PDEs) for real-time prediction, deep neural operators such as deep operator networks (DeepONets) provide a new simulation paradigm in science and engineering. Pure data-driven neural operators and deep learning models, in general, are usually limited to interpolation scenarios, where new predictions utilize inputs within the support of the training set. However, in the inference stage of real-world applications, the input may lie outside the support, i.e., extrapolation is required, which may result to large errors and unavoidable failure of deep learning models. Here, we address this challenge of extrapolation for deep neural operators. First, we systematically investigate the extrapolation behavior of DeepONets by quantifying the extrapolation complexity, via the 2-Wasserstein distance between two function spaces and propose a new strategy of bias–variance trade-off for extrapolation with respect to model capacity. Subsequently, we develop a complete workflow, including extrapolation determination, and we propose five reliable learning methods that guarantee a safe prediction under extrapolation by requiring additional information—the governing PDEs of the system or sparse new observations. The proposed methods are based on either fine-tuning a pre-trained DeepONet or multifidelity learning. We demonstrate the effectiveness of the proposed framework for various types of parametric PDEs. Furthermore, our systematic comparisons provide practical guidelines for selecting a proper extrapolation method depending on the available information, desired accuracy, and required inference speed.

42 ENGINEERING↗
DOE OSTI2023

A multifidelity deep operator network approach to closure for multiscale systems

Projection-based reduced order models (PROMs) have shown promise in representing the behavior of multiscale systems using a small set of generalized (or latent) variables. Despite their success, PROMs can be susceptible to inaccuracies, even instabilities, due to the improper accounting of the interaction between the resolved and unresolved scales of the multiscale system (known as the closure problem). In the current work, we interpret closure as a multifidelity problem and use a multifidelity deep operator network (DeepONet) framework to address it. In addition, to enhance the stability and/or accuracy of the multifidelity-based closure, we employ the recently developed ``in-the-loop" training approach from the literature on coupling physics and machine learning models. The resulting approach is tested on shock advection for the one-dimensional viscous Burgers equation and vortex merging for the two-dimensional Navier-Stokes equations. Furthermore, the numerical experiments show significant improvement of the predictive ability of the closure-corrected PROM over the un-corrected one both in the interpolative and the extrapolative regimes.

42 ENGINEERING↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Data-driven reduced-order models for port-Hamiltonian systems with operator inference

Hamiltonian operator inference has been developed in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. The method constructs a low-dimensional model using only data and knowledge of the functional form of the Hamiltonian. The resulting ROMs preserve the intrinsic structure of the system, ensuring that the mechanical and physical properties of the system are maintained. In this work, we extend this approach to port-Hamiltonian systems, which generalize Hamiltonian systems by including energy dissipation, external input, and output. Based on snapshots of the system’s state and output, together with the information about the functional form of the Hamiltonian, reduced operators are inferred through optimization and are then used to construct data-driven ROMs. To further alleviate the complexity of evaluating nonlinear terms in the ROMs, a hyper-reduction method via discrete empirical interpolation is applied. Accordingly, we derive error estimates for the ROM approximations of the state and output. Lastly, we demonstrate the structure preservation, as well as the accuracy of the proposed port-Hamiltonian operator inference framework, through numerical experiments on a linear mass–spring-damper problem and a nonlinear Toda lattice problem.

97 MATHEMATICS AND COMPUTING↗