Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bifurcation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
DOE OSTI2024

Crystal structure of indacaterol hydrogen maleate (C 24 H 29 N 2 O 3 )(HC 4 H 2 O 4 )

The crystal structure of indacaterol hydrogen maleate has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional techniques. Indacaterol hydrogen maleate crystallizes in space groupP-1 (#24) witha= 8.86616(9),b= 9.75866(21),c= 16.67848(36) Å,α= 102.6301(10), β = 94.1736(6),γ= 113.2644(2)°,V= 1273.095(7) Å 3 , andZ= 2 at 295 K. The crystal structure consists of layers of cations and anions parallel to theab-plane. Traditional N–H⋯O and O–H⋯O hydrogen bonds link the cations and anions into chains along thea-axis. There is a strong intramolecular charge-assisted O–H⋯O hydrogen bond in the non-planar hydrogen maleate anion. There are also two C–H⋯O hydrogen bonds between the anion and cation. The cation makes a strong N–H⋯O hydrogen bond to the anion, but also acts as a hydrogen bond donor to an aromatic C in another cation. The amino group makes bifurcated N–H⋯O hydrogen bonds, one intramolecular and the other intermolecular. The hydroxyl group acts as a donor to another cation. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

Materials Science↗
DOE OSTI2020

Controls on Iron Reduction and Biomineralization over Broad Environmental Conditions as Suggested by the Firmicutes Orenia metallireducens Strain Z6

Microbial iron reduction is a ubiquitous biogeochemical process driven by diverse microorganisms in a variety of environments. However, it is often difficult to separate the biological from the geochemical controls on bioreduction of Fe(III) oxides. In this study, we investigated the primary driving factor(s) that mediate secondary iron mineral formation over a broad range of environmental conditions using a single dissimilatory iron reducer, Orenia metallireducens strain Z6. A total of 17 distinct geochemical conditions were tested with differing pH (6.5–8.5), temperature (22–50 °C), salinity (2–20% NaCl), anions (phosphate and sulfate), electron shuttle (anthraquinone-2,6-disulfonate), and Fe(III) oxide mineralogy (ferrihydrite, lepidocrocite, goethite, hematite, and magnetite). The observed rates and extent of iron reduction differed significantly with k int between 0.186 and 1.702 mmol L –1 day –1 and Fe(II) production ranging from 6.3% to 83.7% of the initial Fe(III). Using X-ray absorption and scattering techniques (EXAFS and XRD), we identified and assessed the relationship between secondary minerals and the specific environmental conditions. It was inferred that the observed bifurcation of the mineralization pathways may be mediated by differing extents of Fe(II) sorption on the remaining Fe(III) minerals. These results expand our understanding of the controls on biomineralization during microbial iron reduction and aid the development of practical applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Photoelectron Spectroscopy and Theoretical Study on Monosolvated Cyanate Analogue Clusters ECX-· Sol (ECX- = NCSe-, AsCSe- and AsCS-; Sol = H2O, CH3CN)

Six monosolvated cyanate analogue clusters ECX-?Sol (ECX– = NCSe-, AsCSe-, and AsCS-; Sol = H2O and CH3CN) were investigated using negative ion photoelectron spectroscopy (NIPES). NIPES experiments show that these clusters possess overall similar spectra compared to their respective isolated ECX- anions but shift to higher electron binding energy with solvent CH3CN stabilizing the excess electrons slightly more than H2O. For the ECX-·H2O series, vertical detachment energies and their increments relative to the bare species are measured to be 3.700 / 0.370, 3.085 / 0.425, 3.085 / 0.430 eV for NCSe-, AsCSe- and AsCS-, respectively, while the corresponding values in ECX-·CH3CN series are 3.835 / 0.505, 3.145 / 0.475, and 3.135 / 0.480 eV. Ab initio electronic structure calculations indicate the excess charges located at the terminal N and Se atoms in NCSe- and migrated to the central C atom in AsCSe- and AsCS-. For NCSe-, the solvation is driven by the interactions with the two negatively charged terminal ends, while for AsCSe- and AsCS-, the solvation revolves around the interactions with the central C atom, where all the excess negative charge is concentrated. Two nearly degenerate isomers for NCSe-·H2O are identified, one forming a single strong N…H-O hydrogen bond (HB) and the other featuring a bidentate HB with two hydroxyl H atoms pointing to N and Se ends. In contrast, the negative central C atom in AsCSe-/AsCS- allows the formation of a bifurcated HB with H2O. Similar effects are observed for the acetonitrile case, in which the three H atoms of methyl group interact with the two negatively charged terminal ends in NCSe-, while prefer binding to the central negative carbon atom in AsCSe-/AsCS-. The different binding motifs derived in this work may suggest distinctly different solvation properties in NCSe- versus AsCSe-/AsCS-, with the former anion leading to asymmetric solvation at the N end of the solute, while the latter species creating more ‘isotropic’ solvation around the central C equatorial plane. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer.

Yuan, Qinqin↗
DOE OSTI2022

Scaling of Shear Rheology of Concentrated Charged Colloidal Suspensions across Glass Transition

Electrostatic interparticle interactions are a key component in controlling and designing rheological characteristics of concentrated charged colloidal suspensions. Herein, we investigate electroviscous effects on shear rheology using highly charged silica particles. By fixing the volume fraction but varying the salinity, the system undergoes a glass transition as evidenced by the evolution of the yield stress and zero-shear viscosity. We show that the steady shear viscosities obey a critical scaling relation that scales the flow curves into a super- and a sub-critical branch with glass transition salinity serving as the bifurcation point; we also demonstrate an isoviscosity scaling that collapses all isoviscosity lines into a single master curve that exhibits no singularity. Based on each scaling relation, in conjunction with common modeling equations, the quantitative relationships between the shear viscosity, stress, and salinity are established. This study demonstrates a new framework to model the steady shear rheology of concentrated charged colloids.

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

Amorphous CeO 2 –Cu Heterostructure Enhances CO 2 Electroreduction to Multicarbon Alcohols

Electrochemical conversion of carbon dioxide (CO 2 ) gas to value-added chemicals such as multicarbon (C 2+ ) alcohols is a promising and attractive decarbonization strategy. However, there are tremendous challenges in tuning the intrinsic activity and selectivity of the catalysts to produce C 2+ alcohols. In this work, we prepared a CeO 2 –Cu composite catalyst via a combination of metallurgy and dealloying method. The interfacial sites of amorphous CeO 2 –Cu heterostructure improve the adsorption of key reaction intermediates *CO and promote the C–C coupling. Significantly, they also stabilize *CH 2 CHO at the bifurcation step, steering the reaction pathway toward the formation of C 2+ alcohols over ethylene. The CeO 2 –Cu catalyst achieves a remarkable faradaic efficiency of 32.9% ± 2.6% for C 2+ alcohols at -0.6 V vs RHE. Overall, this work demonstrates an effective strategy of improving the intrinsic activity and selectivity of the Cu-based catalysts for the generation of C 2+ alcohols.

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

Divergent Bimetallic Mechanisms in Copper(II)-Mediated C–C, N–N, and O–O Oxidative Coupling Reactions

Copper-catalyzed aerobic oxidative coupling of diaryl imines provides a route for conversion of ammonia to hydrazine. Here the present study uses experimental and density functional theory computational methods to investigate the mechanism of N–N bond formation, and the data support a mechanism involving bimolecular coupling of Cu-coordinated iminyl radicals. Computational analysis is extended to Cu II -mediated C–C, N–N, and O–O coupling reactions involved in the formation of cyanogen (NC–CN) from HCN, 1,3-butadiyne from ethyne (i.e., Glaser coupling), hydrazine from ammonia, and hydrogen peroxide from water. The results reveal two different mechanistic pathways. Heteroatom ligands with an uncoordinated lone pair (iminyl, NH 2 , OH) undergo charge transfer to Cu II , generating ligand-centered radicals that undergo facile bimolecular radical–radical coupling. Ligands lacking a lone pair (CN and CCH) form bridged binuclear diamond-core structures that undergo C–C coupling. This mechanistic bifurcation is rationalized by analysis of spin densities in key intermediates and transition states, as well as multiconfigurational calculations. Radical–radical coupling is especially favorable for N–N coupling owing to energetically favorable charge transfer in the intermediate and thermodynamically favorable product formation.

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

Tunable Multisite Proton-Coupled Electron Transfer Mediators: Distinct Pathways for Substrate Reduction Versus Competing Hydrogen Evolution

Proton-coupled electron transfer (PCET) reagents have emerged as powerful tools for transferring net H atoms to organic substrates from relatively weak X–H bonds. One advantage of employing PCET reagents is the tunability of the X–H bond strength by independently varying their redox potential and/or p K a for selective substrate reductions; however, the rational development of modular catalytic PCET reagents based on these features remains underdeveloped. In this work, we address important mechanistic questions relevant to a dimethylaniline-appended cobaltocene PCET mediator that our lab has previously introduced. Specifically, we examine where protonation occurs within the reactive Co(II, NH) + intermediate of a Brønsted-base modified cobaltocene mediator, whether substrate reduction and hydrogen evolution reaction (HER) proceed by a common or bifurcated mechanistic pathway, and how the redox, acid–base, and structural properties of PCET mediators can dictate their reactivity and selectivity. We show that substrate compatibility can be tuned and, via a model study with N -aryl imine substrates, provide data pointing to a multisite PCET (MS-PCET) pathway. Moreover, we rigorously characterize the site of protonation in the reactive reduced, protonated form of the mediator, and through kinetic analysis establish that the pathway for undesired competing HER is fundamentally different and involves Cp-ring protonation. Our findings point to a high degree of flexibility in the design of reductive PCET mediators.

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

Selectivity-Determining Steps in O 2 Reduction Catalyzed by Iron(tetramesitylporphyrin)

The oxygen reduction reaction (ORR) is the cathode reaction in fuel cells and its selectivity for water over hydrogen peroxide production is important for these technologies. Iron porphyrin catalysts have long been studied for the ORR, but the origins of their selectivity are not well understood because the selectivity-determining step(s) usually occur after the rate-determining step. We report here the effects of acid concentration, as well as other solution conditions such as acid pK a , on the H 2 O 2 /H 2 O selectivity in electrocatalytic ORR by iron(tetramesitylporphyrin) (Fe(TMP)) in DMF. The results show that selectivity reflects a kinetic competition in which the dependence on [HX] is one order greater for the production of H 2 O than H 2 O 2 . Based on such experimental results and computational studies, we propose that the selectivity is governed by competition between protonation of the hydroperoxo intermediate, Fe III (TMP)(OOH), to produce water versus dissociation of the HOO – ligand to yield H 2 O 2 . The data rule out a bifurcation based on the regioselectivity of protonation of the hydroperoxide, as suggested in the enzymatic systems. Furthermore, the analysis developed in this paper should be generally valuable to the study of selectivity in other multi-proton/multi-electron electrocatalytic reactions.

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

Evolving CO 2 rather than SST leads to a factor of ten decrease in GCM convergence time

The high computational cost of Global Climate Models (GCMs) is a problem that limits their use in many areas. Recently an inverse climate modeling (InvCM) method, which fixes the global mean sea surface temperature (SST) and evolves the CO 2 mixing ratio to equilibrate climate, has been implemented in a cloud resolving model. In this paper, we apply InvCM to ExoCAM GCM aquaplanet simulations, allowing the SST pattern to evolve while maintaining a fixed global-mean SST. We find that InvCM produces the same climate as normal slab-ocean simulations but converges an order of magnitude faster. We then use InvCM to calculate the equilibrium CO 2 for SSTs ranging from 290 K to 340 K at 1 K intervals and reproduce the large increase in climate sensitivity at an SST of about 315 K at much higher temperature resolution. The speedup provided by InvCM could be used to equilibrate GCMs at higher spatial resolution or to perform broader parameter space exploration in order to gain new insight into the climate system. Additionally, InvCM could be used to find unstable and hidden climate states, and to find climate states close to bifurcations such as the runaway greenhouse transition.

54 ENVIRONMENTAL SCIENCES↗