SEARCH · Engineering Papers

Results for “Multiphase”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Study of retained austenite and nano-scale precipitation and their effects on properties of a low alloyed multi-phase steel by the two-step intercritical treatment

Microstructure evolution and properties were studied in a low carbon low alloyed hot-rolled bainitic steel by annealing and annealing plus tempering. Microstructure of the hot-rolled steel consists of lath bainite and martensite. By annealing at 720 °C for 30 min and water quenching, multi-phase microstructure consisting of intercritical ferrite, tempered bainite/martensite, retained austenite and fresh martensite was obtained. With increasing annealing temperature to 760 °C, microstructure of the steel consisted of intercritical ferrite, fresh martensite without retained austenite. After the second step of tempering at 680 °C for samples annealed both at 720 °C and 760 °C, ~ 8–9% volume fraction of retained austenite was obtained in the multi-phase microstructure. Moreover, fine precipitates of VC with size smaller than 10 nm and copper precipitates with size of ~ 10–50 nm were obtained after tempering. Results from scanning transmission electron microscopy (STEM) give evidence to support that the partitioning of Mn, Ni and Cu is of significance for retained austenite stabilization. Due to the combined contribution of multiphase microstructure, the transformation-induced-plasticity effect of retained austenite and strengthening effect of nanometer-sized precipitates, yield strength greater than 800 MPa, yield to tensile ratio of 0.9, uniform elongation of ~ 9% and good low temperature impact toughness of 147 J at − 40 °C were achieved. - Highlights: • Stable retained austenite was produced in a low alloyed steel. • Partition of Mn, Ni and Cu was confirmed by STEM for austenite stabilization. • Nano-sized VC and Cu precipitates were achieved by second tempering. • High strength–high toughness with low Y/T ratio was obtained.

36 MATERIALS SCIENCE↗

DOE OSTI2017

The effect of Si on the microstructure and mechanical properties of spark plasma sintered boron carbide

Highlights: • Si addition leads to a dense final product at a low sintering temperature. • The added Si preferentially formed SiC and borosilicate glass instead of replacing the chain atoms in boron carbide. • The hardness of our Si-containing boron carbide is 30–35% greater than currently in-use commercial boron carbide plates. - Abstract: Fully dense boron carbide discs were achieved by spark plasma sintering boron carbide powders with 10 wt% silicon. The silicon did not diffuse into boron carbide grains to produce a solid solution of Si-doped boron carbide; instead the silicon reacted with impurities in the starting powder to form β-SiC and borosilicate glass. The resultant new phases facilitated densification of the multiphase ceramic through liquid phase-assisted sintering. The resultant material exhibits improved hardness (34.3 GPa Vikers hardness under 1 kg load) with toughness comparable to both Si-free and commercially available boron carbide.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Thermodynamic, structural and magnetic studies of phase transformations in MnAl nanocomposite alloys

