Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase velocity”

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

Rock Valley Accelerated Weight Drop Seismic Data Processing and Picking of P-wave and S-wave Arrival Times

Rock Valley, in the southern end of the Nevada National Security Site, hosts a fault system that was responsible for a shallow (< 3 km below surface ) magnitude 3.7 earthquake in May 1993. In order to better understand this system, seismic properties of the shallow subsurface need to be better constrained. In April and May of 2021, accelerated weight drop (AWD) active-source seismic data were recorded in order to measure P- and S-wave travel-times for the area. This report describes the processing and phase picking of the recorded seismic waveforms. In total, we picked 7,982 P-wave arrivals at offsets up to ~2500 m, and 4,369 S-wave arrivals at offsets up to ~2200 m. These travel-time picks can be inverted for shallow P-wave and S-wave velocity structure in future studies.

58 GEOSCIENCES↗
DOE OSTI2020

Microstructure Experiments-Enabled MARMOT Simulations of SiC/SiC-based Accident Tolerant Nuclear Fuel System

We have undertaken an experimental-computational project that addresses a few key technology gaps associated with the use of SiC/SiC composites for light water reactor fuel cladding. Two principal endeavors in this project are to assess the irradiation-induced microstructural changes and swelling in SiC/SiC composites, and to characterize and model the porous oxide surface layer (when exposed to steam) along SiC recession that is detrimental to the clad integrity under accident conditions. The project tasks are: (i) ion irradiation and characterization (University of Tennessee, Knoxville), (ii) mechanical tests and analysis (University of South Carolina), (iii) steam exposure tests (Oak Ridge National Laboratory), (iv) electron microscopy and spectroscopic characterization (NC State University), (v) x-ray microscopy and reconstruction (University of South Carolina and NC State University), (vi) phase field modeling and simulations (NC State University and Idaho National Laboratory). The effects of 10 MeV Au ion irradiation at 350°C on the microstructure evolution in SiC/SiC composites are investigated at doses up to 400 displacements per atom (dpa) at the University of Tennessee, Knoxville. Atomic force microscopy and optical profilometry reveal irradiation induced axial and radial shrinkage of the fibers for doses greater than 10 dpa. Based on detailed electron microscopy characterization, the primary cause of the fiber shrinkage is attributed to irradiation-induced loss of carbon packets. Additionally, the multilayer PyC interface is observed to portray high resistance to irradiation damage. The mechanical response and failure mechanisms of un-irradiated samples is also assessed through loading tests and X-ray imaging at the University of South Carolina. Steam exposure tests are performed at the Oak Ridge National Laboratory in a facility that represents a reactor pressure vessel under a loss-of-coolant-accident scenario. The samples that are analyzed methodically in this report are tested for 32/31 hours at 1200°C with a velocity of 0.25 cm/s for the pressures: 0.1 MPa, 0.45 MPa, 0.92 MPa and 1.38 MPa. Scanning/transmission electron microscopy analysis conducted at the NC State University (NCSU) shows that the oxide layer thickness increases with the steam pressure. While the oxide layer is crystalline (α- cristobalite) for the pressures 0.45 MPa, 0.92 MPa and 1.38 MPa, the layer is amorphous at 0.1 MPa. Results from Raman spectroscopy have also confirmed the formation of α-cristobalite phase of SiO₂. The abrupt increase in the integrated Raman intensity ratio between steam pressures 0.1 and 0.45 MPa suggests the onset of accelerated crystallization. Non-destructive three dimensional X-ray microscopy/tomography (XCT) and computational image processing techniques are employed by the University of South Carolina to probe the porous oxidation features on SiC samples at varying pressures. Interestingly, most pores are observed to be located away from the surface as well as the oxide-SiC interface. The average oxide layer thickness assessed from the XCT analysis is seen to be in excellent agreement with the values determined through scanning electron microscopy analysis for the highest pressures where the oxide layer is relatively more uniform. A phase-field model developed by the NC State University and Idaho National Laboratory for simulating oxidation of SiC by steam captures the paralinear kinetics of SiC oxidation with a high degree of fidelity. Results from quasi-one-dimensional and two dimensional simulations show that the pores with oxidizing species lead to a higher volatilization rate. These simulations indicate that the enhanced apparent volatilization rates at higher pressures observed in experiments can be rationalized by the increased volatilization from the pores. The project team has also successfully developed the capability to import images from experiments for realistic evolution of the oxidizing microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Observational constraints on the multiphase nature of outflows using large spectroscopic surveys at z ∼ 0

