Search NASA⌕ Search

SEARCH · Search NASA

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Spatially quasi-periodic bifurcations from periodic traveling water waves and a method for detecting bifurcations using signed singular values

We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

2D analysis of tokamak divertor-plasma detachment-bifurcation with operational parameters and geometries

UEDGE simulations with density scans for various input power, transport coefficients and outer poloidal leg length are performed to study the conditions for the existence of a bifurcation-like drop of T e at the outer strike point, commonly referred to as a detachment cliff, when transitioning to a detached plasma from an attached plasma in the outer divertor as the upstream density increases (McLean et al., 2015). The simulation results show that a detachment cliff tends to occur with a higher power input regardless of diffusivities and leg length. Further analysis of change of plasma profiles at a cliff indicate that, in addition to the sharp reduction of the E x B drift fluxes in the outer divertor studied in Jaervinen et al., (2018), the substantial change of the Mach number in the outer divertor and the decrease of the outer mid-plane T e due to the radiation front moving across the separatrix into the confinement region above the X-point consistently occur for all UEDGE density scans that have a detachment cliff. UEDGE time-dependent simulation of the evolution of a detachment cliff shows that the rapid increase of radiation above the X-point occurs in a time scale of ~0.3–0.5, which could possibly be the trigger for the formation of a detachment cliff, quicker than the Mach number change in a time scale of ~1 ms and the drop of T e in a time scale of ~2–3 ms in the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Electron Bifurcating Hydrogenases

The importance of electron-bifurcating enzymes is manifest by their ability to maximize energy efficiency. Specifically, they couple a downhill oxidation-reduction (redox) reaction with an uphill redox reaction. Since the rapid increase in the discovery of bifurcating enzymes starting in 2008, there has been interest in incorporating their mechanistic principles into artificial/semiartificial systems to drive chemically challenging reactions. This has yet to be achieved, partly because the details of electron bifurcation, i.e. mechanisms, are largely elusive. Nevertheless, much progress has been made in understanding reactivities, structures, and some mechanistic aspects of these enzymes. Notable examples are electron-bifurcating hydrogenases, which are the focus of this chapter. The chapter is organized as follows. Section 11.1 provides an overview of hydrogenases and electron bifurcation. In Section 11.2, some physiological roles of electron-bifurcating hydrogenases are highlighted. Additionally, electron-bifurcating subunit compositions and biochemical reactivities are comprehensively tabulated, and some key points/considerations about these are noted. In Section 11.3, we discuss the known structures of these enzymes, which provide insight into their complex arrangements of redox cofactors, such as iron-sulfur clusters. Also provided are tabulations and discussions of some biophysical properties of the cofactors. In Section 11.4, we discuss the mechanistic proposals of these enzymes, which are primarily based on structural information. Areas of research that are much needed are outlined in Section 11.5. We conclude on the note that what is learned from electron-bifurcating hydrogenases has applicability to other bifurcating enzymes, nonbifurcating analogs, and mechanistic enzymology at large.

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

The anaerobic fungus Caecomyces churrovis produces H 2 via a non-bifurcating NADH-dependent enzyme complex

ABSTRACT Hydrogenosomes are mitochondria-derived organelles that produce ATP and H 2 to support energy metabolism in anaerobic eukaryotes. H 2 production allows reoxidation of reduced cofactors generated during fermentative metabolism; however, the metabolic mechanisms for H 2 production in anaerobic eukaryotes remains incompletely understood. In particular, it remains unclear whether anaerobic fungi (AF) hydrogenosomes use a ferredoxin-dependent pathway or a distinct mechanism to regenerate NAD(P) + and link electron transfer to H 2 formation. Here, by combining genomic search, proteomic analysis, and enzymology, we reveal the molecular mechanism for H 2 production in the AF strain Caecomyces churrovis . Our enzyme assays on the organelle fraction of C. churrovis revealed the activity of H 2 :NAD + oxidoreductase but not pyruvate:ferredoxin oxidoreductase, which is usually linked to H 2 formation. We identified genes encoding [FeFe] hydrogenase (Hyd) and NADH dehydrogenase subunits E and F (NuoE, NuoF) in C. churrovis , and confirmed their expression in the isolated hydrogenosomal fractions by proteomic analysis. Combining the individually purified enzymes, we found Hyd and NuoEF proteins formed H 2 directly from NADH independently of ferredoxin, functioning as a non-bifurcating NADH-dependent enzyme rather than an electron-bifurcating enzyme. We identified homologs of hydrogenosomal NuoE, NuoF, and Hyd in many other AF, indicating this pathway is commonly shared among the AF. This work demonstrates the existence of a non-bifurcating NADH-dependent enzyme complex in eukaryotes. Moreover, this complex could potentially be exploited as a target for controlling AF H 2 production and altering fungal metabolism. IMPORTANCE H 2 production is a prominent feature of anaerobic energy metabolism, yet our understanding of eukaryotic mechanisms remains limited. Anaerobic fungi (AF) are key decomposers of lignocellulose and contribute to hydrogen flux in anaerobic environments. Although it has been more than 40 years since the H 2 production from AF was first reported, the molecular mechanism for hydrogenosomal H 2 production and redox balance remains unclear. We demonstrate that AF produce H 2 from NADH utilizing a non-bifurcating NADH-dependent enzyme complex rather than an electron-bifurcating, ferredoxin-dependent variant. We show that this enzyme complex is conserved across multiple AF lineages and thus demonstrate the occurrence of a non-bifurcating NADH-dependent enzyme in eukaryotes. This discovery expands our understanding of eukaryotic hydrogenosomal metabolism, reveals a previously unknown strategy for redox balancing, and highlights potential targets for manipulating H 2 production. These insights have broad implications for microbial energy metabolism, anaerobic ecosystems, and bioengineering of H 2 -producing systems.