Highlights: • Small changes in MnAl stoichiometry (5 at.%) produce different phase structure. • In Mn{sub 60}Al{sub 40} as-cast alloy only hcp ε phase can be found • In Mn{sub 55}Al{sub 45} as-cast alloy, there is a mixture of γ{sub 2} (Al{sub 8}Mn{sub 5}) and ε as well as ferromagnetic MnAl τ-phase. • Energies of phase transformation of about 165 kJ/mol for Mn{sub 60}Al{sub 40} and 290 kJ/mol for Mn{sub 55}Al{sub 45} alloys are found. • hcp ε phase decomposes through the migration of interphase interfaces. - Abstract: We have undertaken a temperature-dependent, investigation of the thermodynamics, structure, morphology and magnetism of two MnAl nanocomposite alloys (Mn{sub 60}Al{sub 40} and Mn{sub 55}Al{sub 45}). Differential scanning calorimetry (DSC) studies allowed the determination of the exo-effects occurring in the structural phase transformation of MnAl, while by using temperature-dependent X-ray diffraction of synchrotron radiation we were able to monitor the phase transformation effects and the evolution with temperature of various structural phases occurring in the samples. We have shown that slight changes in stoichiometry (5 at.%) give rise to different phase structure in the as-cast state. While in Mn{sub 60}Al{sub 40} only hcp ε phase can be found, in Mn{sub 55}Al{sub 45} as-cast alloy, there is a quite complex phase structure with a mixture of γ{sub 2} (Al{sub 8}Mn{sub 5}) and ε as well as ferromagnetic MnAl τ-phase. Activation energies of about 165 kJ/mol and 290 kJ/mol have been calculated from the Ozawa-Flynn-Wall analysis of DSC experimental data. These findings are in good agreement with the results obtained from XRD. For the Mn{sub 55}Al{sub 45} as-cast alloy, we have confirmed by temperature-dependent synchrotron XRD that the hcp ε phase decomposes through the migration of interphase interfaces with the transformation rate controlled by boundary diffusion processes. High-resolution transmission electron microscopy have confirmed the phase structure obtained by XRD, while magnetic properties, obtained for the as-cast alloys are consistent with the multiphase character of the samples and in good agreement with previously reported results. The magnetization does not saturate for the maximum applied field of about 4 × 10{sup 6} A/m and a marked coercivity of about 160 kA/m is obtained for the Mn{sub 55}Al{sub 45} as-cast alloy.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2018

High-resolution characterization of the martensite-austenite constituent in a carbide-free bainitic steel

Highlights: • Microstructural characterization of martensite-austenite (M-A) areas by TKD and TEM • Auto-tempering effects ranging from C clustering to carbide formation in martensite • Signs of structural modulation of austenite in M-A areas by electron diffraction • Heterogeneous C distribution within M-A areas due to auto-tempering during cooling - Abstract: The multiphase microstructure of carbide-free bainitic steels comprises bainitic ferrite laths, retained austenite with different morphologies, a minor fraction of carbides and so-called martensite-austenite areas, which partially transform during the last cooling step. While the other constituent received much attention, little is known about the structure of the martensite-austenite constituent in carbide-free bainitic steels. Thus, in this study, it was structurally and chemically investigated by high-resolution techniques such as transmission electron microscopy and atom probe tomography after preceded unambiguous identification by electron backscatter diffraction in conventional as well as transmission mode. The results, ranging from carbon segregation to cementite precipitation in the martensitic part, indicate strong auto-tempering during final cooling which is followed by aging. Also, some kind of structural modulation in the austenite belonging to the martensite-austenite areas was observed. Atom probe tomography revealed a heterogeneous carbon distribution, further supporting the findings by transmission electron microscopy.

36 MATERIALS SCIENCE↗

DOE OSTI2018

A construction of Ag-modified raspberry-like AgCl/Ag2WO4 with excellent visible-light photocatalytic property and stability

