Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling method”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An extended numerical manifold method for unsaturated soil-water interaction analysis at micro-scale

To investigate unsaturated soil-water interaction at micro-scale, this work extends the numerical manifold method (NMM) by incorporating a soil-water coupling model considering specific capillary water distribution and capillary force calculation. The soil skeleton is constructed by a soil skeleton generation algorithm with random polygons. To more realistically capture the interaction between soil grains and capillary water, a capillary mechanics-based geometric algorithm is proposed to iteratively calculate the capillary water distribution. The capillary forces corresponding to the capillary water distribution are calculated based on the Young-Laplace equation. The proposed capillary water solving framework is first verified by reproducing the soil-water characteristic curve and the capillary water distribution of an ideal contact-disk model against analytical solutions. To further validate the ability of the capillary water solving framework to predict hydraulic behavior of the real soil, a laboratory test on the Toyoura sand is reproduced numerically. Then an ideal direct shear test is performed to further validate the two-way soil-water coupling procedure, in which a comparison between the numerical and analytical results regarding the shear strength and matric suction is presented. Finally, microscopic hydraulic and compression tests are conducted on two soil specimens with the same porosity and mean grain diameter but different uniformity coefficients. The results elucidate that the extended method is a potential tool to explore unsaturated soil behaviors at micro-scale.

58 GEOSCIENCES↗
DOE OSTI2023

Tensile behavior of diffusion bonded AA6061 - AA6061 with variation in cooling method

Hot isostatically pressed AA6061 cladding is an important structural component of the high performance, Zr-laminated U-10Mo monolithic fuel system for the application in research and test reactors. In this study, the mechanical behavior of two diffusion bonded aluminum alloy, AA6061, was examined using tensile testing. Solid-to-solid diffusion bonding between two pieces of AA6061 was performed by isothermal annealing at 560 °C for 1.5 h, and diffusion couples were subsequently cooled via three different cooling methods: furnace cooling, air cooling, and water quenching. Dog-bone shaped tensile specimens, with 10 mm in gauge length (with diffusion bonded interface in the middle), and 1.5 × 1.5 mm 2 gauge cross-sections, were fabricated from the diffusion bonded AA6061 by electro-discharge machining. Yield strength (% EL at failure) of furnace cooled, air cooled and water quenched tensile specimens determined was 82 ~ 89 MPa (10 ~ 30%), 112 ~ 116 MPa (10 ~ 14%), and 149 ~ 164 MPa (10 ~ 17%), respectively. This variation in mechanical behavior was examined with cooling-rate dependent, concentrated precipitation of Mg 2 Si at the diffusion bonded interface, with due respect for mechanical properties of the AA6061 alloy that inherently vary as a function of cooling rate from 560 °C. Finite element analysis using ABAQUS was employed to augment experimental findings with the appropriate microstructural constituents and alloy properties. Results suggest that the strength is dominated by matrix/bulk properties of AA6061, while ductility is strongly influenced by the cooling method dependent presence of Mg 2 Si precipitates at the interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Higgs boson width and couplings at high energy muon colliders with forward muon detection

We propose a novel method using the ZZ-fusion channel and forward muon detection at high energy muon colliders to address the challenge of the Higgs coupling-width degeneracy. Our approach enables inclusive Higgs rate measurement to 0.75% at 10 TeV muon collider, breaking the coupling-width degeneracy. Results indicate the potential to refine Higgs coupling to subpercent levels and estimate its total width within (–0.41%, +2.1%). Key insights include the effectiveness of forward muon tagging in signal-background separation despite broad recoil mass distribution due to muon energy reconstruction and beam energy spread. The study emphasizes the significance of muon rapidity coverage up to |η⁡(μ)| < 6, enhancing measurement precision. Our findings highlight the unique capabilities of high energy lepton colliders for model-independent Higgs coupling determination and lay the groundwork for future advancements in muon collider technology and Higgs physics research.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Experimental and numerical investigation of sands and Geldart A biomass co-fluidization

