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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Axisymmetric shapes and stability of isolated charged drops

Axisymmetric equilibrium shapes and stability of isolated charged drops are found by solving simultaneously the Young-Laplace equation for surface shape and the Laplace equation for the electric field. Families of two-, three-, and four-lobed shapes that branch from the trunk family of spheres are treated systematically by means of the Galerkin/finite element method and a tessellation that deforms with the free surface. The results show that at the limit found by Rayleigh in 1882 the spherical family exchanges stability with a family of two-lobed shapes, a transcritically bifurcating family, one arm of which proves to consist of stable shapes. The results are reinforced by those of approximating the stable drop shapes as oblate spheroids. Thus oblate drops carrying charge in excess of the Rayleigh limit ought to be seen in experiments, though none have yet been reported.

Basaran, O. A.↗
NASA NTRS1989

Shallow cells in directional solidification

The existing theory on two-dimensional transitions (appropriate to thin parallel-plate geometries) is presented in such a way that it is possible to identify easily conditions for the onset of shallow cells. Conditions are given under which succinonitrile-acetone mixtures should undergo supercritical bifurcation in experimentally accessible ranges. These results suggest a means for the quantitative test of the Mullins and Sekerka (1964) model and its weakly nonlinear extensions.

Merchant, G. J.↗
NASA NTRS1989

Effects of mistuning and matrix structure on the topology of frequency response curves

The stability of a frequency response curve under mild perturbations of the system's matrix is investigated. Using recent developments in the theory of singularities of differentiable maps, it is shown that the stability of a response curve depends on the structure of the system's matrix. In particular, the frequency response curves of a cylic system are shown to be unstable. Consequently, slight parameter variations engendered by mistuning will induce a significant difference in the topology of the forced response curves, if the mistuning transformation crosses the bifurcation set.

Afolabi, Dare↗
NASA NTRS1989

Computation of Navier-Stokes equations for three-dimensional flow separation

Supersonic flows over a sharp and a flat-faced blunt fin mounted on a flat plate are simulated numerically. Several basic issues involved in the resultant three-dimensional steady flow separation are studied. Using the same number of grid points, different grid spacings are employed to investigate the effects of a grid resolution on the origin of the line of separation. Various shock strengths are used to study the so-called separated and unseparated boundary layer and to establish the existence or absence of secondary separation. The length of separation ahead of the flat-faced blunt fin, bifurcation of a horseshoe vortex, and the accessibility of a closed-type separation are investigated. The usual interpretation of the flow field from previous studies and new interpretations arising from the present simulation are discussed.

Hung, Ching-Mao↗
NASA NTRS1990

Primary arm spacing in directionally solidified Pb-10 wt percent Sn alloys

The dependence of primary arm spacings on growth speed was investigated for cellular and dendritic arrays in Pb-10 wt percent Sn samples directionally solidified under a constant positive thermal gradient in the melt. The gradient of constitutional supercooling was varied from almost zero (near the break-down of the planar liquid-solid interface at small growth speeds, cellular morphology) to near unity (large growth speeds, dendritic morphology). The spatial arrangements of cells and dendrites, as given by their coordination number, are not very different from each other. It appears that primary arm spacing maxima and the cell to dendrite transition are strongly influenced by the magnitude of the solute partition coefficient. The planar to cellular bifurcation is supercritical in Pb-Sn which has a high partition coefficient, as compared to the subcritical behavior reported in Al-Cu and succinonitrile-acetone, both of which have low partition coefficients. The primary arm spacing model due to Hunt agrees with the experimentally observed trend for the whole growth regime. There is a good quantitative agreement at higher gradients of supercooling. However, the model overpredicts the primary arm spacings at low gradients of constitutional supercooling.

Chopra, M. A.↗
NASA NTRS1990

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗
NASA NTRS1990

On spurious steady-state solutions of explicit Runge-Kutta schemes

The bifurcation diagram associated with the logistic equation v sup n+1 = av sup n (1-v sup n) is by now well known, as is its equivalence to solving the ordinary differential equation u prime = alpha u (1-u) by the explicit Euler difference scheme. It has also been noted by Iserles that other popular difference schemes may not only exhibit period doubling and chaotic phenomena but also possess spurious fixed points. Runge-Kutta schemes applied to both the equation u prime = alpha u (1-u) and the cubic equation u prime = alpha u (1-u)(b-u) were studied computationally and analytically and their behavior was contrasted with the explicit Euler scheme. Their spurious fixed points and periodic orbits were noted. In particular, it was observed that these may appear below the linearized stability limits of the scheme and, consequently, computation may lead to erroneous results.

