Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution Coefficient”

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 145 records · Page 8

Developing an Accurate CFD Based Gust Model for the Truss Braced Wing Aircraft

The increased flexibility of long endurance aircraft having high aspect ratio wings necessitates attention to gust response and perhaps the incorporation of gust load alleviation. The design of civil transport aircraft with a strut or truss-braced high aspect ratio wing furthermore requires gust response analysis in the transonic cruise range. This requirement motivates the use of high fidelity nonlinear computational fluid dynamics (CFD) for gust response analysis. This paper presents the development of a CFD based gust model for the truss braced wing aircraft. A sharp-edged gust provides the gust system identification. The result of the system identification is several thousand time steps of instantaneous pressure coefficients over the entire vehicle. This data is filtered and downsampled to provide the snapshot data set from which a reduced order model is developed. A stochastic singular value decomposition algorithm is used to obtain a proper orthogonal decomposition (POD). The POD model is combined with a convolution integral to predict the time varying pressure coefficient distribution due to a novel gust profile. Finally the unsteady surface pressure response of the truss braced wing vehicle to a one-minus-cosine gust, simulated using the reduced order model, is compared with the full CFD.

Bartels, Robert E.↗

A Non-Cut Cell Immersed Boundary Method for Use in Icing Simulations

This paper describes a computational fluid dynamic method used for modelling changes in aircraft geometry due to icing. While an aircraft undergoes icing, the accumulated ice results in a geometric alteration of the aerodynamic surfaces. In computational simulations for icing, it is necessary that the corresponding geometric change is taken into consideration. The method used, herein, for the representation of the geometric change due to icing is a non-cut cell Immersed Boundary Method (IBM). Computational cells that are in a body fitted grid of a clean aerodynamic geometry that are inside a predicted ice formation are identified. An IBM is then used to change these cells from being active computational cells to having properties of viscous solid bodies. This method has been implemented in the NASA developed node centered, finite volume computational fluid dynamics code, FUN3D. The presented capability is tested for two-dimensional airfoils including a clean airfoil, an iced airfoil, and an airfoil in harmonic pitching motion about its quarter chord. For these simulations velocity contours, pressure distributions, coefficients of lift, coefficients of drag, and coefficients of pitching moment about the airfoil's quarter chord are computed and used for comparison against experimental results, a higher order panel method code with viscous effects, XFOIL, and the results from FUN3D's original solution process. The results of the IBM simulations show that the accuracy of the IBM compares satisfactorily with the experimental results, XFOIL results, and the results from FUN3D's original solution process.

Sarofeen, Christian M.↗

Assessment of Using Ideal Gas for Predicting Boattail Flow at Cryogenic Temperatures

The applicability of using ideal gas assumptions to simulate high Reynolds number experimental data that was obtained at cryogenic temperatures is examined. Flow over an axisymmetric nozzle boattail model was calculated using reference temperatures of 117 K and 300 K and at Reynolds numbers from 50 to 200 million per meter. From the testing perspective, pressure, compression factor, and isentropic coefficients calculated using one-dimensional real gas equations are used to examine the departure of cryogenic flow from ideal gas flow across the range of temperatures and potential impacts on measured aerodynamic data. Solutions developed using ideal gas assumptions in a three-dimensional Navier-Stokes code are compared with experimental data obtained at cryogenic temperatures at two unit Reynolds numbers at freestream Mach numbers of 0.6 and 0.9. Results for several one- and two-equation turbulence models are shown. Predicted pressure coefficient distributions along the nozzle boattail differed from experimental data between 8% to less than 0.5% depending on the turbulence model and Mach number. The greatest discrepancy occurred in the level of static pressure recovery in the recompression region where the flow was separated. Solutions using warm and cryogenic freestream temperatures predicted similar boattail pressure distributions at the same unit Reynolds number.

Nozzle↗

Assessment of Using Ideal Gas for Predicting Boattail Flow at Cryogenic Temperatures

