Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “representative volume elements”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Loading Capacity and Dilute Nitric Acid Rinse of Diglycolamide Resin for the Recovery of Transplutonium and Rare Earth Elements from Mark-18A Targets

N,N,N′,N′-tetraoctyldiglycolamide (TODGA) as a resin (“DGA resin”) produced by Eichrom Technologies will be used by Savannah River National Laboratory for the indiscriminate extraction of trivalent actinides and rare earth elements with the intent of recovering Cm and Am from dissolved irradiated 242 Pu Mark-18A targets in 7 to 9 M nitric acid. The extracted constituents will be recovered as an oxide by direct thermal decomposition of the loaded resin followed by calcination of the resultant residue. The characteristics of DGA resin with non-radiological feed simulant representative of the anticipated feed in the Mark-18A process, with Sm and Nd as surrogates for Cm and Am, respectively, including breakthrough point and saturation capacity were evaluated in this work. Additionally, this work examined the losses from the loaded resin by rinsing the resin bed with dilute acid to reduce the nitrate concentration in the resin bed prior to thermal decomposition operations to improve the safety posture of the process. A resin loading profile was developed, and the resin was determined to have a trivalent metal saturation capacity of 74 μmol/mL resin under the experimental conditions evaluated. Following a wash of the loaded resin bed with fresh 8 M HNO 3 , the trivalent metals retained was reduced to 68 μmol/mL resin, which represents the practical capacity of the resin for the Mark-18A process. Lighter lanthanides breakthrough the resin well before the saturation capacity is reached. Rinsing the saturated resin bed with 0.26 M HNO 3 was found to result in a rapid reduction in retention of rare earth elements by the resin. After 2.8 bed volumes of dilute acid rinse, the mean resin bed free acid concentration was reduced to 0.28 M and 3.1 bed volumes of dilute acid rinse resulted in a reduction of the cumulative rare earth element retention to 52 μmol/mL of resin.

Transplutonium separations↗
DOE OSTI2025

High-Q Superconducting Lumped-Element Resonators for Low-Mass Axion Searches

Low-frequency superconducting lumped-element resonators have recently attracted significant attention in the context of axion dark matter searches. Here we present the design and implementation of a fixed-frequency superconducting resonator operating near $250~\mathrm{kHz}$, possessing an inductor volume of $\sim 1$ liter and achieving an unloaded quality factor $Q \approx 2.1\times10^{6}$. This resonator represents a significant improvement over the state of the art and informs the design of searches for low-mass axions.

FOS: Physical sciences↗
DOE OSTI2022

An energy-based study of the embedded element method for explicit dynamics

Abstract The embedded finite element technique provides a unique approach for modeling of fiber-reinforced composites. Meshing fibers as distinct bundles represented by truss elements embedded in a matrix material mesh allows for the assignment of more specific material properties for each component rather than homogenization of all of the properties. However, the implementations of the embedded element technique available in commercial software do not replace the material of the matrix elements with the material of the embedded elements. This causes a redundancy in the volume calculation of the overlapping meshes leading to artificially increased stiffness and mass. This paper investigates the consequences in the energy calculations of an explicit dynamic model due to this redundancy. A method for the correction of the edundancy within a finite element code is suggested which removes extra energy and is shown to be effective at correcting the energy calculations for large amounts of redundant volume.

Martin, Valerie A.↗
DOE OSTI2022

An Immersed Finite Element Lagrangian-Eulerian Code-Coupling Framework

This report presents an assessment of immersed Eulerian-Lagrangian code-coupling techniques suitable for use in a broad range of mechanics applications. The coupling algorithm is based on an immersed finite element method that considers the Lagrangian and Eulerian overlap regions in the overall variational formulation. In this report the basic formulation details are presented followed by various aspects of the code-coupling algorithm using OpenIFEM as the Lagrangian/coupling framework. A series of representative test cases that illustrate the code-coupling algorithm are discussed. The current work provides an in-depth investigation into the immersed finite element method for the purposes of providing a rigorous coupling technique that is minimally invasive in the respective Eulerian and Lagrangian codes. A number of extensions to the base immersed finite element method have been examined. These extension include nodal and quadrature-based indicator functions, a Lagrangian volume-fraction calculation in regions of overlap, and the use of penalty constraints between the Lagrangian and Eulerian domains. A unique MPI-based coupling strategy that retains the independent MPI structure of each code has been demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

Here the flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear-stress ratio μ are relatively insensitive to system size for a wide range of inertial numbers I. Second-order effects in strain rate, such as normal stress differences, require large system sizes to accurately extract meaningful results, notably a nonmonotonic dependence in the first normal stress difference with strain rate. The rheological response represented by the μ(I) scalar model works well at describing the lower-order aspects of the rheology, except near the quasistatic limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems N≤1000 near the quasistatic limit can lead to the arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasistatic and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number, and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasistatic flow.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A new hybrid gadolinium nanoparticles-loaded polymeric material for neutron detection in rare event searches

