Search NASA⌕ Search

SEARCH · Search NASA

Results for “WINGS - ASPECT RATIO, LARGE”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗
NASA NTRS1979

On the stability of the boundary layer on a transonic swept wing

Both incompressible and compressible linear stability theory are applied to the three-dimensional compressible boundary layer on a particular transonic sweptback wing of infinite span. A spatial stability theory is used which identifies the growth direction with the real part of the complex angle of the group velocity. It is found that in the forward, but not the rear, crossflow instability region, the maximum amplification rates of the steady disturbances may be calculated to within about 10% by the incompressible stability theory. There is little difference between the sixth and eighth-order compressible theories. The maximum amplification rate of the steady disturbances at any chordwise station is closely related to the maximum crossflow at that station independent of the Reynolds number. For other than crossflow instability, there can be large differences between the incompressible and compressible theories, both as to the amplification rate and the angle of the wavenumber vector for maximum instability. Amplitude ratios of individual wave components are obtained by integrating the spatial amplification rate along the growth direction subject to the constraint that the wavenumber vector is irrotational. This procedure yields steady disturbances aligned with the local potential flow direction whose wavelengths are nearly independent of downstream distance.

Mack, L. M.↗
NASA NTRSDate not supplied

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m 3 . This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion, Large drops, Wind tunnel↗
NASA NTRSDate not supplied

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m3. This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion↗
NASA NTRS2019

Bowles-Tatnall Wake Vortex Encounter Hazard Metric

Wake vortex spacing standards constrict the terminal area throughput and impose severe constraints on the overall capacity and efficiency of the National Airspace System. For more than two decades starting in the early 1990s, the National Aeronautics and Space Administration conducted extensive research on characterizing the formation and evolution of aircraft wakes. This multidisciplinary work included comprehensive field experiments (Pruis et al. 2016), flight tests (Vicroy et al. 1998), and wind tunnel tests (Rossow 1994; Chow et al. 1997). Parametric studies using large eddy simulations (Proctor 1998; Proctor et al. 2006) were conducted in order to develop fast-time models for the prediction of wake transport and decay (Ahmad et al. 2016). Substantial effort was spent on the formulation of acceptable vortex hazard metrics (Tatnall 1995; Hinton and Tatnall 1997). Several wake encounter severity metrics have been suggested in the past, which include the wake circulation strength, vortex-induced rolling moment coefficient (Clv), bank angle, and the roll control ratio (Tatnall 1995; Hinton and Tatnall 1997; Van der Geest 2012). The vortex-induced rolling moment coefficient introduced by Bowles and Tatnall (Tatnall 1995; Gloudemans et al. 2016) has been used extensively for risk and safety analysis of newly proposed air traffic management concepts and procedures. The original method of Bowles and Tatnall assumed a constant wing loading (the wing lift-curve slope, CL is constant), which resulted in an overestimation of the vortexinduced rolling moment coefficient. Bowles (2014) suggested a correction to the original method that provides more accurate values of Clv and which is also consistent with the underlying physics of the problem. The overestimation of Clv in the original method can be corrected by assuming an elliptical lift distribution. Figure 1.1 illustrates the correction in Clv achieved by the modified method.

Joel Malissa↗
NASA NTRS2011

Axisymmetric Calculations of a Low-Boom Inlet in a Supersonic Wind Tunnel

This paper describes axisymmetric CFD predictions made of a supersonic low-boom inlet with a facility diffuser, cold pipe, and mass flow plug within wind tunnel walls, and compares the CFD calculations with the experimental data. The inlet was designed for use on a small supersonic aircraft that would cruise at Mach 1.6, with a Mach number over the wing of 1.7. The inlet was tested in the 8-ft by 6-ft Supersonic Wind Tunnel at NASA Glenn Research Center in the fall of 2010 to demonstrate the performance and stability of a practical flight design that included a novel bypass duct. The inlet design is discussed here briefly. Prior to the test, CFD calculations were made to predict the performance of the inlet and its associated wind tunnel hardware, and to estimate flow areas needed to throttle the inlet. The calculations were done with the Wind-US CFD code and are described in detail. After the test, comparisons were made between computed and measured shock patterns, total pressure recoveries, and centerline pressures. The results showed that the dual-stream inlet had excellent performance, with capture ratios near one, a peak core total pressure recovery of 96 percent, and a large stable operating range. Predicted core recovery agreed well with the experiment but predicted bypass recovery and maximum capture ratio were high. Calculations of offdesign performance of the inlet along a flight profile agreed well with measurements and previous calculations.

