SEARCH · Engineering Papers

Results for “Finite Volume”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2025

2-loop scattering on superstring and supermembrane in flat space

We consider the S-matrix of transverse scalar excitations on an infinite D = 10 GS superstring and D = 11 supermembrane in flat target space. We compute the 4-particle scattering amplitude in the 2-loop approximation and demonstrate that, like in the Nambu string case, the D = 10 GS string S-matrix does not contain non-trivial 2d UV divergences (UV pole not accompanied by terms with logarithms of momenta is an artifact of dimensional regularization). This is consistent with underlying integrability of this model which is maintaining by adding appropriate local counterterms. In the supermembrane case there are no 1-loop divergences but we find a genuine 2-loop UV pole. This demonstrates non-finiteness of the world-volume S-matrix of the M2 brane theory.

AdS-CFT Correspondence↗

DOE OSTI2020

Thermodynamics of spin-orbit-coupled bosons in two dimensions from the complex Langevin method

We investigate the thermal properties of interacting spin-orbit-coupled bosons with contact interactions in two spatial dimensions. To that end, we implement the complex Langevin method, motivated by the appearance of a sign problem, on a square lattice with periodic boundary conditions. We calculate the density equation of state nonperturbatively in a range of spin-orbit couplings and chemical potentials. In this work, our results show that mean-field solutions tend to underestimate the average density, especially for stronger values of the spin-orbit coupling. Additionally, the finite nature of the simulation volume induces the formation of pseudocondensates. These have been observed to be destroyed by the spin-orbit interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

An $\bar{F}$ Meshfree Treatment for Nearly Incompressible Materials

The computational modeling of nearly incompressible materials is a difficult task for many numerical methods, and even after several decades of investigation, it is still an active research area. This report seeks to address the treatment of incompressible materials in meshfree methods using a synergistic combination of two treatments. The first treatment is an $\bar{F}$ method, where the decomposed dilatational and deviatoric parts are calculated over different smoothing domains. The second treatment “activates” additional nodes throughout the domain to increase the flexibility of the model. We implement this synergistic combination in the context of the reproducing kernel particle method (RKPM) and present results for the Cook’s membrane benchmark problem. The results are compared with those using the composite tet10 finite element with a volume-averaged J formulation. We show that the combined treatment is an effective way to deal with nearly incompressible materials in a meshfree framework and compares well with other highly-effective treatments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

In situ effective pumping speed measurements as a tool for non-evaporable getter activation and saturation monitoring

We have developed method of in situ effective pumping speed measurements based on injecting finite gas pulse into a vacuum volume and subsequent recording of pressure response using a Residual Gas Analyzer (RGA). The pressure burst caused by injected gas pulse falls exponentially in time and the exponential coefficient of such pressure decay is proportional to total effective pumping speed of all vacuum pumps acting on vacuum chamber volume. The effective pumping speed can be extracted from the dependence of injected gas species pressure vs time recorded by an RGA. Non-Evaporable Getters (NEGs) are widely used in multiple ultra and extremely high vacuum devices and applications during the last several decades. Areas of NEG applications were recently expanded into vacuum devices with relatively high-pressure levels (up to 10 -7 Torr) by the invention of high capacity ZAO NEG which can be operated at elevated temperatures. The intrinsic property of all NEG-based vacuum pumps is the reduction of their pumping speed, called “NEG saturation”, after an extended period of operation. Typical time interval for such process can span from about a month and up to a few years depending on NEG type, vacuum level of NEG pump operation, and residual gas content. After that, the NEG pump will require a re-activation cycle to restore its pumping speed close to the initial value. For many applications it is difficult to realize what is the current pumping speed of NEG pump is and when the pump should be re-activated. A method for in situ effective pumping speed measurements developed in this study can be used as a tool for NEG-based pump activation and saturation monitoring. Application of such an approach to monitor pumping status of NEG pumps was utilized and will be continually used at Extended EBIS which was recently developed and commissioned for Relativistic Heavy Ion Collider (RHIC) and future Electron Ion Collider (EIC). The results obtained during commissioning are presented and discussed.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Recovering Gardner restacking with purely diffusive operations

The maximum particle kinetic energy that can be extracted from an initial six-dimensional phase space distribution motivates the concept of free or available energy. The free energy depends on the allowed operations that can be performed. A key concept underlying the theoretical treatment of plasmas is the Gardner free energy, where the exchange of the contents of equal phase volumes is allowed. Further, a second free energy concept is the diffusive free energy, in which the contents of volumes are instead averaged. For any finite discretization of phase space, the diffusive free energy is known to be less than the Gardner free energy. However, despite the apparent fundamental differences between these free energies, it is demonstrated here that the Gardner free energy may be recovered from the continuous limit of the diffusive free energy, leading to the surprise that macroscopic phase-space conservation can be achieved by ostensibly entropy-producing microscopic operations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

High Resolution Polar Kerr Effect Studies of CsV 3 Sb 5 : Tests for Time-Reversal Symmetry Breaking below the Charge-Order Transition