This study investigated the fluidization of sands and small Geldart A biomass mixtures. The mixture fluidized like Geldart A type particles with a uniform bed expansion regime before bubbling. The video recorded color distance between pure sands and sands–biomass mixtures was used to estimate the sands–biomass mixing. The coarse-grained computational fluid dynamics–discrete element method with a hybrid drag model which couples the Syamlal–O'Brien drag and a filtered drag can capture the mixing while the simulation with Gidaspow drag predicted a segregated bed. The simulations were further validated with experimental measured pressure drops. The time averaged pressure drop equals the weight of the bed material, however, its fluctuation is about three times of the bed material fluctuation.

42 ENGINEERING↗
DOE OSTI2026

Simultaneous Determination of Halogens and Metals in Waste Plastic Pyrolysis Oil by Inductively Coupled Plasma Mass Spectrometry

A major barrier to integrating pyrolysis-derived oil into conventional refinery technology is the presence of impurities, particularly halogens and metals, that can deactivate catalysts. This study presents a novel, cost-effective approach for the simultaneous analysis of a subset of halogens, metals, nonmetals, and metalloids in complex, industrially relevant distilled pyrolysis oil samples, this was previously achievable only through multiple techniques. Excellent results were obtained using a widely accessible inductively coupled plasma mass spectrometry (ICP-MS) method with helium gas mode and standard laboratory consumables, enabling high-throughput analysis and efficient evaluation of adsorbent performance. Sample preparation is straightforward, requiring only dilution in a compatible matrix, and provides accurate and precise quantification, with 75%–137% spike recovery for Be, Ti, V, Cr, Fe, Ni, Co, Cu, As, Se, Mo, Cd, Sb, Tl, and Pb, and 62%–65% spike recovery for Cl and 109%–133% spike recovery for Br. Additionally, an enhanced version of the method using ICP-MS/MS and hydrogen gas is described, which has higher accuracy with 69%–112% spike recovery for Be, B, Ti, V, Fe, Co, Ni, Cu, As, Se, Mo, Cd, Sb, Tl, and Pb, with 85%–102% spike recovery for Cl and 63%–93% spike recovery for Br. Helium mode detection limits for industrially relevant elements (V, Fe, Ni, Cu, As, and Pb) are less than 1.7 µg/kg, and less than 0.2 mg/kg for Br and Cl. Furthermore, this methodology facilitates rapid systematic evaluation of adsorption capacities of materials under time-on-stream conditions and supports robust comparisons across diverse operating environments.