Highlights: • The Ag-modified raspberry-like AgCl/Ag{sub 2}WO{sub 4} photocatalysts were prepared by a facile solvent-free method. • Ag and AgCl were easily deposited on the Ag{sub 2}WO{sub 4} by one-step solid-sate reaction between Ag{sub 2}WO{sub 4} and HONH{sub 3}Cl. • The morphology of raspberry-like can facilitate the contact of dye molecules and composite photocatalysts. • The photocatalysts exhibit an effective separation of the photoinduced electron-hole pairs. • The photocatalysts show excellent photocatalytic activity and photo-stability under visible light irradiation. - Abstract: A novel multiphase photocatalyst of Ag-modified raspberry-like AgCl/Ag{sub 2}WO{sub 4} composite was successfully fabricated via a facile solvent-free reaction technique, in which the Ag and AgCl were subtlely generated from the one-step reaction between Ag{sub 2}WO{sub 4} particles and HONH{sub 3}Cl as reductant and chlorine source at room temperature. The as-obtained ternary Ag-modified raspberry-like AgCl/Ag{sub 2}WO{sub 4} composites can facilitate the contact of dye molecules and photocatalysts, exhibiting superior photocatalytic property and good recycling ability towards the decomposition of organic dye under the visible-light irradiation. The reductant addition amounts greatly affect the performance of Ag{sub 2}WO{sub 4} photocatalyst. The optimum composite ratio of Ag-modified AgCl/Ag{sub 2}WO{sub 4} samples show significant enhancement of visible-light photocatalytic activities and the degradation rate constant is about 16 times greater than that of pure Ag{sub 2}WO{sub 4}. The highly improved performance of Ag-modified raspberry-like AgCl/Ag{sub 2}WO{sub 4} structure can be ascribed to the surface plasma resonance (SPR) effects of metal Ag and the synergistic effect among composites, accordingly effectively promote the separation of photoinduced electron-hole pairs, widen the absorbance range in the visible light and increase surface area. The h{sup +} and ·OH degraded mechanism for the improved photocatalytic activity of Ag-modified AgCl/Ag{sub 2}WO{sub 4} is proposed based on radical-trapping experiments and ESR spectra records.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Facile synthesis of various α-Fe2O3 micro/nanostructures: Highlighting on the enhanced catalysis activities by formation of bowl-like α-Fe2O3/Au composites

Highlights: • Various α-Fe{sub 2}O{sub 3} micro/nanostructures are fabricated by a hydrothermal method. • The morphologies are tunable by changing the content of NaHCO{sub 3} and reaction time. • Bowl-like α-Fe{sub 2}O{sub 3} has the highest Fenton catalytic efficiency in dark condition. • Fast photo-Fenton reaction of α-Fe{sub 2}O{sub 3}/Au can be tuned by multiphase structures. • Bowl-like α-Fe{sub 2}O{sub 3}/Au is an excellent photo-Fenton catalyst for water treatment. - Abstract: Various α-Fe{sub 2}O{sub 3} micro/nanostructures referring to spherical-like, bowl-like, cubic, and discoid shapes have been fabricated by a facile hydrothermal method using sodium bicarbonate as surfactant and alkaline reagents, displaying superior Fenton catalytic performances for degradation of Rhodamine 6G. Because of abundant exposed crystal surfaces such as (110) and (104) on the surfaces of bowl-like α-Fe{sub 2}O{sub 3}, the rapid degradation rate of 90.0% in 5 min is observed under the dark condition as the adding amount of H{sub 2}O{sub 2} limiting to 300 μL. Significantly, the enhanced catalytic performances of bowl-like α-Fe{sub 2}O{sub 3}/Au catalysts under natural light condition have been investigated and carried out by in situ reduction of well-distributed Au nanoparticles with average sizes of 3–6 nm. Due to the appropriate phase structures and heterogenous interface properties, the highest photo-Fenton catalytic activity with the degradation rate of 90.0% in 2 min can be obtained for α-Fe{sub 2}O{sub 3}/Au{sub 0.01} composites.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Direct, efficient and selective capture of low concentration of CO 2 from natural gas flue gas using a high temperature tubular carbon capture membrane

Natural gas (NG) fired power plants emit low concentration (4–5%) of CO 2 , which presents additional technical and economic challenges to the current benchmark amine absorption technology. The newly emerged high-temperature multiphase membranes operated on molten carbonate (MC) chemistry for CO 2 capture/separation/conversion have been demonstrated with great potential to meet this challenge. In this study, we report on the CO 2 capture performance of such a membrane in tubular geometry from a mockup NG flue gas. The membrane is comprised of a mixture of Gd 0.20 Ce 0.80 O 1.95 (GDC) and MC, in which GDC forms a porous skeleton to contain MC. Here, we show that the membrane with a dimension of 6.1 mm in outer diameter, 5.1 mm in inner diameter and 5 cm in effective length (resulting in 4cm 2 effective surface area) can achieve a CO 2 flux density of 0.46–0.55 mL/min·cm 2 at 650°C, capturing 97% pure CO 2 at a rate of 37–42% from 5%CO 2 –N 2 using moistened Ar as the sweep gas. The level of performance demonstrated by this study suites the membrane well for stationary CO 2 capture from NG power plants.