We report high resolution polar Kerr effect measurements on CsV 3 Sb 5 single crystals in search for signatures of spontaneous time reversal symmetry breaking below the charge order transition at T * = 94 K. Utilizing two different versions of zero-area loop Sagnac interferometers operating at 1550 nm wavelength, each with the fundamental attribute that without a time reversal symmetry breaking sample at its path, the interferometer is perfectly reciprocal, we find no observable Kerr effect to within the noise floor limit of the apparatus at 30 nanoradians. Simultaneous coherent reflection ratio measurements confirm the sharpness of the charge order transition in the same optical volume as the Kerr measurements. At finite magnetic field we observe a sharp onset of a diamagnetic shift in the Kerr signal at T *, which persists down to the lowest temperature without change in trend. Since 1550 nm is an energy that was shown to capture all features of the optical properties of the material that interact with the charge order transition, we are led to conclude that it is highly unlikely that time reversal symmetry is broken in the charge ordered state in CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Mind the crosscap: $τ$-scaling in non-orientable gravity and time-reversal-invariant systems

Spectral statistics of quantum chaotic systems are governed by random matrix universality. In many cases of interest, time-reversal symmetry selects the Gaussian Orthogonal Ensemble (GOE) as the relevant universality class. In holographic CFTs, this is mirrored by the presence of non-orientable geometries in the dual gravitational path integral. In this work, we analyze general properties of these matrix models and their gravitational counterparts. First, we develop a formalism to express the universal level statistics in the canonical ensemble for arbitrary spectral curves, leading to a topological expansion with finite radius of convergence in the late-time $τ$-scaling limit. Then, we focus on topological gravity and study topological recursion on the moduli space of non-orientable surfaces. We find that the Weil-Petersson volumes display non-analytic behaviour multiplying polynomials in the boundary lengths. The volumes give rise to wormholes with late-time divergences, in contrast with the orientable case, which is finite. We identify systematic cancellations among WP volumes implied by the consistency and finiteness of the $τ$-scaling limit. In particular, the cancellation of late-time divergences requires a nontrivial genus resummation. Working in the gravitational microcanonical ensemble, we derive and resum all orders of the topological expansion matching the GOE matrix model in the high-energy regime.

Chaotic Dynamics (nlin.CD)↗

DOE OSTI2020

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Prelude to Simulation Diagnostics in MOOSE: Diagnosing Meshes

As finite element simulations increase in complexity, so does the need for high quality meshes. Without proper diagnostics, issues pertaining to poor mesh quality or improper mesh construction may go undetected, causing simulations to fail unexpectedly. This paper details the new mesh diagnostic system available in MOOSE, which provides users with the tools to inspect mesh quality for both externally and internally generated meshes. This system can detect issues arising from element volume size, non-conformal elements, intersecting edges, missing sidesets, as well as other common mesh issues. As meshes increase in size and complexity, such issues become more difficult to diagnose through visual inspection alone.

99 - GENERAL AND MISCELLANEOUS↗

DOE OSTI2021

Local wave-number model for inhomogeneous two-fluid mixing

We analyze the local wave-number (LWN) model, a two-point spectral closure model for turbulence, as applied to the Rayleigh-Taylor (RT) instability, the flow induced by the relaxation of a statically-unstable density stratification. Model outcomes are validated against data from 3D simulations of the RT instability. In the first part of the study we consider the minimal model terms required to capture inhomogeneous mixing and show that this version, with suitable model coefficients, is sufficient to capture the evolution of important mean global quantities including mix-width, turbulent mass flux velocity, and Reynolds stress, if the start time is chosen such that the earliest transitions are avoided. However, this simple model does not permit the expected finite asymptote of the density-specific-volume covariance b. In the second part of the study, we investigate two forms for a source term for the evolution of the spectrum of density-specific-volume covariance for the LWN model. The first includes an empirically motivated calibration of the source to achieve the final asymptotic state of constant b. Here, the second form does not require calibration but, in conjunction with enhanced diffusion and drag captures the full evolution of all the dynamical quantities, namely, the mix-layer growth, turbulent mass-flux velocity, Reynolds stress, as well as the desired behavior of b.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Hoop Tensile Failure Stress of 500 µm-Diameter Ceramic Spheres

The tensile failure stress of a brittle-material sphere is important to the survival of the sphere itself but is particularly relevant to systems that employ such spheres for applications in which the sphere exterior is subjected to tension during system fabrication or operation. That maximum hoop or tangential stress can be measured by diametrally compressing a C-sphere specimen directly machined from such spheres. This combination of specimen geometry and test method was conceived at Oak Ridge National Laboratory around 20 years ago. Although the C-sphere geometry has been tested with larger spheres, the ease and efficacy of its testing are unknown for submillimeter sphere diameters. To consider this, C-sphere specimens were machined from yttria-stabilized zirconia (YSZ) spheres of ∼500 µm diameter. The maximum (hoop or tangential) tensile stress was determined as a function of compressive failure force and specimen dimensions using finite element analysis. Companion effective area and volume as a function of Weibull modulus were computed. Hoop tensile failure stresses were then measured and fractography conducted to examine flaw type and location and to estimate the mirror constant of this YSZ grade. This novel study demonstrates that hoop tensile failure stress can be determined with sub-millimeter brittle-material spheres.