Experiments aimed at direct searches for WIMP dark matter require highly effective reduction of backgrounds and control of any residual radioactive contamination. In particular, neutrons interacting with atomic nuclei represent an important class of backgrounds due to the expected similarity of a WIMP-nucleon interaction, so that such experiments often feature a dedicated neutron detector surrounding the active target volume. In the context of the development of DarkSide-20k detector at INFN Gran Sasso National Laboratory (LNGS), several R&D projects were conceived and developed for the creation of a new hybrid material rich in both hydrogen and gadolinium nuclei to be employed as an essential element of the neutron detector. Thanks to its very high cross-section for neutron capture, gadolinium is one of the most widely used elements in neutron detectors, while the hydrogen-rich material is instrumental in efficiently moderating the neutrons. In this paper results from one of the R&Ds are presented. In this effort the new hybrid material was obtained as a poly(methyl methacrylate) (PMMA) matrix, loaded with gadolinium oxide in the form of nanoparticles. We describe its realization, including all phases of design, purification, construction, characterization, and determination of mechanical properties of the new material.

Dark Matter detectors (WIMPs, axions, etc.)↗
DOE OSTI2025

Isomerization kinetics of azobenzene crosslinked to a liquid crystal polymer network

Azobenzene is an important chromophore having numerous applications that are a direct implication of its ability to change shape upon irradiation with a suitable wavelength. Here, in this study, photoisomerization kinetics of a diacrylate azobenzene molecule (A3MA) cross-linked to a liquid crystal polymer network (LCN) is investigated. Isomerization experiments were conducted on planar LCN thin films with varying cross-linking densities for a fixed azobenzene concentration. The experiments were conducted at approximately 26 °C, which is significantly below each film’s glass transition temperature. The forward trans-cis and reverse cis-trans isomerization slows down in crosslinked azobenzene liquid crystal polymer networks (ALCNs) as compared to a solution and an uncrosslinked polymer mixture. The observed slowdown is due to the reduction in the free volume and increased steric hindrances with an increase in crosslinking. A single-term first-order kinetics rate equation is known to describe the isomerization kinetics in solutions or uncrosslinked mixtures while the same form fails to represent the reaction kinetics in moderate to highly cross-linked ALCNs. Our results for the first time indicate that a two-term first-order reaction rate is an accurate way to describe the forward trans-cis and reverse thermal cis-trans isomerization in crosslinked azobenzene liquid crystal polymer networks. In addition, the effect of film thickness and incident intensities on the reaction kinetics is investigated. Finally, the two-term isomerization kinetics is incorporated into a finite element photo-chemo-mechanical model to reveal the consequent slowdown in the actuation response of ALCN films. The current findings will be of particular interest to theoreticians and experimentalists alike, for it will help predict and understand accurately the responses of azobenzene-based coatings, and actuators and may have potential applications in underwater soft robotics.

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

Critical mineral source potential from oil & gas produced waters in the United States

The volume of produced water, a by-product of oil & gas operations and other energy processes, has been growing across the United States (U.S.) along with the need to manage or recycle this wastewater. Produced water contains many naturally occurring elements of varying concentrations, including critical minerals which are essential to the clean energy transition. However, the current understanding of critical mineral concentrations in produced water and the associated volumes across the U.S. is limited. This study has assessed available databases and literature to gain insight into the presence and concentration of five high priority critical minerals, namely cobalt, lithium, magnesium, manganese, and nickel. The U.S. Geological Survey's National Produced Waters Geochemical Database was the main data source used for determining average critical mineral concentrations in produced water from the major oil and gas reservoirs in the U.S. The volumes of produced water for these major reservoirs were coupled with these concentrations to provide insights into where critical minerals are likely to have high abundance and therefore more recovery options. The analysis indicated the highest recovery potential for lithium and magnesium from produced water in the Permian basin and the Marcellus shale region. However, these assessments should be considered conservative due to the limited availability of reliable concentration data. Finally, it is expected more critical mineral recovery options could emerge with comprehensive characterization data from more recent and representative sources of produced water.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2022

Comparison of Fission Product Release Predictions using PARFUME and BISON with Results from the AGR-3/4 Irradiation Experiment

