Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

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

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗
DOE OSTI2022

SSR COST – COmposite Structural Technologies for advanced reactor deployment (Final Scientific/Technical Report)

Moltex is designing a Generation 4 or Advanced nuclear reactor, the molten salt fast spectrum SSR-W. The design intent is to reduce the cost and complexity of the nuclear systems by relying on inherent safety characteristics, instead of engineered systems. One such characteristic as a natural tendency for power to decrease as temperature increases. Another characteristic that improves natural safety is the fact that the working fluids (molten salts) are at atmospheric pressure at operating temperatures. This significantly reduces the need for pressure boundary components and systems. Moltex is highly motivated to develop simple and economic reactors to reduce the carbon intensity and improve energy security of global electricity production. The COST project investigated the opportunities for reducing construction cost of advanced reactors with the adoption of Rigid Elastic Composite Civil works (RECC). The construction design and methods used in currently operational light water power plants would generally be suboptimal for the SSR-W since there is no requirement for civil structures to withstand high pressure accidents such as hydrogen combustion or steam explosion. In parallel to downselecting the most promising RECC-type technology, accident scenarios representing the largest credible challenge to the SSR-W civil and construction design were modelled. The capability of an SSR-W reactor constructed with and without the selected construction technology was evaluated and the consequential licensing challenges identified. A construction schedule and cost estimate of the SSR-W with and without the selected RECC construction technology was developed. Once the above data was collected, then Moltex was able to review the safety performance, license-ability and cost/schedule benefit of applying the preferred construction technology. The CAPEX/NPV benefit was quantified, and the licensing risks and opportunities were documented. Finally, the CAPEX/NPV benefit of the advanced construction technology was modest due to the suitability of traditional stick-built construction for the SSR-W design. However, previously unidentified benefits of the advanced construction technology were revealed during the studies that could credibly de-risk certain aspects of the design to be licensed. Since licensing (duration and delays) are very significant in new build project costs, the cost benefits may be significant. Further studies on the risk reduction of the construction technology are foreseen.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Combined TDLAS and chemiluminescence imaging in a flat flame burner operated with NH3/H2 blends

Ammonia is viewed as a viable hydrogen carrier due to favorable storage and transport characteristics. While it can be re-converted to hydrogen at point-of-use via thermal catalytic cracking, direct utilization in combustion systems can result in reduced costs and improved efficiency. A major barrier to this approach is the low flammability and potential for high nitrogen oxide emissions, driven by fuel-bound nitrogen and complex kinetic pathways. While a number of kinetic mechanisms currently exist for simulating ammonia combustion, a major need continues to be direct information about species profiles in easy-to-model systems capable of isolating chemical kinetics from multi-dimensional fluid dynamic effects. This paper reports on recent species measurements made in a flat flame burner using a combined tunable-diode-laser-absorption-spectroscopy (TDLAS) and chemiluminescence imaging approach. Three flame conditions were included representative of NH3/H2 blends, partially cracked NH3 (inc. N2), and 100% NH3 with enhanced air (30% oxygen). Two NIR distributed feedback (DFB) laser diodes were used to determine H2O concentration and temperature (via ratio thermometry) at various radial and axial positions, after which an inverse Abel transform was used to infer centerline values. Companion images were collected using a Princeton Instrument PI-MAX intensified camera equipped with a 105mm UV lens and multiple filter sets targeting OH*, NH*, and NH2* emission. Results were compared to companion Cantera burner-stabilized flame simulations using various kinetic mechanisms. A recent mechanism including excited species chemistry was also included, to investigate whether excited and ground state profiles exhibit significant differences.

ammonia combustion↗
DOE OSTI2020

Thallium salen fluorescent tracers

Thallium salens can be synthesized by reacting thallium alkoxide with a salen ligand. As examples of the invention, the dinuclear complexes Tl2-(salo-But) and Tl2-(saloPh-But) were synthesized by the reaction of thallium ethoxide with (H2-salo-But) or (H2-saloPh-But). These thallium salens may have applications as fluorescent tracers (or taggants) for subterranean fluid flows.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A matrix-free approach for finite-strain hyperelastic problems using geometric multigrid

