Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Improving and Automating Building Model Data Exchange

There are many instances throughout a project’s lifecycle where there arises a need for quick and accurate risk assessment of building designs. For example, an unexpected design change during construction may necessitate structural engineers to perform a seismic risk assessment on analytical models of the updated building design using high fidelity structural analysis software, such as ANSYS or Abaqus. However, the efficiency of such workflows often depends upon the interoperability of architectural design software and structural analysis software. When the quality of this interoperability is lacking or even non-existent, the efficiency of virtual engineering workflows is hampered, which increases project costs. A McGraw Hill industry survey of professional users of Building Information Modeling (BIM) technologies found that there is high demand for BIM interoperability for structural analysis, but that the value/difficulty ratio is currently too low for practical use. There have been efforts by the academic community to facilitate model data exchange between the architectural design and structural analysis domains, but such solutions have not been widely adopted by industry, face technical challenges, and oftentimes are limited in applicability for users of various BIM software. Therefore, INL is developing capabilities to improve, automate, and generalize model data exchange between architectural BIM software (e.g., Revit) and structural analysis software (e.g., SAP2000, ANSYS). The goal is to help expedite and automate as much of the pre-processing step for creating analytical models in finite element analysis software as reasonably as possible. Such a "BIM-to-FEA" conversion tool should provide direct benefit to end-users through accuracy, automation, quick turn-around, and wide applicability. To generalize the application of this BIM-to-FEA conversion tool and increase its useability among the many different commercial BIM software currently used by industry, the program is being developed with the concept of openBIM. OpenBIM is the application of non-proprietary, open data standards that allow for BIM model data exchange in a format that is accessible, retainable, and useable for all users. The most widely used open, non-proprietary data exchange format for BIM is the Industry Foundation Classes (IFC) schema. IFC is developed by buildingSMART international and is ISO certified (ISO 16739-1:2018). The BIM-to-FEA conversion tool is being developed for compatibility with typical commercial building designs of steel framed structures. The tool is currently capable of importing architectural BIM data of framed building structures, recognizing and extracting the aspects of the model that are required for structural analysis, adjusting the connectivity of frame members, and finally exporting to an analytical model stored in the IFC format. The exported IFC analytical model can then be imported into various openBIM compliant software, such as SAP2000. Such capabilities have already been tested on commercial software, as shown above, and continue to be improved. Work is underway to test the conversion on various commercial BIM software, develop a user-friendly interface, incorporate the program into the broader DeepLynx data warehouse project being developed by INL, and to eventually open-source the tool for the benefit of the community. Future development of the tool envisions the ability for efficient iterative risk assessment of generative building designs, all within a workflow utilizing open-source tools. One such open-source tool will be MOOSE, an advanced finite element analysis tool developed at INL. The conversion tool will also branch out from typical commercial building designs and will aim to incorporate nuclear construction. The aim will be to convert both structural and non-structural components of nuclear facilities, such as curved concrete containment structures and piping systems, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Measurement of grain boundary strength of Inconel X-750 superalloy using in-situ micro-tensile testing techniques in FIB/SEM system

Grain boundaries (GBs), known as two-dimensional defects, are omnipresent in polycrystalline metallic alloys and thus influence a wide range of mechanical properties under different environmental conditions like irradiation and corrosion. Therefore, quantifying the strength of individual GBs is critical for understanding the degradation of mechanical properties of materials under different conditions. In this study we developed an efficient approach for the fabrication of micro-tensile specimens with a GB almost perpendicular to the tensile direction, which is expected to advance the development of individual GB tensile testing at micro or nanoscale in a wide scope of materials. Further, an in-situ cantilever micro-tensile testing method was developed and used to quantify the strength of a Σ3 GB in Inconel X-750 with the combination of finite element modeling. The average ultimate tensile strength (UTS) of a non-irradiated Σ3 GB is estimated at around 1.4 GPa, comparable to that of a neutron-irradiated Σ3 GB with a dose of ~1.5 dpa (1.3 GPa). Moreover, the in-situ push-to-pull micro-tensile testing technique developed in this work provides valuable insights into the high-angle GB deformation and fracture behavior. This method generates qualitatively similar ductility behavior before and after neutron irradiation as the bulk material testing. However, the ductility and UTS values obtained from this method are different from bulk measurements due to vastly different specimen dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Oscillating paramagnetic Meissner effect and Berezinskii-Kosterlitz-Thouless transition in underdoped Bi2Sr2CaCu2O8+δ

