Search NASA⌕ Search

SEARCH · Search NASA

Results for “control flow divergence”

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 109 records · Page 6

Preliminary Analysis of the Effect of Flow Separation Due to Rocket Jet Pluming on Aircraft Dynamic Stability During Atmospheric Exit

A theoretical investigation was conducted to determine the effects of body boundary-layer separation resulting from a highly underexpanded jet on the dynamic stability of a typical rocket aircraft during an atmospheric exit trajectory. The particular flight condition studied on a digital computer for five degrees of freedom was at Mach 6.0 and 150,000 feet. In view of the unknown character of the separated flow field, two estimates of the pressures in the separated region were made to calculate the unbalanced forces and moments. These estimates, based on limited fundamental zero-angle-of-attack studies and observations, are believed to cover what may be the actual case. In addition to a fixed control case, two simulated pilot control inputs were studied: rate-limited and instantaneous responses. The resulting-motions with and without boundary-layer separation were compared for various initial conditions. The lower of the assumed misalinement forces and moments led to a situation whereby a slowly damped motion could be satisfactorily controlled with rate-limited control input. The higher assumption led to larger amplitude, divergent motions when the same control rates were used. These motions were damped only when the instantaneous control responses were assumed.

Dryer, Murray↗

Two-mirrored galvanometer laser light sheet generator

Light sheets generated with either laser or noncoherent sources have found widespread application to flow visualization. Previous light sheet generating systems were usually dedicated to a specific viewing geometry. The technique with the most flexibility is the galvanometer mirror based laser light sheet system. A two-mirrored system was designed and developed to provide flexibility and adaptability to a wide range of applications. The design includes the capability to control the size and location of the laser light sheet in real time, to generate horizontal or vertical sheets, to sweep the sheet repeatedly through a volume, to generate multiple sheets with controllable separation and to rotate single or multiple laser light sheets. The system is capable of producing up to 12 sheets of laser light at an angular divergence of + or - 20 degrees. Maximum scan rate of any one line is 500 Hertz. This system has proven to be uniquely versatile and a patent has been applied for.

Leighty, B. D.↗

Highly-scalable GPU-accelerated compressible reacting flow solver for modeling high-speed flows

Emerging supercomputing systems utilize a combination of central processing units (CPUs) and graphics processing units (GPUs) in an effort to reach exascale capabilities while minimizing the energy footprint of operating such systems. Such heterogeneous machines introduce new challenges for fluids solvers because the hardware architecture and operation of a GPU are fundamentally different from conventional CPUs. In this work, a general approach for efficient implementation of finite-volume based reacting flow solvers on such heterogeneous systems is presented. Three main challenges, namely, data access pattern, thread divergence, and thread safety, are addressed. Since compressible reacting flows require special methods to deal with chemical reactions, hyperbolic and nonlinear convection terms, and the presence of turbulence, specific algorithms that ensure GPU-based efficiency are developed. The approach is demonstrated on the widely available OpenFOAM open source software by modifying core algorithms for GPU accessibility. The scalability of the resulting solver, is demonstrated using practical test cases, including flow through a scramjet engine and the dynamics of a rotating detonation engine. Here, the solver provides near-ideal scaleup on a large number of GPUs (>3000), and extremely efficient use of the GPUs, with throughput nearly a constant even when processing a large number of control volumes.

42 ENGINEERING↗

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

Shadowgraph Study of Gradient Driven Fluctuations

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Cannell, David↗

Internal Performance of a Fixed-Shroud Nonaxisymmetric Nozzle Equipped with an Aft-Hood Exhaust Deflector

An investigation was conducted in the model preparation area of the Langley 16-Foot Transonic Tunnel to determine the internal performance of a fixed-shroud nonaxisymmetric nozzle equipped with an aft-hood exhaust deflector. Model geometric parameters investigated included nozzle power setting, aft-hood deflector angle, throat area control with the aft-hood deflector deployed, and yaw vector angle. Results indicate that cruise configurations produced peak performance in the range consistent with previous investigations of nonaxisymmetric convergent-divergent nozzles. The aft-hood deflector produced resultant pitch vector angles that were always less than the geometric aft-hood deflector angle when the nozzle throat was positioned upstream of the deflector exit. Significant losses in resultant thrust ratio occurred when the aft-hood deflector was deployed with an upstream throat location. At each aft-hood deflector angle, repositioning the throat to the deflector exit improved pitch vectoring performance and, in some cases, substantially improved resultant thrust ratio performance. Transferring the throat to the deflector exit allowed the flow to be turned upstream of the throat at subsonic Mach numbers, thereby eliminating losses associated with turning supersonic flow. Internal throat panel deflections were largely unsuccessful in generating yaw vectoring.