Chima, Rodrick V.↗
DOE OSTI2024

Properties of a fading AGN from SDSS-IV MaNGA

ABSTRACT We identify a fading AGN SDSS J220141.64+115124.3 from the internal Product Launch-11 (MPL-11) in Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The central region with a projected radius of $\sim$2.4 kpc is characterized as LINER-like line ratios while the outskirts extended to $\sim$15 kpc show Seyfert-like line ratios. The ${\rm{[O {\small III}]}}$$\lambda$5007 luminosity of the Seyfert regions is a factor of 37 (2) higher than the LINER regions without (with) dust attenuation correction, suggesting that the AGN activity decreases at least $\sim$8 $\times$ 10$^3$ yr ($\sim$2.4 kpc/light-speed) ago. We model the emission line spectra in the central region with double Gaussian components (a narrow core and a broad wing) and analyse the properties of each component. The narrow core component mostly co-rotates with the stellar disc, whereas the broad wing component with a median of the velocity dispersion $\sim$300 km s$^{-1}$ is related to a wind outflow. The kinematic position angle (PA) of the ionized gas shows an $\sim 20^{\circ }$ twist from the galaxy centre to 1.5 effective radius. The median of the PA difference between the gas and stellar components is as large as $\sim 50^{\circ }$ within 0.4 effective radius. The tidal feature in DESI image and star–gas misalignment suggest this galaxy is a merger remnant. Combining all these observational results as well as public available X-ray and MIR luminosities, we confirm this is a fading AGN, the merger process kick-started the central engine to quasar phase which ionized gas composed of tidal debris, and now the activity of the central black hole decreases. The discontinuity in ${\rm{[O {\small III}]}}$$\lambda$5007 flux and EQW maps is due to multiple AGN outbursts triggered by merger remnant gas inflows.

Mo, Hao (ORCID:0000000234430768)↗
NASA NTRS1960

Flight Measurement of Wall-Pressure Fluctuations and Boundary-Layer Turbulence

The results are presented for a flight test program using a fighter type jet aircraft flying at pressure altitudes of 10,000, 20,000, and 30,000 feet at Mach numbers from 0.3 to 0.8. Specially designed apparatus was used to measure and record the output of microphones and hot-wire anemometers mounted on the forward-fuselage section and wing of the airplane. Mean-velocity profiles in the boundary layers were obtained from total-pressure measurements. The ratio of the root-mean-square fluctuating wall pressure to the free-stream dynamic pressure is presented as a function of Reynolds number and Mach number. The longitudinal component of the turbulent-velocity fluctuations was measured, and the turbulence-intensity profiles are presented for the wing and forward-fuselage section. In general, the results are in agreement with wind-tunnel measurements which have been-reported in the literature. For example, the variation the square root of p(sup 2)/q times the square root of p(sup 2) is the root mean square of the wall-pressure fluctuation, and q is the free-stream dynamic pressure) with Reynolds number was found to be essentially constant for the forward-fuselage-section boundary layer, while variations at the wing station were probably unduly affected by the microphone diameter (5/8 in.), which was large compared with the boundary-layer thickness.

Mull, Harold R.↗
NASA NTRSDate not supplied

Design Under Uncertainty with Design-Dependent Uncertain Variables

