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

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗
DOE OSTI2024

Classical dynamics of the antiferromagnetic Heisenberg spin ladder

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the S = 1/2 ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule. Here, we anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

Here, we develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems, including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input–output and the hierarchical equations of motion formalisms into a double hierarchy of density matrix equations. We show analytically that under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an N-photon Fock state input induces no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon with or without the influence of phonons. For short pulses (or, equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δ and an exciton–light coupling constant Γ. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyze the absorption probability in the opposite long pulse (narrow bandwidth) regime. We then derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective vs independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy equations.

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

Quantitative Determination of Total Cesium in Savannah River Site (SRS) Tank Waste by Inductively Coupled Plasma Emission Spectroscropy

Utilization of the significant and unique wavelength dynamic range (165 to 1100 nm) available on the Leeman Prodigy Inductively Coupled Plasma Emission Spectroscopy (ICP-ES) allows for quantification of cesium using the emission line 894.347 nm without the potential of Barium interference in nuclear waste samples. Samples are added to a Sodium Peroxide Fusion and Lithium Carbonate buffer in 2% nitric acid solution to increase sensitivity of the easily ionizable (<6 eV) cesium for quantitation to 5 pbb. Through direct comparison of analytical data results, the method is proven to be complementary to the Inductively Coupled Plasma Mass Spectrometry (ICP-MS) analysis (LOQ 100 ppt) in support of the Department of Energy’s (DOE) Nuclear Waste Mission.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Two-Component Multireference Restricted Active Space Configuration Interaction for the Computation of L-Edge X-ray Absorption Spectra

X-ray absorption spectroscopy is a powerful probe of local electronic and nuclear structure, providing insights into chemical processes. The theoretical prediction and interpretation of metal L-edge X-ray absorption spectra is complicated by both relativistic effects, including spin-orbit coupling, and by the multiconfigurational nature of the states involved. This work details an exact two-component multireference restricted active space configuration interaction (X2C-MRRASCI) scheme that uses an exact two-component state average complete active space self-consistent-field method (X2C-SA-CASSCF), which includes the spin-orbit coupling in a variational manner, for the accurate description of the electronic structure before using a restricted active space configuration interaction method to describe the core excited states of the X-ray spectrum. Benchmark calculations are presented for a series of iron containing complexes, with results showing key features of the spectrum being reproduced, including ligand to metal charge transfer and shake-up excitations.

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

On the convergence of an IEQ-based first-order semi-discrete scheme for the Beris-Edwards system

We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization Method (IEQ). The model consists of the Navier–Stokes equations for the fluid flow, coupled to the Q-tensor gradient flow describing the liquid crystal molecule alignment. By using the Invariant Energy Quadratization Method, we obtain a linearly implicit scheme, accelerating the computational speed. However, this introduces an auxiliary variable to replace the bulk potential energy and it is a priori unclear whether the reformulated system is equivalent to the Beris-Edward system. In this work, we prove stability properties of the scheme and show its convergence to a weak solution of the coupled liquid crystal system. We also demonstrate the equivalence of the reformulated and original systems in the weak sense.

Weber, Franziska↗
DOE OSTI2021

Photoporomechanics: An Experimental Method to Visualize the Effective Stress Field in Fluid-Filled Granular Media

Effective stress governs the mechanical behavior of porous media. In this study, we use photoelasticimetry to visualize the evolving effective stress field in fluid-filled granular media in processes that couple fluid flow and mechanical deformation. We refer to this experimental method as photoporomechanics. We develop a fabrication process to produce millimeter-scale residual-stress-free photoelastic spheres with high geometric accuracy. We use color to quantify the forces acting on the particles over a wide range of forces, while using light intensity for a small range of forces. We then provide an application of photoporomechanics to illustrate the evolution of effective stress during one-dimensional consolidation: a process by which the stresses caused by a sudden load are gradually transmitted through a fluid-filled granular pack as the fluid drains and excess pore pressures dissipate. We report our technique provides a powerful experimental model system to study the grain-scale underpinning of coupled solid-fluid processes in granular media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Multiphysics Meshfree Degradation Modeling of Energy Storage Materials with Kernel Enrichment

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. Reference: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

25 ENERGY STORAGE↗
DOE OSTI2021

Cyclic (Alkyl)(amino)carbene Ligands Enable Cu‐Catalyzed Markovnikov Protoboration and Protosilylation of Terminal Alkynes: A Versatile Portal to Functionalized Alkenes**

Abstract Regioselective hydrofunctionalization of alkynes represents a straightforward route to access alkenyl boronate and silane building blocks. In previously reported catalytic systems, high selectivity is achieved with a limited scope of substrates and/or reagents, with general solutions lacking. Herein, we describe a selective copper‐catalyzed Markovnikov hydrofunctionalization of terminal alkynes that is facilitated by strongly donating cyclic (alkyl)(amino)carbene (CAAC) ligands. Using this method, both alkyl‐ and aryl‐substituted alkynes are coupled with a variety of boryl and silyl reagents with high α‐selectivity. The reaction is scalable, and the products are versatile intermediates that can participate in various downstream transformations. Preliminary mechanistic experiments shed light on the role of CAAC ligands in this process.

Gao, Yang↗
DOE OSTI2024

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗
DOE OSTI2022

Dissociating the phononic, magnetic and electronic contributions to thermal conductivity: a computational study in alpha-iron

Computational tools to study thermodynamic properties of magnetic materials have, until recently, been limited to phenomenological modeling or to small domain sizes limiting our mechanistic understanding of thermal transport in ferromagnets. Herein, we study the interplay of phonon and magnetic spin contributions to the thermal conductivity in a-iron utilizing non-equilibrium molecular dynamics simulations. It was observed that the magnetic spin contribution to the total thermal conductivity exceeds lattice transport for temperatures up to two-thirds of the Curie temperature after which only strongly coupled magnon-phonon modes become active heat carriers. Characterizations of the phonon and magnon spectra give a detailed insight into the coupling between these heat carriers, and the temperature sensitivity of these coupled systems. Comparisons to both experiments and ab initio data support our inferred electronic thermal conductivity, supporting the coupled molecular dynamics/spin dynamics framework as a viable method to extend the predictive capability for magnetic material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

Friction stir welding (FSW) has found increased applications in automotive and aerospace industries due to its advantages of solid-state bonding, no fusion and melting, and versatility in various working conditions and material combinations. However, the relationship among processing parameters, material properties, and bonding extent and fidelity remains largely empirical, primarily because of the lack of the mechanistic understanding of the tool-workpiece frictional behavior that affects our subsequent understanding of microstructural evolution and interface bonding formation. While the tool-workpiece stick-slip condition is believed to dictate the resulting torque and heat generation rate during the welding process, it remains rare and elusive to conduct a quantitative experimental measurement of such interfacial field. On the other hand, numerical simulations based on Computational Fluid Dynamics (CFD) rely on ad hoc assumptions of interfacial pressure and shear-stress conditions, but predictions can only be validated via the medium- and far-range temperature field which is known to be insensitive to the interfacial frictional behavior. This work first presents a comparison among two CFD-based simulation methodologies and the Coupled Eulerian Lagrangian (CEL) model in finite element method, the last of which uses the Coulomb friction so that the stick-slip is naturally developed. Based on the Hill-Bower similarity relationship in the contact analysis, an analytical model is developed here to prove why a constant stick-slip fraction will be developed in the steady state, to correlate the stick-slip fraction to processing parameters such as the tool spin rate, and further to derive dimensionless functions for torque and heat-generation-rate predictions. Pros and cons of various numerical approaches in predicting stick-slip are discussed, and our analytical model has been found to agree well with our numerical simulation and literature experimental results. These analyses provide the critical strain-rate and temperature fields that are needed for the bonding analysis in our future work.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phonons and phase symmetries in bulk CrCl 3 from scattering measurements and theory

Phonon-derived behaviors are important indicators of novel phenomena in transition metal trihalides, including spin liquid behavior, two-dimensional magnetism, and spin-lattice coupling. However, phonons and their dependence on spin structure and excitations have not been adequately explored. In this work, we probe and critically examine the vibrational properties of the prototype ferromagnetic honeycomb lattice material CrCl 3 using inelastic neutron scattering and density functional theory. We demonstrate that magnetic and van der Waals interactions are essential to describing the structure and phonons in CrCl 3 ; however, the specific spin configuration is unimportant. Further, this provides context for understanding thermal transport measurements as governed by dynamical spin-lattice couplings. More importantly, we introduce an efficient dynamic method that exploits translational symmetries in large conventional unit cells that generates insights into phonon dispersions, interactions, and measured spectra in terms of quantum phase interference conditions. This work opens new avenues for understanding phonons in layered magnets and more generally in conventional cell geometries of a variety of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A dissolution model of alite coupling surface topography and ions transport under different hydrodynamics conditions at microscale

Portland cement is the most produced material in the world. The hydration process of cement consists of a group of complex chemical reactions. In order to investigate the mechanism of cement hydration, it is vital to study the hydration of each phase separately. An integrated model is proposed in this paper to simulate the dissolution of alite under different hydrodynamic conditions at microscale, coupling Kinetic Monte Carlo model (KMC), Lattice Boltzmann method (LBM) and diffusion boundary layer (DBL). The dissolution of alite is initialised with KMC. Two Multiple-relaxation-time (MRT) LB models are used to simulate the fluid flow and transport of ions, respectively. For solid-liquid interface, DBL is adapted to calculate the concentration gradient and dissolution flux. The model is validated with experiment from literature. The simulation results show good agreements with the results published in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-physics modeling of injected nanoparticles effect on remediation of CO 2 leakage through cracks

One potential risk in CO 2 sequestration is the leakage of carbon dioxide, which can result in contamination of underground water, creating potential threats to existing ecosystems. The common leakage pathway is through the pre-existing fractures or discontinuities within cement in the wellbore, incurred by the environmental conditions imposed on the cement. Injecting nanoparticles into pre-existing cracks is one of the most recently proposed ideas for mitigating fracture propagation in cement CO 2 sequestration. To demonstrate the feasibility of this new technology, a numerical approach was taken in this work, as it is challenging to investigate it in a laboratory setting. We proposed a coupled ALE (Arbitrary Lagrangian–Eulerian)–DEM (Discrete Element Method)–peridynamic modeling strategy within the LS-Dyna package to investigate the intertwined interaction among the CO 2 fluid flow, native fluid (e.g., brine), particle clusters, and cracks within the cement. The numerical results demonstrate that injected nanoparticles can effectively reduce the pressure exerted on the crack surface. Accordingly, the potential fracture propagation at the crack tip would be reduced compared to corresponding cases without nanoparticels as pressurized by fluid flow. This result verifies the effectiveness of proposed nanoparticle injection technology. Finally, using this established modeling strategy, the effect of filling particles on the fracture mitigation for different crack geometries (e.g. particle cluster patterns, aspect ratio of crack aperture and length) and CO 2 reservoir pressure are examined. The result shows that injected particles successfully reduce the fracture propagation in these scenarios.

42 ENGINEERING↗
DOE OSTI2024

Recovery of value-added compounds through fast pyrolysis of apple pomace hydrochar

The environmental challenges associated with food production can be addressed via the thermochemical upcycling of agro-industrial biomass. Two such methods, hydrothermal carbonization (HTC) and pyrolysis, can be coupled to first reduce the water content of wet biomass wastes by producing a hydrochar (HC) via HTC and then a bio-oil via pyrolysis of the HC. However, HTC of biomass results in the formation of secondary char (SC), an amorphous tar-like mixture resulting from organic compounds released into the aqueous phase that adsorb, recondense and polymerize on the parent biomass. This study investigated how HTC temperature impacts the formation of SC from apple pomace and the SC’s subsequent impact on fast pyrolysis products. HCs were produced at temperatures of 175°C, 200°C, and 250°C. Lower HTC temperatures favor the formation of biorefinery platform chemicals such as 5-hydroxymethylfurfural and levulinic acid, while higher temperatures result in increased lignin degradation products (i.e., phenolics). HCs were subjected to fast pyrolysis before and after SC extraction in two analytical pyrolysis instruments. Fast pyrolysis of HC produced compounds similar to those found in SC, but with variations in CO and CO 2 emissions. In conclusion, the combination of SC extraction and fast pyrolysis demonstrates promise for recovering value-added compounds from agro-industrial waste biomass while retaining a solid char for fuel and carbon management.

09 BIOMASS FUELS↗