This paper investigates matrix-free algorithms for problems in quasi-static finite-strain hyperelasticity. Iterative solvers with matrix-free operator evaluation have emerged as an attractive alternative to sparse matrices in the fluid dynamics and wave propagation communities because they significantly reduce the memory traffic, the limiting factor in classical finite element solvers. Specifically, we study different matrix-free realizations of the finite element tangent operator and determine whether generalized methods of incorporating complex constitutive behavior might be feasible. Furthermore, in order to improve the convergence behavior of iterative solvers, we also propose a method by which to construct level tangent operators and employ them to define a geometric multigrid preconditioner. Additionally, the performance of the matrix-free operator and the geometric multigrid preconditioner is compared to the matrix-based implementation with an algebraic multigrid (AMG) preconditioner on a single node for a representative numerical example of a heterogeneous hyperelastic material in two and three dimensions. We find that matrix-free methods for finite-strain solid mechanics are very promising, outperforming linear matrix-based schemes by two to five times, and that it is possible to develop numerically efficient implementations that are independent of the hyperelastic constitutive law.

42 ENGINEERING↗
DOE OSTI2023

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Asynchronous GPU-based DEM solver embedded in commercial CFD software with polyhedral mesh support

A novel graphical processing unit-based discrete element method solver is introduced to improve stability, performance, and provide seamless integration into commercial or open-source computational fluid dynamics software. A key innovation is eliminating a need for network communication between solvers, which was previously required for cross-platform coupling. This is accomplished by a direct coupling method that employs dynamic-linked libraries. Furthermore, the solver optimizes memory usage by streamlining the particle-cell search algorithm by eliminating the cells' searching grid. This ensures the solver is compatible with a wide range of mesh types, providing high geometric flexibility. The approach simplifies the simulation process by directly incorporating computational fluid dynamics mesh information into the discrete element method solver. The performance analysis indicates about sixteen times boost in computational speed compared to benchmark central processing unit-based solvers. Finally, the solver's compatibility with polyhedral meshes, a vital advantage for complex geometries, is tested against a referenced study regarding the simulation of an immersed-tube fluidized bed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D 2 O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH ~2. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. In this work, the aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Mesoporous material↗
DOE OSTI2021

Pressure drop in a prototypical 3D magnetohydrodynamic flow across contraction of a fusion blanket manifold

Predicting 3D magnetohydrodynamic (MHD) pressure losses associated with liquid metal flows in complex geometry ducts is required to design breeding blankets for fusion reactors. Of such components, manifolds exhibit major pressure losses. Recent correlations [T. Rhodeset al. Magnetohydrodynamic pressure drop and flow balancing of liquid metal flow in a prototypic fusion blanket manifold. Phys. Fluids. 2018; 30: 057101.] demonstrate promise for predicting pressure drops in electrically insulated inlet manifolds. In the present work, the extend to which these correlations can be applied to outlet manifolds, which feature sudden contractions, is investigated in both viscous-electromagnetic (VE) and inertial-electromagnetic (IE) regimes for Reynolds numbers 50 < Re < 1500, Hartmann numbers 2500 < Ha < 5475, and expansion/contraction ratio 10. Numerical computations have shown that the MHD pressure drops in the eutectic lead-lithium (PbLi) flows in contractions are almost identical to those in the flows featuring expansion with slightly higher magnitudes in the contraction cases. The small discripancy in the MHD pressure drop between contractions and expansions (<8%) suggests that the earlier obtained correlations for the 3D MHD pressure drop in a duct flow with a sudden expansion can also be applied to flows with a sudden contraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An energy-based coupling approach to nonlocal interface problems.

Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, settings in which traditional partial differential equation models fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces, such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a physically consistent interface theory which is needed to support numerical developments and, among other features, reduces to classical models in the limit as the extent of nonlocal interactions vanish. In this paper, we use an energy-based approach to develop a formulation of a nonlocal interface problem which provides a physically consistent extension of the classical perfect interface formulation for partial differential equations. Numerical examples in one and two dimensions validate the proposed framework and demonstrate the scope of our theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonequilibrium Physics of Multiphase Flow in Porous Media: Wettability and Disorder

Carbon capture and geologic storage, dissociation of methane hydrates in permafrost, enhanced oil recovery, and water dropout in low-temperature fuel cells, all have something in common: two or more fluids flow simultaneously through a porous medium; and the displacement of one fluid by another is often unstable (either due to gravity or viscous forces). Yet, our ability to model multiphase flow in porous media has remained a challenge. The traditional equations are unable to predict, explain, or even reproduce, the formation of the complex patterns observed in experiments. The overarching goal of the project is to develop new physical understanding of the role of wettability and disorder in multiphase flow through permeable media, and to develop new mathematical and computational models at the pore scale and at the continuum scale.

58 GEOSCIENCES↗
DOE OSTI2024

Quantifying Trapped Powder in Electron Beam Powder Bed Fusion