ABSTRACT Superconducting phase transitions in two dimensions lie beyond the description of the Ginzburg-Landau symmetry-breaking paradigm for three-dimensional superconductors. They are Berezinskii-Kosterlitz-Thouless (BKT) transitions of paired-electron condensate driven by the unbinding of topological excitations, i.e. vortices. The recently discovered monolayers of layered high-transition-temperature (${{{T}}}_{{\rm C}}$) cuprate superconductor Bi2Sr2CaCu2O8+δ (Bi2212) meant that this 2D superconductor promised to be ideal for the study of unconventional superconductivity. But inhomogeneity posed challenges for distinguishing BKT physics from charge correlations in this material. Here, we utilize the phase sensitivity of scanning superconducting quantum interference device microscopy susceptometry to image the local magnetic response of underdoped Bi2212 from the monolayer to the bulk throughout its phase transition. The monolayer segregates into domains with independent phases at elevated temperatures below ${{{T}}}_{{\rm C}}$. Within a single domain, we find that the susceptibility oscillates with flux between diamagnetism and paramagnetism in a Fraunhofer-like pattern up to ${{{T}}}_{{\rm C}}$. The finite modulation period, as well as the broadening of the peaks when approaching ${{{T}}}_{{\rm C}}$ from below, suggests well-defined vortices that are increasingly screened by the dissociation of vortex-antivortex plasma through a BKT transition. In the multilayers, the susceptibility oscillation differs in a small temperature regime below ${{{T}}}_{{\rm C}}$, consistent with a dimensional crossover led by interlayer coupling. Serving as strong evidence for BKT transition in the bulk, we observe a sharp jump in phase stiffness and paramagnetism at small fields just below ${{{T}}}_{{\rm C}}$. These results unify the superconducting phase transitions from the monolayer to the bulk underdoped Bi2212, and can be collectively referred to as the BKT transition with interlayer coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Finite-range pairing in nuclear density functional theory

Pairing correlations are ubiquitous in low-energy states of atomic nuclei. To incorporate them within nuclear density functional theory, often used for global computations of nuclear properties, pairing functionals that generate nucleonic pair densities and pairing fields are introduced. Many pairing functionals currently used can be traced back to zero-range nucleon-nucleon interactions. Unfortunately, such functionals are plagued by deficiencies that become apparent in large model spaces that contain unbound single-particle (continuum) states. In particular, the underlying computational schemes diverge as the single-particle space increases, and the results depend on how marginally occupied states are incorporated. These problems become more pronounced for pairing functionals that contain gradient-density dependence, such as in the Fayans functional. To remedy this, finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. Here, we show that a folding radius of about 1 fm offers the best compromise between quality and stability, and substantially reduces the pathological behavior in different numerical applications.

Nuclear density functional theory↗
DOE OSTI2021

Flexosensitive polarization vortices in thin ferroelectric films

In this work, we present a theoretical study on the influence of the flexoelectric coupling on the spatial distribution and temperature behavior of the spontaneous polarization for several types of stable domain structure in thin ferroelectric films, such as stripe domains and vortices. Finite-element modeling (FEM) for BaTiO 3 films and analytical calculations within the Landau-Ginzburg-Devonshire approach reveals that an out-of-plane polarization component can be very sensitive to the flexoelectric coupling for periodic quasi-2D stripe domains and 3D vortex-antivortex structures. However, the influence is rather different for these structures. The flexoelectric coupling increases significantly the amplitude of a small out-of-plane polarization component in the stripe domains, and the “up” or “down” direction of the component is defined by the sign of the flexoelectric coefficients. Concerning the vortex-antivortex pairs, their antivortices with in-plane antiсirculation have smooth wide dipolar cores through the entire film, whose shape and other features are almost insensitive to the coupling. The vortices with in-plane vorticity have spikelike cores with an out-of-plane quadrupolar moment induced by the flexoelectric coupling. The cores are located near the film–dead-layer interfaces. FEM results corroborated by analytical calculations prove that a change of the flexoelectric coefficient sign leads to a reorientation of the core axial polarization, making the flexosensitive 3D vortices similar to the recently introduced “flexons” in cylindrical nanoparticles. The relatively wide temperature range (from 200 to 400 K) of the flexosensitive vortices' existence gives us confidence that they can be observed experimentally in thin ferroelectric films by scanning-probe and nonlinear optical microscopy methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗
DOE OSTI2021