The applicability of using ideal gas assumptions to simulate high Reynolds number experimental data that was obtained at cryogenic temperatures is examined. Flow over an axisymmetric nozzle boattail model was calculated using reference temperatures of 117 K and 300 K and at Reynolds numbers from 50 to 200 million per meter. From the testing perspective, pressure, compression factor, and isentropic coefficients calculated using one-dimensional real gas equations are used to examine the departure of cryogenic flow from ideal gas flow across the range of temperatures and potential impacts on measured aerodynamic data. Solutions developed using ideal gas assumptions in a three-dimensional Navier-Stokes code are compared with experimental data obtained at cryogenic temperatures at two unit Reynolds numbers at freestream Mach numbers of 0.6 and 0.9. Results for several one- and two-equation turbulence models are shown. Predicted pressure coefficient distributions along the nozzle boattail differed from experimental data between 8% to less than 0.5% depending on the turbulence model and Mach number. The greatest discrepancy occurred in the level of static pressure recovery in the recompression region where the flow was separated. Solutions using warm and cryogenic freestream temperatures predicted similar boattail pressure distributions at the same unit Reynolds number.

Nozzle↗

New Insights on the 7[superscript]Be Cycle in the Ocean

The cosmogenic radionuclide 7[superscript]Be has been applied as a tracer of dynamical processes in the upper ocean and of atmospheric deposition of trace elements at the sea surface. These applications usually assume that 7[superscript]Be is entirely in the dissolved form, and that scavenging and downward export of 7[superscript]Be by settling particles can be neglected. In this work, we explore these assumptions and more generally assess the significance of the 7[superscript]Be activity in the particulate fraction, through the generation of vertical profiles of particulate 7[superscript]Be in the open ocean. From detailed measurements obtained from low-background gamma spectrometers placed in underground facilities, we report vertical profiles of 7[superscript]Be activity in suspended particles (7[superscript]Bep[subscript]) collected in various oceanic regions: the Mediterranean Sea (DYFAMED station), the Indian Sector of the Southern Ocean (station A3-2 from the KEOPS2 cruise), the Sargasso Sea (OFP station), and the subpolar North Atlantic Ocean (GEOVIDE cruise). We find that, in each oceanic region, 7[superscript]Bep[subscript] activities are generally higher in the mixed layer than in the thermocline. They vary in the mixed layer from 3.0 dpm/m3 at DYFAMED to 33.1dpm/m3 at GEOVIDE cross-over station 51/60, i.e., within a range consistent with previous 7[superscript]Bep[subscript] measurements for the open ocean. For the GEOVIDE cruise, the 7[superscript]Bep[subscript] activities measured on different filter types at different depths are corrected for filter offsets derived from multiple 7[superscript]Bep[subscript] measurements at a near-coastal station in the western Mediterranean Sea. We then combine measurements of total 7[superscript]Be activity (Shelley et al., 2017) with our measurements of 7[superscript]Bep[subscript] activity to estimate the solid/solution partitioning of 7[superscript]Be. On average, the particulate fraction would represent 6% of total 7[superscript]Be activity at 5-m water depth (n = 6), 22% at 20 m (n = 2) and 9% at 70 m (n = 3). At GEOVIDE stations, 7[superscript]Bep[subscript] inventories range from 5% to 19% of the total 7[superscript]Be inventories. In the Labrador Sea, the measured 7[superscript]Bep[subscript] inventories are lower than the dry 7[superscript]Be deposition fluxes estimated from aerosol samples collected during GEOVIDE, suggesting that a significant portion of 7[superscript]Bep[subscript] may be removed by sinking particles. The distribution coefficient Kd for 7Be increases with depth, with log10Kd averaging 5.2 ± 0.1 at 5 m to 6.1 ± 0.1 between 70 and 150 m, suggesting that scavenging of dissolved 7[superscript]Be by particles is more pronounced in the thermocline than in the mixed layer when differences in particle concentrations are taken into account. Overall, our study suggests that, at least in some oceanic regions, the removal of 7[superscript]Be by marine particles may be significant and that it may need to be considered in applications of 7[superscript]Be as a tracer of oceanic processes and atmospheric deposition.

