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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sensitive dependence on initial conditions in a formation of magnetic vortices

The magnetic vortex exhibits promise as a true random number generator for hardware-based encryption and probabilistic computing due to its stochastic formation of energetically equivalent fourfold degenerate states, characterized by two topologies: polarity and chirality. However, a comprehensive understanding of the stochastic formation of magnetic vortices remains elusive. In this work, we show that the magnetization relaxation in asymmetric Permalloy disks evolves along a pitchfork bifurcation, with both bifurcation paths leading to the formation of magnetic vortices with the same chirality. In the bifurcation, one formation path is always chosen under weak in-plane magnetic fields, ultimately determining the final magnetic vortex state. By delaying the in-plane magnetic field, we quantitatively investigate when the final vortex state is determined and find that it is closely associated with the initial conditions rather than the bifurcation point itself. Our findings provide valuable insights into future spintronic-based encryption and probabilistic computing.

Jeong, Suyeong↗
DOE OSTI2022

Zonal shear layer collapse and the power scaling of the density limit: old L-H wine in new bottles

Edge shear layer collapse causes edge cooling and aggravates radiative effects. This paper details on the microscopic dynamics of the emergence of power (Q) scaling of density limit (DL) from the shear layer collapse transport bifurcation scenario. The analysis is based on a novel 4-field model, which evolves turbulence energy, zonal flow energy, temperature gradient and density, including the neoclassical screening of zonal flow response. Bifurcation analysis yields power scaling of critical density for shear layer collapse as $n_{crit}\sim Q^{1/3}$. The favorable Q scaling of the DL emerges from the fact that the shear layer strength increases with Q, thus preventing shear layer collapse. This in turn reduces particle transport and improves particle confinement. RMP induced ambient stochastic fields degrade the shear layer by inducing decoherence in the Reynolds stress. As a result the particle transport increases and particle confinement degrades. This leads to the emergence of unfavorable stochastic field intensity ($b_{st}^{2}$) scaling of the critical density as $n_{crit}\sim(1+b_{st}^{2})^{-5/3}$. All fields, including zonal flow shear, exhibit hysteresis when the power (Q) is ramped cyclically across the bifurcation point. The hysteresis is due to dynamical delay in bifurcation on account of critical slowing down. Furthermore, the dynamical hysteresis here is fundamentally different from the hysteresis associated with the existence of bi-stable states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

CRNT4SBML: a Python package for the detection of bistability in biochemical reaction networks

Motivation: Signaling pathways capable of switching between two states are ubiquitous within living organisms. They provide the cells with the means to produce reversible or irreversible decisions. Switchlike behavior of biological systems is realized through biochemical reaction networks capable of having two or more distinct steady states which are dependent on initial conditions. Investigation of whether a certain signaling pathway can confer bistability involves a substantial amount of hypothesis testing. The cost of direct experimental testing can be prohibitive. Therefore, constraining the hypothesis space is highly bene?cial. One such methodology is based on Chemical Reaction Network Theory which uses computational techniques to rule out pathways that are not capable of bistability regardless of kinetic constant values and molecule concentrations. Although useful, these methods are complicated from both pureandcomputationalmathematicsperspectives.Thus,theiradoptionisverylimitedamongstbiologists. Results: We brought Chemical Reaction Network Theory approaches closer to experimental biologists by automating all the necessary steps in CRNT4SMBL. The input is based on SBML format which is the community standard for biological pathway communication. The tool parses SBML and derives C-graph representations of the biological pathway with mass action kinetics. Next steps involve an ef?cient search for potential saddle-node bifurcation points using an optimization technique. This type of bifurcation is important as it has the potential of acting as a switching point between two steady states. Finally, if any bifurcation points are present, numerical continuation analysis extends the equilibria branches for generating the diagram. Presence of an S-shaped bifurcation diagram indicates that the pathway acts as a bistable switch for the given optimization parameters. Availability: CRNT4SBML is available via the Python Package Index. The documentation can be found at https://crnt4sbml.readthedocs.io. CRNT4SBML is licensed under the Apache Software License 2.0. Contact: vladislav.petyuk@pnnl.gov

