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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Experimental investigation of PDI bifurcation of lower hybrid waves during electron density ramp-up in EAST

Abstract The effect of parametric decay instability (PDI) on the current drive efficiency of 4.6 GHz lower hybrid (LH) waves in EAST is investigated experimentally, showing the PDI channel bifurcation of LH waves for the first time in EAST. First, experiments with three platforms of LH power were performed, achieving the LH power required for the PDI occurrence. Second, PDI bifurcation experiments were further carried out by ramping up the plasma electron density. The loop voltage increases with an increase in density, implying a decrease in the current driven by the LH wave. PDI bifurcation during electron density ramp-up was studied by analyzing the parallel refractive index ( n ∥ ) and the frequency spectrum broadening, which is measured with a radio frequency magnetic probe array recently installed close to the LH antenna. It is observed for the first time that they both first increase with density, then there is not much variation and a clear sideband in the frequency spectrum is also observed when the density is up to 4 × 10 19 m − 3 , suggesting a change in the PDI channel. Calculation of the mode growth rate driven by PDI shows that when the edge electron density is up to 1.9 × 10 18 m − 3 , the growth rate of the ion cyclotron quasi-mode (ICQM) will exceed that of the ion sound quasi-mode (ISQM), quantitively explaining that with an increase in density, the PDI channel partly transits from the ISQM to the ICQM channel. Studies provide a possible way to reduce the power deposition in the edge region and improve drive capability by means of mitigating PDI behavior.

Physics↗
DOE OSTI2026

Effect of controlled magnetic island bifurcation on electron diffusion

Magnetic islands strongly influence cross-field electron transport in magnetized plasmas. In particular, bifurcations of the island topology modify the number and location of O-points, X-points, and separatrix boundaries, thereby altering diffusion pathways. In recent DIII-D experiments, external magnetic perturbations were used to rotate and periodically bifurcate the island on the q = 2 surface, causing a switchback between a q = 2/1-dominated structure and a narrower q = 4/2-dominated structure. To investigate how this topological change affects electron transport, we employ the field line tracing code TRIP3D with an implemented collisional operator. Thermal, tracer electrons launched from O-points, X-points, and outside separatrix boundaries reveal distinct diffusion regimes, including classical, subdiffusive, and superdiffusive behavior, depending on both the dominant island mode and launch location. These results suggest that island bifurcation can alter electron diffusion across rational surfaces, with direct implications for particle confinement. While the present work emphasizes diffusion as a general framework, the findings provide insight into the conditions under which electron trapping into an island or stochastization of the island's separatrix can enable additional mechanisms, such as the generation of energetic electrons.

DIII-D↗
DOE OSTI2026

The anaerobic fungus Caecomyces churrovis produces H2 via a non-3 bifurcating NADH-dependent enzyme complex

Anaerobic fungi (AF) decompose lignocellulose-based biomass into fermentable sugars through the production of powerful biomass-degrading enzymes. AF are unusual among fungi in that they generate energy via hydrogenosomes, which are also associated with the release of H2 though yet unknown metabolic mechanisms. In particular, it remains unclear how NAD(P)+ is regenerated within hydrogenosomes and how H2 is formed. Here, we reveal the molecular mechanism for hydrogenosomal H2 production in the AF strain C. churrovis by combining genomic search, proteomic analysis, and enzymology. Our enzyme assays on the large organelle fraction of C. churrovis revealed the activity of H2:NAD+ oxidoreductase but not pyruvate:ferredoxin oxidoreductase activity. 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 proteins, we found that the assay system consisting of Hyd-Strep and NuoEF-Strep reduced NAD+ with H2. Furthermore, this system formed H2 directly from NADH independent 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 widely conserved among the early-branching AF. This work demonstrates the existence of a non-bifurcating NADH-dependent enzyme complex in eukaryotes. Moreover, this complex could be a target for controlling AF H2 production and altering fungal metabolism.