Beryllium-7↗

Transition Modeling Based on the Dual N-factor Method for the CRM-NLF Wind Tunnel Configuration

The dual N-factor method is used to model the boundary-layer transition over the common research model with natural laminar flow (CRM-NLF) aircraft configuration. The flow conditions match selected test conditions from a wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. The paper presents a systematic methodology for transition prediction in the presence of a dual shock system and extends the prior capability for iteratively coupled computational fluid dynamics (CFD) predictions to incorporate three-dimensional, transonic wings. The method employs stability computations based on the linear parabolized stability equations (PSE), along with a dual N-factor criterion. The iterative process begins with the fully turbulent Reynolds-averaged-Navier-Stokes (RANS) mean flow solution. For the first iteration, a mean flow solution is calculated with an imposed transition front that aligns with the shock front from the fully turbulent solution. Subsequently, stability computations are performed along a set of streamlines across the wing to calculate the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) modes. The transition criterion based on the dual N-factor method is used to infer the updated transition front and the process is successively repeated until convergence of the solution. Within three iterations, the predicted fronts for angles of attack of 1.45, 1.98, 2.46 and 2.94 degrees and a mean-aerodynamic-chord Reynolds number equal to 15 million, approach visual convergence in most regions of the studied cases, and the resulting predictions are in good agreement with the transition fronts deduced from measurements of temperature-sensitive paint. Even though surface pressure measurements based on fully-turbulent flow agree well with the measured pressure coefficient distributions, strong viscous-inviscid interaction effects cause significant shifts in the shock locations based on the imposed transition front, underscoring the intrusive nature of static pressure measurements using surface mounted ports on the CRM-NLF configuration.

Boundary Layer Transition↗

Pressure distribution and aerodynamic coefficients associated with heat addition to supersonic air stream adjacent to two-dimensional supersonic wing

The modifications in the pressure distributions and the aerodynamic coefficients associated with additions of heat to the two-dimensional supersonic in viscid flow field adjacetnt to the lower surface of of a 5-percent-thickness symmetrical circular-arc wing are presented in this report. The pressure distributions are obtained by the use of graphical method which gives the two-dimensional supersonic inviscid flow field obtained with moderate heat addition. The variation is given of the lift-drag ratio and of the aerodynamic coefficients of lift, drag, and moment with free stream Mach number, angle of attack, and parameters defining extent and amount of heat addition. The six graphical solutions used in this study included Mach numbers of 3.0 and 5.0 and angles of attack of 0 degrees and 2 degrees.

Pinkel, I Irving↗

Measurements of Pressure Distributions and Force Coefficients in a Squeeze Film Damper. Part 2: Partially Sealed Configuration

Experimental results from a partially sealed squeeze film damper (SFD) test rig, executing a circular centered orbit are presented and discussed. A serrated piston ring is installed at the damper exit. This device involves a new sealing concept which produces high damping values while allowing for oil flow to cool the damper. In the partially sealed damper, large cavitation regions are observed in the pressure fields at orbit radii epsilon equals 0.5 and epsilon equals 0.8. The cavitated pressure distributions and the corresponding force coefficients are compared with a cavitated bearing solution. The experimental results show the significance of fluid inertia and vapor cavitation in the operation of squeeze film dampers. Squeeze film Reynolds numbers tested reach up to Re equals 50, spanning the range of contemporary applications.

Jung, S. Y.↗

Chemical fractionations in group IIIAB iron meteorites - Origin by dendritic crystallization of an asteroidal core