Asbury, Scott C.↗

NASA Tech Briefs, October 2009

Topics covered include: Light-Driven Polymeric Bimorph Actuators; Guaranteeing Failsafe Operation of Extended-Scene Shack-Hartmann Wavefront Sensor Algorithm; Cloud Water Content Sensor for Sounding Balloons and Small UAVs; Pixelized Device Control Actuators for Large Adaptive Optics; T-Slide Linear Actuators; G4FET Implementations of Some Logic Circuits; Electrically Variable or Programmable Nonvolatile Capacitors; System for Automated Calibration of Vector Modulators; Complementary Paired G4FETs as Voltage-Controlled NDR Device; Three MMIC Amplifiers for the 120-to-200 GHz Frequency Band; Low-Noise MMIC Amplifiers for 120 to 180 GHz; Using Ozone To Clean and Passivate Oxygen-Handling Hardware; Metal Standards for Waveguide Characterization of Materials; Two-Piece Screens for Decontaminating Granular Material; Mercuric Iodide Anticoincidence Shield for Gamma-Ray Spectrometer; Improved Method of Design for Folding Inflatable Shells; Ultra-Large Solar Sail; Cooperative Three-Robot System for Traversing Steep Slopes; Assemblies of Conformal Tanks; Microfluidic Pumps Containing Teflon[Trademark] AF Diaphragms; Transparent Conveyor of Dielectric Liquids or Particles; Multi-Cone Model for Estimating GPS Ionospheric Delays; High-Sensitivity GaN Microchemical Sensors; On the Divergence of the Velocity Vector in Real-Gas Flow; Progress Toward a Compact, Highly Stable Ion Clock; Instruments for Imaging from Far to Near; Reflectors Made from Membranes Stretched Between Beams; Integrated Risk and Knowledge Management Program -- IRKM-P; LDPC Codes with Minimum Distance Proportional to Block Size; Constructing LDPC Codes from Loop-Free Encoding Modules; MMICs with Radial Probe Transitions to Waveguides; Tests of Low-Noise MMIC Amplifier Module at 290 to 340 GHz; and Extending Newtonian Dynamics to Include Stochastic Processes.

Source record↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

Jet Spreading Increase by Passive Control and Associated Performance Penalty

This paper reviews the effects of 'screech', 'asymmetric nozzle shaping', 'tabs' and 'overexpansion' on the spreading of free jets. Corresponding thrust penalty for the tabs and overexpanded condition are also evaluated. The asymmetric shapes include rectangular ones with varying aspect ratio. Tabs investigated are triangular shaped 'delta-tabs' placed at the exit of a convergent circular nozzle. The effect of overexpansion is examined with circular convergent-divergent (C-D) nozzles. Tabs and overexpansion are found to yield the largest increase in jet spreading. Each, however, involves a performance penalty, i.e., a loss in thrust coefficient. Variation of the size of four delta-tabs show that there exists an optimum size for which the gain in jet spreading is the maximum per unit loss in thrust coefficient. With the C-D nozzles, the minimum in thrust coefficient is expected near the beginning of the overexpanded regime based on idealized flow calculations. The maximum increase in jet spreading, however, is found to occur at higher pressure ratios well into the overexpanded regime. The optimum benefit with the overexpanded flow, in terms of gain in spreading for unit penalty, is found to be comparable to the optimum tab case.

Zaman, K. B. M. Q.↗

Artificial Intelligence-Aided Wind Plant Optimization for Nationwide Evaluation of Land Use and Economic Benefits of Wake Steering