03 NATURAL GAS↗

DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2014

Correlations between Nanoindentation Hardness and Macroscopic Mechanical Properties in DP980 Steels

Multiphase advanced high strength steels (AHSS) are being increasingly used in the automotive industry due to their low cost, good availability and excellent combination of strength and ductility. There is a keen interest from the automotive and steel industry for more fundamental understandings on the key microstructure features influencing the macroscopic properties, i.e., tensile properties, hole-expansion ratio and localized formability of AHSS. In this study, the micro- and macro-level properties for eight commercial DP980 steels are first characterized and quantified with various experimental methods. Correlations between macroscopic-level properties and relationships between various micro- and macro- properties for these steels are then established based on the experimental measurements. It is found that, despite their differences in their chemistry, processing parameters and sheet thickness, the eight DP980 steels do have common microstructural level properties governing their specific macroscopic properties in terms of strength, elongation and hole expansion performance.

Taylor, Mark D.↗

DOE OSTI2025

A novel microstructural approach for TRIP/TWIP β-Ti alloys using intermetallic precipitation

Interest in metastable β-Ti has grown, due to the exceptional work hardening caused by TRansformation and TWinning Induced Plasticity (TRIP/TWIP). A current challenge is increasing the yield strength, while maintaining the benefits of TRIP/TWIP. This study presents the first report of TRIP/TWIP in a metastable β-Ti alloy containing intermetallic compounds, specifically Ti silicide precipitates. Silicon (Si) additions enable the formation of silicide particles, which can be precipitated, dissolved, and controlled through tailored heat treatments. In-situ synchrotron x-ray diffraction during tensile testing reveals primary {332}<113> twinning activates first, followed by TRIP, leading to extensive work hardening and uniform elongation in the solutionized state. In conclusion, these findings open a novel microstructural approach to the design of multiphase TRIP/TWIP titanium alloys, relying on intermetallic precipitation rather than allotropes such as α or ω phase, and lay the groundwork for future studies on intermetallic precipitation strengthening.

Materials science↗

DOE OSTI2017

Development of an integrated four-channel fast avalanche-photodiode detector system with nanosecond time resolution

In this paper, a four-channel nanosecond time-resolved avalanche-photodiode (APD) detector system is developed at Beijing Synchrotron Radiation. It uses a single module for signal processing and readout. This integrated system provides better reliability and flexibility for custom improvement. The detector system consists of three parts: (i) four APD sensors, (ii) four fast preamplifiers and (iii) a time-digital-converter (TDC) readout electronics. The C30703FH silicon APD chips fabricated by Excelitas are used as the sensors of the detectors. It has an effective light-sensitive area of 10 × 10 mm 2 and an absorption layer thickness of 110 μ m. A fast preamplifier with a gain of 59 dB and bandwidth of 2 GHz is designed to readout of the weak signal from the C30703FH APD. The TDC is realized by a Spartan-6 field-programmable-gate-array (FPGA) with multiphase method in a resolution of 1ns. The arrival time of all scattering events between two start triggers can be recorded by the TDC. The detector has been used for nuclear resonant scattering study at both Advanced Photon Source and also at Beijing Synchrotron Radiation Facility. Finally, for the X-ray energy of 14.4 keV, the time resolution, the full width of half maximum (FWHM) of the detector (APD sensor + fast amplifier) is 0.86 ns, and the whole detector system (APD sensors + fast amplifiers + TDC readout electronics) achieves a time resolution of 1.4 ns.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2017

Interface tracking simulations of bubbly flows in PWR relevant geometries