Abstract Electron beam powder bed fusion (PBF-EB) shows great potential for manufacturing complex parts including those with internal cavities for heat exchanger, manifold systems, or energy absorption purposes. PBF-EB allows for the manufacture of channel geometries without the need for support structures. Due to the nature of the powder spreading process, powder feedstock is often trapped in intentionally manufactured cavities. This trapped powder can often be difficult to remove and can disturb the intended flow of fluid through the cavity or damage downstream components in its use case. These trapped powder particles present a risk of contamination and component failure if not completely evacuated. Ti6Al4V is a choice material for aerospace applications due to its high strength to weight ratio and its composition as a nonferrous metal; however, in weight sensitive applications excess entrapped powders or powders loosely attached to the surface could cause undesirable weight increases. The inherent spreading process of PBF-EB is different than laser powder bed fusion (PBF-LB) in its operational temperature, sintering. In addition, PBF-EB is less commonly studied in literature compared to its PBF-LB counterpart, and as a result the complexity of the semi-sintered powder and its spreading behavior are not well understood. Prior work has investigated the difficulty in removing trapped powder from PBF-EB, but these studies do not address how to quantify the amount of trapped powder in the cavity. Thus, an accurate method to measure the amount of trapped powder in the cavity must be investigated. In this work, Ti6Al4V coupons were manufactured with horizontal and vertical cavities of three different sizes. Archimedes testing allows for the determination of density differences caused by porosity and trapped powders by measuring mass and volumetric dispersion. Computed tomography (CT) is well suited for segmenting the internal structure and features of a part and has been studied for applications including voids, porosity, and dross. Thus, CT was explored as a method for evaluating trapped powder content in this work. The volumetric representation of the segmentation of the reconstructed CT volume can vary greatly depending on the input filter and thresholding methods. In this study, four different types of segmentation approaches were evaluated to determine the best approach for segmenting the volume as compared to an operator labeled ground truth. The percentage density results from the Archimedes testing were compared to the volumetric percent density from the computed tomography approach. Differences in packing density between two different internal channel features were investigated. Overall, this work sought to validate the use of computed tomography for the detection of trapped powders and present a framework for volumetric segmentation.

Johnstone, Brian↗
DOE OSTI2025

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗
DOE OSTI2020

Operation and Optimization of Microwave-Heated Continuous-Flow Microfluidics

Microwave (MW) technology can be powerful for electrification and process intensification but limited fundamental understanding of scalability and design principles hinders its effective use. In this work, we build a continuous-flow microreactor inside a commercial single-mode MW applicator and the corresponding computational fluid dynamics model to simulate the temperature profile. The model is in good agreement with experiments for various microreactor dimensions and operating conditions. The model indicates that MW heating is greatly influenced by reactor geometry as well as the operating parameters. We observe a strong correlation between parameters and develop a gradient boost regression tree model to predict the outlet temperature accurately. This model is then applied to optimize the dimensions and operating conditions to maximize the outlet temperature and energy efficiency, resulting in a Pareto optimal. We demonstrate computationally and experimentally that it is possible to surpass the Pareto optimal and achieve an energy efficiency of ~90% or greater at temperatures relevant for liquid-phase chemistry via salting of the solvent. The present methodology can be applied to other complex MW reactors. Lastly, the combined numerical and experimental approach provides insights into and a framework for scale-up and optimization.

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

Hydrothermal solution calorimetry in acidic aqueous solutions and revisiting the standard partial molal thermodynamic properties of Nd 3+ from 25 to 300 °C

The mobility of rare earth elements (REE) can be predicted in aqueous fluids using geochemical modeling but the accuracy of these models strongly depends on the availability of robust thermodynamic properties for the REE aqueous species. The REE 3+ aqua ions are important in the derivation of the formation constants of all the major REE complexes including the chloride, sulfate, and fluoride species which predominate in many hydrothermal-magmatic systems. However, the thermodynamic properties of the REE 3+ aqua ions are still commonly derived from the Helgeson-Kirkham-Flowers (HKF) equation of state parameters tabulated several decades ago. The standard state thermodynamic properties at reference conditions (25 °C and 1bar) and their extrapolations to high temperature need to be verified, if not revised, based on hydrothermal experiments. In this study, the enthalpy of solution was measured for synthetic Nd hydroxide from 25 to 150 ºC to retrieve the standard partial molal thermodynamic properties of Nd 3+ as a function of temperature. The experiments were conducted in aqueous perchloric acid based solutions with starting pH of 2 and varying ionic strength (0.01 to 0.09 mol/kg NaClO 4 ). The standard partial molal enthalpy of formation (Δ f H°) of Nd 3+ derived from the experimental study displays differences of up to 10 kJ/mol compared to the enthalpy values derived from the HKF equation of state in the studied temperature range. These inaccuracies are resolved by adjusting the standard partial molal Gibbs energy of formation (Δ f G°) of Nd 3+ at 25 °C and 1 bar from -672.0 to -679.7 kJ/mol. The heat capacity function (C p °) derived between 25 and 150 ºC can be described by: C p ° = a 0 + a 1 ·T + a 2 ·T -2 , with a 0 = 1256, a 1 = -2.68, a 2 = -55.56·10 6 and T in Kelvin. A set of recommended thermodynamic properties is provided for the Nd 3+ aqua ions and corrections are provided for the chloride and fluoride species to remain internally consistent with the experimentally derived properties. These results allow predicting accurately the solubility of monazite between 25 and 300 ºC. Before these corrections, the properties for the Nd 3+ aqua ions derived from the HKF parameters resulted in up to ~1.5 orders of magnitude lower monazite solubility than determined experimentally. Here, a revision of the REE +3 aqua ions properties is necessary to accurately predict the mobility of REE in hydrothermal acidic solutions.

