Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bifurcations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Plasticity-mediated deformation instabilities in thin film-compliant substrate systems: direct three-dimensional simulations

Abstract Surface wrinkles driven by mechanical instability commonly form in thin-film structures attached to a compliant substrate. In this study, a recently developed computational approach is employed to simulate the formation and transformation of wrinkles involving plastic yielding of the thin film. The three-dimensional (3D) finite element models contain an embedded imperfection at the film-substrate interface, serving to trigger the bifurcation modes. Successful application of this technique to allow for film plasticity is demonstrated, including the evolution of 3D surface patterns and their correlation with the overall load–displacement response. The simulations reveal that plastic yielding transforms the surface instability patterns into more localized forms. Under uniaxial loading, the sinusoidal elastic wrinkles undergo the wrinkle-to-fold transition. With equi-biaxial loading, the initial square-checkerboard array turns into continuous tall ridges along the 45° directions. In both loading modes, the plasticity-induced instability patterns are only partially relieved upon unloading, leaving permanent features on the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Consequences and benefits of utilizing continuous vibro-impact representations in constrained pipeline conveying fluid systems

Here, the effectiveness of continuous vibro-impact forcing representations for the cantilevered pipe that conveys fluid is explored and analyzed. The previously accepted forcing model utilizing a smoothened trilinear spring is estimated using three continuous forcing representations, namely, polynomial, rational polynomial, and hyperbolic tangent. The accuracy of the estimated forcing functions is investigated and analyzed by calculating the root mean square error, and bifurcation diagrams are generated and compared to the nominal system. Additionally, the dynamic response of the system is further characterized using Poincare maps, power spectra, and basins of attraction. Once all continuous forcing representations are analyzed and compared to the nominal system, the computational cost of each method is examined, and further limitations of the hyperbolic tangent method are discovered. It is proved that the hyperbolic tangent forcing representation most accurately captures the dynamic response of the pipeline, and the least accurate representation is the rational polynomial representation. Additionally, considerable computational cost is saved when employing the hyperbolic tangent representation compared to the discontinuous representation.

42 ENGINEERING↗
DOE OSTI2021

Transport Barriers in Symplectic Maps

Chaotic transport is a subject of paramount importance in a variety of problems in plasma physics, specially those related to anomalous transport and turbulence. On the other hand, a great deal of information on chaotic transport can be obtained from simple dynamical systems like two-dimensional area-preserving (symplectic) maps, where powerful mathematical results like KAM theory are available. In this work we review recent works on transport barriers in area-preserving maps, focusing on systems which do not obey the so-called twist property. For such systems KAM theory no longer holds everywhere and novel dynamical features show up as nonresistive reconnection, shearless curves and shearless bifurcations. After presenting some general features using a standard nontwist mapping, we consider magnetic field line maps for magnetically confined plasmas in tokamaks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Schmid factor crack propagation and tracking crystallographic texture markers of microstructural condition in direct energy deposition additive manufacturing of Ti-6Al-4V

Metallic additive manufacturing (AM) provides a customizable and tailorable manufacturing process for new engineering designs and technologies. The greatest challenge currently facing metallic AM is maintaining control of microstructural evolution during solidification and any solid state phase transformations during the build process. Ti-6Al-4V has been extensively surveyed in this regard, with the potential solid state and solidification microstructures explored at length. In this study, we evaluate the applicability of previously determined crystallo- graphic markers of microstructural condition observed in electron beam melting powder bed fusion (PBF-EB) builds of Ti-6Al-4V in a directed energy deposition (DED) build process. The aim of this effort is to elucidate whether or not these specific crystallographic textures are useful tools for indicating microstructural conditions in AM variants beyond PBF-EB. Parent β-Ti grain size was determined to be directly related to α-Ti textures in the DED build process, and the solid state microstructural condition could be inferred from the intensity of specific α-Ti orientations. Qualitative trends on the as-solidified β-Ti grain size were also determined to be related to the presence of a fiber texture, and proposed as a marker for as-solidified grain size in any cubic metal melted by AM. Analysis of the DED Ti-6Al-4V build also demonstrated a near complete fracture of the build volume, suspected to originate from accumulated thermal stresses in the solid state. Crack propagation was found to only appreciably occur in regions of slow cooling with large α +β colonies. Schmid factors for the basal and prismatic α-Ti systems explained the observed crack pathway, including slower bifurcation in colonies with lower Schmid factors of both slip systems. Colony morphologies and localized equiaxed β-Ti solidification were also found to originate from build pauses during production and uneven heating of the build edges during deposition. Tailoring of DED Ti-6Al-4V microstructures with the insight gained here is proposed, along with cautionary insight on preventing unplanned build pauses to maintain an informed and controlled thermal environment for microstructural control.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An Unmanned Aerial System (UAS) for concurrent measurements of solar-induced chlorophyll fluorescence and hyperspectral reflectance toward improving crop monitoring

