Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bifurcation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Solutal convection during directional solidification

During directional solidification of a binary alloy at constant velocity, buoyancy-driven fluid flow may occur due to the solute gradients generated by the solidification process. Numerical calculations of the solute and fluid flow fields in the melt have been carried out using finite differences in a two-dimensional, time-dependent model that assumes a planar crystal-melt interface and allows time-dependent gravitational accelerations. The container walls are rigid and perfectly insulating for solute. For constant vertical gravitational accelerations, as the solutal Rayleigh number is varied, multiple steady states and time-dependent states may occur. The bifurcation from the quiescent state may be subcritical or transcritical, depending on the aspect ratio of the container. Calculations have also been performed for a gravitational acceleration that is assumed to be uniform in magnitude with its direction rotating uniformly. Numerical results have been obtained for a Schmidt number of 10 and a gravitational acceleration of 0.0001 G. The maximum variation in the solute concentration at the crystal-melt interface is calculated for various values of the rotation rate of the gravitational acceleration.

Mcfadden, G. B.↗
NASA NTRS1988

On the nonlinear interaction of Goertler vortices and Tollmien-Schlichting waves in curved channel flows at finite Reynolds numbers

The flow in a two-dimensional curved channel driven by an azimuthal pressure gradient can become linearly unstable due to axisymmetric perturbations and/or nonaxisymmetric perturbations depending on the curvature of the channel and the Reynolds number. For a particular small value of curvature, the critical neighborhood of this curvature value and critical Reynolds number, nonlinear interactions occur between these perturbations. The Stuart-Watson approach is used to derive two coupled Landau equations for the amplitudes of these perturbations. The stability of the various possible states of these perturbations is shown through bifurcation diagrams. Emphasis is given to those cases which have relevance to external flows.

Daudpota, Q. Isa↗
NASA NTRS1988

Numerical results for axial flow compressor instability

Using Cornell's supercomputing facilities, an extensive study of the Moore-Greitzer model was carried out, which gives accurate and reliable information about compressor instability. The bifurcation analysis in the companion paper shows the dependence of the mode of compressor response on the shape of the rotating stall characteristic. The numerical results verify and extend this with a more accurate representation of the characteristic. The effect of the parameters on the shape of the rotating stall characteristic is investigated, and it is found that the parameters with the strongest effects are the inlet length, and the shape of the compressor pressure rise vs. mass flow diagram (i.e. tall diagrams vs. shallow diagrams). The effects of inlet guide vane loss on the characteristic are discussed.

Mccaughan, F. E.↗
NASA NTRS1988

The influence of and the identification of nonlinearity in flexible structures

Several models were built at NASA Langley and used to demonstrate the following nonlinear behavior: internal resonance in a free response, principal parametric resonance and subcritical instability in a cantilever beam-lumped mass structure, combination resonance in a parametrically excited flexible beam, autoparametric interaction in a two-degree-of-freedom system, instability of the linear solution, saturation of the excited mode, subharmonic bifurcation, and chaotic responses. A video tape documenting these phenomena was made. An attempt to identify a simple structure consisting of two light-weight beams and two lumped masses using the Eigensystem Realization Algorithm showed the inherent difficulty of using a linear based theory to identify a particular nonlinearity. Preliminary results show the technique requires novel interpretation, and hence may not be useful for structural modes that are coupled by a guadratic nonlinearity. A literature survey was also completed on recent work in parametrically excited nonlinear system. In summary, nonlinear systems may possess unique behaviors that require nonlinear identification techniques based on an understanding of how nonlinearity affects the dynamic response of structures. In this was, the unique behaviors of nonlinear systems may be properly identified. Moreover, more accutate quantifiable estimates can be made once the qualitative model has been determined.

Zavodney, Lawrence D.↗
NASA NTRS1988

Radar investigation of asteroids and planetary satellites

The aim is to make radar reconnaissance of near-Earth asteroids, mainbelt ateroids, the Galilean satellites, the Martian satellites, and the largest Saturnian satellites, using the Arecibo 13-cm and the Goldstone 3.5-cm systems. Measurements of echo strength, polarization, and delay/Doppler distribution of echo power provide information about dimensions, spin vector, large-scale topography, cm-to-m-scale morphology, and surface bulk density. The observations also yield refined estimates of target orbital elements. Radar signatures were measured for 31 mainbelt asteroids and 16 near-Earth asteroids since this task began eight years ago. The dispersion in asteroid radar albedoes and circular polarization ratios is extreme, revealing huge differences in surface morphologies, bulk densities, and metal concentration. For the most part, correction between radar signature and VIS/IR class is not high. Many near-Earth asteroids have extremely irregular, nonconvex shapes, but some have polar silhouettes that appear only slightly noncircular. The signatures of 1627 Ivar, 1986 DA, and the approximately 180-km mainbelt asteroid 216 Kleopatra suggest bifurcated shapes. Observational milestones during 1987 and 1988 are noted.

Ostro, Steven J.↗
NASA NTRS1989

Low Reynolds number numerical solutions of chaotic flow

Numerical computations of two-dimensional flow past an airfoil at low Mach number, large angle of attack, and low Reynolds number are reported which show a sequence of flow states leading from single-period vortex shedding to chaos via the period-doubling mechanism. Analysis of the flow in terms of phase diagrams, Poincare sections, and flowfield variables are used to substantiate these results. The critical Reynolds number for the period-doubling bifurcations is shown to be sensitive to mesh refinement and the influence of large amounts of numerical dissipation. In extreme cases, large amounts of added dissipation can delay or completely eliminate the chaotic response. The effect of artificial dissipation at these low Reynolds numbers is to produce a new effective Reynolds number for the computations.

Pulliam, Thomas H.↗
NASA NTRS1989

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Yung, Chain-Nan↗
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.↗