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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Nonlinear Dynamic Modeling and Controls Development for Supersonic Propulsion System Research

This paper covers the propulsion system component modeling and controls development of an integrated nonlinear dynamic simulation for an inlet and engine that can be used for an overall vehicle (APSE) model. The focus here is on developing a methodology for the propulsion model integration, which allows for controls design that prevents inlet instabilities and minimizes the thrust oscillation experienced by the vehicle. Limiting thrust oscillations will be critical to avoid exciting vehicle aeroelastic modes. Model development includes both inlet normal shock position control and engine rotor speed control for a potential supersonic commercial transport. A loop shaping control design process is used that has previously been developed for the engine and verified on linear models, while a simpler approach is used for the inlet control design. Verification of the modeling approach is conducted by simulating a two-dimensional bifurcated inlet and a representative J-85 jet engine previously used in a NASA supersonics project. Preliminary results are presented for the current supersonics project concept variable cycle turbofan engine design.

Connolly, Joseph W.↗
NASA NTRS2003

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗
NASA NTRS2012

A Cross-Validation Approach to Approximate Basis Function Selection of the Stall Flutter Response of a Rectangular Wing in a Wind Tunnel

The stall flutter response of a rectangular wing in a low speed wind tunnel is modelled using a nonlinear difference equation description. Static and dynamic tests are used to select a suitable model structure and basis function. Bifurcation criteria such as the Hopf condition and vibration amplitude variation with airspeed were used to ensure the model was representative of experimentally measured stall flutter phenomena. Dynamic test data were used to estimate model parameters and estimate an approximate basis function.

Kukreja, Sunil L.↗
NASA NTRS2012

Methods for determining the genetic affinity of microorganisms and viruses

Selecting which sub-sequences in a database of nucleic acid such as 16S rRNA are highly characteristic of particular groupings of bacteria, microorganisms, fungi, etc. on a substantially phylogenetic tree. Also applicable to viruses comprising viral genomic RNA or DNA. A catalogue of highly characteristic sequences identified by this method is assembled to establish the genetic identity of an unknown organism. The characteristic sequences are used to design nucleic acid hybridization probes that include the characteristic sequence or its complement, or are derived from one or more characteristic sequences. A plurality of these characteristic sequences is used in hybridization to determine the phylogenetic tree position of the organism(s) in a sample. Those target organisms represented in the original sequence database and sufficient characteristic sequences can identify to the species or subspecies level. Oligonucleotide arrays of many probes are especially preferred. A hybridization signal can comprise fluorescence, chemiluminescence, or isotopic labeling, etc.; or sequences in a sample can be detected by direct means, e.g. mass spectrometry. The method's characteristic sequences can also be used to design specific PCR primers. The method uniquely identifies the phylogenetic affinity of an unknown organism without requiring prior knowledge of what is present in the sample. Even if the organism has not been previously encountered, the method still provides useful information about which phylogenetic tree bifurcation nodes encompass the organism.

Fox, George E.↗
NASA NTRS2012

Design and Analysis Tool for External-Compression Supersonic Inlets

A computational tool named SUPIN has been developed to design and analyze external-compression supersonic inlets for aircraft at cruise speeds from Mach 1.6 to 2.0. The inlet types available include the axisymmetric outward-turning, two-dimensional single-duct, two-dimensional bifurcated-duct, and streamline-traced Busemann inlets. The aerodynamic performance is characterized by the flow rates, total pressure recovery, and drag. The inlet flowfield is divided into parts to provide a framework for the geometry and aerodynamic modeling and the parts are defined in terms of geometric factors. The low-fidelity aerodynamic analysis and design methods are based on analytic, empirical, and numerical methods which provide for quick analysis. SUPIN provides inlet geometry in the form of coordinates and surface grids useable by grid generation methods for higher-fidelity computational fluid dynamics (CFD) analysis. SUPIN is demonstrated through a series of design studies and CFD analyses were performed to verify some of the analysis results.

Slater, John W.↗
NASA NTRS2013

Buckling Analysis of a Honeycomb-Core Composite Cylinder with Initial Geometric Imperfections

Thin-walled cylindrical shell structures often have buckling as the critical failure mode, and the buckling of such structures can be very sensitive to small geometric imperfections. The buckling analyses of an 8-ft-diameter, 10-ft-long honeycomb-core composite cylinder loaded in pure axial compression is discussed in this document. Two loading configurations are considered configuration 1 uses simple end conditions, and configuration 2 includes additional structure that may more closely approximate experimental loading conditions. Linear eigenvalue buckling analyses and nonlinear analyses with and without initial geometric imperfections were performed on both configurations. The initial imperfections were introduced in the shell by applying a radial load at the midlength of the cylinder to form a single inward dimple. The critical bifurcation buckling loads are predicted to be 924,190 lb and 924,020 lb for configurations 1 and 2, respectively. Nonlinear critical buckling loads of 918,750 lb and 954,900 lb were predicted for geometrically perfect configurations 1 and 2, respectively. Lower-bound critical buckling loads for configurations 1 and 2 with radial perturbations were found to be 33% and 36% lower, respectively, than the unperturbed critical loads. The inclusion of the load introduction cylinders in configuration 2 increased the maximum bending-boundary-layer rotation up to 11%.