Uncertainty quantification (UQ) can provide a more robust understanding of a system, leading to better informed decisions earlier in the design process. The additional information that UQ provides can be leveraged during a design optimization process known as design under uncertainty that, when incorporated with multidisciplinary design and optimization, can become computationally infeasible due to the large number of responses required for meaningful results. Previous work addressed reducing the computational expense in design under uncertainty by incorporating analytic derivatives throughout polynomial chaos expansion. Although this allows design under uncertainty to be feasible for more systems, some multidisciplinary systems have design-dependent uncertain variables. This paper details an implementation of design dependent uncertain variables in a manner than preserves derivatives required for efficient gradient-based optimization throughout the process. Two analytic examples of design-dependent uncertain variables are given: the first transforms a uniform uncertain variable with one design variable and the second transforms a normal uncertain variable with two design variables. The polynomial chaos expansion (PCE) results are comparable to both the Monte Carlo (MC) results and the analytic results for the two examples. A case study that maximizes the lift-to-drag ratio with a design-dependence between the wing leading edge sweep angle and uncertain parameter percentage of laminar flow is compared to a MC and alternative optimization formulations. This paper demonstrates that design-dependent uncertain variables are valid and hold throughout PCE.

Joanna N Schmidt↗
NASA NTRSDate not supplied

Design Under Uncertainty with Design-Dependent Uncertain Variables

Uncertainty quantification (UQ) can provide a more robust understanding of a system, leading to better informed decisions earlier in the design process. The additional information that UQ provides can be leveraged during a design optimization process known as design under uncertainty that, when incorporated with multidisciplinary design and optimization, can become computationally infeasible due to the large number of responses required for meaningful results. Previous work addressed reducing the computational expense in design under uncertainty by incorporating analytic derivatives throughout polynomial chaos expansion. Although this allows design under uncertainty to be feasible for more systems, some multidisciplinary systems have design-dependent uncertain variables. This paper details an implementation of design dependent uncertain variables in a manner than preserves derivatives required for efficient gradient-based optimization throughout the process. Two analytic examples of design-dependent uncertain variables are given: the first transforms a uniform uncertain variable with one design variable and the second transforms a normal uncertain variable with two design variables. The polynomial chaos expansion (PCE) results are comparable to both the Monte Carlo (MC) results and the analytic results for the two examples. A case study that maximizes the lift-to-drag ratio with a design-dependence between the wing leading edge sweep angle and uncertain parameter percentage of laminar flow is compared to a MC and alternative optimization formulations. This paper demonstrates that design-dependent uncertain variables are valid and hold throughout PCE.

robust design↗
NASA NTRS1987

High angle of attack hypersonic aerodynamics

A new aerodynamics force model is presented which is based on modified Newtonian theory and empirical correlations. The algebraic model was developed for complete vehicles from take off to orbital speeds and for large angles of attack. Predictions are compared to results for a wind tunnel model at a Mach number of 20, and the full scale Shuttle Orbiter for Mach numbers from 0.25 to 20 for angles of attack from 0 to 50 deg. The maximum shuttle orbiter lift/drag at Mach 10 and 20 is 1.85 at 20-deg angle-of-attack. Aerodynamic force predictions are made for a transatmospheric vehicle, which is a derivative of the Shuttle Orbiter, for Mach numbers from 4 to 27 at angles of attack from 5 to 40 deg. Predicted aerodynamic force data indicate that lift/drag ratios of 5.2 at Mach number 10 and 3.6 at Mach number 26 are obtainable. Changes in force coefficients with changes in: nose angle, sweep angle, and (volume exp 2/3)/planform area are quantified for Mach numbers of 10 and 26. Lift/drag ratios increase with decreasing nose angle and (volume exp 2/3)/planform area and increasing wing sweep angle. Lift/drag ratios are independent of these variables for angles of attack in excess of 20 deg at Mach 10 and 30 deg at Mach 26.

Harloff, Gary J.↗
NASA NTRS1939

Radiator Design and Installation

The fundamental principles of fluid flow, pressure losses, and heat transfer have been presented and analyzed for the case of a smooth tube with fully developed turbulent flow. These equations apply to tubes with large length-diameter ratios where the f1ow is at a high Reynolds Number. The error introduced by using these equations increases as the magnitude of the tube length and the air-flow Reynolds Number approaches the values encountered in modern radiator designs. Accordingly, heat-transfer tests on radiator sections were made and the results are presented in nondimensional form to facilitate their use and for comparison with other heat-transfer data. In addition, pressure losses were measured along smooth tubes of circular, square, and rectangular cross section and the results were also correlated and are presented in nondimensional form. The problem of a radiator design for a particular installation is solved, the experimental heat-transfer and pressure-loss data being used, on a basis of power chargeable to the radiator for form drag, for propelling the weight, and for forcing the air through the radiator. The case of an installation within a wing or an engine nacelle is considered. An illustration of radiator design is carried through for an arbitrary set of conditions. Sufficient detail is given to enable the reader to reproduce the analysis for any given case.