The PARFUME (PARticle Fuel ModEl) fuel performance modeling code and the BISON nuclear fuel performance application built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library were used to predict the fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the third and fourth irradiation experiment of the Advanced Gas Reactor (AGR-3/4) Fuel Development and Qualification Program. The fuel performance modeling codes PARFUME and BISON modeled the AGR-3/4 irradiation experiment using the fuel compact time-averaged volume averaged (TAVA) daily temperatures for a total irradiation duration of 369.1 effective full power days (EFPD) to predict the release fraction of the fission product silver (Ag-110m) from a representative TRISO-coated fuel particle from AGR-3/4 compacts. Post-irradiation examination (PIE) measurements provided data on the release of these fission products in the compacts outside of the silicon carbide (SIC) layer. The PARFUME and BISON results were then compared to the silver release measured from compact gamma scanning. The results showed good agreement between PARFUME and BISON but both codes under-predicted the silver release fraction for all the compacts. In addition, BISON was used to model and predict the fission product concentration radial profile outside of the compacts in capsules’ inner and outer rings. These rings were either comprised of matrix and/or structural graphite. To obtain the concentration profiles of silver, cesium, and strontium, a sorption isotherm model was developed in BISON to capture the effects of fission product transport across the gaps between the concentric rings. The general shape of the concentration radial profiles as calculated by BISON were similar in the inner ring (IR) but varied in the outer ring (OR) depending on the fission product of interest or capsule temperature. Using this methodology and model, BISON now has the capability to aid in developing new fission product diffusion coefficients for matrix or structural graphite materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Level-set topology optimization with PDE generated conformal meshes

This paper presents a level-set topology optimization approach that uses conformal meshes for the analysis of the displacement field. The structure’s boundary is represented by the iso-contour of a level-set field discretized on a fixed background design mesh. The conformal mesh is updated for each design iteration via a PDE based mesh morphing process that identifies the set of facets in the background mesh that are homeomorphic to the boundary and relaxes the homeomorphic mesh to conform to the structure’s boundary and ensure high element quality. The conformal mesh allows for a more accurate computation of the response versus density and some level-set based methods which interpolate material properties using the volume fraction. Numerical examples illustrate the proposed approach by optimizing linear-elastic two- and three-dimensional structures, wherein insight into the performance of the mesh morphing process is provided. The examples also highlight the scalability of the approach.

42 ENGINEERING↗
DOE OSTI2020

Microstructure Generation via Generative Adversarial Network for Heterogeneous, Topologically Complex 3D Materials

Using a large-scale, experimentally captured 3D microstructure dataset, we implement the generative adversarial network (GAN) framework to learn and generate 3D microstructures of solid oxide fuel cell electrodes. The generated microstructures are visually, statistically, and topologically realistic, with distributions of microstructural parameters, including volume fraction, particle size, surface area, tortuosity, and triple phase boundary density, being highly similar to those of the original microstructure. These results are compared and contrasted with those from an established, grain-based generation algorithm (DREAM.3D). Importantly, simulations of electrochemical performance, using a locally resolved finite element model, demonstrate that the GAN generated microstructures closely match the performance distribution of the original, while DREAM.3D leads to significant differences. Finally, the ability of the generative machine learning model to recreate microstructures with high fidelity suggests that the essence of complex microstructures may be captured and represented in a compact and manipulatable form.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Science Capsule: Towards Sharing and Reproducibility of Scientific Workflows