Cha, Gene↗
NASA NTRS2013

Structural Assessment of Advanced Composite Tow-Steered Shells

The structural performance of two advanced composite tow-steered shells, manufactured using a fiber placement system, is assessed using both experimental and analytical methods. The fiber orientation angles vary continuously around the shell circumference from 10 degrees on the shell crown and keel, to 45 degrees on the shell sides. The two shells differ in that one shell has the full 24-tow course applied during each pass of the fiber placement system, while the second shell uses the fiber placement system s tow drop/add capability to achieve a more uniform shell wall thickness. The shells are tested in axial compression, and estimates of their prebuckling axial stiffnesses and bifurcation buckling loads are predicted using linear finite element analyses. These preliminary predictions compare well with the test results, with an average agreement of approximately 10 percent.

Wu, K. Chauncey↗
NASA NTRS2013

An Approximate Solution and Master Curves for Buckling of Symmetrically Laminated Composite Cylinders

Nondimensional linear-bifurcation buckling equations for balanced, symmetrically laminated cylinders with negligible shell-wall anisotropies and subjected to uniform axial compression loads are presented. These equations are solved exactly for the practical case of simply supported ends. Nondimensional quantities are used to characterize the buckling behavior that consist of a stiffness-weighted length-to-radius parameter, a stiffness-weighted shell-thinness parameter, a shell-wall nonhomogeneity parameter, two orthotropy parameters, and a nondimensional buckling load. Ranges for the nondimensional parameters are established that encompass a wide range of laminated-wall constructions and numerous generic plots of nondimensional buckling load versus a stiffness-weighted length-to-radius ratio are presented for various combinations of the other parameters. These plots are expected to include many practical cases of interest to designers. Additionally, these plots show how the parameter values affect the distribution and size of the festoons forming each response curve and how they affect the attenuation of each response curve to the corresponding solution for an infinitely long cylinder. To aid in preliminary design studies, approximate formulas for the nondimensional buckling load are derived, and validated against the corresponding exact solution, that give the attenuated buckling response of an infinitely long cylinder in terms of the nondimensional parameters presented herein. A relatively small number of "master curves" are identified that give a nondimensional measure of the buckling load of an infinitely long cylinder as a function of the orthotropy and wall inhomogeneity parameters. These curves reduce greatly the complexity of the design-variable space as compared to representations that use dimensional quantities as design variables. As a result of their inherent simplicity, these master curves are anticipated to be useful in the ongoing development of buckling-design technology.

Nemeth, Michael P.↗
NASA NTRS2013

An Indirect System Identification Technique for Stable Estimation of Continuous-Time Parameters of the Vestibulo-Ocular Reflex (VOR)

The vestibulo-ocular reflex (VOR) is a well-known dual mode bifurcating system that consists of slow and fast modes associated with nystagmus and saccade, respectively. Estimation of continuous-time parameters of nystagmus and saccade models are known to be sensitive to estimation methodology, noise and sampling rate. The stable and accurate estimation of these parameters are critical for accurate disease modelling, clinical diagnosis, robotic control strategies, mission planning for space exploration and pilot safety, etc. This paper presents a novel indirect system identification method for the estimation of continuous-time parameters of VOR employing standardised least-squares with dual sampling rates in a sparse structure. This approach permits the stable and simultaneous estimation of both nystagmus and saccade data. The efficacy of this approach is demonstrated via simulation of a continuous-time model of VOR with typical parameters found in clinical studies and in the presence of output additive noise.

Kukreja, Sunil L.↗
NASA NTRS2014

Crack Growth Behavior in the Threshold Region for High Cycle Fatigue Loading