If clean energy pathways are to harness massive increases in wind power, innovations with broad geographic viability will be needed to support buildout in diverse locations. However, geodiversity in impact potential is seldom captured in technology assessment. Here we propose a scalable approach to plant-level optimization using artificial intelligence to evaluate land sparing and economic benefits of wake steering at more than 6,800 plausible onshore wind locations in the USA. This emerging controls strategy optimizes plant energy production by directing turbine wakes. On the basis of estimates from our artificial intelligence model trained on engineering wind flow simulations, co-optimizing plant layouts with wake steering can reduce land requirements by an average of 18% per plant (site-specific benefits range from 2% to 34%), subject to errors and uncertainties in the flow model, wind resource estimates, buildout scenario and geographic factors. According to model estimates, wake steering is predicted to increase power production during high-value (relatively low wind) periods, boosting the annual revenue of individual plants by up to US$3.7 million (equivalent to US$13,000 MW-1 yr-1) but producing negligible gains in some settings. Consideration of wake steering’s geographic potential reveals divergent nationwide prospects for improved economics and siting flexibility.

deployment↗

Impact of gas injection location and divertor surface material on ITER fusion power operation phase divertor performance assessed with SOLPS-ITER *

Abstract The ITER divertor design and performance assessment, primarily based on the SOLPS-4.3 burning plasma database (Pitts R. et al 2019 Nucl. Mater. Energy 20 100696), assumes the use of beryllium (Be) as the divertor surface material and the injection of gas from the main chamber top. However, the current ITER baseline favors gas injection from the more toroidally symmetric sub-divertor region. This paper evaluates the implications of these assumptions for divertor performance in the ITER fusion power operation phase. The impact of the divertor surface material and the gas injection location on the main ions mirrors the hydrogen only low power phase scenario shown in Park J.-S. et al (2020 Nucl. Fusion 61 016021). However, during burning plasma operation, extrinsic impurity seeding will be required. In the case of neon (Ne), studied here, impurity retention is influenced by both the divertor surface material and the fueling location. Neon leakage increases due to more energetic reflection from tungsten than beryllium, but equivalent divertor performance can be achieved by adjusting the neon seeding rate. While the impurity seeding location does not affect the distributions of impurity or radiation, the fueling location does. Top fueling provides local ionization sources mainly in the mid-SOL under detached conditions, enhancing divergences of the flux there (source-driven flow), bringing stagnation points close to the fueling location, and equilibrating flows towards both targets. In contrast, the global flow pattern (in the absence of fluid drifts) in the case of sub-divertor fueling is biased towards the inner target. Impurity flows, driven by force balance, largely mirror those of the main ion flow, including the stagnation point. The case with top fueling enhances Ne retention and corresponding radiation in the outer divertor, effectively reducing the total and peak target heat fluxes by 20%–40%, compared to the case with divertor fueling. Meanwhile, the case with outer target fueling also achieves similar reductions by enhancing plasma-neutral interactions. These results suggest the possibility that the selection of the fueling location and throughput can be used as an actuator to control impurity divertor retention and divertor radiation asymmetry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Experimental, Theoretical, and Computational Investigation of Separated Nozzle Flows

A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.

Hunter, Craig A.↗

Extension of Low Dissipative High Order Hydrodynamics Schemes for MHD Equations

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD (magnetohydrodynamic) equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls to limit the amount and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvi-linear grids. The three features of the present MHD scheme over existing schemes in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion magnetized flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike existing schemes for the conservative MHD equations which suffer from ill-conditioned eigen-decompositions, the present scheme makes use of a well-conditioned eigen-decomposition to solve the conservative form of the MHD equations. This is due to, partly. the fact that the divergence of the magnetic field condition is a different type of constraint from its incompressible Navier-Stokes cousin. Third, a new approach to minimize the numerical error of the divergence free magnetic condition for high order scheme is introduced.

Yee, H. C.↗

Development of a Common Research Model for Applied CFD Validation Studies