[en] Highlights: • Simulations were performed for turbulent bubbly flows in PWR subchannel geometry. • Liquid turbulence is fully resolved by direct numerical simulation approach. • Bubble behavior is captured using level-set interface tracking method. • Time-averaged single- and two-phase turbulent flow statistical quantities are obtained. - Abstract: The advances in high performance computing (HPC) have allowed direct numerical simulation (DNS) approach coupled with interface tracking methods (ITM) to perform high fidelity simulations of turbulent bubbly flows in various complex geometries. In this work, we have chosen the geometry of the pressurized water reactor (PWR) core subchannel to perform a set of interface tracking simulations (ITS) with fully resolved liquid turbulence. The presented research utilizes a massively parallel finite-element based code, PHASTA, for the subchannel geometry simulations of bubbly flow turbulence. The main objective for this research is to demonstrate the ITS capabilities in gaining new insight into bubble/turbulence interactions and assisting the development of improved closure laws for multiphase computational fluid dynamics (M-CFD). Both single- and two-phase turbulent flows were studied within a single PWR subchannel. The analysis of numerical results includes the mean gas and liquid velocity profiles, void fraction distribution and turbulent kinetic energy profiles. Two sets of flow rates and bubble sizes were used in the simulations. The chosen flow rates corresponded to the Reynolds numbers of 29,079 and 80,775 based on channel hydraulic diameter (D_h) and mean velocity. The finite element unstructured grids utilized for these simulations include 53.8 million and 1.11 billion elements, respectively. This has allowed to fully resolve all the turbulence scales and the deformable interfaces of individual bubbles. For the two-phase flow simulations, a 1% bubble volume fraction was used which resulted in 17 bubbles in the smaller case and 262 bubbles in the larger case. In the larger simulation case the size of the resolved bubbles is 0.65 mm in diameter, and the bulk mesh cell size is about 30 microns. Those large-scale simulations provide new level of details previously unavailable and were enabled by the excellent scaling performance of our two-phase flow solver and access to the state-of-the-art supercomputing resources. The presented simulations used up to 256 thousand processing threads on the IBM BG/Q supercomputer “Mira” (Argonne National Laboratory).

Nuclear Science & Technology↗

DOE OSTI2017

Wall-resolved spectral cascade-transport turbulence model

A spectral cascade-transport model has been developed and applied to turbulent channel flows (Reτ= 550, 950, and 2000 based on friction velocity, uτ ; or ReδΜ= 8,500; 14,800 and 31,000, based on the mean velocity and channel half-width). This model is an extension of a spectral model previously developed for homogeneous single and two-phase decay of isotropic turbulence and uniform shear flows; and a spectral turbulence model for wall-bounded flows without resolving the boundary layer. Data from direct numerical simulation (DNS) of turbulent channel flow was used to help develop this model and to assess its performance in the 1D direction across the channel width. The resultant spectral model is capable of predicting the mean velocity, turbulent kinetic energy and energy spectrum distributions for single-phase wall-bounded flows all the way to the wall, where the model source terms have been developed to account for the wall influence. We implemented the model into the 3D multiphase CFD code NPHASE-CMFD and the latest results are within reasonable error of the 1D predictions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Direct numerical simulation of reactor two-phase flows enabled by high-performance computing

Nuclear reactor two-phase flows remain a great engineering challenge, where the high-resolution two-phase flow database which can inform practical model development is still sparse due to the extreme reactor operation conditions and measurement difficulties. Owing to the rapid growth of computing power, the direct numerical simulation (DNS) is enjoying a renewed interest in investigating the related flow problems. A combination between DNS and an interface tracking method can provide a unique opportunity to study two-phase flows based on first principles calculations. More importantly, state-of-the-art high-performance computing (HPC) facilities are helping unlock this great potential. This paper reviews the recent research progress of two-phase flow DNS related to reactor applications. The progress in large-scale bubbly flow DNS has been focused not only on the sheer size of those simulations in terms of resolved Reynolds number, but also on the associated advanced modeling and analysis techniques. Specifically, the current areas of active research include modeling of sub-cooled boiling, bubble coalescence, as well as the advanced post-processing toolkit for bubbly flow simulations in reactor geometries. A novel bubble tracking method has been developed to track the evolution of bubbles in two-phase bubbly flow. Also, spectral analysis of DNS database in different geometries has been performed to investigate the modulation of the energy spectrum slope due to bubble-induced turbulence. In addition, the single-and two-phase analysis results are presented for turbulent flows within the pressurized water reactor (PWR) core geometries. The related simulations are possible to carry out only with the world leading HPC platforms. Furthermore, these simulations are allowing more complex turbulence model development and validation for use in 3D multiphase computational fluid dynamics (M-CFD) codes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