Lazarcik, James [University of Wisconsin-Madison, ↗
DOE OSTI2021

An adaptive granular representative volume element model with an evolutionary periodic boundary for hierarchical multiscale analysis

The hierarchical multiscale analysis normally utilizes a microscopic representative volume element (RVE) model to capture path/history-dependent macroscopic responses instead of using phenomenological constitutive models. However, for problems involving large deformation, the current RVE model used in geomechanics may lose representative properties due to the progressive distortion of the RVE box, unless a proper reinitialization is applied. This work develops an adaptive RVE model in conjunction with an evolutionary periodic boundary (EPB) algorithm for hierarchical multiscale analysis of granular materials undergoing large deformation based on a recent RVE model proposed for coupling molecular dynamics and the material point method. The proposed adaptive RVE model avoids the reinitialization of the RVE box that even undergoes extremely large shear deformation; meanwhile, it accounts for the deformation history of the RVE model and treats the interaction between boundary particles and other image particles in a more efficient way. Numerical examples with extremely large deformation are used to illustrate the adaptive granular RVE model enhanced by the proposed EPB algorithm. Furthermore, some key features of this new methodology are further discussed for clarification.

42 ENGINEERING↗
DOE OSTI2020

Near conformal perturbation theory in SYK type models

We present a systematic procedure to extract the dynamics of the low energy soft mode in SYK type models with a single energy scale J and emergent reparametrization symmetry in the IR. This is given in the framework of the perturbative scheme of arXiv:1608.07567 based on a specific (off-shell) breaking of conformal invariance in the UV, adjusted to yield the exact large-N saddle point. While this breaking term formally vanishes on-shell, it has a non-trivial effect on correlation functions and the effective action. In particular, it leads to the Schwarzian action with a specific coupling to bi-local matter. The method is applied to the evaluation of O(1) corrections to the correlation function of bi-locals. As a byproduct we confirm precise agreement with the explicit, symmetry breaking procedure. We provide a verification in the large q limit (Liouville theory), where the correlators can be calculated exactly at all length scales. In this case, our scheme illuminates how the enhanced O(J) and the subleading O(1) contributions originate from the Schwarzian dynamics of the soft mode and its interaction with h = 2 (bi-local) matter.

1/N expansion↗
DOE OSTI2020

GMRES with embedded ensemble propagation for the efficient solution of parametric linear systems in uncertainty quantification of computational models

In a previous work, embedded ensemble propagation was proposed to improve the efficiency of sampling-based uncertainty quantification methods of computational models on emerging computational architectures. It consists of simultaneously evaluating the model for a subset of samples together, instead of evaluating them individually. A first approach introduced to solve parametric linear systems with ensemble propagation is ensemble reduction. In Krylov methods for example, this reduction consists in coupling the samples together using an inner product that sums the sample contributions. Ensemble reduction has the advantages of being able to use optimized implementations of BLAS functions and having a stopping criterion which involves only one scalar. However, the reduction potentially decreases the rate of convergence due to the gathering of the spectra of the samples. In this paper, we investigate a second approach: ensemble propagation without ensemble reduction in the case of GMRES. This second approach solves each sample simultaneously but independently to improve the convergence compared to ensemble reduction. This raises two new issues which are solved in this paper: the fact that optimized implementations of BLAS functions cannot be used anymore and that ensemble divergence, whereby individual samples within an ensemble must follow different code execution paths, can occur. We tackle those issues by implementing a high-performing ensemble GEMV and by using masks. The proposed ensemble GEMV leads to a similar cost per GMRES iteration for both approaches, i.e. with and without reduction. For illustration, we study the performances of the new linear solver in the context of a mesh tying problem. Furthermore, this example demonstrates improved ensemble propagation speed-up without reduction.

BLAS↗
DOE OSTI2022

Crystal plasticity modeling of strain-induced martensitic transformations to predict strain rate and temperature sensitive behavior of 304 L steels: Applications to tension, compression, torsion, and impact

This paper advances crystallographically-based Olson-Cohen (direct γ → α’) and deformation mechanism (indirect γ→ε→α’) phase transformation models for predicting strain-induced austenite to martensite transformation. Here, the advanced transformation models enable predictions of not only strain-path sensitive, but also of strain-rate and temperature sensitive deformation of polycrystalline stainless steels (SSs). The deformation of constituent grains in SSs is modeled as a combination of anisotropic elasticity, crystallographic slip, and phase transformation, while the hardening is based on the evolution of dislocation density and explicit shifts in phase fractions. Such grain-scale deformation is implemented within the meso-scale elasto-plastic self-consistent (EPSC) homogenization model, which is coupled with the implicit finite element (FE) method to provide a constitutive response at each FE integration point for solving boundary value problems at the macro-scale. Parameters pertaining to the hardening and transformation models within FEEPSC are calibrated and validated on a suite of data including flow curves and phase fractions for monotonic compression, tension, and torsion as a function of strain-rate and temperature for wrought and additively manufactured (AM) SS304L. To illustrate the potential and accuracy of the integrated multi-level FE-EPSC simulation framework, geometry, mechanical response, phase fractions, and texture evolution are simulated during gas-gun impact deformation of a cylinder and quasi-static tension of a notched specimen made of AM SS304L. Details of the simulation framework, comparison between experimental and simulation results, and insights from the results are presented and discussed.

304L steels↗
DOE OSTI2021

CFD-DEM study of bubble properties in a cylindrical fluidized bed of Geldart Group D particles and comparison with prior MRI data

Numerical simulations of a freely bubbling cylindrical fluidized bed are carried out using a coupled computational fluid dynamics and discrete element method (CFD-DEM) model and compared to recent experimental data. The experiments were conducted using high-resolution and high-frequency magnetic resonance imaging providing high-fidelity data of the bubbling within a central 10 mm slice of the bed. Qualitatively, we nd more regular (i.e., less chaotic) structures in the simulated beds than observed experimentally. Quantitatively, however, the bubble diameter and number of bubbles as a function of height within the bed is predicted well by the base model. Unfortunately, the regularity in the simulations manifests as a considerable discrepancy in the speed of the (dense) emulsion phase. The simulated velocity probability distribution functions show an accumulation of low-speed regions and deficiency of high-speed regions. A simple parametric study of the base model is also carried out considering many of the most common CFD-DEM modeling parameters. It is found that the fluid grid size, geometry resolution, transfer kernel and drag law did not have a significant effect on bubble or particle dynamics.

42 ENGINEERING↗
DOE OSTI2023

Estimating Future Surface Water Availability Through an Integrated Climate‐Hydrology‐Management Modeling Framework at a Basin Scale Under CMIP6 Scenarios

Abstract Climate change and increasing water demand due to population growth pose serious threats to surface water availability. The biggest challenge in addressing these threats is the gap between climate science and water management practices. Local water planning often lacks the integration of climate change information, especially with regard to its impacts on surface water storage and evaporation as well as the associated uncertainties. Using Texas as an example, state and regional water planning relies on the use of reservoir “Firm Yield” (FY)—an important metric that quantifies surface water availability. However, this existing planning methodology does not account for the impacts of climate change on future inflows and on reservoir evaporation. To bridge this knowledge gap, an integrated climate‐hydrology‐management (CHM) modeling framework was developed, which is generally applicable to river basins with geographical, hydrological, and water right settings similar to those in Texas. The framework leverages the advantages of two modeling approaches—the Distributed Hydrology Soil Vegetation Model (DHSVM) and Water Availability Modeling (WAM). Additionally, the Double Bias Correction Constructed Analogues method is utilized to downscale and incorporate Coupled Model Intercomparison Project Phase 6 GCMs. Finally, the DHSVM simulated naturalized streamflow and reservoir evaporation rate are input to WAM to simulate reservoir FY. A new term—“Ratio of Firm Yield” (RFY)—is created to compare how much FY changes under different climate scenarios. The results indicate that climate change has a significant impact on surface water availability by increasing reservoir evaporation, altering the seasonal pattern of naturalized streamflow, and reducing FY.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulation of Microcapsule Transport in Fractured Media Using Coupled CFD‐DEM

Geothermal energy is sustainable and gaining momentum as a solution to energy crises and environmental issues. However, challenges like production temperature and thermal breakthrough can impact geothermal project efficiency. One innovative solution to alleviate the thermal breakthrough is to inject polymer-based materials that are encapsulated in microcapsules into fractures to modify fracture permeability and prevent preferential flow. In our study, we utilized a coupled computational fluid dynamics and discrete element method to simulate the transport of microcapsules under various scenarios controlled by microcapsule size, microcapsule concentration, and fracture roughness. For a smooth fracture, the results indicate that small microcapsules can travel through a smooth fracture regardless of their concentrations. Large microcapsules can transport through a smooth fracture when present in lower concentrations. However, medium and mixed-size microcapsules tend to cause the sealing of a smooth fracture, irrespective of their concentrations. For a rough fracture, the transport of microcapsules is complicated by their interactions with the rough fracture walls. The presence of two sealing positions in a rough fracture adds further complexity to this transport phenomenon. The size and concentration of microcapsules control one sealing location, while the rough fracture walls determine the other sealing location. The rough walls substantially affect microcapsule transport, rendering the role of microcapsule size and concentration less significant. The simulation results suggest that complex fracture surfaces significantly elevate the occurrence of sealing behavior. To mitigate sealing behavior within more complex fractures, it would be beneficial to use smaller and lower concentrations of microcapsules.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Multiscale computational understanding and growth of 2D materials: a review

Abstract The successful discovery and isolation of graphene in 2004, and the subsequent synthesis of layered semiconductors and heterostructures beyond graphene have led to the exploding field of two-dimensional (2D) materials that explore their growth, new atomic-scale physics, and potential device applications. This review aims to provide an overview of theoretical, computational, and machine learning methods and tools at multiple length and time scales, and discuss how they can be utilized to assist/guide the design and synthesis of 2D materials beyond graphene. We focus on three methods at different length and time scales as follows: (i) nanoscale atomistic simulations including density functional theory (DFT) calculations and molecular dynamics simulations employing empirical and reactive interatomic potentials; (ii) mesoscale methods such as phase-field method; and (iii) macroscale continuum approaches by coupling thermal and chemical transport equations. We discuss how machine learning can be combined with computation and experiments to understand the correlations between structures and properties of 2D materials, and to guide the discovery of new 2D materials. We will also provide an outlook for the applications of computational approaches to 2D materials synthesis and growth in general.

Momeni, Kasra↗
DOE OSTI2026

X-ray tomography of damage dynamics in advanced materials using a laser wakefield accelerator

Additively manufactured (AM) metals offer the potential for customizable, cost-effective components, but qualification and certification are crucial. Key to this process is understanding pore dynamics under stress, typically analyzed using micro-computed tomography. This study introduces laboratory-scale “betatron” x-rays from laser wakefield acceleration as a high-throughput alternative for x-ray tomography of advanced materials, such as AM AlSi10Mg alloys. Coupled with 3D finite element modeling, this method provides detailed insights into stress-porosity interactions. The approach delivers high-resolution scans, revealing that pore shape and local triaxiality significantly influence fracture dynamics, supporting advanced material characterization. This work also demonstrates the potential and versatility of laser-betatron x-ray μCT for generating large datasets to accelerate our understanding of the stochastic, process-specific nature of pore formation in AM alloys.

Senthilkumaran, Vigneshvar↗
DOE OSTI2022

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

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

Optical telecommunications-band clock based on neutral titanium atoms

Here, we propose an optical clock based on ultranarrow transitions in neutral titanium, which exhibit small blackbody radiation and quadratic Zeeman shifts and have wavelengths in the S-, C-, and L-telecommunications fiber bands, allowing for integration with robust laser technology. We calculate relevant properties using a high-precision relativistic hybrid method that combines configuration interaction and coupled-cluster approaches. To identify magic wavelengths, we have completed the largest-to-date direct dynamical polarizability calculations. Finally, we identify challenges that arise from magnetic dipole-dipole interactions and describe an approach to overcome them. A telecommunications-band atomic frequency standard will aid the deployment of optical clock networks and clock comparisons over long distances.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Field-induced magnetization of δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of magnetization, induced by an external magnetic field, for δ-phase plutonium. The fully relativistic electronic structure accounts for Zeeman splitting effects through a Hamiltonian that couples the magnetic field to both spin and orbital magnetic moments. The electronic-structure model is further improved by an extension to DFT in terms of an orbital-orbital coupling via the conventional orbital-polarization method, as has routinely been done for plutonium. The response to the applied magnetic field is shown to be weak in the DFT model, with induced magnetic moments of the order of in magnetic fields up to 30 T. These results are in accord with recent assessments from x-ray magnetic circular dichroism measurements in magnetic field on δ-plutonium. Somewhat surprisingly, a model assuming δ-plutonium to be absent magnetic moments shows stronger response to an external static magnetic field than models allowing for antiferromagnetism or magnetic disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