The development of a wing/body/nacelle/pylon/horizontal-tail configuration for a common research model is presented, with focus on the aerodynamic design of the wing. Here, a contemporary transonic supercritical wing design is developed with aerodynamic characteristics that are well behaved and of high performance for configurations with and without the nacelle/pylon group. The horizontal tail is robustly designed for dive Mach number conditions and is suitably sized for typical stability and control requirements. The fuselage is representative of a wide/body commercial transport aircraft; it includes a wing-body fairing, as well as a scrubbing seal for the horizontal tail. The nacelle is a single-cowl, high by-pass-ratio, flow-through design with an exit area sized to achieve a natural unforced mass-flow-ratio typical of commercial aircraft engines at cruise. The simplicity of this un-bifurcated nacelle geometry will facilitate grid generation efforts of subsequent CFD validation exercises. Detailed aerodynamic performance data has been generated for this model; however, this information is presented in such a manner as to not bias CFD predictions planned for the fourth AIAA CFD Drag Prediction Workshop, which incorporates this common research model into its blind test cases. The CFD results presented include wing pressure distributions with and without the nacelle/pylon, ML/D trend lines, and drag-divergence curves; the design point for the wing/body configuration is within 1% of its max-ML/D. Plans to test the common research model in the National Transonic Facility and the Ames 11-ft wind tunnels are also discussed.

Vassberg, John C.↗

An Estimation of the Flying Qualities of the Kaiser Fleetwing All-Wing Airplane from Tests of a 1/7-Scale Model, TED No. NACA 2340

An investigation of a 1/7-scale powered model of the Kaiser Fleetwing all-wing airplane was made in the Langley full-scale tunnel to provide data for an estimation of the flying qualities of the airplane. The analysis of the stability and control characteristics of the airplane has been made as closely as possible in accordance with the requirements of the Bureau of Aeronautics, Navy Department's specifications, and a summary of the more significant conclusions is presented as follows. With the normal center of gravity located at 20 percent of the mean aerodynamic chord, the airplane will have adequate static longitudinal stability, elevator fixed, for all flight conditions except for low-power operation at low speeds where the stability will be about neutral. There will not be sufficient down-elevator deflection available for trim above speeds of about 130 miles per hour. It is probable that the reduction in the up-elevator deflections required for trim will be accompanied by reduced elevator hinge moments for low-power operation at low flight speeds. The static directional stability for this airplane will be low for all rudder-fixed or rudder-free flight conditions. The maximum rudder deflection of 30 deg will trim only about 15 deg yaw for most flight conditions and only 10 deg yaw for the condition with low power at low speeds. Also, at low powers and low speeds, it is estimated that the rudders will not trim the total adverse yaw resulting from an abrupt aileron roll using maximum aileron deflection. The airplane will meet the requirements for stability and control for asymmetric power operation with one outboard engine inoperative. The airplane would have no tendency for directional divergence but would probably be spirally unstable, with rudders fixed. The static lateral stability of the airplane will probably be about neutral for the high-speed flight conditions and will be only slightly increased for the low-power operation in low-speed flight. The airplane will not roll against the ailerons in a side-slip maneuver. Although the airplane would probably meet the minimum requirements of pb/2V of 0.07 at all speeds, there will be a loss in rolling ability of the airplane at high aileron deflections and at low flight speeds. It is estimated that the wing stall will be a gradual movement forward from the trailing edge and will be accompanied by no sudden pitching or rolling accelerations. Some stall warning may be indicated by reduction in the elevator and aileron force gradients and by the shaking of the controls caused by unsteady flow over the surfaces near the stall.

Brewer, Gerald W.↗

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗

Nozzle Aerodynamic Stability During a Throat Shift

An experimental investigation was conducted on the internal aerodynamic stability of a family of two-dimensional (2-D) High Speed Civil Transport (HSCT) nozzle concepts. These nozzles function during takeoff as mixer-ejectors to meet acoustic requirements, and then convert to conventional high-performance convergent-divergent (CD) nozzles at cruise. The transition between takeoff mode and cruise mode results in the aerodynamic throat and the minimum cross-sectional area that controls the engine backpressure shifting location within the nozzle. The stability and steadiness of the nozzle aerodynamics during this so called throat shift process can directly affect the engine aerodynamic stability, and the mechanical design of the nozzle. The objective of the study was to determine if pressure spikes or other perturbations occurred during the throat shift process and, if so, identify the caused mechanisms for the perturbations. The two nozzle concepts modeled in the test program were the fixed chute (FC) and downstream mixer (DSM). These 2-D nozzles differ principally in that the FC has a large over-area between the forward throat and aft throat locations, while the DSM has an over-area of only about 10 percent. The conclusions were that engine mass flow and backpressure can be held constant simultaneously during nozzle throat shifts on this class of nozzles, and mode shifts can be accomplished at a constant mass flow and engine backpressure without upstream pressure perturbations.

Kawecki, Edwin J.↗