Sweby, P. K.↗
NASA NTRS1990

H2, fixed architecture, control design for large scale systems

The H2, fixed architecture, control problem is a classic linear quadratic Gaussian (LQG) problem whose solution is constrained to be a linear time invariant compensator with a decentralized processing structure. The compensator can be made of p independent subcontrollers, each of which has a fixed order and connects selected sensors to selected actuators. The H2, fixed architecture, control problem allows the design of simplified feedback systems needed to control large scale systems. Its solution becomes more complicated, however, as more constraints are introduced. This work derives the necessary conditions for optimality for the problem and studies their properties. It is found that the filter and control problems couple when the architecture constraints are introduced, and that the different subcontrollers must be coordinated in order to achieve global system performance. The problem requires the simultaneous solution of highly coupled matrix equations. The use of homotopy is investigated as a numerical tool, and its convergence properties studied. It is found that the general constrained problem may have multiple stabilizing solutions, and that these solutions may be local minima or saddle points for the quadratic cost. The nature of the solution is not invariant when the parameters of the system are changed. Bifurcations occur, and a solution may continuously transform into a nonstabilizing compensator. Using a modified homotopy procedure, fixed architecture compensators are derived for models of large flexible structures to help understand the properties of the constrained solutions and compare them to the corresponding unconstrained ones.

Mercadal, Mathieu↗
NASA NTRS1990

An investigation of chaotic Kolmogorov flows

A two dimensional flow governed by the incompressible Navier-Stokes equations with a steady spatially periodic forcing (known as the Kolmogorov flow) is numerically simulated. The behavior of the flow and its transition states as the Reynolds number (Re) varies is investigated in detail, as well as a number of the flow features. A sequence of bifurcations is shown to take place in the flow as Re varied. Two main regimes of the flow were observed: small and large scale structure regimes corresponding to different ranges of Re. Each of the regimes includes a number of quasiperiodic, chaotic, and relaminarization windows. In addition, each range contains a chaotic window with non-ergodic chaotic attractors. Spatially disordered, but temporally steady states were discovered in large scale structure regime. Features of the diverse cases are displayed in terms of the temporal power spectrum, Poincare sections and, where possible, Lyapunov exponents and Kaplan-Yorke dimension.

Platt, N.↗
NASA NTRS1989

Evolution of the Orszag-Tang vortex system in a compressible medium. I - Initial average subsonic flow

The results of fully compressible, Fourier collocation, numerical simulations of the Orszag-Tang vortex system are presented. The initial conditions for this system consist of a nonrandom, periodic field in which the magnetic and velocity field contain X points but differ in modal structure along one spatial direction. The velocity field is initially solenoidal, with the total initial pressure field consisting of the superposition of the appropriate incompressible pressure distribution upon a flat pressure field corresponding to the initial, average Mach number of the flow. In these numerical simulations, this initial Mach number is varied from 0.2-0.6. These values correspond to average plasma beta values ranging from 30.0 to 3.3, respectively. It is found that compressible effects develop within one or two Alfven transit times, as manifested in the spectra of compressible quantities such as the mass density and the nonsolenoidal flow field. These effects include (1) a retardation of growth of correlation between the magnetic field and the velocity field, (2) the emergence of compressible small-scale structure such as massive jets, and (3) bifurcation of eddies in the compressible flow field. Differences between the incompressible and compressible results tend to increase with increasing initial average Mach number.

Dahlburg, R. B.↗
NASA NTRS1990

Application of dynamical systems theory to the high angle of attack dynamics of the F-14

Dynamical systems theory has been used to study the nonlinear dynamics of the F-14. An eight degree of freedom model that does not include the control system present in operational F-14s has been analyzed. The aerodynamic model, supplied by NASA, includes nonlinearities as functions of the angles of attack and sideslip, the rotation rate, and the elevator deflection. A continuation method has been used to calculate the steady states of the F-14 as continuous functions of the control surface deflections. Bifurcations of these steady states have been used to predict the onset of wing rock, spiral divergence, and jump phenomena which cause the aircraft to enter a spin. A simple feedback control system was designed to eliminate the wing rock and spiral divergence instabilities. The predictions were verified with numerical simulations.

Jahnke, Craig C.↗
NASA NTRS1990

Generation of strong MHD Alfvenic turbulence

Strong Alfvenic turbulence containing a number of solitonlike structures propagating at super-Alfvenic speeds is generated self-consistently and studied by means of computer simulation. A one-dimensional hybrid (kinetic ions, fluid electrons) code is used to investigate the nonlinear evolution of an electromagnetic ion-beam instability that generates low-frequency Alfven-like waves. As the instability develops, the field-aligned hydromagnetic waves steepen, forming a soliton that bifurcates several times, leading to a fully turbulent state.

