Engineering Papers⌕ Search

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.

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

At least 433 records · Page 24

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations

Three-dimensional equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes: VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near axis. The importance of higher-accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier–Zernike basis also yields solutions with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences and allows for exponential convergence. Further, DESC can compute a solution with the same accuracy as VMEC in order-of-magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of geometric and magnetic boundary conditions on magnetic islands in 3D force-free ideal MHD equilibria

Three-dimensional boundary features and magnetic perturbations affect the magnetic topology of the force-free ideal MHD equilibria formed within a given volume. The use of the PSI-Tet 3D finite-element equilibrium code reveals the existence of magnetic islands inside such equilibria within boundaries with three-dimensional features in configurations that directly depend on the natures of those features. Here, the nature and limits of the influence of such features are explored, as is the influence of flux boundary conditions to qualitatively identify the structure of magnetic islands in Taylor state equilibria with three-dimensional boundary conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Properties of the neutron star crust: Quantifying and correlating uncertainties with improved nuclear physics

A compressible liquid-drop model (CLDM) is used to correlate uncertainties associated with the properties of the neutron star (NS) crust with theoretical estimates of the uncertainties associated with the equation of state (EOS) of homogeneous neutron and nuclear matter. For the latter, we employ recent calculations based on Hamiltonians constructed using chiral effective-field theory (χEFT). Fits to experimental nuclear masses are employed to constrain the CLDM further, and we find that they disfavor some of the χEFT Hamiltonians. The CLDM allows us to study the complex interplay between bulk, surface, curvature, and Coulomb contributions, and their impact on the NS crust. It also reveals how the curvature energy alters the correlation between the surface energy and the bulk symmetry energy. Here our analysis quantifies how the uncertainties associated with the EOS of homogeneous matter implies significant uncertainties for the composition of the crust, its proton fraction, and the volume fraction occupied by nuclei. We find that the finite-size effects impact the crust composition but have a negligible effect on the net isospin asymmetry of matter. The isospin asymmetry is largely determined by the bulk properties and the isospin dependence of the surface energy. The most significant uncertainties associated with matter properties in the densest regions of the crust, the precise location of the crust-core transition, are found to be strongly correlated with uncertainties associated with the Hamiltonians. By adopting a unified model to describe the crust and the core of NSs, we tighten the correlation between their global properties such as their mass-radius relationship, moment of inertia, crust thickness, and tidal deformability with uncertainties associated with the nuclear Hamiltonians.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experimental Assessment of Elastic Modulus vs. Relative Density for Stretch- and Bend-Dominated Lattices

This project was designed to study the possibility of using structural properties of lattices to replicate the material properties of certain hard to manufacture designs and use topology optimization to determine the lattice type and density required to mimic these properties. The integration of additively manufactured lattice structures with topology optimization highlights the need for well characterized mechanical properties and uncertainty analyses to insure these optimized structures respond as predicted. Stereolithographically printed octet and rhombic dodecahedron lattices were manufactured at 10%, 25% and 65% density by volume. As a separate task, yet integrated into this work, finite element analysis (FEA) was used to predict the printed lattice’s mechanical properties, which were then compared to our experimental results. After comparing FEA and experimentallymeasured elastic moduli, it was determined that the FEA provides highly reliable predictions for the modulus of these printed lattice structures. These lattices also exhibited greater tensile stiffness than that of the solid material demonstrating the flexibility that lattices provide to designing parts with designer structural properties. The accurate printing and reliable modeling of these lattices will enable topology optimization of complex parts from well-characterized rhombic dodecahedron and octet lattice structures of varying densities.

36 MATERIALS SCIENCE↗

Pore-scale visualization of natural hydrate-bearing sediments

Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.