Martensitic transformation in superlattices of two non-transforming metals

A number of studies have shown that the coherent integration of a second non-transforming phase in martensitic materials results in unprecedented thermo-mechanical properties. We demonstrate the possibility of martensitic transformation by epitaxially combining two non-transforming materials, Mg (hexagonal closed packed) and MgLi (body centered cubic). Density functional theory calculations predict that MgLi/Mg superlattices undergo a transformation crystallographically equivalent to MgSc alloys and can be expected to result in finite temperature shape memory behavior. We predict the zero-temperature energy of the martensite phase to be 9.3 meV/atom lower than austenite for a 50–50 at. % MgLi/Mg nanolaminate, which is expected to result in a martensitic transformation temperature of approximately 150 K. This energy difference can be tuned between 4.9 and 18.1 meV/atom by varying the fraction of the pure Mg phase in the superlattice, potentially enabling the transformation from cryogenic to room and higher temperatures. While Li diffusion might limit the applications of the systems proposed, this work opens the door to a new approach to create martensitic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Distributed fiber-optic sensing in a subscale high-temperature superconducting dipole magnet

High-temperature superconductors, such as REBa2Cu3O7−x (REBCO, RE = rare earth), are becoming pivotal for high-field magnet technology for future circular colliders and compact fusion reactors. The U.S. Magnet Development Program, in collaboration with industry, is developing REBCO magnet technology using round conductors consisting of multiple REBCO tapes. For these multi-tape cables, traditional instrumentation, such as voltage taps and resistive strain gauges, become insufficient to help measure and understand the performance-limiting factors in these model magnets. Distributed fiber-optic sensing (DFOS) is a potential solution to address this challenge. Although DFOS is well established for various applications, measuring temperature and strain in high-temperature superconducting magnets is in its infancy. Here we report the detailed implementation and test results of DFOS based on Rayleigh scattering in a subscale canted cosθ (CCT) dipole magnet using high-temperature superconducting CORC® wires. We co-wound optical fibers in each layer of the CCT magnet and compared different types of commercial fibers and mold-release agents to reduce the power attenuation in the fibers. The DFOS allowed us to measure mechanical deformation and temperature along the conductor during tests at 77 and 4.2 K. The measured strain agreed quantitively with a finite-element mechanical model of the subscale magnet. Our results indicate that DFOS can effectively identify locations of strain and temperature changes, offering unique insight into magnet performance that can advance our understanding and development of the REBCO magnet technology for high-energy physics and fusion applications.

Luo, Linqing↗
DOE OSTI2025

Soil Carbon Saturation: What Do We Really Know?

Managing soils to increase organic carbon storage presents a potential opportunity to mitigate and adapt to global change challenges, while providing numerous co-benefits and ecosystem services. However, soils differ widely in their potential for carbon sequestration, and knowledge of biophysical limits to carbon accumulation may aid in informing priority regions. Consequently, there is great interest in assessing whether soils exhibit a maximum capacity for storing organic carbon, particularly within organo–mineral associations given the finite nature of reactive minerals in a soil. While the concept of soil carbon saturation has existed for over 25 years, recent studies have argued for and against its importance. Here, we summarize the conceptual understanding of soil carbon saturation at both micro- and macro-scales, define key terminology, and address common concerns and misconceptions. We review methods used to quantify soil carbon saturation, highlighting the theory and potential caveats of each approach. Critically, we explore the utility of the principles of soil carbon saturation for informing carbon accumulation, vulnerability to loss, and representations in process-based models. We highlight key knowledge gaps and propose next steps for furthering our mechanistic understanding of soil carbon saturation and its implications for soil management.

Environmental sciences↗
DOE OSTI2021

Shell-structure and asymmetry effects in level densities

Level density [Formula: see text] is derived for a nuclear system with a given energy [Formula: see text], neutron [Formula: see text], and proton [Formula: see text] particle numbers, within the semiclassical extended Thomas–Fermi and periodic-orbit theory beyond the Fermi-gas saddle-point method. We obtain [Formula: see text], where [Formula: see text] is the modified Bessel function of the entropy [Formula: see text], and [Formula: see text] is related to the number of integrals of motion, except for the energy [Formula: see text]. For small shell structure contribution one obtains within the micro–macroscopic approximation (MMA) the value of [Formula: see text] for [Formula: see text]. In the opposite case of much larger shell structure contributions one finds a larger value of [Formula: see text]. The MMA level density [Formula: see text] reaches the well-known Fermi gas asymptote for large excitation energies, and the finite micro-canonical limit for low excitation energies. Fitting the MMA [Formula: see text] to experimental data on a long isotope chain for low excitation energies, due mainly to the shell effects, one obtains results for the inverse level density parameter [Formula: see text], which differs significantly from that of neutron resonances.