Reyes, Brandon C.↗
DOE OSTI2023

Recombination smooths the time-signal disrupted by latency in within-host HIV phylogenies

Within-host HIV evolution involves several features that may disrupt standard phylogenetic reconstruction. One important feature is re-activation of latently integrated provirus, which has the potential to disrupt the temporal signal, leading to variation in the branch lengths and apparent evolutionary rates in a tree. Yet, real within-host HIV phylogenies tend to show clear, ladder-like trees structured by the time of sampling. Another important feature is recombination, which violates the fundamental assumption that evolutionary history can be represented by a single bifurcating tree. Thus, recombination complicates the within-host HIV dynamic by mixing genomes and creating evolutionary loop structures that cannot be represented in bifurcating trees. In this paper, we develop a coalescent-based simulator of within-host HIV evolution that includes latency, recombination, and effective population size dynamics that allows us to study the relationship between the true, complex genealogy of within-host HIV evolution, encoded as an Ancestral Recombination Graph (ARG), and the observed phylogenetic tree. To compare our ARG results to the familiar phylogeny format, we calculate the expected bifurcating tree after decomposing the ARG into all unique site trees, their combined distance matrix, and the overall corresponding bifurcating tree. While latency and recombination separately disrupt the phylogenetic signal, remarkably, we find that recombination recovers the temporal signal of within-host HIV evolution caused by latency by mixing fragments of old, latent genomes into the contemporary population. In effect, recombination averages over extant heterogeneity, whether it stems from mixed time-signals or population bottlenecks. Further, we establish that the signals of latency and recombination can be observed in phylogenetic trees despite being an incorrect representation of the true evolutionary history. Using an Approximate Bayesian Computation method, we develop a set of statistical probes to tune our simulation model to nine longitudinally-sampled within-host HIV phylogenies. Because ARGs are exceedingly difficult to infer from real HIV data, our simulation system allows investigating effects of latency, recombination, and population size bottlenecks by matching decomposed ARGs to real data as observed in standard phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1974

Plastic buckling

The present article is concerned mainly with the post-bifurcation and imperfection-sensitivity aspects of plastic buckling. A simple two-degree-of-freedom model is used to introduce post-bifurcation behavior and a second model illustrates features of the behavior of continuous solids and structures. Hill's bifurcation criterion for a class of three-dimensional solids is applied to the Donnell-Mushtari-Vlasov (DMV) theory of plates and shells. A general treatment of the initial post-bifurcation behavior of plates and shells is given within the context of the DMV theory. This is illustrated by problems involving columns and circular plates under radial compression. Numerical results are given for a column under axial compression, a circular plate under radial compression, and spherical and cylindrical shells.

Hutchinson, J. W.↗
NASA NTRS1983

Multiple buoyancy driven flows in a vertical cylinder heated from below

The structure of axisymmetric buoyancy-driven convection in a vertical cylinder heated from below is probed by finite element solution of the Boussinesq equations coupled with computed-implemented perturbation techniques for detecting and tracking multiple flows and for determining flow stability. Results are reported for fluids with Prandtl number of one and for cylinders with aspect ratio (Lambda) (defined as the height to radius of the cylinder) between 0.5 and 2.25. Extensive calculations of the neutral stability curve for the static solution and of the nonlinear motions along the bifurcating flow families show a continuous evolution of the primary cellular motion from a single toroidal cell to two and three cells nested radially in the cylinder, instead of the sharp transitions found for a cylinder with shear-free sidewalls. The smooth transitions in flow structure with Rayleigh number and lambda are explained by nonlinear connectivity between the first two bifurcating flow families formed either by a secondary bifurcation point for Lambda or = Lambda * approximately 0.80 or by a limit point for Lambda Lambda *. The transition between these two modes may be described by the theory of multiple limit point bifurcation.

Yamaguchi, Y.↗
NASA NTRS1985

Nonlinear global stability analysis of compressor stall phenomena