Liu, Mengwei [Oak Ridge Institute for Science and ↗

Large-volume centimeter-wave cavities for axion searches

The scan rate of an axion haloscope is proportional to the square of the cavity volume. In this paper, a new class of thin-shell cavities are proposed to search for axionic dark matter. These cavities feature active volume much larger (>20×) than that of a conventional cylindrical haloscope, comparable quality factor Q, and a similar frequency tuning range. Furthermore, full 3D numerical finite-element analyses have been used to show that the TM 010 eigenmodes are singly polarized throughout the volume of the cavity and can facilitate axion-photon conversion in uniform magnetic field produced by a superconducting solenoid. To mitigate spurious mode crowding and volume losses due to localization, a pre-amplification binary summing network will be used for coupling. Because of the favorable frequency-scaling, the new cavities are most suitable for centimeter-wavelength (~ 10–100 GHz), corresponding to the promising post-inflation axion production window. In this frequency range, the tight machining tolerances required for high-Q thin-shell cavities are achievable with standard machining techniques for near-infrared mirrors.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards determination of the strong coupling $α_s(m_Z)$ from four-flavor lattice QCD using the continuous $β$-function method

The precise value of the strong coupling $α_s(m_{Z})$ at the $Z$-boson mass $m_{Z}$ is essential for high-energy phenomenology and precision tests of quantum chromodynamics (QCD). We present the status of a program targeting a $\sim 0.3\%$ determination of $α_s(m_{Z})$ using the renormalization group $β$-function in the infinite volume gradient flow scheme based on lattice QCD simulations of degenerate four-flavor highly improved staggered quark (HISQ) ensembles. In particular, we analyze both tree-level cutoff effects and finite-mass effects. We also outline the next steps of the analysis, including the infinite-volume and continuum extrapolations required for a precise determination of $α_s(m_Z)$.

Mandlecha, Yash (ORCID:000000020587962X)↗

Towards determination of the strong coupling $α_s(m_Z)$ from four-flavor lattice QCD using the continuous $β$-function method

The precise value of the strong coupling $α_s(m_{Z})$ at the $Z$-boson mass $m_{Z}$ is essential for high-energy phenomenology and precision tests of quantum chromodynamics (QCD). We present the status of a program targeting a $\sim 0.3\%$ determination of $α_s(m_{Z})$ using the renormalization group $β$-function in the infinite volume gradient flow scheme based on lattice QCD simulations of degenerate four-flavor highly improved staggered quark (HISQ) ensembles. In particular, we analyze both tree-level cutoff effects and finite-mass effects. We also outline the next steps of the analysis, including the infinite-volume and continuum extrapolations required for a precise determination of $α_s(m_Z)$.

Mandlecha, Yash [Michigan State U.; Michigan State↗

An Analytical Solution for the Initiation and Early Progression of Fretting Wear in Spherical Contacts

Abstract This article derives analytical solutions to calculate the wear volume at the initiation of fretting motion and its early progression over the first few oscillation cycles. The Archard-based model considers a deformable hemisphere that is contact with a deformable flat block. The material pairs investigated are special alloys, the Inconel 617/Incoloy 800H, and Inconel 617/Inconel 617. The analytical study begins with a unidirectional frictional sliding contact, where the local interfacial sliding distance and the nominal sliding distance at the initiation of gross slip are derived. The obtained analytical expressions for unidirectional sliding are then used to derive the corresponding wear volume for the initiation and early progression of gross slip and the wear volume for a general fretting cycle under elastic conditions. These analytical derivations are all verified by the finite element analysis (FEA). The FEA method and the analytical solutions render virtually identical results for both similar and dissimilar material pairs. The effects of plasticity on the wear volume under elastic–plastic conditions are also investigated. It is found that the fretting wear volumes obtained from the FEA simulations, which include plasticity, are close to those obtained from the analytical expressions for purely elastic regimes. All the results are presented in normalized forms, which can be easily generalized and applied to three-dimensional fretting wear of other material pairs.

Engineering↗

Optimization of quasi-axisymmetric stellarators with varied elongation

An optimization study of Quasi-Axisymmetric (QA) stellarators with varied elongation has been carried out using the optimization code STELLOPT. The starting point of our optimization is a previously obtained QA stellarator with 3 field periods and aspect ratio of 6. Here, a series of QA stellarators are obtained at zero plasma beta with varied elongation value ranging from 2.5 to 3.7. Good quasi-symmetry is kept when the elongation value is reduced from the original value of 3.7. The rotational transform profile and aspect ratio are kept fixed. The plasma volume is ether kept fixed or varied linearly with elongation. Furthermore, finite beta QA stellarators are considered. The corresponding bootstrap currents are calculated using the kinetic code SFINCS. A series of kink-stable QA stellarators are obtained via optimization with varied plasma beta up to 5% and self-consistent bootstrap current. This work demonstrates that good QA stellarators with finite beta and varied elongation exist that are stable to external kink modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coalescence, the thermal model and multi-fragmentation: the energy and volume dependence of light nuclei production in heavy ion collisions

Here, we present results of a phase space coalescence approach within the UrQMD transport and -hybrid model for a very wide range of beam energies from SIS to LHC. The coalescence model is able to qualitatively describe the whole range of experimental data with a fixed set of parameters. Some systematic deviations are observed for very low beam energies where the role of feed down from heavier nuclei and multi-fragmentation becomes relevant. The coalescence results are mostly very close to the thermal model fits. However, both the coalescence approach as well as thermal fits are struggling to simultaneously describe the triton multiplicities measured with the STAR and ALICE experiment. The double ratio of tp/d 2 , in the coalescence approach, is found to be essentially energy and centrality independent for collisions of heavy nuclei at beam energies of E lab > 10A GeV. On the other hand the clear scaling of the d/p 2 and t/p 3 ratios with the systems volume is broken for peripheral collisions, where a canonical treatment and finite size effects become more important.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NO x emissions under low-load gasoline compression ignition conditions

Here, in the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NO x ), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NO x sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NO x emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NO x and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NO x emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NO x emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NO x decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NO x emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Computational fluid dynamics↗

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite temperature properties of uranium mononitride

Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys, such as high melting temperature and high uranium density, and thermal conductivity, respectively. A better understanding of UN behavior at operating temperatures can be obtained from finite temperature data, such as elastic properties. To get this information, ab initio molecular dynamics (AIMD) simulations were performed at five different temperatures using constant volume (NVT) and constant pressure (NPT) ensembles. Initially, the performance of PBE functional in reproducing experimental crystallographic properties and magnetic ordering is assessed. The finite temperature phonon dispersions are calculated using NVT simulation results, which show a softening of the phonon modes with increasing temperature. The NPT results are used to obtain the thermal expansion of UN and finite temperature electronic properties. The calculated thermal expansion is compared with our measurements using neutron diffraction. Additionally, the temperature dependent elastic properties of UN are evaluated using the strain-stress method in AIMD simulations, indicating that UN becomes softer and more compressible with increasing temperature. Also, the calculated Young’s modulus slope is in very good agreement with the experiment. The finite temperature heat capacity and electronic thermal conductivity are calculated from AIMD simulations, which are in better agreement with the experiment than the heat capacity and thermal conductivity calculated using the structures relaxed at 0K. Finally, the thermal diffusivity from AIMD has opposite temperature dependence compared to experimental results, which we argued comes from the underestimated electronic thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Finite element modeling of biomass hopper flow

Hoppers are widely used biomass handling devices that channel bulk biomass from storage to subsequent handling equipment. Jenike’s longstanding approach, based on the Mohr-Coulomb model, has been successfully used to design hoppers handling cohesionless granular materials such as grains and other agricultural produces. However, designing a hopper to ensure reliable biomass flow is found to be challenging due to cohesion, irregular particle shape, and bulk material elastoplasticity. This study aims to address the biomass handling engineering challenge with alternative constitutive material models concerning the flow behavior of bulk solids. Finite element modeling is an approach that allows for implementation of different material models, whose underlying constitutive theories assist in investigating the origin and manifestation of bulk mechanical behavior of granular materials. This study focuses on the incipient gravity hopper flow of two types of biomass feedstocks, i.e., ground corn stover and Douglas fir wood. Three widely used constitutive material models, i.e., Mohr-Coulomb model, modified Cam-Clay model, and Drucker-Prager/Cap model, are implemented. Using the flow pattern represented by the volume of biomass exhibiting more than 7% of axial strain (Kamath and Puri, 1999), the finite element model predicts that the bulk corn stover particulate material forms an arch, which represents a hampered transition from the static state to the dynamic flow-state out of the hopper, whereas bulk Douglas fir wood particulate material develops a reliable mass flow pattern. A laboratory scale hopper was used to experimentally determine the biomass flow conditions, which were subsequently compared with the predicted onset of flow by a finite element model (FEM). The developed FEM was found to correctly predict the initiation of mass flow for the milled Douglas fir wood, whereas corn stover was predicted to establish a strong core flow suggesting an unreliable handling characteristic. This observation aligns with the reported poor handling of milled corn stover.

09 BIOMASS FUELS↗

Stringy canonical forms

Canonical forms of positive geometries play an important role in revealing hidden structures of scattering amplitudes, from amplituhedra to associahedra. In this paper, we introduce “stringy canonical forms”, which provide a natural definition and extension of canonical forms for general polytopes, deformed by a parameter α'. They are defined by real or complex integrals regulated with polynomials with exponents, and are meromorphic functions of the exponents, sharing various properties of string amplitudes. As α'→ 0, they reduce to the usual canonical form of a polytope given by the Minkowski sum of the Newton polytopes of the regulating polynomials, or equivalently the volume of the dual of this polytope, naturally determined by tropical functions. At finite α', they have simple poles corresponding to the facets of the polytope, with the residue on the pole given by the stringy canonical form of the facet. There is the remarkable connection between the α'→ 0 limit of tree-level string amplitudes, and scattering equations that appear when studying the α'→ ∞ limit. We show that there is a simple conceptual understanding of this phenomenon for any stringy canonical form: the saddle-point equations provide a diffeomorphism from the integration domain to the interior of the polytope, and thus the canonical form can be obtained as a pushforward via summing over saddle points. When the stringy canonical form is applied to the ABHY associahedron in kinematic space, it produces the usual Koba-Nielsen string integral, giving a direct path from particle to string amplitudes without an a priori reference to the string worldsheet. We also discuss a number of other examples, including stringy canonical forms for finite-type cluster algebras (with type A corresponding to usual string amplitudes), and other natural integrals over the positive Grassmannian.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