Physics↗
DOE OSTI2025

H420 Imager Calibration Procedure: Field of View Scan

Gamma-ray imagers with coded apertures have a finite field of view (FOV) within which an image of a source can be generated. A “FOV scan” is a facet of imager calibration procedure that involves collecting specialized data for a range of source locations filling the FOV in order to understand and correct for differences in imager performance. Systematic distortions reveal relative displacements between internal imager components, which can be accounted for in data analysis. The success of this scan requires a stable imager position during the full measurement duration (~hours). This document summarizes the hardware and software tools required to conduct a FOV calibration scan of an H3D H420 Coded Aperture Gamma-Ray Imager.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Study on Fracture of Fiber-Reinforced Polymeric Composites Using Spiral Notch Torsion Test

Fiber-reinforced polymer-matrix composites have been widely applied in various areas due to their excellent mechanical properties. In this study, carbon fiber-reinforced polymer-matrix composites and glass fiber-reinforced polymer-matrix composites failure tests under quasi static torsion have been demonstrated by using spiral notch torsion test methodology. Three types of cylindrical specimens with different fiber orientations were prepared and tested. For both composites, after increase stage, A type exhibited limited stress variations before the sudden failure, whereas B type showed gradual failure with torque drops. Displacement contours, stress contours and energy release rates were acquired using the anisotropic elastic finite element models. Microscopic characterization exhibited that the failure of A type occurred along interfaces between fibers and matrices, whereas circular cracks in B type propagated alternatively between polymeric matrices and interfaces. Type C failure is similar to type A whereas with a curvature observed on the fracture surface.

Du, Fen↗
DOE OSTI2024

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Generalized Korn’s inequalities for piecewise $H¹$ and $H²$ vector fields

The purpose of this paper is to construct a new class of discrete generalized Korn’s inequalities for piecewise $H^1$ vector fields and piecewise $H^2$ vector fields in three-dimensional space. The resulting Korn’s inequalities are different from the standard Korn’s inequalities, as they involve the trace-free symmetric gradient operator, in place of the usual symmetric gradient operator. Furthermore, it is anticipated that the new generalized Korn’s inequalities will be useful for the analysis of a broad range of finite element methods, including mixed finite element methods and discontinuous Galerkin methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Anomalous thermodynamic properties of quantum critical superconductors

Recent high-precision measurements employing different experimental techniques have unveiled an anomalous peak in the doping dependence of the London penetration depth which is accompanied by anomalies in the heat capacity in iron-pnictide superconductors at the optimal composition associated with the hidden antiferromagnetic quantum critical point. Here, we argue that finite temperature effects can be a cause of observed features. Specifically we show that quantum critical magnetic fluctuations under superconducting dome can give rise to a nodal-like temperature dependence of both specific heat and magnetic penetration depth in a fully gapped superconductor. In the presence of line nodes in the superconducting gap fluctuations can lead to the significant renormalization of the relative slope of T-linear penetration depth which is steepest at the quantum critical point. The results we obtain are general and can be applied beyond the model we use.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Shear Is Not Always Simple: Rate-Dependent Effects of Flow Type on Granular Rheology

we report that despite there being an infinite variety of types of flow, most rheological studies focus on a single type such as simple shear. Using discrete element simulations, we explore bulk granular systems in a wide range of flow types at large strains and characterize invariants of the stress tensor for different inertial numbers and interparticle friction coefficients. We identify a strong dependence on the type of flow, which grows with increasing inertial number or friction. Standard models of yielding, repurposed to describe the dependence of the stress on flow type in steady-state flow and at finite rates, are compared with data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Validation and Parametric Investigations of an Internal Permanent Magnet Motor Using a Lumped Parameter Thermal Model

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multiphysics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a lumped parameter thermal network (LPTN) model of an Internal Permanent Magnet synchronous motor within the commercially available MOTOR-CAD modeling environment. The validation is based on temperature comparison with experimental data and with more detailed finite element analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations, and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% - for the highest temperature in the slot-winding and the end-winding - was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot-winding and end-winding. As for the LPTN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance.

42 ENGINEERING↗