Compressor stall phenomena are analyzed from the point of view of nonlinear control theory, based on bifurcation-catastrophe techniques. This new approach appears promising and offers insight into such well-known compressor instability problems as surge and rotating stall and suggests strategies for recovery. Three interlocking dynamic nonlinear state space models are developed. It is shown that the problem of rotating stall can be viewed as an induced bifurcation of solution of the unstalled model. Hysteresis effects are shown to exist in the stall/recovery process. Surge cycles are observed for some critical parameter values. The oscillatory behavior is seen to be due to development of limit cycles, generated by Hopf bifurcation of solutions. More specifically, it is observed that at certain critical values of parameters, a family of stable limit cycles with growning and then diminishing amplitudes is generated, then giving rise to an unstable family of limit cycles. This unstable family in turn bifurcates into other unstable families. To further illustrate the utility of the methodology, some partial computation of domains is carried out, and parameter sensitivity analysis is performed.

Razavi, H.↗
NASA NTRS1991

Nonlinear response of infinitely long circular cylindrical shells to subharmonic radial loads

The nonlinear response of infinitely long circular cylindrical shells (thin circular rings) in the presence of a two-to-one internal (autoparametric) resonance to a subharmonic excitation of order one-half of the higher mode is analyzed with the multiple-scale method. Four autonomous first-order ordinary differential equations are derived for the modulation of the amplitudes and phases of the interacting models. These modulation equations are used to determine the fixed points and their stability. The fixed points correspond to periodic oscillations of the shell, whereas the limit-cycle solutions of the modulation equations correspond to amplitude and phase-modulated oscillations of the shell. The force response curves exhibit saturation, jumps, and Hopf bifurcation. As excitation frequency changes, all limit cycles deform and lose stability through either pitchfork or cyclic-fold (saddle-node) bifurcations. Some of these saddle-node bifurcations cause a transition to chaos. The pitchfork bifurcations break the symmetry of the limit cycles.

Nayfeh, Ali H.↗
NASA NTRS1992

On the structure of cellular solutions in Rayleigh-Benard-Marangoni flows in small-aspect-ratio containers

Multiple steady flow patterns occur in surface-tension/buoyancy-driven convection in a liquid layer heated from below (Rayleigh-Benard-Marangoni flows). Techniques of numerical bifurcation theory are used to study the multiplicity and stability of two-dimensional steady flow patterns (rolls) in rectangular small-aspect-ratio containers as the aspect ratio is varied. For pure Marangoni flows at moderate Biot and Prandtl number, the transitions occurring when paths of codimension 1 singularities intersect determine to a large extent the multiplicity of stable patterns. These transitions also lead, for example, to Hopf bifurcations and stable periodic flows for a small range in aspect ratio. The influence of the type of lateral walls on the multiplicity of steady states is considered. 'No-slip' lateral walls lead to hysteresis effects and typically restrict the number of stable flow patterns (with respect to 'slippery' sidewalls) through the occurrence of saddle node bifurcations. In this way 'no-slip' sidewalls induce a selection of certain patterns, which typically have the largest Nusselt number, through secondary bifurcation.

Dijkstra, Henk A.↗
NASA NTRS2004

Flow field topology of transient mixing driven by buoyancy

Transient mixing driven by buoyancy occurs through the birth of a symmetric Rayleigh-Taylor morphology (RTM) structure for large length scales. Beyond its critical bifurcation the RTM structure exhibits self-similarity and occurs on smaller and smaller length scales. The dynamics of the RTM structure, its nonlinear growth and internal collision, show that its genesis occurs from an explosive bifurcation which leads to the overlap of resonance regions in phase space. This event shows the coexistence of regular and chaotic regions in phase space which is corroborated with the existence of horseshoe maps. A measure of local chaos given by the topological entropy indicates that as the system evolves there is growth of uncertainty. Breakdown of the dissipative RTM structure occurs during the transition from explosive to catastrophic bifurcation; this event gives rise to annihilation of the separatrices which drives overlap of resonance regions. The global bifurcation of explosive and catastrophic events in phase space for the large length scale of the RTM structure serves as a template for which mixing occurs on smaller and smaller length scales. Copyright 2004 American Institute of Physics.

Nonlinear Dynamics↗
NASA NTRS2005

Correlation Between Low Frequency Auroral Kilometric Radiation (AKR) and Auroral Structures