Brevoort, M.J.↗
NASA NTRS1961

Experimental Investigation of a Hypersonic Glider Configuration at a Mach Number of 6 and at Full-Scale Reynolds Numbers

The aerodynamic characteristics of a hypersonic glider configuration, consisting of a slender ogive cylinder with three highly swept wings, spaced 120 apart, with the wing chord equal to the body length, were investigated experimentally at a Mach number of 6 and at Reynolds numbers from 6 to 16 million. The objectives were to evaluate the theoretical procedures which had been used to estimate the performance of the glider, and also to evaluate the characteristics of the glider itself. A principal question concerned the viscous drag at full-scale Reynolds number, there being a large difference between the total drags for laminar and turbulent boundary layers. It was found that the procedures which had been applied for estimating minimum drag, drag due to lift, lift curve slope, and center of pressure were generally accurate within 10 percent. An important exception was the non-linear contribution to the lift coefficient which had been represented by a Newtonian term. Experimentally, the lift curve was nearly linear within the angle-of-attack range up to 10 deg. This error affected the estimated lift-drag ratio. The minimum drag measurements indicated that substantial amounts of turbulent boundary layer were present on all models tested, over a range of surface roughness from 5 microinches maximum to 200 microinches maximum. In fact, the minimum drag coefficients were nearly independent of the surface smoothness and fell between the estimated values for turbulent and laminar boundary layers, but closer to the turbulent value. At the highest test Reynolds numbers and at large angles of attack, there was some indication that the skin friction of the rough models was being increased by the surface roughness. At full-scale Reynolds number, the maximum lift-drag ratio with a leading edge of practical diameter (from the standpoint of leading-edge heating) was 4.0. The configuration was statically and dynamically stable in pitch and yaw, and the center of pressure was less than 2-percent length ahead of the centroid of plan-form area.

Seiff, Alvin↗
DOE OSTI2022

Active galactic nuclei signatures in Red Geyser galaxies from Gemini GMOS-IFU observations

ABSTRACT Red Geysers are quiescent galaxies with galactic scale ionized outflows, likely due to low-luminosity active galactic nuclei (AGN). We used Gemini GMOS-IFU observations of the inner ∼1–3 kpc of nine Red Geysers selected from the MaNGA survey to study the gas ionization and kinematics. The emission-line ratios suggest the presence of Seyfert/LINER (Low Ionization Nuclear Emission Region) nuclei in all sources. Two galaxies show Hα equivalent width ( Hα EW) larger than 3 Å (indicative of AGN ionization) within an aperture 2 ${_{.}^{\prime\prime}}$5 of diameter (1.3–3.7 kpc at the distance of galaxies) for MaNGA data, while with the higher resolution GMOS data, four galaxies present Hα EW>3 Å within an aperture equal to the angular resolution (0.3–0.9 kpc). For two objects with GMOS-IFU data, the Hα EW is lower than 3 Å but larger than 1.5 Å, most probably due to a faint AGN. The spatially resolved electron density maps show values between 100 and 3000 cm −3 and are consistent with those determined in other studies. The large (MaNGA) and the nuclear scale (GMOS-IFU) gas velocity fields are misaligned, with a kinematic position angle difference between 12° and 60°. The [N ii] λ6583 emission-line profiles are asymmetrical, with blue wings on the redshifted side of the velocity field and red wings on the blueshifted side. Our results support previous indications that the gas in Red Geysers is ionized by an AGN, at least in their central region, with the presence of outflows, likely originating in a precessing accretion disc.

Ilha, Gabriele S. (ORCID:0000000191006106)↗
NASA NTRS1991

Static performance of a multiaxis thrust vectoring cruciform nozzle