Unmanned aerial system (UAS)-based remote sensing can provide high spatial- and temporal-resolution crop monitoring for precision management. Existing crop monitoring UAS primarily use conventional techniques such as multi-spectral broadband vegetation indices (VIs) to track canopy structure/biomass. Recently developed lightweight hyperspectral sensors can track crop physiology and performance via complementary signals such as solar-induced chlorophyll fluorescence (SIF) and photochemical reflectance index (PRI) that can be derived from hyperspectral reflectance. However, few existing UAS can acquire high-quality SIF and hyperspectral reflectance. Here, we designed a novel UAS that simultaneously captures far-red SIF and hyperspectral reflectance, using a single bifurcated fiber and motorized arm to measure downwelling and upwelling irradiance. The UAS was tested over (1) a heterogeneous corn field with spatially varying crop yield potential and (2) a set of soybean and corn plots under differing nutrient treatments. Our UAS can maintain stability under low-to-moderate winds (2–3 m/s) with <0.4 m geolocation accuracy, <0.1 m altitude drift, and <1.5° error of the fiber optic from nadir position during scanning. Weekly seasonal campaigns reveal that SIF outperforms conventional VIs such as the normalized difference vegetation index (NDVI) for distinguishing plots with different crop yield potential, especially after canopy closure when NDVI tends to saturate. The UAS can capture the diurnal dynamics of SIF and PRI for both corn and soybean. At the seasonal scale, SIF acquired from the UAS correlates well with that from a fixed SIF tower system ( R 2 = 0.81), which measured the same target albeit with different footprints, demonstrating the capability of UAS in characterizing the seasonal progression of crop activity. In conclusion, our newly developed UAS provides high-quality SIF and hyperspectral reflectance, facilitating mechanistic understanding of the physiological control on photosynthesis dynamics. Our findings also imply that this UAS could enable extension of fixed tower-based systems to monitor multiple targets from diurnal to seasonal scales and improve monitoring of heterogeneous fields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Machine learning informed rational design of high entropy double perovskite oxide universal air/steam electrodes for solid oxide electrochemical cells

Due to their high efficiency and versatility, solid oxide electrochemical cells (SOCs) are poised to play a significant role in future energy conversion and storage applications. In recent years, SOCs have bifurcated into two distinct categories: traditional oxygen-ion conducting SOCs that typically operate from ∼650—850 °C and the more recent proton-conducting ceramic (PCC) SOCs that typically operate from ∼400—650 °C. Current performance and lifetime of both oxygen-ion conducting SOCs and PCCs is primarily limited by the air/steam electrode, which facilitates the oxygen reduction reaction (ORR) during fuel cell operation and must also facilitate the oxygen evolution reaction (OER) during electrolysis operation. Here, we present a newly designed high-entropy double perovskite oxide suitable as a universal ORR/OER electrode for both oxygen-ion conducting SOCs and PCCs. Machine learning methods are applied to identify chemical descriptors for highly catalytic high-entropy double perovskite oxides (AA’B 2 O 6 ) across a large compositional space. Based on the machine-learning guidance, we ultimately converge on Ba 0.9 Cs 0.1 (Ca 0.2 Gd 0.2 La 0.2 Pr 0.2 Sr 0.2 )Co 1.5 Fe 0.5 O 6 (CsBaHEO) as a universal air/steam electrode. Structure stabilization is accomplished by an equimolar five-cation high-entropy composition on the A’-site, while cesium substitution on the A-site enhances the electrical conductivity and leads to a higher oxygen vacancy concentration. This material exhibits versatility and high performance in reversible oxygen-ion SOCs, reversible PCCs, and also large-scale tubular PCCs. For example, the CsBaHEO-based PCC reaches 1018 mW∙cm −2 at 600°C, while a large-scale tubular PCC using CsBaHEO for electrolysis achieves a hydrogen production rate of 21.314 ML∙min −1 at 600 °C.

Cell↗
DOE OSTI2023

Vibro-impact analysis and characterization of pipeline conveying fluids with multi-segmented motion-limiting constraints