High-fidelity pool boiling simulations on multiple nucleation sites using interface capturing method

Boiling has proved to be one of the most efficient means for heat transfer and is a very important phenomenon during severe accident scenarios in light water reactors. High-fidelity pool boiling simulations can provide a numerical database for improving mechanistic boiling models by allowing for specific evaluation of interactions among bubbles. Previously published pool boiling simulations investigated two nucleation sites in which bubble growth at one site suppressed nucleation at the other site. Based on previous study results, more complicated interface-capturing simulations on pool boiling were conducted using PHASTA code with locally refined unstructured mesh. First, different boundary conditions (BCs) were assessed to support robustness and reproducibility of the boiling model. Then, a scale study was conducted at a larger domain with nine nucleation sites where either nine or four nucleation sites are activated. Involving more nucleation sites increased the complexity of bubble interactions from surrounding sites. Finally, bubble departure behavior influenced by wall heat flux was investigated. In this work, when heat flux was increased, the order of bubble departure changed, but diagonal bubbles always departed after one another. The departure time interval between the first and second bubble reduced as heat flux increased. The corresponding frequency was almost linearly proportional to the heat flux. In addition, bubble departure behavior was found to be greatly influenced by the nucleation site pattern. Multiple nucleation sites resulted in superimposed inhibitive effects from surrounding sites to each bubble, which extensively delayed the departure. This new observation was not discussed in previously published works. The work presented here provides new insight on the fundamental understanding of boiling phenomena, contributes to the development of a 3D multiphase computational fluid dynamics (M-CFD) model, and provides a more comprehensive database for data-driven pool boiling studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗

DOE OSTI2018

Integrated simulation of vertical fracture propagation induced by water injection and its borehole electromagnetic responses in shale gas systems

We analyze fracture propagation induced by hydraulic fracturing with water injection and examine their detectability with crosswell electromagnetic (EM) geophysical methods. For rigorous 3D coupled flow-geomechanical modeling, we employ a numerical method that can model failure by tensile and shear stresses, dynamic nonlinear permeability, dual continuum approach, and thermo-poro-mechanical effects. From numerical simulation, we find that the fracture propagation is not the same as propagation of the water front, because fracturing is governed by geomechanics whereas water saturation is determined by multiphase flow. At early times, the water front is almost identical to the fracture tip, suggesting that the fracture is mostly filled with the injected water. However, at late times, movement of the water front is retarded compared to fracture propagation, yielding a significant gap between the water front and the fracture top, filled with reservoir gas. During fracture propagation, the coupled flow-geomechanical models are transformed via a rock-physics model into electrical conductivity models. We employ a full 3D finite-element EM geophysical simulator to evaluate the sensitivity of the crosswell EM method to fracture propagation. It is shown that anomalous distribution of electrical conductivity is closely related to the injected water saturation, but not closely related to newly created unsaturated fractures. Our numerical modeling experiments demonstrate that the crosswell EM method can be highly sensitive to electrical conductivity changes that directly indicate the migration pathways of the injected fluid. Accordingly, the EM method can serve as an effective monitoring tool for monitoring the injected fluids during hydraulic fracturing operations.

58 GEOSCIENCES↗