fungi↗
DOE OSTI2023

Spectral Deconvolution of Electron-Bifurcating Flavoproteins

Electron-bifurcating flavoproteins catalyze the tightly coupled reduction of high- and low-potential acceptors using a median-potential electron donor, and are invariably complex systems with multiple redox-active centers in two or more subunits. Methods are described that permit, in favorable cases, the deconvolution of spectral changes associated with reduction of specific centers, making it possible to dissect the overall process of electron bifurcation into individual, discrete steps.

electron bifurcation↗
DOE OSTI2022

Unusual reactivity of a flavin in a bifurcating electron-transferring flavoprotein leads to flavin modification and a charge-transfer complex

From the outset, canonical electron transferring flavoproteins (ETFs) earned a reputation for containing modified flavin. We now show that modification occurs in the recently recognized bifurcating (Bf) ETFs as well. In Bf ETFs, the ’electron transfer’ (ET) flavin mediates single electron transfer via a stable anionic semiquinone state, akin to the FAD of canonical ETFs, whereas a second flavin mediates bifurcation (the Bf FAD). We demonstrate that the ET FAD undergoes trans- formation to two different modified flavins by a sequence of protein-catalyzed reactions that occurs specifically in the ET site, when the enzyme is maintained at pH 9 in an amine-based buffer. Our optical and mass spectrometric characterizations identify 8-formyl flavin early in the process and 8-amino flavins (8AFs) at later times. The latter have not previously been documented in an ETF to our knowledge. Mass spectrometry of flavin products formed in Tris or bis-tris-aminopropane solutions demonstrates that the source of the amine adduct is the buffer. Stepwise reduction of the 8AF demonstrates that it can explain a charge transfer band observed near 726 nm in Bf ETF, as a complex involving the hydroquinone state of the 8AF in the ET site with the oxidized state of unmodified flavin in the Bf site. This supports the possibility that Bf ETF can populate a conformation enabling direct electron transfer between its two flavins, as has been proposed for cofactors brought together in complexes between ETF and its partner proteins.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Bifurcation of equilibrium positions for ellipsoidal particles in inertial shear flows between two walls

We conducted a systematic numerical investigation of spherical, prolate and oblate particles in an inertial shear flow between two parallel walls, using smoothed particle hydrodynamics (SPH). It was previously shown that above a critical Reynolds number, spherical particles experience a supercritical pitchfork bifurcation of the equilibrium position in shear flow between two parallel walls, namely that the central equilibrium position becomes unstable, leading to the emergence of two new off-centre stable positions (Foxet al.,J. Fluid Mech., vol. 915, 2021). This phenomenon was unexpected given the symmetry of the system. In addition to confirming this finding, we found, surprisingly, that ellipsoidal particles can also return to the centre position from the off-centre positions when the particle Reynolds number is further increased, while spherical particles become unstable under this increased Reynolds number. By utilizing both SPH and the finite element method for flow visualization, we explained the underlining mechanism of this reverse of bifurcation by altered streamwise vorticity and symmetry breaking of pressure. Furthermore, we expanded our investigation to include asymmetric particles, a novel aspect that had not been previously modelled, and we observed similar trends in particle dynamics for both symmetric and asymmetric ellipsoidal particles. While further validation through laboratory experiments is necessary, our research paves the road for development of new focusing and separation methods for shaped particles.

Mechanics↗
DOE OSTI2022

Bifurcation-driven vertical plasma displacement

This paper considers vertical plasma motion resulting from plasma current decay during the disruption event. The presented filament-based model describes the motion in the ideal wall limit as an adiabatically slow evolution of the plasma equilibrium. Here, the equilibrium exhibits a pitchfork bifurcation when the decaying plasma current passes a critical value determined by the external magnetic field. This bifurcation affects the disruption-induced mechanical loads on the first wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic-ballooning-bifurcation in tokamak pedestals across shaping and aspect-ratio