ABSTRACT Mass outflow rates and loading factors are typically used to infer the quenching potential of galactic-scale outflows. However, these generally rely on observations of a single gas phase that can severely underestimate the total ejected gas mass. To address this, we use observations of high mass (≥1010 M⊙), normal star-forming galaxies at z ∼ 0 from the MaNGA, xCOLD GASS, xGASS, and ALFALFA surveys and a stacking of Na d, Hα, CO(1–0), and H i 21 cm tracers with the aim of placing constraints on an average, total mass outflow rate, and loading factor. We find detections of outflows in both neutral and ionized gas tracers, with no detections in stacks of molecular or atomic gas emission. Modelling of the outflow components reveals velocities of |vNa d| = 131 km s−1 and |vHα| = 439 km s−1 and outflow rates of $\dot{M}_{\rm {Na\,\small{D}}}$ = 7.55 M⊙ yr−1 and $\dot{M}_{\text{H}\alpha }$ = 0.10 M⊙ yr−1 for neutral and ionized gas, respectively. Assuming a molecular/atomic outflow velocity of 200 km s−1, we derive upper limits of $\dot{M}_{\text{CO}}\lt 19.43$ M⊙ yr−1 and $\dot{M}_{\rm {H\,\small {I}}}\lt $ 26.72 M⊙ yr−1 for the molecular and atomic gas, respectively. Combining the detections and upper limits, we find average total outflow rates of $\dot{M}_{\text{tot}}\lesssim$27 M⊙ yr−1 and a loading factor of ηtot ≲ 6.39, with molecular gas likely contributing ≲72 per cent of the total mass outflow rate, and neutral and ionized gas contributing ∼28 and <1 per cent, respectively. Our results suggest that, to first order, a degree of quenching via ejective feedback could occur in normal galaxies when considering all gas phases, even in the absence of an active galactic nucleus.

Roberts-Borsani, G. W.↗
DOE OSTI2026

Ultrafast x-ray scattering of photodissociation dynamics in 2-iodothiophene

Here, a time-resolved x-ray scattering (TRXS) investigation of the photodissociation dynamics of gas-phase 2-iodothiophene molecules following 252 nm excitation is presented. Structural evolution of the molecule and dynamical information on the resulting photofragments were captured using femtosecond x-ray free-electron laser pulses. Two dissociation pathways were identified, arising via excitation to ππ* and (n/π)σ* states, respectively, yielding distinct interfragment recoil velocities of ∼6.4 Åps −1 and 17.0 Åps −1 . A comparison of asymptotic scattering data with simulated patterns indicates that the thiophene ring remains closed following dissociation at this wavelength. Modeling the experimental data yields a branching ratio of ∼3:1 in favor of the high velocity channel. These findings demonstrate the capability of TRXS to disentangle concurrent ultrafast reaction pathways and provide detailed structural insight into energy redistribution during photoinduced bond fission in complex molecular systems.