Results of an investigation of the crystallization history of the asteroidal core that produced nearly 200 iron meteorites of group IIIAB are presented. The accuracy of the published distribution of coefficients between solid and liquid metal for seven elements and their dependence on the concentrations of S and P are critically assessed to develop a nonideal fractional crystallization model that reproduces all of the major features of the overall chemical trends in group IIIAB. It is suggested that the apparent distribution coefficient of S was was much higher than its equilibrium value (not greater than 0.01) and that it increased during crystallization. In the present preferred model the apparent distribution coefficient for S increases from 0.6 to 0.8, while the S concentration of the liquid increases from 6 to 13 wt pct. It is inferred that light S-rich liquid accumulates preferentially at the top of the core and in structural traps formed by the advancing solid.

Haack, Henning↗

Flow distributions and discharge coefficient effects for jet array impingement with initial crossflow

To model the impingement cooled mid-chord region of gas turbine airfoils in cases where an initial crossflow is present, the paper presents experimentally determined flow distributions for jet arrays with ten spanwise rows of holes, which range from uniform to highly nonuniform. The jet flow after impingement is constrained to exit in a single direction along the channel formed by the jet orifice plate and the impingement surface. The streamwise distributions and crossflow velocities are presented for ratios of the initial crossflow rate to the total jet flow rate ranging from zero to unity. For crossflow to jet velocity ratios greater than a value somewhat less than unity, jet orifice discharge coefficients do not remain constant, but decrease significantly, showing a secondary dependence on z/d, where z is the channel height and d is the jet hole diameter.

Florschuetz, L. W.↗

Higher-order skewness and excess coefficients of some probability distributions applicable to optical propagation phenomena

Expressions are derived for higher-order skewness and excess coefficients using central moments and cumulants up to 8th order. These coefficients are then calculated for three probability distributions: (1) Log-normal, (2) Rice-Nakagami, and (3) Gamma distributions. Curves are given to shown the variation of skewness with excess coefficients for these distributions. These curves are independent of the particular distribution parameters. This method is useful for studying fluctuating phenomena, which obey non-Gaussian statistics.

Majumdar, A. K.↗

Heat Transfer and Pressure Distributions on a Gas Turbine Blade Tip

Heat transfer coefficient and static pressure distributions are experimentally investigated on a gas turbine blade tip in a five-bladed stationary linear cascade. The blade is a 2-dimensional model of a first stage gas turbine rotor blade with a blade tip profile of a GE-E(sup 3) aircraft gas turbine engine rotor blade. The flow condition in the test cascade corresponds to an overall pressure ratio of 1.32 and exit Reynolds number based on axial chord of 1.1 x 10(exp 6). The middle 3-blade has a variable tip gap clearance. All measurements are made at three different tip gap clearances of about 1%, 1.5%, and 2.5% of the blade span. Heat transfer measurements are also made at two different turbulence intensity levels of 6.1 % and 9.7% at the cascade inlet. Static pressure measurements are made in the mid-span and the near-tip regions as well as on the shroud surface, opposite the blade tip surface. Detailed heat transfer coefficient distributions on the plane tip surface are measured using a transient liquid crystal technique. Results show various regions of high and low heat transfer coefficient on the tip surface. Tip clearance has a significant influence on local tip beat transfer coefficient distribution. Heat transfer coefficient also increases about 15-20% along the leakage flow path at higher turbulence intensity level of 9.7% over 6.1 %.

Azad, Gm S.↗

Satellite drag coefficients calculated from measured distributions of reflected helium atoms

The primary objectives of this study were to obtain the necessary data and develop a calculation procedure that would facilitate predicting the atmospheric-helium contribution to the drag of a satellite having a predominantly convex exterior. Molecular-beam techniques were used to measure, for several incidence angles, the spatial and energy distributions of 7000 m/s helium atoms scattered from a 6061-T6 aluminum plate and an anodized 1235-0 aluminum surface. From these measured distributions, tangential and normal momentum accommodation coefficients were calculated as functions of incidence angle. Using these calculated accommodation coefficients, one can predict drag coefficients for satellites having predominantly convex exteriors. For spherical satellites, drag coefficients of 2.64 and 2.62 were predicted for the subject surfaces.