Previous studies of the cantilevered pipeline conveying fluid system have included motion-limiting constraints in the form of trilinear springs. While this is desirable in experimental scenarios, it may not be representative of real-world applications. Therefore, here, this study focuses on multi-segmented motion-limiting constraints. As this type of motion-limiting constraint has not been investigated with a cantilevered pipeline system, a wide variety of outer and inner constraint stiffness and constraint gap sizes are investigated in this study to gain a comprehensive understanding of how the multi-segmented constraints affect the dynamics of the cantilevered pipeline. In this effort, bifurcation diagrams, phase portraits, Poincare maps, time histories, and power spectra are used to investigate the dynamics of the system, and the fluid flow speeds where dynamic characteristics are considered. In general, it is found that critical flow speeds like when the pipe sticks in the constraints are reduced as the constraint stiffnesses are increased. Additionally, the sticking flow speed occurred at lower flow speeds as the gap sizes of the inner and outer constraints decrease, and a larger constraint offset results in a smaller inner gap size leading to critical behaviors occurring at earlier flow speeds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

VRC34-Antibody Lineage Development Reveals How a Required Rare Mutation Shapes the Maturation of a Broad HIV-Neutralizing Lineage

Rare mutations have been proposed to restrict the development of broadly neutralizing antibodies against HIV-1, but this has not been explicitly demonstrated. We hypothesized that such rare mutations might be identified by comparing broadly neutralizing and non-broadly neutralizing branches of an antibody-developmental tree. Because sequences of antibodies isolated from the fusion peptide (FP)-targeting VRC34-antibody lineage suggested it might be suitable for such rare mutation analysis, we carried out next-generation sequencing (NGS) on B cell transcripts from donor N123, the source of the VRC34 lineage, and functionally and structurally characterized inferred intermediates along broadly neutralizing and poorly neutralizing developmental branches. The broadly neutralizing VRC34.01 branch required the rare heavy-chain mutation Y33P to bind FP, whereas the early bifurcated VRC34.05 branch did not require this rare mutation and evolved less breadth. Our results demonstrate how a required rare mutation can restrict development and shape the maturation of a broad HIV-1-neutralizing antibody lineage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Reflected shock-initiated ignition probed via simultaneous lateral and endwall high-speed imaging with a transparent, cylindrical test-section

A TRAnsparent cylindrical (TRAC) shock tube test section has been developed to explore the ignition structure behind reflected shockwaves via simultaneous lateral and endwall high-speed imaging using three fuel/oxidizer mixtures: 0.1% n-C 7 H 16 /1.10% O 2 /4.14% N 2 /94.66% Ar, 0.1% iC 8 H 18 /1.25% O 2 /4.69% N 2 /93.96% Ar, and 5% CH 4 /10% O 2 /85% CO 2 . During mild ignition, the first exothermic centers were located outside of the observable test section. It is shown that with increasing temperature the far-wall ignition events move closer to the endwall, until a strong ignition is achieved. The temperature gradients that promote mild ignition were quantified and were found to have a slight dependence on fuel with the RMS temperature gradients being 0.80% for n-heptane and 1.07% for isooctane. The impact of temperature gradients on metrics used to improve kinetic models, such as species time histories, is noteworthy. From this temperature gradient, an experimental correlation is provided to estimate the transition temperature from mild to strong ignition for a mixture. The effect of diagnostic location in non-homogeneous reactions is quantified in reference to ignition delay. In conclusion, the unique and insightful perspective on the ignition process under heavy shock bifurcation is also discussed.

02 PETROLEUM↗
DOE OSTI2022

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗
DOE OSTI2023

Dynamics of coupled thermoacoustic modes: Noise and frequency spacing effects

Practical combustion systems consist of multiple thermoacoustic modes that are linearly stable or unstable at a wide range of frequencies. When coupled, the modal interactions are a strong function of the frequency spacing between them – having very significant implications on interpretation of linear stability or instability of modes from experimental data. Frequency spacing between the modes particularly influences the stability and existence of potential limit cycle solutions. In this paper, we extend earlier deterministic studies to include the effects of noise, studying how deterministic dynamics change with the addition of noise and the simultaneous impact of frequency spacing (i.e., closely, or widely spaced). Specifically, this paper characterizes three key effects of noise, which are (1) distribution of limit cycle amplitudes (“diffusion”), (2) shift in the average limit cycle amplitudes (“noise induced drift”), and (3) introduction of new bifurcations that are not present in the deterministic system. Key conclusions from this work are that in highly noisy practical environments, a deterministic/ low noise analysis will not be sufficient to explain even critical qualitative features of the system – such as existence of stable limit cycles or their amplitudes. Furthermore, the nature of noise-induced effects on nonlinear modal interactions can be fundamentally different for degenerate, closely spaced, and widely spaced modes.