Delia, Daniel [ORNL] (ORCID:0000000209023239)↗

DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗

DOE OSTI2024

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

DOE OSTI2022

Extreme variation in fatigue: Fatigue life prediction and dependence on build volume location in laser powder bed fusion of 17-4 stainless steel

Laser powder bed fusion (LPBF), a metal additive manufacturing technology, is well-suited for design optimization but fatigue life is limited by manufacturing defects. In this work, 17-4 stainless steel components were manufactured in densely populated build volumes, simulating at-scale LPBF production. Tests revealed extreme variability in fatigue life data, analyzed via rigorous statistical tools. The El-Haddad model, modified for finite-life, enabled defect-based life prediction. Specimen location within the build volume correlated to life, which was heteroscedastic. Investigating defect concentration over the build volume explained typical life and scatter. Finally, these findings argue for qualification approaches which acknowledge high material lot variability.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effect of copper content on the tensile elongation of Al–Cu–Mn–Zr alloys: Experiments and finite element simulations

Microstructures of cast aluminum alloys used in automotive engine applications often consist of intermetallic particles that can impact the tensile elongation of these alloys. Here in this paper, we investigate the effect of intermetallic grain boundary particles on tensile elongation by fabricating a series of alloys with Cu content varying between 6.0 - 9.0 wt% in cast Al–Cu–Mn–Zr (ACMZ) type compositions. The tensile elongation of as-aged ACMZ alloys decreases monotonically with increase in Cu content. While the microstructure within the grains and yield stress of the alloy remains invariant with Cu content, the decrease in tensile elongation correlates well with increase in the size and volume fraction of grain boundary particles. Microstructural observations are combined with finite element simulations to explain the trend in tensile elongation with changing Cu content. Crack initiation is found to occur by brittle fracture of the grain boundary particles. Increase in particle size promotes crack initiation by reduction in size dependent particle fracture strength. Lower inter-particle spacing at higher particle volume fraction further facilitates crack initiation by increasing stress within the particles caused by the interaction between stress fields of neighboring particles. Increase in particle volume fraction also accelerates crack propagation through the formation of macro shear zones in the microstructure. The increase in Cu content of cast ACMZ alloys, therefore, decreases tensile elongation by promoting both crack initiation and crack propagation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

DOE OSTI2020

Evaluation of Detritiation Strategies for Mitigating Tritium Releases

Tritium confinement is performed using different barriers to minimize releases to the environment. Primary confinement is usually performed by process piping and components, secondary confinement typically by inerted gloveboxes connected to a continuously operating tritium stripper system. Some facility designs employ tertiary tritium confinement by using process room confinement with a tritium stripper system activated after an accident to further mitigate tritium releases to the environment. Secondary confinement of some tritium systems such as radiation-hardened or shielded-cell enclosures are not easily achieved using designs for inerted glovebox. These enclosures can be challenging to seal for secondary tritium confinement due to heating, ventilation, and air conditioning (HVAC) needs for temperature control, as well as feedthroughs for equipment and instrumentation, are potential pathways for tritium leaks or loses. Additional tritium leak paths can be created in secondary confinement enclosures due to a Design Basis Accident (DBA) seismic event. These new, larger tritium leak rates from the DBA event are usually not part of the stripper system design basis so all tritium released to the secondary confinement enclosure is assumed released to the environment. Detritiation strategies for shielded enclosures and/or gloveboxes with significant leak rates are necessary to minimize off-site radiological dose consequences. To address this, the relationship between the confinement system leak rate and stripper system (recirculating) flow rate was evaluated to meet prescribed maximum allowable environmental tritium emissions from the event. The simple analysis described assumes a constant leak rate from the confinement system (i.e. shielded enclosure or glovebox) while a recirculating stripper system strips tritium from the system– a competition between tritium release and tritium recovery. At a high level, the resulting expression summarizes the relationship between system leak rate (F L ) and stripper flow rate (F S ) relative to allowable tritium release (Q A ) and initial tritium release (Q 0 ): $\frac{F_L}{F_s}$ ≈ $\frac{Q_A}{Q_0}$. This report derives the relationship between these parameters, examines the impact of finite stripping times followed by purging of the stripped volume. The analysis provides example results for up to 30 gram tritium releases: the maximum tritium inventory of a Hazard Category III (Department of Energy) Nuclear Facility. The detritiation model presented represents a high-level analysis of tritium recovery versus losses for “leaky” confinement enclosures such as shielded cells after large tritium releases. Key parameters for doing the analyses are the system leak rate, stripper flow rate, initial tritium release, and allowable tritium release values. The analyses show a proportional decrease in releases with reduction in initial tritium release but a non-linear increase in tritium recovery with increased stripper flow rate or reduced leak-rate. The analyses also apply to glovebox confinement systems which could have significant increases in leak rates after a DBA. The analysis recognizes but does not include tritium absorption followed by re-emission from the walls of the confinement volume – an analysis which could be pursued in future studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