The research results described in this paper presents a new understanding of the behavior of fatigue crack growth in the threshold region. It is believed by some crack growth experts that the ASTM load shedding test method does not produce true or valid threshold properties. The concern involves the observed fanning of threshold region da/dN data plots for some materials in which the low R-ratio data fans out or away from the high R-ratio data. This data fanning or elevation of threshold values is obviously caused by an increase in crack closure in the low R-ratio tested specimens. This increase in crack closure is assumed by some investigators to be caused by a plastic wake on the crack surfaces that was created during the load shedding test phase. This study shows that the increase in crack closure is the result of an extensive occurrence of crack bifurcation behavior in some materials, particularly in aluminum alloys, when the crack tip cyclic yield zone size becomes less than the grain size of the alloy. This behavior is related to the high stacking fault energy (SFE) property of aluminum alloys which results in easier slip characteristics. Therefore, the particular fanning behavior in aluminum alloys is a function of intrinsic dislocation property of the materials and that the fanned data represents valid material properties. However, for corrosion sensitive steel alloys used in this study the fanning was caused by a build-up of iron oxide at the crack tip from fretting corrosion.

Forman, Royce G.↗
NASA NTRS2014

Crack Growth Behavior in the Threshold Region for High Cycle Fatigue Loading

This paper describes the results of a research program conducted to improve the understanding of fatigue crack growth rate behavior in the threshold growth rate region and to answer a question on the validity of threshold region test data. The validity question relates to the view held by some experimentalists that using the ASTM load shedding test method does not produce valid threshold test results and material properties. The question involves the fanning behavior observed in threshold region of da/dN plots for some materials in which the low R-ratio data fans out from the high R-ratio data. This fanning behavior or elevation of threshold values in the low R-ratio tests is generally assumed to be caused by an increase in crack closure in the low R-ratio tests. Also, the increase in crack closure is assumed by some experimentalists to result from using the ASTM load shedding test procedure. The belief is that this procedure induces load history effects which cause remote closure from plasticity and/or roughness changes in the surface morphology. However, experimental studies performed by the authors have shown that the increase in crack closure is a result of extensive crack tip bifurcations that can occur in some materials, particularly in aluminum alloys, when the crack tip cyclic yield zone size becomes less than the grain size of the alloy. This behavior is related to the high stacking fault energy (SFE) property of aluminum alloys which results in easier slip characteristics. Therefore, the fanning behavior which occurs in aluminum alloys is a function of intrinsic dislocation property of the alloy, and therefore, the fanned data does represent the true threshold properties of the material. However, for the corrosion sensitive steel alloys tested in laboratory air, the occurrence of fanning results from fretting corrosion at the crack tips, and these results should not be considered to be representative of valid threshold properties because the fanning is eliminated when testing is performed in dry air.

Forman, R. G.↗
NASA NTRS2014

Structural Performance of Advanced Composite Tow-Steered Shells With Cutouts

The structural performance of two advanced composite tow-steered shells with cutouts, manufactured using an automated fiber placement system, is assessed using both experimental and analytical methods. The shells' fiber orientation angles vary continuously around their circumference from +/-10 degrees on the crown and keel, to +/-45 degrees on the sides. The raised surface features on one shell result from application of all 24 tows during each fiber placement system pass, while the second shell uses the system's tow drop/add capability to achieve a more uniform wall thickness. These unstiffened shells were previously tested in axial compression and buckled elastically. A single cutout, scaled to represent a passenger door on a commercial aircraft, is then machined into one side of each shell. The prebuckling axial stiffnesses and bifurcation buckling loads of the shells with cutouts are also computed using linear finite element structural analyses for initial comparisons with test data. When retested, large deflections were observed around the cutouts, but the shells carried an average of 92 percent of the axial stiffness, and 86 percent of the buckling loads, of the shells without cutouts. These relatively small reductions in performance demonstrate the potential for using tow steering to mitigate the adverse effects of typical design features on the overall structural performance.

Wu, K. Chauncey↗
NASA NTRS2014

Buckling Analysis for Stiffened Anisotropic Circular Cylinders Based on Sanders Nonlinear Shell Theory

Nonlinear and bifurcation buckling equations for elastic, stiffened, geometrically perfect, right-circular cylindrical, anisotropic shells subjected to combined loads are presented that are based on Sanders' shell theory. Based on these equations, a three-parameter approximate Rayleigh-Ritz solution and a classical solution to the buckling problem are presented for cylinders with simply supported edges. Extensive comparisons of results obtained from these solutions with published results are also presented for a wide range of cylinder constructions. These comparisons include laminated-composite cylinders with a wide variety of shell-wall orthotropies and anisotropies. Numerous results are also given that show the discrepancies between the results obtained by using Donnell's equations and variants of Sanders' equations. For some cases, nondimensional parameters are identified and "master" curves are presented that facilitate the concise representation of results.

Nemeth, Michael P.↗
NASA NTRS2014

Computational Simulations of Convergent Nozzles for the AIAA 1st Propulsion Aerodynamics Workshop