Auroral Kilometric Radiation (AKR) is a radio wave emission that has long been associated with auroral activity. AKR is normally observed in the frequency range from -60 - 600 kHz. Low frequency AKR (or LF-AKR) events are characterized as a rapid extension of AKR related emissions to 30 kHz or lower in frequency for typically much less than 10 minutes. LF-AKR emissions predominantly occur within a frequency range of 20 kHz - 30 kHz, but there are LF-AKR related emissions that reach to a frequency of 5 kHz. This study correlates all instances of LF-AKR events during the first four years of observations from the IMAGE spacecraft's Radio Plasma Imager (WI) instrument with auroral observations from the wideband imaging camera (WIC) onboard IMAGE. The correlation between LF-AKR occurrence and WIC auroral observations shows that in the 295 confirmed cases of LF-AKR emissions, bifurcation of the aurora is seen in 74% of the cases. The bifurcation is seen in the dusk and midnight sectors of the auroral oval, where AKR is believed to be generated. The polarization of these LF-AKR emissions has yet to be identified. Although LF-AKR may not be the only phenomena correlated with bifurcated auroral structures, bifurcation will occur in most instances when LF-AKR is observed. The LF-AKR emissions may be an indicator of specific auroral processes sometimes occurring during storm-time conditions in which field-aligned density cavities extend a distance of perhaps 5-6 RE tailward from the Earth for a period of 10 minutes or less.

Paxamickas, Katherine A.↗
DOE OSTI2022

How a Formate Dehydrogenase Responds to Oxygen: Unexpected O 2 Insensitivity of an Enzyme Harboring Tungstopterin, Selenocysteine, and [4Fe–4S] Clusters

The reversible two-electron interconversion of formate and CO 2 is catalyzed by both nonmetallo- and metallo-formate dehydrogenases (FDHs). The latter group comprises molybdenum- or tungsten-containing enzymes with the metal coordinated by two equivalents of a pyranopterin cofactor, a cysteinyl or selenocysteinyl (Sec) ligand supplied by the polypeptide, and a catalytically essential terminal sulfido ligand. In addition, these biocatalysts incorporate one or more [4Fe–4S] clusters for facilitating long-distance electron transfer. However, an interesting dichotomy arises when attempting to understand how the metallo-FDHs react with O 2 . Whereas existing scholarship portrays these enzymes as being unable to perform in air due to extreme O 2 lability of their metal centers, studies dating as far back as the 1930s emphasize that some of these systems exhibit formate oxidase (FOX) activity, coupling formate oxidation to O 2 reduction. Therefore, to reconcile these conflicting views, we explored context-dependent functional linkages between metallo-FDHs and their cognate electron acceptors within the same organism vis-à-vis catalysis under atmospheric O 2 . Here, we report the discovery and characterization of an O 2 -insensitive FDH2 from the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough (DvH) that ligates tungsten, Sec, and four [4Fe–4S] clusters. By advancing a robust expression platform for its recombinant production, we eliminate both the requirement of nitrate or azide during purification and reductive activation with thiols and/or formate prior to catalysis. Because the distinctive spectral signatures of formate-reduced DvH-FDH2 remain invariant under anaerobic and aerobic conditions, we benchmarked the enzyme activity in air, identifying CO 2 as the catalytic product. Full reaction progress curve analysis discloses a high catalytic efficiency when probed with a high-potential artificial electron acceptor. Furthermore, we show that DvH-FDH2 enables near-stoichiometric hydrogen peroxide production without superoxide release to achieve O 2 insensitivity. Notably, simultaneous electron transfer to cytochrome c and O 2 reveals that metal-based electron bifurcation is operational in this system. Taken together, our work proves the co-occurrence of redox bifurcated FDH and FOX activities within a metalloenzyme scaffold. These findings set the stage for uncovering previously unknown O 2 -insensitive flavin-based electron bifurcation mechanisms, as well as for developing authentic formate/air biofuel cells, engineering O 2 -stable FDHs and biohybrid metallocatalysts, and discerning formate bioenergetics of gut microbiota.