We use a new gyrokinetic threshold model to predict a bifurcation in tokamak pedestal width-height scalings that depends strongly on plasma shaping and aspect-ratio. The bifurcation arises from the first and second stability properties of kinetic-ballooning-modes that yields wide and narrow pedestal branches, expanding the space of accessible pedestal widths and heights. The wide branch offers potential for edge-localized-mode-free pedestals with high core pressure. For negative triangularity, low-aspect-ratio configurations are predicted to give steeper pedestals than conventional-aspect-ratio. Both wide and narrow branches have been attained in tokamak experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations

Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Bifurcation from L-mode to internal transport barrier triggered by a magnetic island in tokamak plasmas

Magnetic islands are believed to be responsible for triggering the internal transport barrier (ITB) around the rational surface and hence improves confinement in magnetic confinement fusion plasma, for more than two decades. Although some recent theoretical or modelling works support this hypothesis, the direct experimental demonstration on the ITB triggered by a magnetic island is still missing. Here, we report the first experimental observation that the increase in the magnetic island width triggers the bifurcation of the plasma transport state from a low confinement mode (L-mode) to an ITB mode in tokamak plasma. At medium island widths, the bifurcation appears with a dithering phase, i.e. plasma transitions from L-mode to ITB or vice versa quasi-periodically. As the island width further increases, the dithering ITB transitions into a steady ITB. The electron temperature gradient inside the ITB increases with island width. The trigger of the ITB and further enhancement of the ITB performance are also demonstrated by locking the island via applying a resonant magnetic perturbation field. The reduction of density fluctuations in the ITB region is observed and indicates turbulence suppression during ITB formation. These findings offer new insights for understanding ITB formation and robust ITB control.

internal transport barrier↗
DOE OSTI2022

Bifurcation Analysis Reveals Solution Structures of Phase Field Models

The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems. Here, in this study, we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models. Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions. To elucidate the idea, we apply this analytical approach to three representative phase field equations: the Allen-Cahn equation, the Cahn-Hilliard equation, and the Allen-Cahn-Ohta-Kawasaki system. The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Experimental observations of bifurcated power decay lengths in the near Scrape-Off Layer of ST40 High Field Spherical Tokamak

The scrape-off layer parallel heat flux decay lengths measured at ST40, a high field, low aspect ratio spherical tokamak, have been observed to bifurcate into two groups. The wide group follows established H-mode scalings (ranging between 2 to 8 mm) while the narrow group falls up to 10 times below these scalings (between 0.2 and 0.8 mm), being comparable to the ion total Larmor radius rather than the ion poloidal Larmor radius. The heat flux profiles of the latter group can only be described by a multi-exponential function, rather than the single exponential function convoluted with a Gaussian. The onset of the narrow scrape-off layer width is observed to be associated with suppressed magnetic fluctuations, suggesting reduced electromagnetic turbulence levels in the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗
DOE OSTI2023

Modeling the Effects of Drift Orbit Bifurcation on the Magnetopause Shadowing Loss of Radiation Belt Electrons

Abstract Last closed drift shell (LCDS) has been identified as a crucial parameter for investigating the magnetopause shadowing loss of radiation belt electrons. However, drift orbit bifurcation (DOB) effects have not been physically incorporated into the LCDS calculation. Here we calculate event‐specific LCDS using different approaches to dealing with the DOB effects, that is, tracing field lines ignoring DOB, tracing test particles rejecting field lines with DOB, and tracing particles including field lines with DOB, and then incorporate them into a radial diffusion model to simulate the fast electron dropout observed by Van Allen Probes in May 2017. The model effectively captures the fast dropout at high L* and exhibits the best agreement with data when LCDS is calculated by tracing test particles with DOB more physically included. This study represents the first quantitative modeling of the DOB effects on radiation belt magnetopause shadowing loss via a more physical specification of LCDS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