Computational Fluid Dynamics (CFD) simulations were completed for a series of convergent nozzles in participation of the American Institute of Aeronautics and Astronautics (AIAA) 1st Propulsion Aerodynamics Workshop. The simulations were performed using the Wind-US flow solver. Discharge and thrust coefficients were computed for four axisymmetric nozzles with nozzle pressure ratios (NPR) ranging from 1.4 to 7.0. The computed discharge coefficients showed excellent agreement with available experimental data; the computed thrust coefficients captured trends observed in the experimental data, but over-predicted the thrust coefficient by 0.25 to 1.0 percent. Sonic lines were computed for cases with NPR >= 2.0 and agreed well with experimental data for NPR >= 2.5. Simulations were also performed for a 25 deg. conic nozzle bifurcated by a flat plate at NPR = 4.0. The jet plume shock structure was compared with and without the splitter plate to the experimental data. The Wind-US simulations predicted the shock structure well, though lack of grid resolution in the plume reduced the sharpness of the shock waves. Unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations and Detached Eddy Simulations (DES) were performed at NPR = 1.6 for the 25 deg conic nozzle with splitter plate. The simulations predicted vortex shedding from the trailing edge of the splitter plate. However, the vortices of URANS and DES solutions appeared to dissipate earlier than observed experimentally. It is believed that a lack of grid resolution in the region of the vortex shedding may have caused the vortices to break down too soon

DES↗
NASA NTRS2014

Wind-US Results for the AIAA 2nd Propulsion Aerodynamics Workshop

This presentation contains Wind-US results presented at the 2nd Propulsion Aerodynamics Workshop. The workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion Systems Integration Technical Committee with the purpose of assessing the accuracy of computational fluid dynamics for air breathing propulsion applications. Attendees included representatives from government, industry, academia, and commercial software companies. Participants were encouraged to explore and discuss all aspects of the simulation process including the effects of mesh type and refinement, solver numerical schemes, and turbulence modeling. The first set of challenge cases involved computing the thrust and discharge coefficients for a 25deg conical nozzle for a range of nozzle pressure ratios between 1.4 and 7.0. Participants were also asked to simulate two cases in which the 25deg conical nozzle was bifurcated by a solid plate, resulting in vortex shedding (NPR=1.6) and shifted plume shock (NPR=4.0). A second set of nozzle cases involved computing the discharge and thrust coefficients for a convergent dual stream nozzle for a range of subsonic nozzle pressure ratios. The workshop committee also compared the plume mixing of these cases across various codes and models. The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 and an offset of 1.34 times the inlet diameter. Boundary layer profiles, wall static pressure, and total pressure at downstream rake locations were examined.

Wind-US↗
NASA NTRS2014

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗
NASA NTRS2015

Structural Characterization of Advanced Composite Tow-Steered Shells with Large Cutouts

The structural performance of two advanced composite tow-steered shells with large cutouts, manufactured using an automated fiber placement system, is assessed using both experimental and analytical methods. The fiber orientation angles of the shells vary continuously around their circumference from +/- 10 degrees on the crown and keel, to +/- 45 degrees on the sides. The raised surface features on one shell result from application of all 24 tows during each fiber placement system pass, while the second shell uses the tow drop/add capability of the system to achieve a more uniform wall thickness. These unstiffened shells, both without and with small cutouts, were previously tested in axial compression and buckled elastically. In this study, a single unreinforced cutout, scaled to represent a cargo door on a commercial aircraft, is machined into one side of each shell. The prebuckling axial stiffnesses and bifurcation buckling loads of these shells with large cutouts are also computed using linear finite element structural analyses for preliminary comparisons with test data. During testing, large displacements are observed around the large cutouts, but the shells maintain an average of 91 percent of the axial stiffness, and also carry 85 percent of the buckling loads, when compared to the pristine shells without cutouts. These relatively small reductions indicate that there is great potential for using tow steering to mitigate the adverse effects of large cutouts on the overall structural performance.

Wu, K. Chauncey↗
NASA NTRS2015

Kodiak: An Implementation Framework for Branch and Bound Algorithms

Recursive branch and bound algorithms are often used to refine and isolate solutions to several classes of global optimization problems. A rigorous computation framework for the solution of systems of equations and inequalities involving nonlinear real arithmetic over hyper-rectangular variable and parameter domains is presented. It is derived from a generic branch and bound algorithm that has been formally verified, and utilizes self-validating enclosure methods, namely interval arithmetic and, for polynomials and rational functions, Bernstein expansion. Since bounds computed by these enclosure methods are sound, this approach may be used reliably in software verification tools. Advantage is taken of the partial derivatives of the constraint functions involved in the system, firstly to reduce the branching factor by the use of bisection heuristics and secondly to permit the computation of bifurcation sets for systems of ordinary differential equations. The associated software development, Kodiak, is presented, along with examples of three different branch and bound problem types it implements.

Smith, Andrew P.↗