13C NMR↗
DOE OSTI2025

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Experimental signatures of electron cyclotron wave energy condensation in magnetic islands

Here, bifurcation of the electron temperature perturbation (ΔT e ) in magnetic islands due to electron cyclotron current drive (ECCD) is reported for the first time in toroidal plasmas. In this DIIID experiment ΔT e spontaneously increases by 50% in a growing 2,1 island which in turn rapidly shrinks, consistent with improved stabilization by the ECCD due to the temperature bifurcation. This non-linearity of ΔT e maps out a hysteresis in the island evolution cycle, in qualitative agreement with the theoretically predicted radio frequency current condensation effect. Therefore, this observed bifurcation can be critical for the avoidance of disruptions caused by tearing modes in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Experiments on stiffened conical shell structures using cast epoxy models

Description of a casting technique for fabricating high-quality plastic structural models, and review of results regarding the use of such specimens to parametrically study the effect of base ring stiffness on the critical buckling pressure of a ring-stiffened conical shell. The fabrication technique involves machining a metal mold to the desired configuration and vacuum-drawing the plastic material into the mold. A room-temperature curing translucent thermoset epoxy was the casting material selected. A shell of revolution computer program which employs a nonlinear axisymmetric prebuckling strain field to obtain a bifurcation buckling solution was used to guide the selection of congifurations tested. The shell experimentally exhibited asymmetric collapse behavior, and the ultimate load was considerably higher than the analytical bifurcation prediction. The asymmetric buckling mode shape, however, initially appeared at a pressure near the analysis bifurcation solution.

Williams, J. G.↗
NASA NTRS1975

Sustained small oscillations in nonlinear control systems

Some results of bifurcation theory were used to study the existence of small-amplitude periodic behavior in launch vehicle dynamics, assuming that nonlinearity exists as a cubic term in the rudder response. The analysis follows closely Sattinger's (1973) approach to the theory of periodic bifurcations. The conditions under which a bifurcating branch of orbitally stable periodic solutions will exist are determined. It is shown that in more complicated cases, the conditions under which the system matrix has a pair of simple purely imaginary eigenvalues can be determined with the aid of linear stability techniques.

George, J. H.↗
NASA NTRS1979

Some potential blood flow experiments for space

Blood is a colloidal suspension of cells, predominantly erythrocytes, (red cells) in an aqueous solution called plasma. Because the red cells are more dense than the plasma, and because they tend to aggregate, erythrocyte sedimentation can be significant when the shear stresses in flowing blood are small. This behavior, coupled with equipment restrictions, has prevented certain definitive fluid mechanical studies from being performed with blood in ground-based experiments. Among such experiments, which could be satisfactorily performed in a microgravity environment, are the following: (1) studies of blood flow in small tubes, to obtain pressure-flow rate relationships, to determine if increased red cell aggregation can be an aid to blood circulation, and to determine vessel entrance lengths, and (2) studies of blood flow through vessel junctions (bifurcations), to obtain information on cell distribution in downstream vessels of (arterial) bifurcations, and to test flow models of stratified convergent blood flows downstream from (venous) bifurcations.

Cokelet, G. R.↗
NASA NTRS1981

Chaotic /strange/ and periodic behavior in instability saturation by the oscillating two-stream instability

The nonlinear Schroedinger equation with linear growth and damping is truncated to three waves. The resulting system of nonlinear ordinary differential equations describes the excitation of linearly damped waves by the oscillating two-stream instability driven by a linearly unstable pump wave. This system represents a simple model for the nonlinear saturation of a linearly unstable wave. The model is examined analytically and numerically as a function of the dimensionless parameters of the system. It is found that the model can exhibit a wealth of characteristic dynamical behavior including stationary equilibria, Hopf bifurcations to periodic orbits, period doubling bifurcations, chaotic solutions characteristic of a strange attractor, tangent bifurcations from chaotic to periodic solutions, transient chaos, and hysteresis. Many of these features are shown to be explainable on the basis of one-dimensional maps. In the case of chaotic solutions, evidence for the presence of a strange attractor is provided by demonstrating Cantor set-like structure (i.e., scale invariance) in the surface of section.

Russell, D. A.↗