58 GEOSCIENCES↗
DOE OSTI2025

Evidence for metal sources, fluid-mixing processes, and S isotope recycling within the feeder zone of an Irish type Zn-Pb deposit

The origin and evolution of fluids in Irish-type Zn-Pb deposits remains debated, particularly regarding the mobility of metals such as Cu and Ni, sources of sulphur, and the role of fluid mixing and replacement. The Lisheen Zn-Pb deposit, Ireland, offers a well-defined natural laboratory to investigate these questions. While most studies have focused on the Waulsortian Limestone Formation, the primary sulphide host, less is known about mineralisation in underlying units, such as the Lisduff Oolite Member (LOM). The LOM displays enrichment in Cu and Ni and displays intense replacement textures compared to other hosts at Lisheen, making it an ideal target for studying metal mobility and sulphur recycling in carbonate-hosted systems. Through characterising and studying LOM-hosted sulphides, valuable insights into mineralisation processes, especially related to Cu-Ni metals, can be defined. This study integrates petrography, EMPA, and in situ sulphur isotope (δ 34 S) analysis to investigate sulphide paragenesis, mineral chemistry, and fluid evolution across LOM ore zones. Results reveal a multistage mineralising system involving extensive replacement of early pyrite (Py0, δ 34 S = −28.4 to −21.9 ‰) by sphalerite and galena, with zoned pyrite (Py1) enriched in As-Cu-Ni-Tl. The δ 34 S values and trace element trends indicate mixing between hydrothermal and bacteriogenic sulphur-rich fluids, with evidence for sulphur recycling during replacement. Pyrite textures and compositions capture this evolving fluid regime, with trace element enrichment linked to paragenetic stage. The steel ore region, adjacent to major fault intersections, records intense hydrothermal fluid interaction, hosting Ni- and As-rich phases such as nickeline, gersdorffite, and arsenopyrite. These findings highlight the importance of structural controls and fluid mixing in metal transport and deposition, positioning the LOM as a key stratigraphic unit for understanding ore-forming processes in Irish-type systems. These results have implications for targeting similar carbonate-hosted systems globally, especially where deeper or structurally complex ore zones remain underexplored.

58 - GEOSCIENCES↗
DOE OSTI2021

How Equivalent Are Equivalent Porous Media?

Abstract Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

58 GEOSCIENCES↗
DOE OSTI2025

The influence of laser power modulation on melt pool dynamics in laser powder bed fusion

While the majority of laser powder bed fusion (LPBF) metal additive manufacturing uses a continuous wave (CW) laser heat source, some commercial applications of LPBF additive manufacturing instead involve the modulation of the laser power on tens-of-microsecond timescales as an adjustable process variable. This article reports the use of in situ, high speed x-ray and optical imaging to probe melt pool fluid flow, defect formation, and nearby powder motion during LPBF with both modulated and CW laser heat sources. We observe melt pool dynamics unique to modulated laser melting even at very high duty cycles that are related to fluctuations in vapor depression depth, complex pore formation mechanisms, and changes to denudation physics when compared to CW melting. These behaviors are present in Ti–6Al–4V, 316L stainless steel, and AL1100 alloys but vary slightly as a function of material, indicating a substantial dependence on the viscosity and surface tension of the liquid metal. While high duty cycles produce weld tracks of comparable quality to CW melting, lower duty cycles introduce substantial defect concentrations. At intermediate duty cycles, careful control of modulation parameters can repeatably and precisely yield one pore per laser pulse, suggesting a method for intentionally inserting engineered porosity at specific sites during an LPBF build.

3D printing↗