42 ENGINEERING↗
DOE OSTI2023

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗
DOE OSTI2024

Structure, dynamics, and redox reactivity of an all-purpose flavodoxin

The flavodoxin of Rhodopseudomonas palustris CGA009 (Rp9Fld) supplies highly reducing equivalents to crucial en- zymes such as hydrogenase, especially when the organism is iron-restricted. By acquiring those electrons from photodriven electron flow via the bifurcating electron transfer flavoprotein, Rp9Fld provides solar power to vital metabolic processes. To understand Rp9Fld’s ability to work with diverse partners, we solved its crystal structure. We observed the canonical fla- vodoxin (Fld) fold and features common to other long-chain Flds but not all the surface loops thought to recognize part- ner proteins. Moreover, some of the loops display alternative structures and dynamics. To advance studies of protein– protein associations and conformational consequences, we assigned the 19F NMR signals of all five tyrosines (Tyrs). Our electrochemical measurements show that incorporation of 3- 19 F-Tyr in place of Tyr has only a modest effect on Rp9Fld’s redox properties even though Tyrs flank the flavin on both sides. Meanwhile, the 19 F probes demonstrate the expected paramagnetic effect, with signals from nearby Tyrs becoming broadened beyond detection when the flavin semiquinone is formed. However, the temperature dependencies of chemical shifts and linewidths reveal dynamics affecting loops close to the flavin and regions that bind to partners in a variety of systems. These coincide with patterns of amino acid type conservation but not retention of specific residues, arguing against detailed specificity with respect to partners. We pro- pose that the loops surrounding the flavin adopt altered con- formations upon binding to partners and may even participate actively in electron transfer.

19F NMR↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

In situ high-temperature 3D imaging of the damage evolution in a SiC nuclear fuel cladding material

Silicon carbide (SiC)-based nuclear fission fuel rod cladding has been considered as one of the possible designs for accident tolerant fuels. It is in the form of a SiC fibre reinforced SiC matrix composite tube (SiC f -SiC m ) with monolithic SiC outer and/or inner coating layers. This study focuses on the deformation and fracture processes in this material using in situ X-ray micro-computed tomography (XCT) at room temperature (RT) and 1200 °C in an inert gas environment in a C-ring compression loading configuration. Prior to testing, local properties and residual stresses were characterised using nanoindentation and Raman spectroscopy since they can impact the mechanical behaviour of the material. The 3D strain distribution, crack formation and propagation processes including the toughening mechanisms (e.g., crack deflection, micro-cracking, crack bridging and bifurcation) are investigated in the coating and underlying composites at RT and 1200 °C. There is no particular sequence which toughening mechanism occurs first – this is very different from the conventional toughening theory in ceramic-matrix composites under uniaxial tension loading. Indeed, no evidence of fibre pull-out or fibre fracture was observed in this SiC f -SiC m nuclear cladding material in the current C-ring compression configuration. The correlation between local measurements and bulk mechanical behaviour are discussed.

3D imaging↗
DOE OSTI2025

Improved damage tolerance of SiC-based nuclear fuel cladding with novel multi-layered SiC coating design at 1200 °C

Continuous SiC fibre reinforced SiC matrix composites (SiC f -SiC m ) with monolithic SiC outer coatings are considered as a damage-tolerant cladding design for loss of coolant accident (LOCA) conditions in light water reactors. However, monolithic SiC coatings are brittle and prone to catastrophic failure. In this study, a SiC f -SiC m cladding with a novel multi-layer SiC outer coating (11 sub-layers, ∼260 µm in total thickness) was investigated under C-ring compression at room temperature and 1200 °C in argon environment. Real-time synchrotron X-ray computed tomography (XCT) was employed to capture crack initiation and propagation processes. Compared to conventional monolithic SiC outer coatings, the multi-layer coating structure facilitated crack deflection and bifurcation enhancing its damage tolerance at both temperatures. Despite pre-existing surface cracks, claddings exhibited stable mechanical-performance at both temperatures. These initial cracks did not critically affect the failure processes as they were not aligned with the maximum stress direction. Furthermore, the microstructure, distribution of residual stresses, and local properties of individual components in the material were thoroughly characterized, and compared with open literature on conventional claddings with monolithic outer coatings. These results provide new insights into the failure mechanisms of multi-layer SiC coatings and offer guidance for the future design of accident-tolerant nuclear fuel claddings.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