Razmus, Weronika O. [Univ. of Southampton (United ↗
DOE OSTI2020

Burning of particles undergoing stripping breakup

A simple model for condensed-phase particle burning has been presented. Particles are assumed to start burning in quiescent air. Particles are then accelerated, adding convective effects. Convective effect due to the relative motion between the particle and air on heat and mass transfer is included by introducing the Nusselt number. Stripping of liquid is included at higher particle velocity. Evaporation and burning from the stripped tail are also included.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Individual Wave Detection and Tracking within a Rotating Detonation Engine through Computer Vision Object Detection applied to High-Speed Images

Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.

Johnson, Kristyn↗
DOE OSTI2022

Evidence for X-Ray Emission in Excess to the Jet-afterglow Decay 3.5 yr after the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first ~3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm -3 . We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 10 38 erg s -1 at 1234 days, the recently detected X-ray emission represents a ≥ 3.2σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥ 3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥ 3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Individual Wave Detection and Tracking within a Rotating Detonation Engine through Computer Vision Object Detection applied to High-Speed Images

Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.

Johnson, Kristyn↗
DOE OSTI2021

Manipulation of Magnetic Skyrmion in a 2D van der Waals Heterostructure via Both Electric and Magnetic Fields

As a promising candidate for the much-desired low power consumption spintronic devices, 2D magnetic van der Waals material also provides a versatile platform for the design and control of topological spin textures. In this work on WTe 2 /CrCl 3 bilayer van der Waals heterostructures, a complete Néel-type skyrmion–bimeron–ferromagnet phase transition is demonstrated, accompanied by the evolution of the topological number. This cyclic transition, mediated by a perpendicular magnetic field, is largely driven by the competition between the out-of-plane magnetocrystalline anisotropy and magnetic dipole–dipole interaction. In the presence of a driving current, the Néel-type skyrmion gains a higher velocity yet larger skyrmion Hall angle, in comparison to the bimeron. By incorporating a ferroelectric CuInP 2 S 6 monolayer as a substrate, writing and erasing of skyrmions may be regulated using a ferroelectric polarization. Here this work sheds light on a novel approach to the design and control of magnetic skyrmions on 2D van der Waals materials.

multiferroic↗
DOE OSTI2025

Is There a Scalar Atmospheric Surface Layer Within a Convective Boundary Layer? Implications for Flux Measurements

Top‐down entrainment shapes the vertical gradients of sensible heat, latent heat, and CO 2 fluxes, influencing the interpretation of eddy covariance (EC) measurements in the unstable atmospheric surface layer (ASL). Using large eddy simulations for convective boundary layer flows, we demonstrate that decreased temperature gradients across the entrainment zone increase entrainment fluxes by enhancing the entrainment velocity, amplifying the asymmetry between top‐down and bottom‐up flux contributions. These changes alter scalar flux profiles, causing flux divergence or convergence and leading to the breakdown of the constant flux layer assumption (CFLA) in the ASL. As a result, EC‐measured fluxes either underestimate or overestimate “true” surface fluxes during divergence or convergence phases, contributing to energy balance non‐closure. The varying degrees of the CFLA breakdown are a fundamental cause for the non‐closure issue. These findings highlight the underappreciated role of entrainment in interpreting EC fluxes, addressing non‐closure, and understanding site‐to‐site variability in flux measurements.

eddy covariance fluxes of scalars↗
DOE OSTI2025

First Order Transport of a Cold Uniform Ellipsoid of Charge

A single particle has a 6-element phase space vector $s = (x, γv)$ where $γv =\frac{p}{m}$ s proportional to the relativistic momentum. The ellipsoid is represented by a 21-element vector $(s, C, V)$ where $s$ is the phase space location of the centroid. $C$ is the 3 × 3 covariance matrix defined by $C_{ij} = ⟨δx_iδx_j⟩$, where $δx$ is the distance of a particle from the centroid $x$. $C$ is symmetric so only counts as 6 elements of the overall vector. $V$ is a general 3 × 3 matrix that determines the velocity distribution in the bunch by $δ(γv) = Vδx$. In total this makes 6 + 6 + 9 = 21 elements.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

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

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Atomic oxygen density determination in the effluent of the COST reference source using in situ effective lifetime measurements in the presence of a liquid interface

Abstract Spatially resolved, absolute densities of atomic oxygen are measured for several helium-based admixtures in the effluent of the COST Reference Microplasma Jet using two-photon absorption laser induced fluorescence (TALIF). Admixtures investigated include a helium-only admixture, four helium/oxygen admixtures (0.1 % , 0.5 % , 0.6 % , and 1.0 % oxygen), and a helium/water admixture (2500 ppm water), chosen to coincide with previously published characterizations of plasma-treated liquid. The atomic oxygen TALIF signal is calibrated for density using the noble gas xenon, which possesses a very similar two-photon excitation and fluorescence scheme. Measurements are conducted for the jet operating in ambient air with both an open effluent and a liquid surface present, allowing for comparison with liquid phase measurements conducted under similar conditions. The presence of a water surface does not appear to alter the background chemistry in the effluent but reduces O densities close to the liquid interface when compared to a similar distance from the nozzle in an open effluent case. This may be the result of a reduction in flow velocity caused by the liquid obstructing the gas flow. Additionally, measurements near the liquid surface revealed a region of atomic oxygen well outside of where the core of the effluent impinges on the liquid. This is likely relevant for applications as it considerably expands the surface area subject to O absorption. Critically, in situ measurements of the effective lifetimes of the laser-excited 3p 3 P J state of atomic oxygen were recorded in the effluent by employing a picosecond (ps) laser and a nanosecond (ns) ICCD. By experimentally determining the contribution from collisional quenching via the in situ effective lifetime measurements, significant improvements in the accuracy of the atomic oxygen density calibration were made, with differences of approximately 30 % from existing methods of estimating quenching rates at atmospheric pressure. Finally, spatially resolved atomic oxygen densities allow for an investigation of O formation and extinction pathways in the effluent and a comparison between admixtures.

laser diagnostics↗
DOE OSTI2021

Self-organization of nanoparticles and molecules in periodic Liesegang-type structures

Chemical organization in reaction-diffusion systems offers a strategy for the generation of materials with ordered morphologies and structural hierarchy. Periodic structures are formed by either molecules or nanoparticles. On the premise of new directing factors and materials, an emerging frontier is the design of systems in which the precipitation partners are nanoparticles and molecules. We show that solvent evaporation from a suspension of cellulose nanocrystals (CNCs) and L -(+)-tartaric acid [ L -(+)-TA] causes phase separation and precipitation, which, being coupled with a reaction/diffusion, results in rhythmic alternation of CNC-rich and L -(+)-TA–rich rings. The CNC-rich regions have a cholesteric structure, while the L -(+)-TA–rich bands are formed by radially aligned elongated bundles. The moving edge of the pattern propagates with a finite constant velocity, which enables control of periodicity by varying film preparation conditions. This work expands knowledge about self-organizing reaction-diffusion systems and offers a strategy for the design of self-organizing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Motion of Air Bubbles in a Cement Slurry

The dynamics of air (gas) bubbles in a column of cement slurry is examined numerically. The air injected at the bottom of a laboratory-scale column through a porous distributor plate spatially distributes and migrates as a swarm of bubbles throughout the slurry toward the freeboard. The two-phase system of the cement slurry and the air bubbles is modeled using the conservation equations of mass and linear momentum in the framework of the volume-of-fluid (VOF) approach. The cement slurry is modeled using the Herschel–Bulkley and Bingham fluid models. Results show that the mean Sauter diameter and the mean rise velocity of the bubbles decrease with the gas flow rate. Meanwhile, it is found that the rising of the bubbles is controlled by breakup events, along with relatively weak path instabilities of the bubbles resulting in relatively straight trajectories, independent of the gas flow rate. The extent of the yielded region appears larger for the Herschel–Bulkley model compared to the Bingham fluid model (by approximately 10%).

36 MATERIALS SCIENCE↗