Akimoto, K.↗
NASA NTRS1989

Origin of the asteroid belt

Earlier and current concepts relevant to the origin of the asteroid belt are discussed and are considered in the framework of the solar system origin. Numerical and analytical solutions of the dynamical theory of planetesimal accumulation are characterized by bifurcations into runaway and nonrunaway solutions, and it is emphasized that the differences in time scales resulting from runaway and nonrunaway growth can be more important than conventional time scale differences determined by heliocentric distances. It is concluded that, in principle, it is possible to combine new calculations with previous work to formulate a theory of the asteroidal accumulation consistent with the meteoritic record and with work on the formation of terrestrial planets. Problems remaining to be addressed before a mature theory can be formulated are discussed.

Wetherill, George W.↗
NASA NTRS1990

Are earthquakes an example of deterministic chaos?

A simple mass-spring model is used to systematically examine the dynamical behavior introduced by fault zone heterogeneities. The model consists of two sliding blocks coupled to each other and to a constant velocity driver by elastic springs. The state of this system can be characterized by the positions of the two blocks relative to the driver. A simple static/dynamic friction law is used. When the system is symmetric, cyclic behavior is observed. For an asymmetric system, where the frictional forces for the two blocks are not equal, the solutions exhibit deterministic chaos. Chaotic windows occur repeatedly between regions of limit cycles on bifurcation diagrams. The model behavior is similar to that of the one-dimensional logistic map. The results provide substantial evidence that earthquakes are an example of deterministic chaos.

Huang, Jie↗
NASA NTRS1989

Transition to two-dimensional turbulent convection in a rapidly rotating annulus

Convection in a rapidly rotating, differentially heated annulus with sloping top and bottom lids is investigated using a semianalytical model. A relatively simple two-dimensional structure is preserved in the experimentally observed flow under rapid rotation, with temporal complexity increasing with the Rayleigh number. The model is, therefore, two-dimensional and exibits a sequence of bifurcations from steadily drifting, azimuthally periodic convection columns (thermal Rossby waves) through vacillation and a period-doubling cascade, to aperiodic weakly turbulent solutions. The results obtained here agree to within a few percent with an earlier limited-resolution two-dimensional model.

Lin, R.-Q.↗
NASA NTRS1988

Statistical detection of the mid-Pleistocene transition

Statistical methods have been used to show quantitatively that the transition in mean and variance observed in delta O-18 records during the middle of the Pleistocene was abrupt. By applying these methods to all of the available records spanning the entire Pleistocene, it appears that this jump was global and primarily represents an increase in ice mass. At roughly the same time an abrupt decrease in sea surface temperature also occurred, indicative of sudden global cooling. This kind of evidence suggests a possible bifurcation of the climate system that must be accounted for in a complete explanation of the ice ages. Theoretical models including internal dynamics are capable of exhibiting this kind of rapid transition.

Maasch, K. A.↗
NASA NTRS1990

Transitions from Hadley to Rossby flows in internally forced rotating spherical systems

The possibility of global-scale transitions between the Hadley and Rossby atmospheric regimes is investiated using a simple three-dimensional rotating spherical model without a boundary layer structure, bottom topography, or cumulus friction, and the expected occurrence of the Hadley to Rossby transition is demonstrated. It is shown that a transition from Hadley flow to wavenumber-5 Rossby flow is preferred, in agreement with standard baroclinic instability results. This result gives a reasonable Rossby wave bifurcation from the Hadley solution. For the cases examined, it was found that the upper symmetric Hadley regime does not exist and that the Hadley to Rossby transition depends on the values of the eddy viscosities.

Higgins, Wayne R.↗
NASA NTRS1990

Eigenfrequencies and horizontal structure of divergent barotropic instability originating in tropical latitudes

A systematic investigation of eigenfunctions of the generalized Laplace tidal equation, for the case of tropical-latitude monotonic mean zonal wind profiles with a single narrow reversed vorticity gradient region, has revealed that, in the limit of low planetary zonal wavenumber, the barotropic instability modes bifurcate into weakly divergent modes of hemispheric scale, on the one hand, and on the other strongly divergent 'internal' modes trapped about the source region. These latter disturbances penetrate into the deep tropics as a Kelvin wave-type phenomenon not previously seen in this context. These results suggest that hemispheric barotropic instability need not be purely nondivergent, and that the growth of weakly divergent modes is preferred.

Dunkerton, Timothy J.↗