Zhang, Bo [Department of Chemical Engineering, Uni↗
DOE OSTI2025

Electrochemical Observation and pH Dependence of All Three Expected Redox Couples in an Extremophilic Bifurcating Electron Transfer Flavoprotein with Fused Subunits

Bifurcating enzymes employ energy from a favorable electron transfer to drive unfavorable transfer of a second electron, thereby generating a more reactive product. They are therefore highly desirable in catalytic systems, for example, to drive challenging reactions such as nitrogen fixation. While most bifurcating enzymes contain air-sensitive metal centers, bifurcating electron transfer flavoproteins (bETFs) employ flavins. However, they have not been successfully deployed on electrodes. Herein, we demonstrate immobilization and expected thermodynamic reactivity of a bETF from a hyperthermophilic archaeon, Sulfolobus acidocaldarius (SaETF). SaETF differs from previously biochemically characterized bETFs in being a single protein, representing a concatenation of the two subunits of known ETFs. However, SaETF retains the chemical properties of heterodimeric bETFs, including possession of two FADs: one that undergoes sequential 1-electron (1e) reductions at high E° and forms an anionic semiquinone, and another that is amenable to lower-E° 2e reduction, including by NADH. We found homologous monomeric ETF genes in archaeal and bacterial genomes, accompanied by genes that also commonly flank heterodimeric ETFs, and SaETF’s sequence conservation is 50% higher with bETFs than with canonical ETFs. Thus, SaETF is best described as a bETF. Our direct electrochemical trials capture reversible redox couples for all three thermodynamically expected redox events. We document electrochemical activity over a range of pH values and reveal a conformational change coupled to proton acquisition that affects the electrochemical activity of the higher-E° FAD. Thus, this well-behaved monomeric bETF opens the door to bioinspired bifurcating devices or bifurcation on a chip.

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

Insight Into the High-Potential Branch of the Alternate Nad+-Dependent NADPH:Ferredoxin Oxidoreductase II (NfnII) from Pyrococcus Furiosus. Evidence for the Gating Step in Electron Bifurcation

We have investigated the rapid-reaction kinetics of the NAD+-dependent NADPH:ferredoxin oxidoreductase II (NfnII) from Pyrococcus furiosus, permitting a comparison with recent work done with the paralog NfnI from the same organism. The half-potentials of the electron-bifurcating L-FAD are highly crossed in both NfnI, meaning the potential for the quinone/semiquinone couple is significantly lower than that for the semiquinone/hydroquinone couple so that the semiquinone oxidation state is thermodynamically unstable. The same appears to be the case with NfnII on the basis of its similar behavior in transient absorption spectroscopy experiments and the absence of any evidence for FAD*- accumulation in the course of reductive titrations (which would be manifested as a transient increase in absorbance at ~380 nm)1. Reductive titrations with both the one-electron donor sodium dithionite and the obligate two-electron donor NADPH demonstrate that, in contrast to NfnI, little FADH* accumulates in the high-potential pathway of NfnII in the course of reduction, as reflected in the absence of a transient increase in absorbance in the 500-600 nm region. This indicates that the half-potentials of the S-FAD are crossed in NfnII by a minimum of 120 mV. Rapid-reaction experiments mixing oxidized NfnII with NADPH also show no evidence of S-FADH* accumulation. Furthermore, the enzyme is only partially reduced at the end of the reaction with NADPH, indicating that there is little electron transfer into the high-potential pathway of NfnII. When the reaction is carried out in the presence of the one-electron carrier ferredoxin, only approximately one equivalent of ferredoxin becomes reduced, in contrast to the maximum of three equivalents seen with NfnI. This is consistent with only a single electron bifurcation event taking place under these conditions with NfnII, with only a single high-potential electron passing into the [2Fe-2S] cluster of the high-potential pathway of NfnII but not progressing further to the S-FAD. This accounts for the observation that NfnII has greatly reduced bifurcating activity compared to NfnI. The lack of electron transfer into the S-FAD of NfnII due to its crossed half-potentials prevents successive bifurcation activity and indicates that electron transfer within the high-potential branch gates bifurcation activity in NfnI.

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

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗
DOE OSTI2025

Differential Ligation Alters Electronic State and Coupling Signals of Iron-Sulfur Clusters in Flavin-Based Electron Bifurcation

Flavin-based electron bifurcation (FBEB) is employed by microorganisms for controlling pools of redox equivalents by reversibly splitting electron pairs into high- and low-energy levels from an initial midpoint potential. Our ability to harness this phenomenon is crucial for biocatalytic design which is limited by our understanding of energy coupling in the bifurcation system. In Pyrococcus furiosus, FBEB is carried out by the NADH-dependent ferredoxin:NADP+-oxidoreductase (NfnSL), coupling the uphill reduction of ferredoxin in NfnL to the downhill reduction of NAD+ in NfnS from oxidation of NADPH. Flanking the bifurcating flavin are two site-differentiated iron-sulfur clusters; the nearest is a glutamate-ligated [4Fe-4S] cluster in NfnL. Recent biochemical experiments substituting the native glutamate with cysteine led to loss of coupling between the uphill and downhill pathways, in contrast to the tight thermodynamic coupling in the native system. To understand how this decoupling is biochemically manifested by the cysteine-substituted [4Fe-4S] in NfnL, we employed electron paramagnetic resonance (EPR) spectroscopy to identify changes in electronic architecture and square wave voltammetry (SWV) to probe thermodynamic shifts produced by the substitution. We observed notable g-value shifts in the EPR for the cysteine-substituted iron-sulfur cluster in addition to significant downward shifts in the redox potential, as well as the disappearance of several low-field signals observed in the native NfnSL complex. These results suggest the site-differentiated glutamate residue facilitates higher spin states in the [4Fesingle bond4S] cluster to bridge energetic gaps in electron transfer to the bifurcating flavin in the native complex, preventing unwanted short-circuiting seen in the cysteine-substituted complex.

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

The Sagittarius stream in Gaia Early Data Release 3 and the origin of the bifurcations

The Sagittarius dwarf spheroidal (Sgr) is a dissolving galaxy being tidally disrupted by the Milky Way (MW). Its stellar stream still poses serious modelling challenges, which hinders our ability to use it effectively as a prospective probe of the MW gravitational potential at large radii. Our goal is to construct the largest and most stringent sample of stars in the stream with which we can advance our understanding of the Sgr-MW interaction, focusing on the characterization of the bifurcations. We improved on previous methods based on the use of the wavelet transform to systematically search for the kinematic signature of the Sgr stream throughout the whole sky in the Gaia data. We then refined our selection via the use of a clustering algorithm on the statistical properties of the colour-magnitude diagrams. Our final sample contains more than 700 000 candidate stars and is three times larger than previous Gaia samples. With it, we have been able to detect the bifurcation of the stream in both the northern and southern hemispheres, which requires four branches (two bright and two faint) to fully describe the system. We present the detailed proper motion distribution of the trailing arm as a function of the angular coordinate along the stream, showing, for the first time, the presence of a sharp edge (on the side of the small proper motions) beyond which there are no Sgr stars. We also characterize the correlation between kinematics and distance. Finally, the chemical analysis of our sample shows that the faint branch of the bifurcation is more metal poor than the bright. We provide analytical descriptions for the proper motion trends as well as for the sky distribution of the four branches of the stream. Based on our analysis, we interpret the bifurcations as a misaligned overlap of the material stripped at the antepenultimate pericentre (faint branches) with the stars ejected at the penultimate pericentre (bright branch), given that Sgr just went through its perigalacticon. The source of this misalignment is still unknown, but we argue that models with some internal rotation in the progenitor – at least during the time of stripping of the stars that are now in the faint branches – are worth exploring.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An electron-bifurcating “plug” to a protein nanowire in tungsten-dependent aldehyde detoxification

Members of the tungsten-containing oxidoreductase (WOR) family, which contain a tungstopyranopterin (Tuco) cofactor, are typically either monomeric (WorL) or heterodimeric (WorLS). These enzymes oxidize aldehydes to the corresponding acids while reducing the redox protein ferredoxin. They have been structurally characterized mainly using WORs from hyperthermophilic archaea. The WORs of some bacteria contain three additional subunits of the BfuABC family and these chimeric WorABCSL enzymes catalyze an electron-bifurcating reaction in which aldehyde oxidation is coupled to the simultaneous reduction of ferredoxin and nicotinamide adenine dinucleotide. In human gut microbes, electron bifurcation by WorABSL is proposed to enable the detoxification of aldehydes generated from cooked foods and in the tungstocentric production of beneficial short chain fatty acids from lactate, potentially impacting health. Herein we present the high-resolution cryogenic electron microscopy (cryo-EM) structure of the WorABCSL purified from the bacteriumAcetomicrobium mobile.The structure reveals a surprising 1:3 stoichiometry between WorABC and WorSL, with the WorSL units forming a nanowire-like architecture leading from three Tuco-containing catalytic sites in WorL via strings of multiple iron-sulfur clusters in WorS to a single bifurcating WorABC core. Our structure uncovers a distinct domain arrangement that links three Tuco-dependent aldehyde oxidation sites with the bifurcation process and potentially facilitates environmental aldehyde oxidation.

Science & Technology - Other Topics↗
DOE OSTI2025

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Photoinduced Electron Transfer Informs on Pathway Coupling in Flavin-Based Electron Bifurcation

Flavin-based electron bifurcation (FBEB) is an enzymatic mechanism that generates extremely high-energy electrons to drive unfavorable chemical reactions. It is utilized by the NADHdependent ferredoxin:NADP + -oxidoreductase (Nfn) enzyme in hyperthermophile Pyrococcus furiosus to bifurcate electrons from NADPH into the coupled low-potential (endergonic) and highpotential (exergonic) pathways. This process enables P. furiosus to live in harsh and uninviting environments. Despite its biological importance, the mechanisms used by Nfn to facilitate exceptional directional control over short-lived, high-energy electrons and to prevent undesired transfer, particularly along the low-potential pathway, are still not well understood. To elucidate how the protein environment contributes to electronic control in the lowpotential pathway, new techniques must be utilized to probe these unstable intermediates. In this study, we have adapted lowtemperature photoexcitation combined with electron paramagnetic resonance (EPR) to accumulate the short-lived intermediate and place it in the context of the other cofactors involved in the low-potential pathway of Nfn. We observed coincident growth of both the radical intermediate and its nearby [4Fe-4S] cluster over 4.5 h of illumination with NADPH at cryogenic temperatures. The photogenerated paramagnetic species were stable in LN 2 storage indefinitely and recombined when warmed to higher temperatures. The results provide insights into the electron transfer steps and cofactor interactions along the low potential pathway, facilitating a more robust mechanistic understanding of the high-energy events of electron bifurcation. Furthermore, through comparison of cryogenic and room temperature experiments, a potential gating step involving the movement of key residues important for the reversibility of electron flow along this pathway is suggested.

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

Mechanistic implications of excited high-spin states, spin–spin coupling, and differential [2Fe–2S] + cluster temperature relaxations in the electron-bifurcating NfnABC from Thermococcus sibiricus

Electron bifurcation (EB) is a mechanism of biological energy transduction in which multiple oxidation–reduction (redox) reactions are thermodynamically coupled within a single enzyme, enabling the enzyme to harness the excess free energy from an exergonic process to drive an endergonic process. Because of this unprecedented chemistry, there is interest to translate EB principles to artificial and bioengineered systems, but a hurdle is that knowledge pertaining to the fundamental design principles of EB enzymes remains scarce. Here, we investigated the fundamental physical and electronic properties of electron transfer sites in a spectroscopically uncharacterized member of the BfuABC family of EB enzymes, the NADH-dependent reduced-ferredoxin:NADP + oxidoreductase from Thermococcus sibiricus (Tsi NfnABC). Cryo-EM structures of Tsi NfnABC previously demonstrated that it contains twelve redox cofactors: two flavins (one FAD and one FMN), eight [4Fe–4S] clusters, and two [2Fe–2S] clusters. The FMN, one [4Fe–4S] cluster, and one [2Fe–2S] cluster comprise the bifurcating active site termed the electron-bifurcating flavobicluster (BF-FBC), which is found in all BfuABC family members. By using electron paramagnetic resonance spectroscopy, we identified spectral signatures originating from interactions between the FMN radical and [4Fe–4S] + cluster in the BF-FBC and observed temperature dependent behavior of the BF-FBC's [2Fe–2S] + cluster indicative of moderately slow spin–lattice relaxation. Additionally, we uncovered numerous spectral features corresponding to half-integer, S > ½ spin states of [4Fe–4S] + clusters, including one attributable to the consequences of lysine-ligation of a [4Fe–4S] cluster unique to NfnABC. We contextualize these findings to electron transfer theory and NfnABC's structure. Our insights further the understanding of how enzymes are designed to exert control over electron transfer to conduct thermodynamically challenging reactions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Nonlinear dynamics, bifurcations, and multi-stability in a vibro-impact system with geometric and multi-segmented freeplay nonlinearities

Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. Here, the effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.

42 ENGINEERING↗
DOE OSTI2025

Global River Topology (GRIT): A Bifurcating River Hydrography

Existing global river networks underpin a wide range of hydrological applications but do not represent channels with divergent river flows (bifurcations, multi‐threaded channels, canals), as these features defy the convergent flow assumption that elevation‐derived networks (e.g., HydroSHEDS, MERIT Hydro) are based on. Yet, bifurcations are important features of the global river drainage system, especially on large floodplains and river deltas, and are also often found in densely populated regions. Here we developed the first raster and vector‐based Global RIver Topology that not only represents the tributaries of the global drainage network but also the distributaries, including multi‐threaded rivers, canals and deltas. We achieve this by merging a 30 m Landsat‐based river mask with elevation‐generated streams to ensure a homogeneous drainage density outside of the river mask for rivers narrower than approximately 30 m. Crucially, we employ the new 30 m digital terrain model, FABDEM, based on TanDEM‐X, which shows greater accuracy over the traditionally used SRTM derivatives. After vectorization and pruning, directionality is assigned by a series of elevation, flow angle and continuity approaches. The new global network and its attributes are validated using gauging stations, comparison with existing networks, and randomized manual checks. The new network represents 19.6 million km of streams and rivers with drainage areas greater than 50 km 2 and includes 67,495 bifurcations. With the advent of hyper‐resolution modeling and artificial intelligence, GRIT is expected to greatly improve the accuracy of many river‐based applications such as flood forecasting, water availability and quality simulations, or riverine habitat mapping.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

On the frequency bifurcations of the MHD startup modes in NSTX

The observed bifurcations of the low frequency (<50 kHz) and low toroidal periodicity (n < 5) magnetohydrodynamic (MHD) activity often present in the initial part of the National Spherical Tokamak Experiment (NSTX) discharges can be explained by the evolution of the radial profile of the safety factor (q=rB φ /RB θ ) crossing multiple rational surfaces in the core. Important performance limiting instability mechanisms in the NSTX spherical tokamak are often linked to low frequency and low-n MHD activity. They are quite common in long-pulse NSTX plasmas. They can be present at the beginning of the plasma current flat-top, at the end of the discharge or during the whole duration, and they have been observed to deleteriously impact performance over a wide range of q 95 . An interesting feature observed in some NSTX discharges is the presence of a bifurcation in the frequency of the low n modes, as low as n = 1, that have frequencies comparable to the plasma core rotation divided by n. Equilibrium reconstructions constrained by magnetic diagnostics data and motional stark effect pitch angle radial profiles suggest that the observed bifurcations are linked to a fast evolving minimum value of q. Finally, 3D non-linear resistive MHD simulations show that these modes are ideal and exist as non-resonant before the correspondent rational surface enters the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

E × B flow driven electron temperature bifurcation in a closed slot divertor with ion B × ∇ B away from the X-point in the DIII-D tokamak

An electron temperature bifurcation is observed in the small angle slot divertor, which has been developed to enhance neutral cooling across the divertor target by coupling a closed slot structure with appropriate target shaping. Experiments in the DIII-D tokamak and associated SOLPS-ITER modeling with full drifts find a strong interplay between drifts and divertor geometry on divertor dissipation. The coupling of divertor geometry and drift flows can strongly affect the path towards divertor detachment onset as the plasma density is raised. With the strike point on the inner slanted surface and ion B × ∇B away from the magnetic X-point, bifurcative transitions were observed with sharp decrease of T e towards detachment onset both experimentally and computationally. This differs from the situation for the open divertor where the T e cliff was only observed for ion B × ∇B towards the X-point. SOLPS-ITER modeling with full drifts demonstrates that the magnitude of the E × B drift flow is comparable with the main plasma flow. The reversal of both the poloidal and radial E × B flows near the strike point leads to rapid density accumulation right near the separatrix, which results in bifurcative step transition of divertor conditions with cold plasma across the entire divertor target plate. Furthermore, these results indicate that the interplay between geometry and drifts should be fully taken into account in future fusion reactor divertor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