Liu, S.-M.↗

A comparison of experimental and theoretical results for leakage, pressure gradients, and rotordynamic coefficients for tapered annular gas seal

A brief review of current annular seal theory and a discussion of the predicted effect on stiffness of tapering the seal stator are presented. An outline of Nelson's analytical-computational method for determining rotordynamic coefficients for annular compressible-flow seals is included. Modifications to increase the maximum rotor speed of an existing air-seal test apparatus at Texas A&M University are described. Experimental results, including leakage, entrance-loss coefficients, pressure distributions, and normalized rotordynamic coefficients, are presented for four convergent-tapered, smooth-rotor, smooth-stator seals. A comparison of the test results shows that an inlet-to-exit clearance ratio of 1.5 to 2.0 provides the maximum direct stiffness, a clearance ratio of 2.5 provides the greatest stability, and a clearance ratio of 1.0 provides the least stability. The experimental results are compared to theoretical results from Nelson's analysis with good agreement. Test results for cross-coupled stiffness show less sensitivity of fluid prerotation than predicted.

Elrod, D. A.↗

Application of program LAURA to three-dimensional AOTV flowfields

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is a robust, finite volume, single-level storage, implicit upwind differencing algorithm which has been documented and tested on several three-dimensional blunt-body flows. The algorithm can run at unlimited Courant numbers (relaxing the steady-state equations) but requires the inversion of only a 5 x 5 matrix per computational cell. An alternating directional sweep Gauss-Seidel substitution strategy is used to relax the governing equations. At present, the Euler and thin-layer Navier-Stokes equations using Sutherland's law for viscosity have been modeled for a perfect gas, equilibrium air, and nonequilibrium air chemistry neglecting diffusion. The equilibrium and nonequilibrium air chemistry options have been described in a companion paper. Good comparisons with experimental data and another calculation method for pressure distributions, aerodynamic coefficients, and heat-transfer distributions have been demonstrated for three-dimensional blunt-body flows.

Gnoffo, P. A.↗

An experimental investigation of control surface effectiveness and real-gas simulation for biconics

Aerodynamic coefficients, pressure distributions, and oil-flow patterns measured on spherically blunted, 13 degrees/7 degrees on-axis and bent biconics at Mach 6 and 10 in air are compared. Angle of attack, referenced to the axis of the aft cone, was varied from zero to 40 degrees. The effect of deflection of base-mounted flaps on aerodynamic characteristics was examined at Mach 10. Real-gas effects on aerodynamic coefficients, pressure distributions and shock detachment distance were simulated by testing the biconic models in Mach 6 air and Mach 6 CF4 flows. (Density ratio equal to 5.3 for air and 12 for CF4.) The on-axis biconic is stable but cannot be trimmed at angles of attack above zero with the assumed center of gravity and flap configuration. Flaps on the bent biconic produced a stable and trim configuration for angles of attack from zero to 23 degrees. A significant effect of density ratio on aerodynamic coefficients, pressure distribution, and shock detachment distance was observed.

Miller, C. G.↗

Calculated spanwise lift distributions, influence functions, and influence coefficients for unswept wings in subsonic flow

Spanwise lift distributions have been calculated for nineteen unswept wings with various aspect ratios and taper ratios and with a variety of angle-of-attack or twist distributions, including flap and aileron deflections, by means of the Weissinger method with eight control points on the semispan. Also calculated were aerodynamic influence coefficients which pertain to a certain definite set of stations along the span, and several methods are presented for calculating aerodynamic influence functions and coefficients for stations other than those stipulated. The information presented in this report can be used in the analysis of untwisted wings or wings with known twist distributions, as well as in aeroelastic calculations involving initially unknown twist distributions.

Diederich, Franklin W↗