Workflows are increasingly processing large volumes of data from scientific instruments, experiments and sensors. These workflows often consist of complex data processing and analysis steps that might include a diverse ecosystem of tools and also often involve human-in-the-loop steps. Sharing and reproducing these workflows with collaborators and the larger community is critical but hard to do without the entire context of the workflow including user notes and execution environment. In this paper, we describe Science Capsule, which is a framework to capture, share, and reproduce scientific workflows. Science Capsule captures, manages and represents both computational and human elements of a workflow. It automatically captures and processes events associated with the execution and data life cycle of workflows, and lets users add other types and forms of scientific artifacts. Science Capsule also allows users to create `workflow snapshots' that keep track of the different versions of a workflow and their lineage, allowing scientists to incrementally share and extend workflows between users. Our results show that Science Capsule is capable of processing and organizing events in near real-time for high-throughput experimental and data analysis workflows without incurring any significant performance overheads.

Ghoshal, Devarshi↗
DOE OSTI2023

Permanent Magnets Featuring Heavy Main Group Elements for Magnetic Anisotropy

Permanent magnets are the functional component of electric motors and generators found in numerous renewable energy applications. To improve energy conversion in such applications, we require fundamentally new magnets that generate higher magnetic flux per volume while retaining the properties conferred by rare-earth elements incorporated into current technologies. We hypothesize that by engendering a covalent interaction between two elements, we can access a new regime of magnetic materials where the two components of a magnetic moment—spin and orbital angular momentum—come from two separate atoms to form a complete magnetic moment. Our previous research utilized high-pressure conditions to discover two new candidate materials ideal for assessing this hypothesis. The first material, FeBi 2 , enables the study of an unprecedented solid state metal-metal bonding interaction. The second, MnBi 2 represents the second member of the promising Mn–Bi family known for its magnetic properties. Importantly, MnBi 2 is isostructural to FeBi 2 . Together these chemically simple but magnetically rich materials provide an elegant platform for elucidating fundamental design principals of magnetic anisotropy while inspiring the synthesis of new magnetic materials. More generally, solid-state chemistry remains a synthetic black box. To fully harness the potential of such new materials, it is vital to create a window into both the structure and the properties of our high-pressure materials. One missing area of importance to energy science is high-pressure magnetometry. To this end, we assessed the magnetic structure of these materials through two objectives: (1) we performed magnetometry measurements at high-pressures, including performing high-pressure X-ray Magnetic Circular Dichroism (XMCD) experiments on these materials that are only synthesized at high-pressure and cannot be recovered to ambient conditions and (2) we targeted the recovery of these materials using dynamic compression approaches, starting with laser shock techniques. Future directions include expanding our focus to ternary phase space, employing high-throughput computational approaches to direct our search for magnetic materials across the Periodic Table. This award was originated at Northwestern, and reopened as a new grant at MIT. This is the close out for Northwestern.

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

Impacts of Lanthanide Ion Complexation on the Radiation Robustness of Diglycolamide Extractants

Nuclear power reactors represent a clean and reliable source of baseload energy. However, nuclear reactors produce used nuclear fuel (UNF), which, if not recycled/reprocessed, must be disposed of as high-level radioactive waste, necessitating long-term storage options. UNF contains approximately a third of the periodic table, including the radioactive minor actinides (MA) americium and curium. Removal of these MA from UNF would greatly reduce the radiotoxic burden, volume, and cost of a storage facility or geological repository. Thus, a significant global effort has been devoted to the design of reprocessing strategies for the clean separation of the MA from UNF. As MA separations are difficult, due to the presence of lanthanide elements in UNF and their remarkably similar chemical properties, many extractants have been studied. The diglycolamide (DGA) class of extractants are promising for this separation as they exhibit high distribution coefficients for the MAs, high metal loading capacities, and are robust in highly acidic radiation environments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Characterization of Incremental Markings in the Sagittal Otolith of the Pacific Sardine ( Sardinops sagax ) Using Different Imaging Modalities

Teleost fish possess calcium carbonate otoliths located in separate chambers (utriculus, sacculus, and lagena) of their membranous labyrinth. This study analyzed the surface topography of the sagittal otolith of the Pacific sardine (Sardinops sagax) and the daily and annual increments in these otoliths. The otolith surface, characterized by laser scanning confocal microscopy for the first time, consisted of a system of prominent ridges and valleys (grooves), but it is unclear whether these structures are functional or represent time-resolving markings reflecting growth periodicity. Within the first-year volume, daily increments, each consisting of an incremental (more mineralized) and a discontinuous (less mineralized) zone, were resolved by optical microscopy and backscattered electron (BSE) imaging in the scanning electron microscope (SEM). Daily growth increments could, however, not be resolved in volumes formed after the first year, presumably because otolith growth markedly slows down and spacing of incremental markings narrows in older fish. Throughout otolith growth, the crystalline network continues across the discontinuous zones. Fluorochrome labeling provides additional information on growth after the first year. Compared with optical and BSE imaging, synchrotron microComputed Tomography of intact otoliths (with 0.69 µm volume elements) was less able to resolve daily increments; X-ray phase contrast reconstructions provided more detail than reconstructions with absorption contrast. Future research directions are proposed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Influence of cementation on the yield surface of rocks numerically determined from digital microstructures

Digital Rock Physics has reached a level of maturity on the characterisation of primary properties that depend on the microstructure – such as porosity, permeability or elastic moduli – by numerically solving field equations on μCT scan images of rock. After small deformations or at depth though, most rocks eventually reach their limit of elasticity and the complementary plastic properties are needed to describe the full mechanical behaviour. Currently, determination of a rock’s yield surface from its microstructure is often restricted to semi-analytical criteria derived by limit analysis or numerical simulations performed on idealised geometries. Such simplification lacks representativeness, particularly for processes that affect directly the pore-grain interface such as the cementation phenomenon, happening during diagenesis. Eventually, only direct numerical simulation of elasto-plasticity performed on digitalised microstructures can be used to assess the strength of different cemented materials and its evolution with the alteration of the microstructure. In this study, we provide a comprehensive parametric study on the impact of cementation on rock strength for real microstructures of cemented granular materials. Compared to most previous studies, the whole yield surface is determined numerically (using Finite Element Method) in order to assess the influence of cementation for different stress-paths. The previously known tendency of rock to strengthen with increasing cementation volume is verified. New results on the influence of cement property namely Young’s modulus, friction and cohesion on the rock’s yield surface are explored. The envelopes obtained are compared to the ones obtained by experimental data and existing models. Finally, the framework presented in this study showcases the wider possibility of determining any rock’s or porous material’s yield surface from its microstructure.

36 MATERIALS SCIENCE↗