A multiaxis thrust vectoring nozzle designed to have equal flow turning capability in pitch and yaw is investigated, and methods for obtaining large multiaxis thrust vector angles without physical vectoring flap interference are studied. These methods include restricting the pitch flaps from the path of the yaw flaps, and shifting the flow path at the throat off the nozzle centerline in order to permit larger pitch flap deflections without interfering with operation of the yaw flaps. The results obtained show that unvectored performance of the cruciform nozzle is only slightly lower (0.5 to 1.0 percent) than that of previously tested axisymmetric and nonaxisymmetric nozzles, despite the complex internal geometry. The shifted-throat nozzle design has larger thrust vector angles at the nozzle pressure ratio of peak resultant thrust efficiency, but has 1 to 2 percent lower peak thrust performance than the restricted-flap nozzle design.

Wing, David J.↗
NASA NTRS1952

Longitudinal Stability and Control Characteristics from a Flight Investigation of a Cruciform Canard Missile Configuration Having an Exposed Wing-canard Area Ratio of 16:1

A flight investigation has been made to determine the longitudinal stability and control characteristics of a 60 0 delta-wing-canard missile configuration with an exposed wing-canard area ratio of 16:1. The results presented include the longitudinal stability derivatives, control effectiveness, and drag characteristics for a Mach number range of 0.75 to 1.80 and are compared with the results of a similar configuration having larger 6ontrols. Stability characteristics are also presented from the flights of an interdigitated canard configuration at a Mach number of 2.08 and a wing-body configuration at Mach numbers of 1.25 to 1.45. The stability derivatives varied gradually with Mach number with the exception of the damping-in-pitch derivative. Aerodynamic damping in pitch decreased to a minimum at a Mach number of 1.0 3, then increased to a peak value at a Mach number of 1.26 followed by a gradual decrease at higher Mach numbers. The aerodynamic-center location of the in-line canard configuration shifted rearward 13 percent of the mean aerodynamic chord at transonic speeds. The pitching-moment curve slope was 25 percent greater for the model having no canards than for the in-line configuration. No large effects of interdigitation were noted in the stability derivatives. Pitching effectiveness of the in-line configuration was maintained throughout the Mach number range. A comparison of the stability and control characteristics of two canard configurations having different area controls showed that decreasing the control area 44 percent decreased the pitching effectiveness proportionally, shifted the aerodynamic-center location rearward 9 to 14 percent of the mean aerodynamic chord, and reduced the total hinge moments required for 10 trimmed flight about 50 percent at transonic speeds.

Moul, Martin T↗
NASA NTRS1959

Transonic Aerodynamic Characteristics of a 45 deg Swept Wing Fuselage Model with a Finned and Unfinned Body Pylon Mounted Beneath the Fuselage or Wing, Including Measurements of Body Loads

An investigation of a model of a standard size body in combination with a representative 45 deg swept-wing-fuselage model has been conducted in the Langley 8-foot transonic pressure tunnel over a Mach number range from 0.80 to 1.43. The body, with a fineness ratio of 8.5, was tested with and without fins, and was pylon-mounted beneath the fuselage or wing. Force measurements were obtained on the wing-fuselage model with and without the body, for an angle-of-attack range from -2 deg to approximately 12 deg and an angle-of-sideslip range from -8 deg to 8 deg. In addition, body loads were measured over the same angle-of-attack and angle-of-sideslip range. The Reynolds number for the investigation, based on the wing mean aerodynamic chord, varied from 1.85 x 10(exp 6) to 2.85 x 10(exp 6). The addition of the body beneath the fuselage or the wing increased the drag coefficient of the complete model over the Mach number range tested. On the basis of the drag increase per body, the under-fuselage position was the more favorable. Furthermore, the bodies tended to increase the lateral stability of the complete model. The variation of body loads with angle of attack for the unfinned bodies was generally small and linear over the Mach number range tested with the addition of fins causing large increases in the rates of change of normal-force coefficient and nose-down pitching-moment coefficient. The variation of body side-force coefficient with sideslip for the unfinned body beneath the fuselage was at least twice as large as the variation of this load for the unfinned body beneath the wing. The addition of fins to the body beneath either the fuselage or the wing approximately doubled the rate of change of body side-force coefficient with sideslip. Furthermore, the variation of body side-force coefficient with sideslip for the body beneath the wing was at least twice as large as the variation of this load with angle of attack.

Wornom, Dewey E.↗
DOE OSTI2020

Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Li, Sibo↗